### Re: [R] R Anova

Keizer_61 wrote: I am really struggling with this question. Three students take part of an experiment. student smoking non-smokingcancer 110.5 7.5 6.5 2 9.5 6.5 8.4 3 8.5 7.2

### [R] How to solve difficult equations?

This below is not solvable with uniroot to find a: fn=function(a){ b=(0.7/a)-a (1/(a+b+1))-0.0025 } uniroot(fn,c(-500,500)) gives Error in uniroot(fn, c(-500, 500)) : f() values at end points not of opposite sign I read R-help posts and someone wrote a function:

### Re: [R] simulate values

Soare Marcian-Alin wrote: Hello, I want to simulate 100 values of the ARMA Process with this function: x[i] = 0.5 * x[i-1] + 0.2 * x[i-2] + x[i] + 0.9 * x[i-1] + 0.2 * x[i-2] + 0.3 * x[i-3] There is no kind of noise in your model, hence no need to do any simulation so far ... Uwe

### [R] How to solve difficult equations?

I don't see the problem, except that you might want to think about what the error message is telling you. A little exploration of your function always helps, too. ss - seq(-2, 2, len = 100) plot(ss, fn(ss), type = l) uniroot(fn, c(-1, 1)) Erreur dans uniroot(fn, c(-1, 1)) : f() values at

### Re: [R] Analysis of Variance

CrazyJoe wrote: Hello Blind toy_Car toy_truck toy_boat 1 6.3 7.5 5.4 2 3.4 8.1 6.1 3 2.2 4.4 5.1 How do we calculate the F-statistic in R. Any help is really appreciated. You may want to try to

### Re: [R] How to solve difficult equations?

plot(fn,-1,1) could be helpful, hth, Ingmar On 25 Apr 2007, at 09:15, francogrex wrote: This below is not solvable with uniroot to find a: fn=function(a){ b=(0.7/a)-a (1/(a+b+1))-0.0025 } uniroot(fn,c(-500,500)) gives Error in uniroot(fn, c(-500, 500)) : f() values at end points not of

### Re: [R] simulate values

On Wed, 25 Apr 2007, Uwe Ligges wrote: Soare Marcian-Alin wrote: Hello, I want to simulate 100 values of the ARMA Process with this function: x[i] = 0.5 * x[i-1] + 0.2 * x[i-2] + x[i] + 0.9 * x[i-1] + 0.2 * x[i-2] + 0.3 * x[i-3] There is no kind of noise in your model, hence no need to

### [R] Unsubscription Confirmation

Thank you for subscribing. You have now unsubscribed and no more messages will be sent. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

### [R] for loops

Hello everybody I'm very new at using R so probably this is a very stupid question. I have a matrix of p columns and I have to calculate for each of them the two sample t-statistic and p-value and to save the results into two different vectors. I have divided my matrix into two submatrices:

### Re: [R] How to solve difficult equations?

On 25-Apr-07 07:15:55, francogrex wrote: This below is not solvable with uniroot to find a: fn=function(a){ b=(0.7/a)-a (1/(a+b+1))-0.0025 } uniroot(fn,c(-500,500)) gives Error in uniroot(fn, c(-500, 500)) : f() values at end points not of opposite sign I read R-help posts and

### [R] Permutations with samr

I have a question regarding the samr package. For a 2 class unpaired problem, with sample 1 of size N1 and sample 2 of size N2, samr computes at most (N1+N2)! permutations of the two samples (if the user-supplied parameter nperms allows it). However, there are only (N1+N2)/(N1!*N2!) DISTINCT

### [R] omit y=zero line in histogram

Dear all, hist ( ) plots a horizontal line at y=0 when the respective bin is empty. I can deal with this by modifying the hist object before plotting it (x$density[x$density == 0] - NA), but I'm sure I've seen a more elegant way. Perhaps this was in truehist (MASS). I have looked but can't find

### [R] RE : for loops

You can see with this simple example. matrix.t.test-function(mx){ p-dim(mx)[2] #number of column in the matrix n-dim(mx)[1] #number of row n.tests- p*(p-1)/2 #Number of tests to be done tested.var -rep(,n.tests) #Keep rang of tested column r.t.stat-rep(0, n.tests)#contain t.stat r.p.val

### [R] correlation table

hello, is it possible to create a correlation table between factors? ___ [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

### [R] unit testing frameworks for R

Greetings! After a quick look at current programming tools, especially with regards to unit-testing frameworks, I've started looking at both butler and RUnit. I would be grateful to receieve real world development experience and opinions with either/both.Please send to me directly (yes,

### [R] attributable risk

Hi everybody, Does anyone know a function to compute the attributable risk of a factor in a logistic regression or a proportional hazard cox model (both with confounding variables)? I need also obtain the confidence interval. Thanks in advance. Isaac Subirana. La informaciÃ³

### [R] regular expressions with grep() and negative indexing

Dear R-helpers, Does anyone know how to use regular expressions to return vector elements that don't contain a word? For instance, if I have a vector x - c(seal.0,seal.1-exclude) I'd like to get back the elements which do not contain the word exclude, using something like (I know this doesn't

### [R] Plotting minimum spanning tree in graph/RBGL

I wonder if anyone could tell me how I can plot a Minimum Spanning Tree using the functions provided in the graph and RBGL packages? I am able to build the MST using the following set of commands: library(graph) library(RBGL) x - read.table(h:/pole.tab,header=T,row.names=1) y - dist(x)

### [R] Help with saptial analysis (cluster)

Hi R-users I'm a beginner with R and statistics, so I need some help to start my data analysis. I've been reading some docs and tutorials on R and cluster analysis. I've got a large dataset (102000 points) with values of longitude, latitude and temperature and want to see if I can find groups

### Re: [R] Problem installing Rmpi with lam on SGI SLES9

On 24/04/07, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Tue, 24 Apr 2007, Hendrik Fuß wrote: Hi, I've been trying here to install Rmpi on an SGI IA-64 machine with 64 processors, running SuSE Linux Enterprise Server 9, R 2.4.0 and lam-mpi 7.1.3. While I've read of similar problems

### Re: [R] regular expressions with grep() and negative indexing

Stephen Tucker wrote: Dear R-helpers, Does anyone know how to use regular expressions to return vector elements that don't contain a word? For instance, if I have a vector x - c(seal.0,seal.1-exclude) I'd like to get back the elements which do not contain the word exclude, using something

### Re: [R] regarding 3d Bar Plot

On 4/24/2007 9:38 AM, [EMAIL PROTECTED] wrote: [EMAIL PROTECTED] wrote: I have data in a two dimensional table. each row of the data adds upto 100 ( hence they are percentages ). it can be interpreted as like this A - I are the matches and P - X are the players. Thus Player P scored 20% of

### Re: [R] regular expressions with grep() and negative indexing

Find the ones that match and then remove them from the full set with 'setdiff'. x - c(seal.0,seal.1-exclude) x.match - grep(exclude, x) # find matches x.match [1] 2 setdiff(seq_along(x), x.match) # exclude the matches [1] 1 On 4/25/07, Peter Dalgaard [EMAIL PROTECTED] wrote: Stephen

### Re: [R] regarding 3d Bar Plot

Hi Duncan I am restating the problem and thanks you for sending me such a good function histogram in 3d. Thanks for that but i think my problem has been misinterpreted. I just wanted a simple 3d bar Plot. Although I have not written anything for R but i will surely like to contribute to R and

### Re: [R] correlation table

Possible: yes - just calcuate correlation of as.numeric(your.factors) Meaningful: no It will depend on the coding for your factors, which may be absolutely arbitrally... Petr elyakhlifi mustapha napsal(a): hello, is it possible to create a correlation table between factors?

### Re: [R] regarding 3d Bar Plot --- correction.

Hi Rolf, If it is possible then please share the code as i am not able to locate any pointers. Thanks in advance :-) cheers and chiao regards -gaurav [EMAIL PROTECTED] 24-04-07 07:27 PM To [EMAIL PROTECTED], r-help@stat.math.ethz.ch cc Subject Re: [R] regarding 3d Bar Plot ---

### Re: [R] Coercing data types for use in model.frame

You could try to follow the code in the dyn package. It intercepts model.frame calls involving time series objects so that it can align lagged objects: e.g. z - ts(seq(10)^2) library(dyn) dyn$lm(z ~ lag(z, -1)) It transforms the last line above to: dyn(lm(dyn(z ~ lag(z, -1))) and the inner

### Re: [R] RE : for loops

I know that a looping solution was requested, but this is exactly what apply() should be used for... Petr justin bem napsal(a): You can see with this simple example. matrix.t.test-function(mx){ p-dim(mx)[2] #number of column in the matrix n-dim(mx)[1] #number of row n.tests-

### Re: [R] A very simple question

If by objects you mean the packages you installed and assuming you are on Windows then there are two batch files movedir.bat and copydir.bat in the batchfiles distribution that will either move them from an older version (they won't be available to the older version any more but its much faster

### Re: [R] Help with saptial analysis (cluster)

Dear Fransico, The distance matrix would be 102000 x 102000. So it would contain 1040400 values. If you need one bit for each value, this would requier 9,7 GB. So the distance matrix won't fit in the RAM of your computer. Cheers, Thierry

### Re: [R] prelim.norm() function not working

Thank you very much, that was indeed the problem. (And now that I read more carefully the help page, it did in fact say that the input was a data matrix and not a data frame.) Brant -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley Sent:

### [R] program avail. for simulating spatial patterns?

Hi all, I am wondering if there is a function available in R for simulating spatial distribution of objects (plants in this case) in order to simulate sampling of a population . Specifically, I would like to be able to change the spatial correlation of individuals. I don't want to reinvent the

### [R] assigning two conditions to grep()

Hi, i have a problem in assigning 2 conditions to grep() , my data look like this: DA 24 N7 Rad= 3.4 20 Sac= 0.93 Acc= 4.76 DA 24 N7 Rad= 3.4 14 Sac= 0.65 Acc= 3.33 DA 24 N7 Rad= 3.4 3 Sac= 0.14 Acc= 0.71 DA 24 N7 Rad= 3.4 11 Sac= 0.51 Acc= 2.62 DG 23 N7 Rad= 3.4 8 Sac= 0.37 Acc= 1.91 DG 23 N7

### Re: [R] regarding 3d Bar Plot

On 4/25/2007 7:56 AM, [EMAIL PROTECTED] wrote: Hi Duncan I am restating the problem and thanks you for sending me such a good function histogram in 3d. Thanks for that but i think my problem has been misinterpreted. I just wanted a simple 3d bar Plot. Although I have not written anything

### Re: [R] program avail. for simulating spatial patterns?

Have a look at packages in the spatial taskview (spatstat, splancs). Cheers, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Reseach Institute for Nature and Forest Cel biometrie, methodologie en

### Re: [R] regarding 3d Bar Plot

I am restating the problem and thanks you for sending me such a good function histogram in 3d. Thanks for that but i think my problem has been misinterpreted. I just wanted a simple 3d bar Plot. Although I have not written anything for R but i will surely like to contribute to R and if i can

### [R] R CMD CHECK and require() / library() methods

Hi, I have a piece of code that decides at runtime whether to load a data package (and which package to load). This is then done with a call to: library(x) (where x is a character variable containing the package name). This causes R CMD check to throw out a warning: 'library' or 'required'

### [R] How to identify and exclude the outliers with R?

Hello, everyone, I want to ask a simple question. If I have a set of data,and I want to identify how many outliers there are in the data.Which packages and functions can I use? Thanks. Shao chunxuan. [[alternative HTML version deleted]] __

### [R] creating random numbers

l want to create a column of 1 and 2 randomly what command should l use eg treatment strata 1 1 2 0 1 1 2 1 2 0 2 1 2

### Re: [R] regarding 3d Bar Plot

On 4/25/2007 9:26 AM, hadley wickham wrote: I am restating the problem and thanks you for sending me such a good function histogram in 3d. Thanks for that but i think my problem has been misinterpreted. I just wanted a simple 3d bar Plot. Although I have not written anything for R but i will

### Re: [R] NA and NaN randomForest

Hi Clayton, If you use the formula interface, then it should do what you want: R library(randomForest) randomForest 4.5-18 Type rfNews() to see new features/changes/bug fixes. R iris1 - iris[-(1:5),] R iris2 - iris[1:5,] R iris2[1, 3] - NA R iris2[3, 1] - NA R iris.rf - randomForest(Species ~

### Re: [R] creating random numbers

sample(1:2, 10, replace=TRUE) -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of raymond chiruka Sent: Wednesday, April 25, 2007 9:45 AM To: r-help@stat.math.ethz.ch Subject: [R] creating random numbers l want to create a column of 1 and 2

### Re: [R] Random Number Generator of Park and Miller

On Tue, 24 Apr 2007, gracezhang wrote: Hi, I failed to search for R package providing random number generator of Park and Miller. Anyone know any R package supporting this kind of function? rng.lcg-function(x,p1=16807,p2=0,N=2147483647){(x*p1+p2)%%N} Dave -- Dr. David Forrest [EMAIL

### [R] print format - fixed number of digits

Hi, is it possible to format real (double) numbers in a data frame to an exact format? I need something like format(..., digits=5) but this is a suggestion only and I need exactly 5 digits. Example where there are more than 5 digits printed follows. I can do it via

### [R] Split column of concatenated data

Hi I have a column of concatenated information stored in an RG object in the limma package and I need to split this information and then paste the first two pieces of data in each case back into two columns of the RG object. This is how I am currently doing this

### [R] barchart producing incorrect number of barcharts when columns renamed

Hi everybody, I'm having problems with the barchart command in the lattice package. I'm creating barcharts from matrices with with anything from 20 to 71 columns. When I leave the column names alone, that is they are set in the read.table command (and inherited by subsequent commands) the

### Re: [R] simulate values

Hello, Yes I tried arima.sim and everything works fine. Thanks for the help! Alin Soare 2007/4/25, Leeds, Mark (IED) [EMAIL PROTECTED]: he's just being a wise guy bcause I'm sure you meant to have an epsilon_t on the end. And he Surely knows that also. Did you Check out ?arima.sim.

### [R] levelplot and unequal cell sizes

I am using levelplot() from lattice with grids that have unequal cell sizes. This means that the boundary between two cells is not always half-way between nodes, as levelplot() assumes. The result is that some cell sizes are rendered incorrectly, which can be painfully obvious if using

### Re: [R] assigning two conditions to grep()

janek0 wrote: Dnia 25-04-2007, śro o godzinie 15:04 +0200, Abi Ghanem josephine napisał(a): use data[grep(^24|N7, data$V2)] and see ?regexp thanks for replying but actually my problem is that the column containing 23 is in data$V2 and the one containing N7 is in data$V3, so the line doen't

### Re: [R] R CMD CHECK and require() / library() methods

On Wed, 25 Apr 2007, Crispin Miller wrote: Hi, I have a piece of code that decides at runtime whether to load a data package (and which package to load). This is then done with a call to: library(x) (where x is a character variable containing the package name). This causes R CMD check

### Re: [R] How to solve difficult equations?

At 03:15 AM 4/25/2007, francogrex wrote: This below is not solvable with uniroot to find a: fn=function(a){ b=(0.7/a)-a (1/(a+b+1))-0.0025 } uniroot(fn,c(-500,500)) gives Error in uniroot(fn, c(-500, 500)) : f() values at end points not of opposite sign I read R-help posts and someone wrote a

### Re: [R] creating random numbers

?sample sample(c(1,2),10,replace=TRUE) P. raymond chiruka napsal(a): l want to create a column of 1 and 2 randomly what command should l use eg treatment strata 1 1 2 0 1 1 2

### Re: [R] R CMD CHECK and require() / library() methods

Many thanks! -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley Sent: 25 April 2007 16:03 To: Crispin Miller Cc: R-help@stat.math.ethz.ch Subject: Re: [R] R CMD CHECK and require() / library() methods On Wed, 25 Apr 2007,

### Re: [R] regular expressions with grep() and negative indexing

Peter Dalgaard wrote: Stephen Tucker wrote: Dear R-helpers, Does anyone know how to use regular expressions to return vector elements that don't contain a word? For instance, if I have a vector x - c(seal.0,seal.1-exclude) I'd like to get back the elements which do not contain the

### Re: [R] assigning two conditions to grep()

On 25/04/07, Abi Ghanem josephine [EMAIL PROTECTED] wrote: Hi, i have a problem in assigning 2 conditions to grep() , my data look like this: DA 24 N7 Rad= 3.4 20 Sac= 0.93 Acc= 4.76 DA 24 N7 Rad= 3.4 14 Sac= 0.65 Acc= 3.33 DA 24 N7 Rad= 3.4 3 Sac= 0.14 Acc= 0.71 DA 24 N7 Rad= 3.4 11

### [R] Getting Confused

Hello, I'm getting confused with my experience of R installing. I had R installed on January without any trouble. (I just had to install gcc4.1.1) Now I'd like to install a packages which requires tcl/tk. So basically I need to reconfigure and re install R right after having installed

### [R] Box Ljung Statistics

Hi All R Experts, I met with below mentioned statistics in paper Stock Index Volatility Forecasting with High Frequency Data by Eugenie Hol, Siem Jan Koopman http://ideas.repec.org/p/dgr/uvatin/20020068.html I would like to ask that what is Box-Ljung portmantacau statistic based on N squared

### [R] Unsubscription Confirmation

Thank you for subscribing. You have now unsubscribed and no more messages will be sent. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

### Re: [R] NA and NaN randomForest

Hi Andy, It worked for classification, but not regression. For example: iris1 - iris[-(1:5),] iris2 - iris[(1:5),] iris2[1,3] - NA iris2[3,1] - NA iris_sum - sum (iris$Sepal.Length + iris$Sepa.Width + iris$Petal.Length + iris$Petal.Width) iris_sum1 - iris_sum[-(1:5)] iris_sum2 -

### Re: [R] regular expressions with grep() and negative indexing

I use regexpr() instead of grep() in cases like this, e.g.: x2[regexpr(exclude,x2)==-1] (regexpr returns a vector of the same length as character vector given it, so there's no problem with it returning a zero length vector) -- Tony Plate Peter Dalgaard wrote: Stephen Tucker wrote: Dear

### [R] R News, volume 7, issue 1 is now available

Hi The October 2006 issue of R News is now available on CRAN under the Documentation/Newsletter link. Torsten (on behalf of the R News Editorial Board) ___ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-announce

### Re: [R] How to identify and exclude the outliers with R?

It depends on the nature of your data set. There is a package simply called 'outliers', which has the Grubbs/Dixon/Cochran tests. There is also the Bonferroni outlier test in 'car' package. I'm sure there are more in the hundreds of packages on CRAN. HTH Horace Shao [EMAIL PROTECTED]

### Re: [R] Getting Confused [Broadcast]

If you are serious in getting useful help, please do try to follow suggestions in the Posting Guide. You have not told us anything about your OS, the versions of R you tried to install, and exactly what you typed to build/install them. Many Linux distro by default do not install the Fortran part

### Re: [R] R News, volume 7, issue 1 is now available

Torsten Hothorn [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] The October 2006 issue of R News is now available on CRAN under the Documentation/Newsletter link. Direct links are useful: R News http://cran.r-project.org/doc/Rnews/ April 2007 Issue:

### Re: [R] Problem installing Rmpi with lam on SGI SLES9

Hendrik, Are you starting the lam daemons before starting R? % lamboot You might need to specify a 'hosts' argument to lamboot. The default way Rmpi calls lamboot is with no arguments, and this might simply create a single lam daemon. I don't usually use papply, but glancing at it's code

### Re: [R] levelplot and unequal cell sizes

On 4/25/07, Waichler, Scott R [EMAIL PROTECTED] wrote: I am using levelplot() from lattice with grids that have unequal cell sizes. This means that the boundary between two cells is not always half-way between nodes, as levelplot() assumes. The result is that some cell sizes are rendered

### Re: [R] R News, volume 7, issue 1 is now available

Earl F. Glynn wrote: Torsten Hothorn [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] The October 2006 issue of R News is now available on CRAN under the Documentation/Newsletter link. Direct links are useful: R News http://cran.r-project.org/doc/Rnews/ April 2007

### Re: [R] Coercing data types for use in model.frame

Prof Brian Ripley wrote: Moved to R-devel What is the 'data class'? In particular what is its underlying type? And where in model.frame[.default] are you trying to use it (in the formula, data, in ..., etc). This is an example of where some reproducible code and the error

### Re: [R] Box Ljung Statistics

Gaurav, I met with below mentioned statistics in paper Stock Index Volatility Forecasting with High Frequency Data by Eugenie Hol, Siem Jan Koopman http://ideas.repec.org/p/dgr/uvatin/20020068.html I would like to ask that what is Box-Ljung portmantacau statistic based on N squared

### [R] Scripting graph generation

Hi, I'm looking to automate the generation of some graphs in R. I can't seem to figure out how to script R, and change the output device of hist() or plot() to create a .gif or .png file. This seems like something that is probably really simple, and I've just overlooked the call do do it. Can

### [R] GLS terminology question not related to R

This is a terminology question not related to R. The literature often says that OLS is inefficient relative to GLS if the residuals in the system are correlated ( and the RHS sides of each are not identical ). Does this mean that OLS overestimates residual and coefficient variances ,

### [R] [R-pkgs] new package adegenet

The new package *adegenet* (linked to the ade4 package for multivariate analysis) has been released on CRAN. Its main focus is on molecular marker data handling for multivariate analysis. Adegenet offers data import/export functions (from GENETIX, Genepop, Fstat and to the packages genetics

### Re: [R] GLS terminology question not related to R

On Wed, 25 Apr 2007, Leeds, Mark (IED) wrote: This is a terminology question not related to R. The literature often says that OLS is inefficient relative to GLS if the residuals in the system are correlated ( and the RHS sides of each are not identical ). Does this mean that OLS overestimates

### Re: [R] levelplot and unequal cell sizes

On 4/25/07, Waichler, Scott R [EMAIL PROTECTED] wrote: I am using levelplot() from lattice with grids that have unequal cell sizes. This means that the boundary between two cells is not always half-way between nodes, as levelplot() assumes. levelplot() is not supposed to make any such

### [R] help

Hi all, I have 2 questions: 1)How do I calculate the mean on an imported txt file? I've imported the file below and that's what it looks like imported. How do I then calcuate the mean, median, or mode on the column LeafArea using the desktop R package? Any help would be greatly appreciated!!

### Re: [R] Scripting graph generation

On Wed, 2007-04-25 at 13:21 -0400, Mike Huber wrote: Hi, I'm looking to automate the generation of some graphs in R. I can't seem to figure out how to script R, and change the output device of hist() or plot() to create a .gif or .png file. This seems like something that is probably

### [R] Самый доступный и недорогой метод похудения

«ÇÎËÎÒÀß ÑÅÐÜÃÀ» - ÈÇÁÀÂËÅÍÈÅ ÎÒ ËÈØÍÈÕ ÊÈËÎÃÐÀÌÌÎÂ (4 ÊÃ Â ÌÅÑ.)! ÌÈËËÈÎÍÛ ËÞÄÅÉ ÈÇÁÀÂÈËÈÑÜ ÎÒ ËÈØÍÈÕ ÊÈËÎÃÐÀÌÌÎÂ ÁËÀÃÎÄÀÐß ÇÎËÎÒÎÉ_ÑÅÐÜÃÅ ÁÅÇ ÔÈÇÈ×ÅÑÊÈÕ ÓÏÐÀÆÍÅÍÈÉ, ÈÑÒßÇÀÞÙÈÕ ÄÈÅÒ, ÒÀÁËÅÒÎÊ!! ÑÒÎÈÌÎÑÒÜ «ÇÎËÎÒÎÉ ÑÅÐÜÃÈ» ÂÑÅÃÎ 3000 ÐÓÁ.!! Ïðèíöèï ìåòîäà Çîëîòàÿ Ñåðüãà ñîçäàí

### [R] dmnorm not meant for 1024-dimensional data?

Hello, I have some data generated by a simple mixture of Gaussians (more like K-means) and (as a test) am using dmnorm to calculate the probability of each data point coming from each Gaussian. However, I get only zero probabilities. This code works in low dimensions (tried 2 and 3 already). I

### Re: [R] negative number to positive number

Quoting H. Paul Benton [EMAIL PROTECTED]: Hello all, I know this is a pretty easy question but I can't find it in S poetry or R help. How can I make a negative number positive. Such as -5 to be +5 I tried +(-5), but that didn't work. So no, I don't mean taking a -5^2 just to get a

### [R] Help on 'find.BIB' function

Hello everyone, I am trying to use the 'find.BIB' function to construct a balanced incomplete block design. When I ran the example given in the help file (find.BIB(10,30,4)), I obtained the following error message: Error in optBlock(~., withinData = factor(1:trt), blocksize = rep(k,

### Re: [R] barchart producing incorrect number of barcharts when columns renamed

On 4/25/07, Colm G. Connolly [EMAIL PROTECTED] wrote: Hi everybody, I'm having problems with the barchart command in the lattice package. [... ... ...] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

### Re: [R] Help on 'find.BIB' function

Jason Parcon wrote: Hello everyone, I am trying to use the 'find.BIB' function to construct a balanced incomplete block design. When I ran the example given in the help file (find.BIB(10,30,4)), I obtained the following error message: Error in optBlock(~., withinData =

### [R] help

Hi, Would anyone know how to calculate the modal value of LeafArea? Thank-you very much!! Nat __ Hi all, I have 2 questions: 1)How do I calculate the mean on an imported txt file? I've imported the file below and that's what it looks like imported. How do I then calcuate

### Re: [R] negative number to positive number

?abs Clint BowmanINTERNET: [EMAIL PROTECTED] Air Dispersion Modeler INTERNET: [EMAIL PROTECTED] Air Quality Program VOICE: (360) 407-6815 Department of Ecology FAX:(360) 407-7534 USPS: PO Box

### Re: [R] Problem installing Rmpi with lam on SGI SLES9

On 25/04/07, Martin Morgan [EMAIL PROTECTED] wrote: Hendrik Fuß [EMAIL PROTECTED] writes: I'm trying to use the papply package. However, when I do: library(papply) papply(list(1:10, 1:15, 1:20), sum) 1 slaves are spawned successfully. 0 failed. master (rank 0, comm 1) of size

### Re: [R] GLS terminology question not related to R

Thomas Lumley wrote: On Wed, 25 Apr 2007, Leeds, Mark (IED) wrote: This is a terminology question not related to R. The literature often says that OLS is inefficient relative to GLS if the residuals in the system are correlated ( and the RHS sides of each are not identical ). Does this

### Re: [R] Fwd: dmnorm not meant for 1024-dimensional data?

On Wed, Apr 25, 2007 at 01:25:09PM -0500, Daniel Elliott wrote: Hello. Thank you for your mnormt package. Below is an email I sent to R-help regarding the function dmnorm. Any help you can provide would be greatly appreciated. Hello. You have not provided information on what

### [R] applying rbind to list elements

Hi, I have a list of n data.frames (or matrices) which I would like to convert to a single data.frame using rbind: x - rbind( l[[1]], l[[2]], l[[3]], l[[4]], ..., l[[n]] ) Is there a simple way to do this? thanks Hendrik -- Hendrik Fuß PhD student Systems Biology Research Group

### Re: [R] applying rbind to list elements

On 4/25/2007 4:09 PM, Hendrik Fuß wrote: Hi, I have a list of n data.frames (or matrices) which I would like to convert to a single data.frame using rbind: x - rbind( l[[1]], l[[2]], l[[3]], l[[4]], ..., l[[n]] ) Is there a simple way to do this? do.call(rbind, l). Duncan Murdoch

### Re: [R] applying rbind to list elements

I knew there would be a simple solution. thanks everybody. On 25/04/07, Duncan Murdoch [EMAIL PROTECTED] wrote: On 4/25/2007 4:09 PM, Hendrik Fuß wrote: Hi, I have a list of n data.frames (or matrices) which I would like to convert to a single data.frame using rbind: x - rbind(

### Re: [R] levelplot and unequal cell sizes

Hadley and Deepayan, Thank you for responding. Here is a simple example of what I'm talking about. It is a grid that is 5 cells wide by 2 cells tall. The width of the cells in the x-direction is variable; the cells at either end have width = 4 units, and the three cells in the middle have

### [R] dist label names

Hello, I am trying to do a multi-dimensional scaling of the World Bank's quality of governance indicators for the Balkan region. I am having trouble labelling my plot. Could some kind person help me out. How do I set the attribute Label by a variable (say, Code)? At present I get this:

### Re: [R] applying rbind to list elements

do.call(rbind, l) or, in the case of matrices, using the abind package: abind(l, along=1) library(abind) l - list(matrix(1:6, ncol=2), matrix(11:14, ncol=2)) abind(l, along=1) [,1] [,2] [1,]14 [2,]25 [3,]36 [4,] 11 13 [5,] 12 14 Hendrik Fuß wrote:

### [R] Sum of specific column

I have a data set that I have imported (not sure if that makes a difference) and I would like to calculate the sum of only specific columns. I have tried colSums(dataset, by=list(dataset$col5), dims=1) and I get an error of unused arguments I have also tried aggregate(dataset,

### [R] aggregate similar to SPSS

Hi, Does anyone know if: with R can you take a set of numbers and aggregate them like you can in SPSS? For example, could you calculate the percentage of people who smoke based on a dataset like the following: smoke = 1 non-smoke = 2 variable 1 1 1 2 2 1 1 1 2 2 2 2 2 2 When aggregated, SPSS

### Re: [R] levelplot and unequal cell sizes

On 4/25/07, Waichler, Scott R [EMAIL PROTECTED] wrote: Hadley and Deepayan, Thank you for responding. Here is a simple example of what I'm talking about. It is a grid that is 5 cells wide by 2 cells tall. The width of the cells in the x-direction is variable; the cells at either end have

### [R] Biostatistician Opportunities at Vanderbilt

The Department of Biostatistics at Vanderbilt University's School of Medicine has openings for biostatisticians at all levels. Details and application procedures may be found at http://biostat.mc.vanderbilt.edu/JobOpenings . For M.S. and B.S. biostatisticians we are especially interested in

### Re: [R] aggregate similar to SPSS

?table On Wednesday 25 April 2007 14:32, Natalie O'Toole wrote: Hi, Does anyone know if: with R can you take a set of numbers and aggregate them like you can in SPSS? For example, could you calculate the percentage of people who smoke based on a dataset like the following: smoke = 1

### Re: [R] aggregate similar to SPSS

Hi Nat, can I suggest, without offending, that you purchase and read Peter Dalgaard's Introductory Statistics with R or Michael Crawley's Statistics: An Introduction using R or Venables and Ripley's Modern Applied Statistics with S or Maindonald and Braun's Data Analysis and Graphics Using R: An

### [R] ANOVA results in R conflicting with results in other software packages

Hi, I'm wrestling with an analysis of a dataset, which I previously analyzed in SYSTAT, but am now converting to R and was doing a re-analysis. I noticed, however, that the same model yields different results (different sums of squares) from the two programs. I first thought this might be