Re: [R] dumpClass, hasSlot in R?

Marc Schwartz MedAnalytics.com> writes: : : On Wed, 2004-07-21 at 15:53, hadley wickham wrote: : > There are a few notes about difference between the R implementation : > and the book at http://developer.r-project.org/methodsPackage.html : > : > I found the hardest thing to get to grips in R wa

Re: [R] questions about principle component analysis (princomp)

Yunfeng Hu wrote: Hi, I am a new R user and am currently using princomp to conduct a PCA. I have read the help(princomp) and still do not quite understand everything in the help. Basically I want to get the covariance matrix, and eigenvector/eigenvalues (loadings()?) so that I can find the prin

[R] Programmation pour MLE

Bonjour, Je veux cherché l’estimateur de vraisemblance maximal (MLE)d’une fonction à 3 paramètre inconue étant donné une échantillon de taille 50 (les observations des valeurs de x) alors comment je peux procédé La fonction de densité est définie par : f(x)= 1/3(g(a1)+g(a2)+g(a3))

Re: [R] Cutting heatmap dendrogram

> "paul" == paul boutros <[EMAIL PROTECTED]> > on Wed, 21 Jul 2004 12:03:39 -0400 writes: paul> Hi Sean, Martin >>> "Sean" == Sean Davis >>> on Wed, 21 Jul 2004 06:01:33 -0400 writes: Sean> Paul, You can certainly get a heatmap of a subset of Sean> your da

[R] Regression Modeling Strategies Short Course

I will be giving a one-day short course related to my book Regression Modeling Strategies in Toronto as part of the Joint Statistical Meetings on August 8. For more information visit the American Statistical Association web site amstat.org and biostat.mc.vanderbilt.edu/rms. The course applies

RE: [R] add more data to a data file, problem with .Rprofile

> From: Francisco J Molina > > Hi, > > Is there a way to incorporate an R object x to a file already > containing R objects? > It seems that 'save' is not capable to do this. If I save x > to a file containing > previously saved data, then I will lose this data. See if the following helps: >

[R] questions about principle component analysis (princomp)

Hi, I am a new R user and am currently using princomp to conduct a PCA. I have read the help(princomp) and still do not quite understand everything in the help. Basically I want to get the covariance matrix, and eigenvector/eigenvalues (loadings()?) so that I can find the principle components.

[R] CCA code

Hello, Is there a fortran or C program to do the canonical Coefficient Analysis (CCA)? Could someone give it to me? Thanks for any help in advance. Qiaozhen __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE

[R] add more data to a data file, problem with .Rprofile

Hi, Is there a way to incorporate an R object x to a file already containing R objects? It seems that 'save' is not capable to do this. If I save x to a file containing previously saved data, then I will lose this data. I use the funcion hsv in .Rprofile but it is recognized. Couriously, if I so

Re: [R] How to do a "go to " in a loop in R

Well, there ain't no such thing in R And your code doesn't really make sense anyway. You talk about ``a[i]'' where a is apparently a ***matrix*** (with the same number of columns as the matrix b, but with fewer --- possibly 0 --- rows). And in such as setting a[i] has a meaning, but probabl

[R] Building problem: leftover from old OS

Hello, I have updated my OS from Solaris 2.7 to Solaris 2.8, and I am trying to build R 1.9.1. Configure seems to run fine, resulting in the message: R is now configured for sparc-sun-solaris2.8 Source directory: . Installation directory:/usr/local C compiler:

Re: [R] dumpClass, hasSlot in R?

On Wed, 2004-07-21 at 15:53, hadley wickham wrote: > There are a few notes about difference between the R implementation > and the book at http://developer.r-project.org/methodsPackage.html > > I found the hardest thing to get to grips in R was method calling - > using multiple dispatch (totally d

Re: [R] dumpClass, hasSlot in R?

There are a few notes about difference between the R implementation and the book at http://developer.r-project.org/methodsPackage.html I found the hardest thing to get to grips in R was method calling - using multiple dispatch (totally different to what I'm used to from Java, Python etc.). I foun

Re: [R] How to do a "go to " in a loop in R

Lisa Wang uhnres.utoronto.ca> writes: : : Hi there, : : I'm writing a function which involves a loop. What to write in the "?" : place would allow it skips the "for loop" and goes to "a[i]<-0". : : a<-b[s>3,] : : if (nrow(a)==0) ? : : for (i in 1:nrow(a)){ :

[R] How to do a "go to " in a loop in R

Hi there, I'm writing a function which involves a loop. What to write in the "?" place would allow it skips the "for loop" and goes to "a[i]<-0". a<-b[s>3,] if (nrow(a)==0) ? for (i in 1:nrow(a)){ a[i]<-1 } a[i]<-0 Lisa Wang

RE: [R] function ms

There's an optimizer in S-PLUS called ms(), but not R. Are you sure you're using the version of MICE for R? Andy > From: Marie-Pierre Sylvestre > > Dear R users, > > I am using the MICE package. Specifically, at some point in > my code I have > > imp2=mice(PoptotalMICE,imputationMethod="logr

[R] function ms

Dear R users, I am using the MICE package. Specifically, at some point in my code I have imp2=mice(PoptotalMICE,imputationMethod="logreg2") And R returns... iter imp variable 1 1 MICEYError in logitreg(xobs, yobs, intercept=F) : couldn't find function "ms" I have been looking

Re: [R] Problem using xfig()

Christophe Grova wrote: Hello ... I tried to generate .fig figures with R, using the xfig() function ...When I open the > figures using xfig software under linux ... the foreground color change strangely ... > whereas when I display the same figure in R or when I saved it in using postscript()

Re: [R] 2 images with 2 color scales on one graphic device. How to get it?

On Wednesday 21 July 2004 09:10 am, Wolski wrote: > Hi! > > Would like to plot two image plots with color scale (like levelplot > (lattice) implements) in one graphic device. Tried to plot 2 levelplot's > (lattice) (the data has not much in common) in one graphics device by > setting par(mfrwo=c(2,

[R] R/S-plus Course***In Houston, TX***R/Splus Fundamentals and Programming Techniques, August 19-20, 2004

XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce June 2004 2-day "R/S-plus Fundamentals and Programming Techniques". Houston, TX --> August, 19-20 Interested in our R/Splus Advanced Programming course? Please email us! Reserve yo

[R] Problem using xfig()

Hello ... I tried to generate .fig figures with R, using the xfig() function ...When I open the figures using xfig software under linux ... the foreground color change strangely ... whereas when I display the same figure in R or when I saved it in using postscript() function there is no problem

Re: [R] RE: Comparison of correlation coefficients - Details

That sounds very close to a meta-analytic comparison of two statistics. As a matter of fact, the Rosenthal & Rubin approach transforms all primary statistics into Pearson r and then to Fisher's z and then follows with comparisons. More, comparisons can take into account sample sizes, or the value

Re: [R] Cutting heatmap dendrogram

Hi Sean, Martin >> "Sean" == Sean Davis >> on Wed, 21 Jul 2004 06:01:33 -0400 writes: Sean>> Paul, You can certainly get a heatmap of a subset of Sean>> your data by simply subsetting. If you have a group Sean>> of genes obtained from cutree, simply do a heatmap on S

Re: [R] 2 images with 2 color scales on one graphic device. How to get it?

Wolski wrote: Hi! Would like to plot two image plots with color scale (like levelplot (lattice) implements) in one graphic device. Tried to plot 2 levelplot's (lattice) (the data has not much in common) in one graphics device by setting par(mfrwo=c(2,1)) which does not work. Is there a way to for

RE: [R] nonparametetric bivariate regression

Ximing: Hope you don't mind me cc'ing R-help. > From: ximing wu > > I am looking for something that can fit > f(x,y)=g(z_1, z_2, z_3,...). > > It seems that the packages you mentioned work for: > z=f(x,y). > > Correct me if I am wrong. Thanks anyway for your help. You're right. My mistake.

Re: [R] Can R work on very large of data?

On Wed, 21 Jul 2004 09:05:18 -0500, Hayashi Soichi - shayas <[EMAIL PROTECTED]> wrote : >Is there anyway I can tweak R around so that it will start processing as >data comes and not load everthing on memory at once? See the ?connections help topic. You can open a file, and then read lines from i

Re: [R] Can R work on very large of data?

Hi On Wed, 21 Jul 2004, Hayashi Soichi - shayas wrote: > Is there anyway I can tweak R around so that it will start processing as > data comes and not load everthing on memory at once? The reason for this is Have you read through R Data Import/Export? There are several other ways, to name a few:

[R] 2 images with 2 color scales on one graphic device. How to get it?

Hi! Would like to plot two image plots with color scale (like levelplot (lattice) implements) in one graphic device. Tried to plot 2 levelplot's (lattice) (the data has not much in common) in one graphics device by setting par(mfrwo=c(2,1)) which does not work. Is there a way to force levelplot

[R] dumpClass, hasSlot in R?

Venables and Ripley (2000) R Programming (Springer, sec. 5.1) and Chambers (1998) Programming with Data (Springer) describe functions such as dumpClass and hasSlot that I can not find in R 1.9.1 for Windows nor in the "R site search". What do people use for similar functionality? Also, w

[R] Can R work on very large of data?

Hello.. I am pretty new for R, but if I understand it correctly, when R read data by doing something like "d <- read.table("/dev/stdin")", it read entire data first then start processing.. Is there anyway I can tweak R around so that it will start processing as data comes and not load evert

Re: [R] MySQL mismatch

On Wed, 2004-07-21 at 13:52, Joerg Schaber wrote: > Hi, > > library(RMySQL) works, but connecting to the database results in > > Warning message: > RS-DBI driver warning: (MySQL mismatch between compiled version 3.23.49 > and runtime version 3.23.51) > > > Any idea? > > Thanks, > > joerg Jo

RE: [R] Rose Diagrams

Lars, Look at diag.rose in CircStats and circular packages. I haven't used these, but I may need to in the future. David Germain Las Vegas -Original Message- From: Lars Peters [mailto:[EMAIL PROTECTED] Sent: Wednesday, July 21, 2004 5:49 AM To: [EMAIL PROTECTED] Subject: [R] Rose Diag

Re: [R] Rose Diagrams

Hi, library(CircStats) ?rose.diag Hope this helps, Matthias > Hi, > > Is it possible to create Rose Diagrams of wind data (speed & > direction) with R?? > > Best regards, > > Lars Peters > > > - > Lars Peters > > University of Konstanz > Limnological Institute > D-78457 Konstanz > Ge

[R] MySQL mismatch

Hi, library(RMySQL) works, but connecting to the database results in Warning message: RS-DBI driver warning: (MySQL mismatch between compiled version 3.23.49 and runtime version 3.23.51) Any idea? Thanks, joerg __ [EMAIL PROTECTED] mailing list https://

[R] Rose Diagrams

Hi, Is it possible to create Rose Diagrams of wind data (speed & direction) with R?? Best regards, Lars Peters - Lars Peters University of Konstanz Limnological Institute D-78457 Konstanz Germany phone: +49 (0)7531 88-2930 fax: +49 (0)7531 88-3533 e-mail: [EMAIL PROTECTED] web: Lars Pe

[R] PCA of vectors

Hello, Is it possible to perform PCA using princomp or prcomp on two-component wind data - ie wind vectors? If not, is anyone able to advise on a suitable methodology for acheiving this? My data is a time series of observations taken at different locations. Thank you in advance for any hints! Lau

Re: [R] Precision in R

Bjørn-Helge Mevik wrote: > Since you didn't say anything about _what_ you did, either in SAS or > R, my first thought was: Have you checked that you use the same > parametrization of the models in R and SAS? Good point! cheers,

Re: [R] Precision in R

Since you didn't say anything about _what_ you did, either in SAS or R, my first thought was: Have you checked that you use the same parametrization of the models in R and SAS? -- Bjørn-Helge Mevik __ [EMAIL PROTECTED] mailing list https://www.stat.ma

[R] RE: Comparison of correlation coefficients - Details

Dear all I apologize for cross-posting, but first it is accepted custom to thank the repliers and give a summary, and second I have still the feeling that this problem might be a general statistical problem and not necessarily related to microarrays only, but I might be wrong. First, I want to

Re: [R] regression slope

On Tue, 2004-07-20 at 17:02, Avril Coghlan wrote: Hello, I'm a newcomer to R so please forgive me if this is a silly question. It's that I have a linear regression: fm <- lm (x ~ y) and I want to test whether the slope of the regression is significantly less than 1. How can I do this in R? Another

Re: [R] Cutting heatmap dendrogram

> "Sean" == Sean Davis <[EMAIL PROTECTED]> > on Wed, 21 Jul 2004 06:01:33 -0400 writes: Sean> Paul, You can certainly get a heatmap of a subset of Sean> your data by simply subsetting. If you have a group Sean> of genes obtained from cutree, simply do a heatmap on Sean

RE: [R] nonparametetric bivariate regression

There may very well be others, but these packages ought to get you started: locfit, sm, and if I'm not mistaken, mgcv (which is shipped with R). Andy > From: ximing wu > > Hi there, > Does R has built-in codes for nonpara. bivariate regression > so that I can > estimate the joint distribution

Re: [R] Cutting heatmap dendrogram

Paul, You can certainly get a heatmap of a subset of your data by simply subsetting. If you have a group of genes obtained from cutree, simply do a heatmap on that set of genes. If you obtain a set of genes, say A, and want to do a heatmap on that subset, simply do heatmap(as.matrix(data[A,])

Re: [R] Testing autocorrelation & heteroskedasticity of residuals in ts

On Wed, 21 Jul 2004 10:28:41 +0200 Pfaff, Bernhard wrote: > > > > Hi, > > > > I'm dealing with time series. I usually use stl() to > > estimate trend, stagionality and residuals. I test for > > normality of residuals using shapiro.test(), but I > > can't test for autocorrelation and heteroskedas

RE: [R] Testing autocorrelation & heteroskedasticity of residuals in ts

Hi Vito, I would treat the residuals as a time series. You can fit an arima model to the residuals. For example: data(nottem) temp<-stl(nottem, "per") my.arima<-arima(temp$time.series[,3], order=c(1,0,0), include.mean = FALSE) my.arima tsdiag(my.arima) If the ar1 term is significant then the

Re: [R] error "evaluation nested too deeply" {was "Heatmap"}

Hi, Trying to plot a large dendrogram, I get the recursion limit error Error in match.fun(FUN) : evaluation nested too deeply: infinite recursion / options(expression=)? but setting options(expressions = 1) or something large enough to solve the infinite recursion limit

RE: [R] Testing autocorrelation & heteroskedasticity of residuals in ts

> > Hi, > > I'm dealing with time series. I usually use stl() to > estimate trend, stagionality and residuals. I test for > normality of residuals using shapiro.test(), but I > can't test for autocorrelation and heteroskedasticity. > Is there a way to perform Durbin-Watson test and > Breusch-Paga

RE: [R] regression slope

see also the contributed document by John Verzani, Simple R, page 87f. > Adaikalavan Ramasamy wrote: > > > I would try to construct the confidence intervals and > compare them to > > the value that you want > > > >>x <- rnorm(20) > >>y <- 2*x + rnorm(20) > >>summary( m1 <- lm(y~x) ) > > > > >

Re: [R] How to sort TWO columns ?

Jacques VESLOT wrote: Sorry sorry, I didn't pay enough attention to order() help pages indeed... and it costs me much useless indexing. As concerns how to put legend in image() plot, is there a piece of information about somewhere ? thanks a lot for helping, Jacques VESLOT In the same way as in usu

[R] Testing autocorrelation & heteroskedasticity of residuals in ts

Hi, I'm dealing with time series. I usually use stl() to estimate trend, stagionality and residuals. I test for normality of residuals using shapiro.test(), but I can't test for autocorrelation and heteroskedasticity. Is there a way to perform Durbin-Watson test and Breusch-Pagan test (or other si

RE: [R] How to sort TWO columns ?

Sorry sorry, I didn't pay enough attention to order() help pages indeed... and it costs me much useless indexing. As concerns how to put legend in image() plot, is there a piece of information about somewhere ? thanks a lot for helping, Jacques VESLOT -Message d'origine- De : Uwe L

Re: [R] How to sort TWO columns ?

Hi, It would be discussed yesterday. See e.g. in Google for ("sorting a data frame" R) http://www.r-project.org/nocvs/mail/r-help/2002/0088.html Or one entry yesterday: Best, Matthias >Hi >Lets assign your data frame to the variable yourdf >then: >yourdf[ , order( yourdf$year ) ] >should so

Re: [R] How to sort TWO columns ?

Hoi Jacques, --On woensdag 21 juli 2004 10:35 +0400 Jacques VESLOT <[EMAIL PROTECTED]> wrote: Could somedy please let me know how to sort two columns of a dataframe, with priority to one of them, just like in access ? Would this help? kind regards,