[R] How to sort TWO columns ?

Dear ALL, I fear my question has already been answered many times before, but I haven't fund that in archives... I am working on spatial datasets and, in most arrays I'm handling, there are two columns dedicated to (x,y)-coords. For different reasons - notably to draw image() plots, I need to

Re: [R] LaTeX errors

Jim Lemon wrote: Hello again, Having successfully mastered the creation of fairly unsophisticated print methods, I have encountered the dreaded LaTeX error, i.e.: * checking concord-manual.tex ... ERROR LaTeX errors when creating DVI version. This typically indicates Rd problems. As I have no

Re: [R] How to sort TWO columns ?

Jacques VESLOT wrote: Dear ALL, I fear my question has already been answered many times before, but I haven't fund that in archives... I am working on spatial datasets and, in most arrays I'm handling, there are two columns dedicated to (x,y)-coords. For different reasons - notably to draw image()

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 http://www.ku.edu/~pauljohn/R/Rtips.html#2.12 help? kind regards,

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 sort it.

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

[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

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

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) ) snip Coefficients:

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-Pagan test

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 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] 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 heteroskedasticity. Is

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

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 of

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 that set of

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

[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

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

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,

[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!

[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

[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

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 Germany phone:

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

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 Joerg, Have

[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 everthing

[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,

[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

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:

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 it

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. I don't

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

Re: [R] Cutting heatmap dendrogram

Hi Sean, Martin Sean == Sean Davis sdavis2 at mail.nih.gov 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

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

[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

[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 your

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,1))

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()

[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 for

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=logreg2)

[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[s3,] if (nrow(a)==0) ? for (i in 1:nrow(a)){ a[i]-1 } a[i]-0 Lisa Wang Cancer

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

Lisa Wang lisawang at 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[s3,] : : if (nrow(a)==0) ? : : for (i in 1:nrow(a)){ :

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

[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] 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

[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

[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] 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.

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: x - 2

[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] 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 sdavis2 at mail.nih.gov 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

[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))