Re: [R] R Anova

2007-04-25 Thread Uwe Ligges
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?

2007-04-25 Thread francogrex
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

2007-04-25 Thread Uwe Ligges
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?

2007-04-25 Thread ken knoblauch
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

2007-04-25 Thread Uwe Ligges
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?

2007-04-25 Thread Ingmar Visser
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

2007-04-25 Thread Prof Brian Ripley
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

2007-04-25 Thread Egoldsystem
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

2007-04-25 Thread [EMAIL PROTECTED]
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?

2007-04-25 Thread Ted Harding
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

2007-04-25 Thread Isabelle Rivals
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

2007-04-25 Thread Paul Artes
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

2007-04-25 Thread justin bem
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

2007-04-25 Thread elyakhlifi mustapha
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

2007-04-25 Thread anthony . rossini
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

2007-04-25 Thread SUBIRANA CACHINERO, ISAAC
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

2007-04-25 Thread Stephen Tucker
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

2007-04-25 Thread Wilson, Andrew
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)

2007-04-25 Thread Francisco Pastor
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

2007-04-25 Thread Hendrik Fuß
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

2007-04-25 Thread Peter Dalgaard
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

2007-04-25 Thread Duncan Murdoch
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

2007-04-25 Thread jim holtman
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

2007-04-25 Thread gyadav
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

2007-04-25 Thread Petr Klasterecky
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.

2007-04-25 Thread gyadav
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

2007-04-25 Thread Gabor Grothendieck
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

2007-04-25 Thread Petr Klasterecky
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

2007-04-25 Thread Gabor Grothendieck
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)

2007-04-25 Thread ONKELINX, Thierry
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

2007-04-25 Thread Inman, Brant A. M.D.
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?

2007-04-25 Thread Wade Wall
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()

2007-04-25 Thread Abi Ghanem josephine
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

2007-04-25 Thread Duncan Murdoch
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?

2007-04-25 Thread ONKELINX, Thierry
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

2007-04-25 Thread hadley wickham
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

2007-04-25 Thread Crispin Miller
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?

2007-04-25 Thread Shao
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

2007-04-25 Thread raymond chiruka
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

2007-04-25 Thread Duncan Murdoch
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

2007-04-25 Thread Liaw, Andy
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

2007-04-25 Thread Doran, Harold
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

2007-04-25 Thread David Forrest
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

2007-04-25 Thread Petr Klasterecky
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

2007-04-25 Thread Street N.R.
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

2007-04-25 Thread Colm G. Connolly
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

2007-04-25 Thread Soare Marcian-Alin
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

2007-04-25 Thread Waichler, Scott R
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()

2007-04-25 Thread Abi Ghanem josephine
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

2007-04-25 Thread Prof Brian Ripley
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?

2007-04-25 Thread Robert A LaBudde
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

2007-04-25 Thread Petr Klasterecky
?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

2007-04-25 Thread Crispin Miller
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

2007-04-25 Thread Peter Dalgaard
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()

2007-04-25 Thread chris
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

2007-04-25 Thread Steiner, Julien
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

2007-04-25 Thread gyadav
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

2007-04-25 Thread Egoldsystem
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

2007-04-25 Thread clayton . springer
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

2007-04-25 Thread Tony Plate
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

2007-04-25 Thread Torsten Hothorn
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?

2007-04-25 Thread Horace Tso
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]

2007-04-25 Thread Liaw, Andy
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

2007-04-25 Thread Earl F. Glynn
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

2007-04-25 Thread Martin Morgan
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

2007-04-25 Thread hadley wickham
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

2007-04-25 Thread Romain Francois
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

2007-04-25 Thread Frank E Harrell Jr
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

2007-04-25 Thread Rogerio Porto
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

2007-04-25 Thread Mike Huber
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

2007-04-25 Thread Leeds, Mark \(IED\)
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

2007-04-25 Thread Thibaut Jombart
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

2007-04-25 Thread Thomas Lumley
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

2007-04-25 Thread Deepayan Sarkar
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

2007-04-25 Thread Natalie O'Toole
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

2007-04-25 Thread Rajarshi Guha
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] Самый доступный и недорогой метод похудения

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

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

2007-04-25 Thread Daniel Elliott
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

2007-04-25 Thread J . delasHeras
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

2007-04-25 Thread Jason Parcon
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

2007-04-25 Thread Deepayan Sarkar
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

2007-04-25 Thread Chuck Cleland
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

2007-04-25 Thread Natalie O'Toole
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

2007-04-25 Thread Clint Bowman
?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

2007-04-25 Thread Hendrik Fuß
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

2007-04-25 Thread Peter Dalgaard
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?

2007-04-25 Thread Adelchi Azzalini
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

2007-04-25 Thread Hendrik Fuß
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

2007-04-25 Thread Duncan Murdoch
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

2007-04-25 Thread Hendrik Fuß
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

2007-04-25 Thread Waichler, Scott R
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

2007-04-25 Thread L . Stirton
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

2007-04-25 Thread Tony Plate
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

2007-04-25 Thread Spilak,Jacqueline [Edm]
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

2007-04-25 Thread Natalie O'Toole
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

2007-04-25 Thread Deepayan Sarkar
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

2007-04-25 Thread Frank E Harrell Jr
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

2007-04-25 Thread Dylan Beaudette
?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

2007-04-25 Thread Andrew Robinson
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

2007-04-25 Thread Daniel Bolnick
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

  1   2   >