[R] Adding a We think R rocks page

Hi, there is a reference given for R. It should be used to prove its value to donators. OK, I quoted R but probably nobody will ever recognize that. A web page where dummies and no name users like me were pointed to and could leave a short statement of use and usefulness might help in

[R] R accuracy

Hello, I am trying to test the precision of R on datasets from The Statistical Reference Datasets Project http://www.itl.nist.gov/div898/strd/index.html and I don't manage to understand how R is storing its results. For example, I calculate a mean on the michelso dataset (100 values) and

Re: [R] R accuracy

Anthony Landrevie wrote: Hello, I am trying to test the precision of R on datasets from The Statistical Reference Datasets Project http://www.itl.nist.gov/div898/strd/index.html and I don't manage to understand how R is storing its results. For example, I calculate a mean on the michelso

[R] Question on statistics

Hi, Can anyone help me with the following (although not directly correlated to R functionality)? I have been looking on the internet but can not find the answer. My question: what is the variation on the mean of a limited distribution (total N points normally distributed), when I have a small

Re: [R] Tool for update

Hello, Am Mittwoch, 23. März 2005 06:12 schrieb Yuandan Zhang: Hi, Is there any tool to check if there is update version of a package available? I look for things alike YUM for linux? Start R --no-save on a root console and launch update.packages() from within R-Enviroment. Surely, this is

Re: [R] Problem encounter during graphics device driver

On Wed, 23 Mar 2005 [EMAIL PROTECTED] wrote: Hello I am facing the following problem using the R-version 1.9.1 The PDF or PS none of these device drivers are opening while I am using R-1.9.1, the following error message is coming Error in PS(file, old$paper, old$family, old$encoding, old$bg,

Re: [R] more classes to one class in one dataset

try this: split(dat, dat$class) where 'dat' is your data.frame I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web:

Re: [R] more classes to one class in one dataset

Thanks. It's work. Jan Sabee On Wed, 23 Mar 2005 11:30:22 +0100, Dimitris Rizopoulos [EMAIL PROTECTED] wrote: try this: split(dat, dat$class) where 'dat' is your data.frame I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School

[R] smallest/biggest number

Hi, I'm running monte carlo and i wonder what is the biggest/smallest number that can reliably be represented in R? Thanks, Chris __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] smallest/biggest number

On Wed, 23 Mar 2005, chris desimpelaere wrote: I'm running monte carlo and i wonder what is the biggest/smallest number that can reliably be represented in R? Well, -Inf and Inf, of course. But if you meant a finite number, see ?.Machine : the values are OS-specific. -- Brian D. Ripley,

[R] R on SuSe 9.2 AMD 64 bit: Solved

Hi everybody, I just downloaded the file R-2.0.1.tar.gz and followed the instruction written in doc/R-admin.html. In particular we installed for SuSe 9.2 AMD64 the following packages: gcc, gcc++, gcc-g77, Perl, te-latex, te-pdf, libpng, libbz, PCRE, Tcl/Tk, BLAS, LAPACK. After we just made:

[R] Browser to replace the internal browser?

Hello, I see that the more I work with R and the more the code gets larger I would like to have some graphic support in my quellcode. Is there a browser that could be easily implemened in R? And how do I call it from R? It would be nice if the browser replaces the fix() function. Carsten

Re: [R] Browser to replace the internal browser?

On Mar 23, 2005, at 6:44 AM, Carsten Steinhoff wrote: Hello, I see that the more I work with R and the more the code gets larger I would like to have some graphic support in my quellcode. Is there a browser that could be easily implemened in R? And how do I call it from R? It would be nice if

[R] manova and contrasts, again

Hi R-people, To determine contrasts after MANOVA I've found a piece of R-code provided by Yves Rosseel (http://tolstoy.newcastle.edu.au/R/help/04/06/0134.html), which has been very helpful. Now I Would like to determine contrasts for a model which has a main effect and an interaction effect,

RE: [R] Question on statistics

Ehh, by limited distribution, I meant to say a population of N points. ... -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Roy Werkman Sent: Wednesday, March 23, 2005 10:22 AM To: r-help@stat.math.ethz.ch Subject: [R] Question on statistics Hi, Can

RE: [R] Tool for update

If you're talking about R itself, I believe the answer is no. However, the release schedule for R is rather predictable (two major releases per year, one in Spring and another in Fall, with patch releases in between as needed), so the need is not that great, IMHO. Andy From: Yuandan Zhang

RE: [R] Question on statistics

If the sample is drawn with replacement from the finite population, then the usual formula applies (assuming iid samples); i.e., var(sample mean) = var(population) / n. There's some problem in your description: A finite population, I believe, is necessarily discrete (since there are only N

Re: [R] Package vignette and build

Yes, you need the package installed first. Something like: R CMD build --no-vignettes DLM R CMD install DLM. R CMD build DLM R CMD install DLM. At least you had to do this with 1.9.1, can't recall looking again since then. On Tue, 22 Mar 2005 11:23:36 -0600 (CST), Giovanni Petris

Re: [R] sampling from a mixture distribution

you have also to sample the mixture compoment membership; check this for a mixtrue of two normals: rnorm.mixture - function(n, prob=0.5, mu1=0, sigma1=1, mu2=0, sigma2=1){ u - runif(n) out - numeric(n) for(i in 1:n) out[i] - if(u[i] prob) rnorm(1, mu1, sigma1) else rnorm(1, mu2,

[R] replace values in a matrix subject to boolean condition

Hi everybody! I am sorry to bother you with a question so simple but I think there might be a better solution: I have a matrix of size 360x501 where I want to check the value of each 5th column of each row and replace it (and the 6th, 7th, 8th column) by zero if the value is less than 1000. I

Re: [R] sampling from a mixture distribution

For each variate, generate it from f1() with probability p1, and from f2() with probability p2. In other words, flip a p1-biased coin to decide which distribution, f1 or f2, to generate from. HTH, Giovanni Date: Wed, 23 Mar 2005 13:53:10 + From: Vumani Dlamini [EMAIL PROTECTED] Sender:

RE: [R] sampling from a mixture distribution

Here's one possible way: rmix2 - function(n, p1, rF1, rF2, argF1=NULL, argF2=NULL) { ## n is the number of deviates to simulate ## p1 is the probability of a point coming from the 1st component ## rF1, rF2 are functions for generating random deviates ## from the two components

[R] Negative binomial GLMMs in R

Dear R-users, A recent post (Feb 16) to R-help inquired about fitting a glmm with a negative binomial distribution. Professor Ripley responded that this was a difficult problem with the simpler Poisson model already being a difficult case:

Re: [R] sampling from a mixture distribution

I would like to sample from a mixture distribution p1*f(x1)+p2*f(x2). ***Surely*** you mean ``p1*f1(x)+p2*f2(x)'' !!! I usually sample variates from both distributions and weight them with their respective probabilities, but someone told me that was wrong. What is the correct way?

[R] Question on class 1, 2 output for RandomForest

Hi All, I read the R-newsletter Volum 2/3, December 2002 on page 18. I tried the example there, too. Then, I used a different data set with random Forest from the UCI respository. The results for the credit data generated 2 additional columns, column 1 and a column 2 that the example given in

Re: [R] Will R work on this 64 bit machine?...

I believe R will run out of the box on your setup. I personally haven't tried the RPMs but you can always build R from the sources (fairly straightforward on a Linux box). -roger dsandif wrote: Hello, Will R work on this 64 bit machine?, Here are the specs. of our linux box: *Red Hat

Re: [R] Will R work on this 64 bit machine?...

On Wed, 23 Mar 2005, dsandif wrote: Hello, Will R work on this 64 bit machine?, Yes. Here are the specs. of our linux box: *Red Hat Enterprise Linux WS (v.3 Standard for AMD64 and Intel EM64T) *OS: redhat-release Release: 3WS CPU Arch: ia32e-redhat-linux That's not at all clear: what is `ia32e'?

[R] nl regression with 8 parameters, help!

I'm doing a non linear regression with 8 parameters to be fitted: J.Tl.nls-nls(Gw~(a1/(1+exp(-a2*Tl+a3))+a4)*(b1/(1+exp(b2*Tl-b3))+b4),data=Enveloppe, start=list(a1=0.88957,a2=0.36298,a3=10.59241,a4=0.26308,

RE: [R] Question on class 1, 2 output for RandomForest

The `1' and `2' columns are the error rates within those classes. E.g., the last row of the `1' column should correspond to the class.error for -, and the last row of the `2' column to the class.error for +. (I would have thought that that should be fairly obvious, but I guess not. It

Re: [R] R accuracy

Try signif(m,8) At 1:03 AM -0800 3/23/05, Anthony Landrevie wrote: Hello, I am trying to test the precision of R on datasets from The Statistical Reference Datasets Project http://www.itl.nist.gov/div898/strd/index.html and I don't manage to understand how R is storing its results. For

Re: [R] alternative to 'groups' for lattice bwplot()

Mulholland, Tom [EMAIL PROTECTED] wrote: I'm afraid you have lost me. What is it that you want that reordering the formula does not achieve. bwplot(yield ~ year | site, data = barley) has sites next to each other. Yes, they are next to each other, but in different panels, as expected when

Re: [R] alternative to 'groups' for lattice bwplot()

On Wednesday 23 March 2005 00:10, Sebastian Luque wrote: Hi, Is there some alternative to the 'groups' argument in lattice's bwplot function for boxplots? Say in the example below: bwplot(yield ~ site | year, data = barley) you want to have two side by side boxplots per site, corresponding

Re: [R] nl regression with 8 parameters, help!

Does this error always occur independently of the starting values that you provide? I guess so, because I think that the parameters in your equation are not identifiable, since the first term (a1 to a4) is identical to the second term (b1 to b4) with a1 = b1, -a2 = b2, a3 = -b3, and a4 = b4 .

[R] How to do such MDS in R

i know cmdscale and isoMDS inR can do classical and non-metric MDS.but i want to konw if there is packages can carry on individual differences scaling and multidimensional analysis og preference?both method are important one,but i can not find any clue on how to do it using R. anyone can help?

[R] how to test for equality of covariance Matrices in lda

when using the two-group discriminant analysis,we need to test for equality of covariance Matrices in lda.as whenm we formed our estimate of the within-group covariance matrix by pooling across groups,we implicitly assumed that the covariance structure was the same across groups.so it seems

Re: [R] mixtures as outcome variables

Jason W. Martinez wrote: Dear R-users, I have an outcome variable and I'm unsure about how to treat it. Any advice? I have spending data for each county in the state of California (N=58). Each county has been allocated money to spend on any one of the following four categories: A, B, C, and D.

[R] nested random effects

Hi I am struggling with nested random effects and hope someone can help. I have individuals (ID) who are nested within families (FAM). I want to model an outcome variable, and take account of the intercorrelation of individuals within each family. I think this amounts to two random

Re: [R] Question on statistics

Roy Werkman wrote: Yes, it is discrete, but the underlying distribution is Gaussian. / I guess you mean what somebody calls the superpopulation distribution. Kjetil / Just got the following from a college: Var(mean of finite population) = ((N - n)/(N - 1)) * var(population) / n This should be

Re: [R] nl regression with 8 parameters, help!

Arne Henningsen [EMAIL PROTECTED] writes: Does this error always occur independently of the starting values that you provide? I guess so, because I think that the parameters in your equation are not identifiable, since the first term (a1 to a4) is identical to the second term (b1 to b4)

Re: [R] How to do such MDS in R

On 21 Mar 2005, at 13:29, ronggui wrote: i know cmdscale and isoMDS inR can do classical and non-metric MDS.but i want to konw if there is packages can carry on individual differences scaling and multidimensional analysis og preference?both method are important one,but i can not find any clue

[R] Does R work in 64 bit on apple G5?

Hi, I am working with R on 2xG5 1.8Ghz from Apple under 10.3.8 The G5 chip is 64 bits but does R run in 64 bit or 32 under OS X? How can know? I think it run in 32 bits... but not sure... anyway thanks for this fabulous soft... ;-) David __

Re: [R] nested random effects

On Wed, 2005-03-23 at 11:58 -0500, Shaw, Philip (NIH/NIMH) wrote: Hi I am struggling with nested random effects and hope someone can help. I have individuals (ID) who are nested within families (FAM). I want to model an outcome variable, and take account of the intercorrelation of

Re: [R] Does R work in 64 bit on apple G5?

On Wed, 23 Mar 2005, David Ruau wrote: Hi, I am working with R on 2xG5 1.8Ghz from Apple under 10.3.8 The G5 chip is 64 bits but does R run in 64 bit or 32 under OS X? How can know? I think it run in 32 bits... but not sure... Under the current OS X it runs 32bit. You can tell by looking at

Re: [R] Does R work in 64 bit on apple G5?

On Wed, 23 Mar 2005, David Ruau wrote: I am working with R on 2xG5 1.8Ghz from Apple under 10.3.8 The G5 chip is 64 bits but does R run in 64 bit or 32 under OS X? How can know? From the size of the ncells! 32-bit machine: gc() used (Mb) gc trigger (Mb) Ncells 144907 3.9 35 9.4

Re: [R] mixtures as outcome variables

Jason W. Martinez [EMAIL PROTECTED] 03/22/05 04:11PM Dear R-users, I have an outcome variable and I'm unsure about how to treat it. Any advice? Below are a couple of ideas/suggestions of things to think about I have spending data for each county in the state of California

RE: [R] nested random effects

An interaction random effect/fixed effect is noted as random ~1|random/fixed in your case random =~1|ID/FAM (but I don't uderstand why indiviuals withing families are fixed and and families are random, but there you go). 1. Fixed effects cannot be nested within random effects. 2.

RE: [R] nested random effects

It should be random=~1|FAM/ID indicating individuals are nested within families. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Federico Calboli Sent: Wednesday, March 23, 2005 12:34 PM To: Shaw, Philip (NIH/NIMH) Cc: r-help Subject: Re: [R] nested

Re: [R] Does R work in 64 bit on apple G5?

Thanks, I was sure the pre-compile version was 32 bit but not if you compile it your self... It give the same infos when you run gc() or .Machine$sizeof.pointer either on OS X client with a pre-compiled version or on OS X Server with a home compile version. .Machine$sizeof.pointer [1] 4 gc()

RE: [R] nested random effects

On Wed, 2005-03-23 at 10:04 -0800, Berton Gunter wrote: An interaction random effect/fixed effect is noted as random ~1|random/fixed in your case random =~1|ID/FAM (but I don't uderstand why indiviuals withing families are fixed and and families are random, but there you go).

RE: [R] nested random effects

I should have added that if you have only one Y observation per ID (within family), then the ID variance component is residual error and the model becomes (without any covariates) Y~1, rand=~1|FAM -- Bert On Wed, 2005-03-23 at 11:58 -0500, Shaw, Philip (NIH/NIMH) wrote: Hi I am

[R] go to msn!!!!

Some people think that this server is like msn messenger! I see that lots of people talk about some uninteresting things like G5 stuff or whatever, but nobody is able to think about real useRs' problems! Guillaume Storchi __ R-help@stat.math.ethz.ch

[R] Looking for function for Double to raw to double conversions

Hello, I am trying to implement function for reading/writing some XML file format. One feature of that XML format is that a lot of binary data is stored in Base64 format, and since R's XML package does not seem to support it, I just wrote my own converter from raw format to Base64, and back.

[R] R on red hat 2.1 problem while trying to generate image

Running R 1.9.1 under red hat 2.1 version When I try to generate an image, we get an error as in the following plot(rnorm(100)) Error in PS(file, old$paper, old$family, old$encoding, old$bg, old$fg, : unable to start device PostScript In addition: Warning message: cannot open

[R] Courses***April R/Splus Advanced and Intermediate level courses by XLSolutions

Here are our April courses: R/Splus Advanced Programming: March 31st - April 1st, San Francisco http://www.xlsolutions-corp.com/Radv.htm R/Splus Programming Techniques: April 14th - April 15th, New York City http://www.xlsolutions-corp.com/Rfund.htm Microarrays Data Analysis with R/S+ and

RE: [R] Gini's Importance Value Variable = Inf

That result looks fishy: Not only there shouldn't be Inf, but there shouldn't be negative values in that measure (look at V6). I will look into it. I hope by now you realize that there's not much point in asking such package-specific questions on R-help... Not all package maintainers are on

[R] Error in unitrootTest (fSeries)

Hello, I am getting the following error message from unitrootTest. Do you have any clue of what could be wrong. Details: AMD64 (x86_64) Gentoo Linux system. library(fSeries) kmodel - list(ar=c(.3,0,0,0,0.7,-.4*.7),d=1) x=armaSim(nobs,model=kmodel)

Re: [R] Error in unitrootTest (fSeries)

On Wed, 23 Mar 2005, Roberto Bertolusso wrote: Hello, I am getting the following error message from unitrootTest. Do you have any clue of what could be wrong. A bug in the package: please contact the maintainer. This *may* work if you run in R_HOME. Hint to Diethelm: use system.file(libs,

[R] Complete Linkage Clustering techniques

Dear R I recently asked for a cluster analysis Using * cluster.results - hclust(iris.dist, method=complete) * but nothing happened i.e the previous scatterplot matrix still showed up whereas I was expecting a dendogram. Could it be that because I had used cutree before on the scatter

[R] non-derivative based optimization and standard errors.

Hi AlL, I ahve this problem that my objective function is discontinous in the paramaters and I need to use methods such as nelder-mead to get around this. My question is: How do i compute standard errors to a problem that does not have a gradient? Any literature on this is greatly appreciated.

[R] Mapping actual to expected columns for princomp object

I am working with data sets in which the number and order of columns may vary, but each column is uniquely identified by its name. E.g., one data set might have columns MW logP Num_Rings Num_H_Donors while another has columns Num_Rings Num_Atoms Num_H_Donors logP MW I would like

[R] parallel r job on sun gridengine

greetings all, this may be the wrong forum for my problem - if so please advise. i am addressing this list because of an error i am getting from the snow library rmpi (i think) after lam has booted the mpi nodes i have a script (provided by a faculty member - i am not an R user but have the task

[R] summing values as image

I'm trying to summarize irregularly spaced data (in data.frame with x,y,z) and need to sum (not average as the as.image() function in fields does) and I'm not sure if there is a function in on of the packages or if I'm going to need to string a few functions together like fields::as.image() and

Re: [R] Complete Linkage Clustering techniques

- Original Message - From: Brett Stansfield [EMAIL PROTECTED] To: R help (E-mail) R-help@stat.math.ethz.ch Sent: Wednesday, March 23, 2005 5:25 PM Subject: [R] Complete Linkage Clustering techniques Dear R I recently asked for a cluster analysis Using * cluster.results -

[R] font sizes for row.names of dendograms

Dear R I recently performed a cluster analysis. It produced the dendogram no problem but unfortunately the font size of the row.names were all cluttered due to their large size So I tried to change the font size using plclust(cluster.results, labels=iris$specie, cex=0.8) and R came back to me

Re: [R] extracting numerical data from text field

Luis Tercero luis.tercero at ebi-wasser.uni-karlsruhe.de writes: : : I have imported a data frame that looks like this: : :Measurement.Date.and.Time Z.Average..nm. PDI : 572 Dienstag, 22. Mrz 2005 11:05:59 366,4 0,468 : 573 Dienstag, 22. Mrz 2005 11:09:30 353,4

Re: [R] extracting numerical data from text field

Gabor Grothendieck ggrothendieck at myway.com writes: Luis Tercero luis.tercero at ebi-wasser.uni-karlsruhe.de writes: : : I have imported a data frame that looks like this: : :Measurement.Date.and.Time Z.Average..nm. PDI : 572 Dienstag, 22. Mrz 2005 11:05:59

[R] Prediction using GAM

Recently I was using GAM and couldn't help noticing the following incoherence in prediction: data(gam.data) data(gam.newdata) gam.object - gam(y ~ s(x,6) + z, data=gam.data) predict(gam.object)[1] 1 0.8017407 predict(gam.object,data.frame(x=gam.data$x[1],z=gam.data$z[1]))

Re: [R] non-derivative based optimization and standard errors.

Have you considered bootstrap or Monte Carlo? spencer graves Jean Eid wrote: Hi AlL, I ahve this problem that my objective function is discontinous in the paramaters and I need to use methods such as nelder-mead to get around this. My question is: How do i compute standard errors to a

RE: [R] font sizes for row.names of dendograms

The error message states that you are passing a parameter called cex which has not been used. If you look at ?plclust more closely you will see it does not have cex parameter. However the S3 method for class hclust, plot, does? So does this help? hc - hclust(dist(USArrests), ave) plot(hc,cex =

[R] Robust multivariate regression with rlm

Dear Group, I am having trouble with using rlm on multivariate data sets. When I call rlm I get Error in lm.wfit(x, y, w, method = qr) : incompatible dimensions lm on the same data sets seem to work well (see code example). Am I doing something wrong? I have already browsed through

Re: [R] R on red hat 2.1 problem while trying to generate image

Ghosh, Sandeep wrote: Running R 1.9.1 under red hat 2.1 version Please upgrade. When I try to generate an image, we get an error as in the following plot(rnorm(100)) Error in PS(file, old$paper, old$family, old$encoding, old$bg, old$fg, : unable to start device PostScript In addition:

RE: [R] Robust multivariate regression with rlm

lm works for multivariate responses rlm does not - check what the help file says about the response. That's about it, really. Bill Venables. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Markku Mielityinen Sent: Thursday, 24 March 2005 5:20 PM To: