[R] Warning message when Loading package svMisc in R 2.10.0

2009-10-27 Thread Qi Li
Hi, I have a problem when using R 2.10.0, Loading required package: svMisc Warning message: package 'svMisc' was built under R version 2.9.1 and help will not work correctly HOW to fix it? __ R-help@r-project.org mailing list

Re: [R] explalinig the output of my linear model analysis

2009-10-27 Thread Peter Ehlers
Ista, Here's the quote from MASS (the book, 4e, p.142): Terms of the form a/x, where a is a factor, are best thought of as separate regression models of type 1 + x within the levels of a. I'm not answering the OP's question because in my view if one doesn't understand the output of lm(), then

Re: [R] Warning message when Loading package svMisc in R 2.10.0

2009-10-27 Thread Orvalho Augusto
I did not try. But you try to update all the packages with: update.packages() Good luck Caveman On Tue, Oct 27, 2009 at 6:28 AM, Qi Li ericq...@hotmail.com wrote: Hi, I have a problem when using R 2.10.0, Loading required package: svMisc Warning message: package 'svMisc' was built under R

[R] Error in solve.default peforming Competing risk regression

2009-10-27 Thread kende jan
Dear all, I am trying to use the crr function in the cmprsk package version 2.2 to analyse 198 observations.I have receive the error in solve.default. Can anyone give me some insights into where the problem is? Thanks here is my script : cov=cbind(x1,x2) z-crr(ftime,fstatus,cov)) and data

[R] tk_choose.files without multiple selection?

2009-10-27 Thread Eger, Marcus
Dear all, I compared tk_choose.files on linux with choose.files or tk_choose.files on windows and found that tk_choose.files does not allow multiple selection on linux - independent on keyword setting multi. tk_choose.files(default=paste(path.original,,sep=),caption=Select file(s)) On

Re: [R] Output pairwise.t.test to data.frame

2009-10-27 Thread Peter Ehlers
Sam, The relevant part of the output is in test[[p.value]] which is a matrix. Use write.table(test[[p.value]], file=output.csv, sep=,) and perhaps use the arguments row.names=FALSE and col.names=FALSE. -Peter Ehlers Sam Player wrote: # I'm doing a pairwise.t.test on a large dataset and

Re: [R] extract day or month as in Splus

2009-10-27 Thread Karl Ove Hufthammer
On Fri, 23 Oct 2009 12:06:28 -0600 Fang (Betty) Yang fang.y...@ualberta.ca wrote: dates - c(02/27/1992, 02/27/1992, 01/14/1992, 02/28/1992, 02/01/1992) Could anyone give me some R codes to get the same results as above' (extract days from dates), please? If your dates are real dates, and

[R] Installing RPostgreSQL from source on Mac OS X

2009-10-27 Thread christiaan pauw
Hi Everybody I use R2.9.2 on Mac OS X 10.5.2. I tried to install RPostgreSQL from source on CRAN via the package installer. Apparently the program cannot find my C compiler gcc . Here is the output message trying URL ' http://cran.za.r-project.org/src/contrib/RPostgreSQL_0.1-6.tar.gz' Content

[R] confidence interval for return levels

2009-10-27 Thread Madan Sigdel
Dear List   From a sample size of 40 I have calculated retrun level values for 2,5,10,25,50 and 100 years using a exponentail disrtibution as follows;  Now how can I calculate confidence interval for each return levels? seeking your help. Thank you  

[R] User input when running R code in batch mode

2009-10-27 Thread Kaushik Krishnan
Hi I've been stumbling over a simple issue that undoubtedly has an easy solution. I need to have some way for a user to enter some values into a data frame which R will then work on. I know that data entry should ideally be done otherwise and I should use R only for the computation, but R's

[R] boxplot using grid

2009-10-27 Thread Stefan JansevanRensburg
--- *** Disclaimer *** --- This e-mail and its contents are subject to the SA Reserve Bank's Disclaimer and Confidentiality Clause, which can be viewed at: http://www.reservebank.co.za/disclaimer Should you be unable to access the link provided, please send

Re: [R] basic statistics to csv

2009-10-27 Thread Adaikalavan Ramasamy
It would be useful to have a simplified version of the 'nsu' object. I am guessing it is a list of some sort (e.g. mean is single value, quantiles here returns 5 numbers) and not a matrix or dataframe (i.e. regular array). So you can have several choices here: 1) print nsu to a file. e.g.

Re: [R] Severe growing pains with new help system

2009-10-27 Thread Petr PIKAL
Hi I have just installed R 2.10.0 version (I used dev version before) and have similar problems with help files. starting plain version freshly installed R on WXP R version 2.10.0 (2009-10-26) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 ?plot only text

Re: [R] User input when running R code in batch mode

2009-10-27 Thread Barry Rowlingson
On Tue, Oct 27, 2009 at 9:20 AM, Kaushik Krishnan kaushik.s.krish...@gmail.com wrote: $ r --vanilla test.r a - scan(what='character',n=1); a 1: Read 0 items character(0) Now it's not working. Assuming this is a unix environment, the syntax ' test.r' means 'my standard input stream

Re: [R] Severe growing pains with new help system

2009-10-27 Thread Duncan Murdoch
On 27/10/2009 5:08 AM, Petr PIKAL wrote: Hi I have just installed R 2.10.0 version (I used dev version before) and have similar problems with help files. starting plain version freshly installed R on WXP R version 2.10.0 (2009-10-26) Copyright (C) 2009 The R Foundation for Statistical

Re: [R] API question (embedding R)

2009-10-27 Thread Duncan Murdoch
On 26/10/2009 6:56 PM, Bernd Kreuss wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hello list, i want to access R from a Pascal program. Basically i only need to be able to * start the engine / the interpreter (and keep it running) * pass strings to it (R code) that should be parsed and

Re: [R] Severe growing pains with new help system

2009-10-27 Thread Petr PIKAL
Thank you. Duncan Murdoch murd...@stats.uwo.ca napsal dne 27.10.2009 10:21:20: On 27/10/2009 5:08 AM, Petr PIKAL wrote: Hi I have just installed R 2.10.0 version (I used dev version before) and have similar problems with help files. starting plain version freshly installed R on WXP

Re: [R] Severe growing pains with new help system

2009-10-27 Thread Duncan Murdoch
On 27/10/2009 5:47 AM, Petr PIKAL wrote: Thank you. Duncan Murdoch murd...@stats.uwo.ca napsal dne 27.10.2009 10:21:20: On 27/10/2009 5:08 AM, Petr PIKAL wrote: Hi I have just installed R 2.10.0 version (I used dev version before) and have similar problems with help files. starting

Re: [R] Creating a sparse matrix from a file

2009-10-27 Thread Pallavi Palleti
Hi all, I used sparseM package for creating sparse Matrix and followed below commands. The sequence of commands are: ex - read.table('fileName',sep=',') M - as.matrix.csr(0,22638,80914) for (i in 1:nrow(ex)) { M[ex[i,1],ex[i,2]]-ex[i,3]} Even after 4 hours, I can still see the above command

Re: [R] reset par() within plot layout

2009-10-27 Thread Janke ten Holt
Peter Ehlers wrote: Janke ten Holt wrote: Tom Gottfried wrote: Janke, Janke ten Holt schrieb: Dear list, I would like to produce a matrix of plots, where par() is reset after each plot (see below [simplified] example). When I use layout() to do so, I seem to also reset the layout.

Re: [R] reset par() within plot layout

2009-10-27 Thread Janke ten Holt
Tom Gottfried wrote: I think it would be much easier to use grid-graphics for your task. It's made for such things. This is a good introduction: Paul Murrell. The grid graphics package. R News, 2(2):14-19, June 2002 Thank you for your suggestion. I will certainly check it out. Janke Tom

[R] ipredknn - How may I find values?

2009-10-27 Thread Grzes
Hi everybody! I want to find a closer neighbourins observation. This is my code: ## library(klaR) library(ipred) library(mlbench) data(PimaIndiansDiabetes2) dane=na.omit(PimaIndiansDiabetes2)[,c(2,5,9)] dane[,2]=log(dane[,2]) dane[,1:2]=scale(dane[,1:2])

[R] output (p-values) of fastbw in Design package

2009-10-27 Thread North, Bernard V
I am using the validate option in the Design package with the Cox survival model. I am using the bw=T option which, like the fastbw function, performs a backward elimination variable selection The output includes a series of columns (below) giving information on eliminated variables. My

Re: [R] Connecting to winbugs from R

2009-10-27 Thread dolberdinho
Ben Bolker wrote: dolberdinho wrote: Hi, I'm having trouble running a script within winbugs when calling winbugs from R. Here are a few lines from the winbugs log: display(log) check(C:\DOCUME~1\michael\LOCALS~1\Temp\RtmpjKnHMu\mm_TO_bugs9bc21e.txt) cannot open

Re: [R] Creating a sparse matrix from a file

2009-10-27 Thread Martin Maechler
PP Hi all, PP I used sparseM package for creating sparse Matrix and PP followed below commands. I'd strongly recommend to use package 'Matrix' which is part of every R distribution (since R 2.9.0). PP The sequence of commands are: ex - read.table('fileName',sep=',')

Re: [R] reset par() within plot layout

2009-10-27 Thread baptiste auguie
Hi, From ?par, Value When parameters are set, their former values are returned in an invisible named list. Therefore opar - par(col=red) will not contain col=red. HTH, baptiste 2009/10/27 Janke ten Holt j.c.ten.h...@rug.nl: This seems to work indeed. But I don't understand why... I would

Re: [R] reset par() within plot layout

2009-10-27 Thread Janke ten Holt
baptiste auguie wrote: Hi, From ?par, Value When parameters are set, their former values are returned in an invisible named list. Wow, I must have missed that. Thanks, Janke Therefore opar - par(col=red) will not contain col=red. HTH, baptiste 2009/10/27 Janke ten Holt

Re: [R] Warning message when Loading package svMisc in R 2.10.0

2009-10-27 Thread Richardson, Patrick
svMisc was built and submitted to CRAN under 2.9.2. If you are trying to use it within TINN, TINN will still work. The package needs to be built under 2.10 and resubmitted. I'm sure Phillipe will get to it in due time. Best regards, Patrick -Original Message- From:

[R] Stack overflow in R 2.10.0 with sub()

2009-10-27 Thread Kenneth Roy Cabrera Torres
Hi R developers: Congratulations for the new R 2.10.0 version. It is a huge effort! Thank you for your work and dedication. I just want to ask how to make this strip blank function to work again (it works on R.2.9.2). alumnos$AL_NUME_ID-sub((^ +)|( +$),,alumnos$AL_NUME_ID),) alumnos is a data

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread David Winsemius
On Oct 26, 2009, at 11:25 PM, ankush...@yahoo.com wrote: Hi, I am using the fitdistr of MASS to get the MLE for the lambda of a Poisson distribution. When i run the fitdistr command, i get an output that looks like - lambda 3.75 (0.03343) Couple of questions - 1. is the

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Kjetil Halvorsen
What is wrong with using mean(x) to get the MLE of the poisson lambda? Kjetil On Tue, Oct 27, 2009 at 9:17 AM, David Winsemius dwinsem...@comcast.net wrote: On Oct 26, 2009, at 11:25 PM, ankush...@yahoo.com wrote: Hi, I am using the fitdistr of MASS to get the MLE for the lambda of a

Re: [R] Creating a sparse matrix from a file

2009-10-27 Thread Pallavi P
Hi Martin, Thanks for the help. Just to make sure I understand correctly. The below steps are for creating an example table similar to the one that I read from file. n - 22638 m - 80914 nnz - 30 # no idea if this is realistic for you set.seed(101) ex - cbind(i = sample(n,nnz,

Re: [R] ipredknn - How may I find values?

2009-10-27 Thread David Winsemius
On Oct 27, 2009, at 6:02 AM, Grzes wrote: Hi everybody! I want to find a closer neighbourins observation. This is my code: ## library(klaR) library(ipred) library(mlbench) data(PimaIndiansDiabetes2) dane=na.omit(PimaIndiansDiabetes2)[,c(2,5,9)] dane[,2]=log(dane[,2])

Re: [R] Stack overflow in R 2.10.0 with sub()

2009-10-27 Thread Duncan Murdoch
On 10/27/2009 8:15 AM, Kenneth Roy Cabrera Torres wrote: Hi R developers: Congratulations for the new R 2.10.0 version. It is a huge effort! Thank you for your work and dedication. I just want to ask how to make this strip blank function to work again (it works on R.2.9.2).

Re: [R] Error in solve.default peforming Competing risk regression

2009-10-27 Thread Ravi Varadhan
Please send reproducible code when you report problems. I am unable to reproduce your problem. It seems to work just fine. dat - read.table(e:/temp/crrdata.txt, sep=\t, header=TRUE) names(dat) [1] x1 x2 fstatus ftime dim(dat) [1] 198 4 cov=cbind(dat$x1,dat$x2) z -

Re: [R] ipredknn - How may I find values?

2009-10-27 Thread Max Kuhn
I think we are having some difficulty understanding what you are looking for. If you are looking to find which of the training samples were closest to the prediction sample, I don't think that you can get it from this function. If this is what you want, I use the dist function in the proxy

Re: [R] How to reduce key strokes when defining S4 classes?

2009-10-27 Thread Peng Yu
On Mon, Oct 26, 2009 at 11:56 PM, mtmor...@fhcrc.org wrote: Quoting Peng Yu pengyu...@gmail.com: I feel tedious when I define a S4 class and the methods. For each method, I have to call both setMethod and setGeneric (or both setReplaceMethod and setGeneric). I would like a more compact

[R] if there is a package can do blast for gene in R?

2009-10-27 Thread echo_july
hello if there is a package can do blast in R? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] output (p-values) of fastbw in Design package

2009-10-27 Thread Frank E Harrell Jr
North, Bernard V wrote: I am using the validate option in the Design package with the Cox survival model. I am using the bw=T option which, like the fastbw function, performs a backward elimination variable selection The output includes a series of columns (below) giving information on

[R] Why 'return' is needed in R?

2009-10-27 Thread Peng Yu
It seems that 'return' is not necessary when returning a value. If this is the case, I don't understand why 'return' is a keyword in R. Is there a case in which I have to use 'return'? f-function(x) { + x + } g-function(x) { + return(x) + } print(f(2)) [1] 2 print(g(2)) [1] 2

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Peter Ehlers
Kjetil Halvorsen wrote: What is wrong with using mean(x) to get the MLE of the poisson lambda? and mean(x)/length(x) to get its estimated variance. -Peter Ehlers Kjetil On Tue, Oct 27, 2009 at 9:17 AM, David Winsemius dwinsem...@comcast.net wrote: On Oct 26, 2009, at 11:25 PM,

Re: [R] Why 'return' is needed in R?

2009-10-27 Thread Sarah Goslee
Did you read the help? If you use return(), that value is returned regardless of where in the function it appears. This can be used in conditionals, for example. If the function never hits a return value, the value of the last evaluated expression is returned. testfunction - function(x,

Re: [R] Why 'return' is needed in R?

2009-10-27 Thread Duncan Murdoch
On 10/27/2009 10:12 AM, Peng Yu wrote: It seems that 'return' is not necessary when returning a value. If this is the case, I don't understand why 'return' is a keyword in R. Is there a case in which I have to use 'return'? If you want to return early from a function you need it, e.g. f -

Re: [R] if there is a package can do blast for gene in R?

2009-10-27 Thread James W. MacDonald
It depends on what you mean by 'do blast'. If you are looking to align sequences to a reference genome with possible gaps, then you should look at (to start) the Biostrings package of Bioconductor. There are other alignment tools, depending on what exactly you want to do. If you are simply

Re: [R] ipredknn - How may I find values?

2009-10-27 Thread Grzes
Yes, I want to know which points in my picture are in red or green area. For example: .glucose..insulin.diabetes 609 0.95177272 1.139969011 - I want to know that it's for example: black point in red area 253 -1.05724970 -1.158814331 - it's for

Re: [R] Stack overflow in R 2.10.0 with sub()

2009-10-27 Thread Kenneth Roy Cabrera Torres
Dr. Murdoch: I am puzzled! As you adviced me I do this: x - as.character(alumnos$AL_NUME_ID) x - x[-seq_len(length(x)/2)] y - gsub((^ +)|( +$),,x) And it fails, But, trying to locate the problem I do: x - as.character(alumnos$AL_NUME_ID) x - x[-seq_len(length(x)/2)] x -

[R] column names of a correlation matrix

2009-10-27 Thread Lee William
Hi! All, I am working on a correlation matrix of 4217x4217 named 'cor_expN'. I wish to obtain pairs with highest correlation values. So, I did this b=matrix(data=NA,nrow=4217,ncol=1) rownames(b)=rownames(cor_expN) for(i in 1:4217){b[i,]=max(cor_expN[i,])} head(b) [,1]

Re: [R] Cannot activate chm help in R 2.10

2009-10-27 Thread Ulrike Groemping
Duncan Murdoch-2 wrote: Ulrike Groemping wrote: ... My system is Windows XP, German locale computer. On installation, I was asked to decide for text or html help and chose html (there was no radio button for chm help). ... If you set options(help_type=html) you'll get HTML

Re: [R] Stack overflow in R 2.10.0 with sub()

2009-10-27 Thread Duncan Murdoch
On 10/27/2009 10:46 AM, Kenneth Roy Cabrera Torres wrote: Dr. Murdoch: I am puzzled! As you adviced me I do this: x - as.character(alumnos$AL_NUME_ID) x - x[-seq_len(length(x)/2)] Please try the following. After doing the lines above, do save(x, file=x.RData) and exit from R. Then

[R] Basic Sweave

2009-10-27 Thread alexander . weidstam
Hi, I have this basic R code that I would like to rewrite into Sweave code. Ideally I would like Latex to display both the code and the statistics. R and Sweave is completely new to me so I would really appreciate some help. read.table(df1.txt, header=TRUE, sep=;, na.strings=c(na,.), skip=1,

Re: [R] Cannot activate chm help in R 2.10

2009-10-27 Thread Duncan Murdoch
On 10/27/2009 10:56 AM, Ulrike Groemping wrote: Duncan Murdoch-2 wrote: Ulrike Groemping wrote: ... My system is Windows XP, German locale computer. On installation, I was asked to decide for text or html help and chose html (there was no radio button for chm help). ... If you set

Re: [R] Creating a sparse matrix from a file

2009-10-27 Thread Martin Maechler
PP == Pallavi P pallavip...@gmail.com on Tue, 27 Oct 2009 18:13:22 +0530 writes: PP Hi Martin, PP Thanks for the help. Just to make sure I understand correctly. PP The below steps are for creating an example table similar to the one that I PP read from file. yes, exactly

[R] wilcox.exact() problem

2009-10-27 Thread David Croll
Dear R friends, here I write again about the wilcox.exact() problem. I want to compare two sets of categorical data, and in one case it says negative length vectors not allowed, and in the other one I get the error cannot allocate vector of length On

Re: [R] Basic Sweave

2009-10-27 Thread James W. MacDonald
http://www.stat.uni-muenchen.de/~leisch/Sweave/Sweave-manual.pdf alexander.weids...@gmail.com wrote: Hi, I have this basic R code that I would like to rewrite into Sweave code. Ideally I would like Latex to display both the code and the statistics. R and Sweave is completely new to me so I would

Re: [R] wilcox.exact() problem

2009-10-27 Thread Achim Zeileis
On Tue, 27 Oct 2009, David Croll wrote: Dear R friends, here I write again about the wilcox.exact() problem. I want to compare two sets of categorical data, and in one case it says negative length vectors not allowed, and in the other one I get the error cannot allocate vector of length

Re: [R] User input when running R code in batch mode

2009-10-27 Thread Bernd Kreuss
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Kaushik Krishnan wrote: Is there any way to make R stop for the user to enter values when running in batch mode either by changing the way I invoke scan() or readLines() or by using any other function? At least on linux this works:

Re: [R] Cannot activate chm help in R 2.10

2009-10-27 Thread Uwe Ligges
Duncan Murdoch wrote: On 10/27/2009 10:56 AM, Ulrike Groemping wrote: Duncan Murdoch-2 wrote: Ulrike Groemping wrote: ... My system is Windows XP, German locale computer. On installation, I was asked to decide for text or html help and chose html (there was no radio button for chm

Re: [R] Code improvement

2009-10-27 Thread Sebastien Bihorel
Hi Baptisti, Sorry for the late reply. I wanted to thank you for putting me on the right track. I finally got something to work not really by extracting the legend from a xyplot call, but by building my own legend grid.frame and passing it to my low-level function. Understanding Grobs and

[R] Detection Times and Poisson Distribution

2009-10-27 Thread Lorenzo Isella
Dear All, Apologies if my questions are too basic for this list. I am given a set of data corresponding to list of detection times (real, non-integer numbers in general) for some events, let us say nuclear decays to fix the ideas. It is a small dataset, corresponding to about 400 nuclear decay

[R] What is the difference between anova {stats} and aov?

2009-10-27 Thread Peng Yu
There are anova {stats} and aov in R. It seems that anova takes an object returned by a model fitting function. But I don't see any examples for anova. Can somebody give me a simple example on anova? What is the difference between anova and aov? __

[R] How to express a tree?

2009-10-27 Thread Hesen Peng
Hello everybody, I'm write a function whose output represents a tree. Can anyone please recommend me some data structure in R which are good for expressing a tree? The possible idea I'm having in my mind is to represent the whole tree as a list in which a sublist (an element to the main list)

[R] lasso plot using LARS

2009-10-27 Thread Russell Ivory
When plotting a lars object, I cannot find a way to plot solid lines. Even when the arguments breaks=F and lty=solid are used, the vertical lines at the break points do not plot but asterisks indicating the breaks still plot as part of each path leaving solid lines broken up by asterisks at the

[R] New vector based on if/else statement within for loop?

2009-10-27 Thread Douglas M. Hultstrand
Hello, I am trying to create a new vector (w) that is based on comparing two vectors (P and Z). The compaison is simple (I created a for loop that reassigns w based on if statement), all Z values = 24 and P values =1, w=88 else w=77. I am not getting the correct results for w, see

Re: [R] Stack overflow in R 2.10.0 with sub()

2009-10-27 Thread Duncan Murdoch
On 10/27/2009 1:05 PM, Kenneth Roy Cabrera Torres wrote: Thank you very much for your interest. I make this: x - as.character(alumnos$AL_NUME_ID) x - x[-seq_len(length(x)/2)] save(x, file=x.RData) I exit form R, and then restart R and I make this: load(x.RData) y - gsub((^ +)|( +$),,x) It

Re: [R] What is the difference between anova {stats} and aov?

2009-10-27 Thread Chuck Cleland
On 10/27/2009 1:06 PM, Peng Yu wrote: There are anova {stats} and aov in R. It seems that anova takes an object returned by a model fitting function. But I don't see any examples for anova. Can somebody give me a simple example on anova? What is the difference between anova and aov?

Re: [R] New vector based on if/else statement within for loop?

2009-10-27 Thread Duncan Murdoch
On 10/27/2009 1:14 PM, Douglas M. Hultstrand wrote: Hello, I am trying to create a new vector (w) that is based on comparing two vectors (P and Z). The compaison is simple (I created a for loop that reassigns w based on if statement), all Z values = 24 and P values =1, w=88 else w=77. I

Re: [R] Stack overflow in R 2.10.0 with sub()

2009-10-27 Thread Kenneth Roy Cabrera Torres
Thank you very much for your interest. I make this: x - as.character(alumnos$AL_NUME_ID) x - x[-seq_len(length(x)/2)] save(x, file=x.RData) I exit form R, and then restart R and I make this: load(x.RData) y - gsub((^ +)|( +$),,x) It shows me: Error en gsub((^ +)|( +$), , x) : input string

[R] sm.regression

2009-10-27 Thread Val
Hi all, I was looking for a non parametric survival analysis and I came up with the following sample from the web. However, I could not run it. Which library or function does sm.regression require? x - runif(100,-2, 2) y - x^2 + rnorm(50) sm.regression(x, y, h=0.5) Error: could not find

Re: [R] RPostgreSQL package and libpq.dll file

2009-10-27 Thread christiaan pauw
Update to R 2.9 As far as I know the RPostgreSQL package doesn't work in R 2.8 2009/9/14 Lore M tchiba...@hotmail.com It still doesn't work. Now, R is asking me for the SSLEAY32.dll. If I download it from the internet, and then put it in the file \R-2.8.1\library\RPostgreSQL\libs, now they

Re: [R] operations on sparse matrices, and dense intermediary steps

2009-10-27 Thread Martin Maechler
Hi Jose, JQ == Jose Quesada ques...@gmail.com on Sat, 24 Oct 2009 23:49:08 +0200 writes: JQ -BEGIN PGP SIGNED MESSAGE- JQ Hash: SHA1 JQ Hi, JQ I'm doing some basic operations on large sparse matrices, for example JQ getting a row. JQ it takes close to 30

[R] automatically adjusting axis limits

2009-10-27 Thread Servet Ahmet Cizmeli
Dear R users, I am a newbie. Just switched from MATLAB. So thanks a lot for your patience. I have 5 spectra collected in field. Each spectra has two columns : Wavelength (56) and the actual measurement. Each measurement came in a different .txt file on disk (5 files in total). I wrote

Re: [R] Stack overflow in R 2.10.0 with sub()

2009-10-27 Thread Kenneth Roy Cabrera Torres
El mar, 27-10-2009 a las 10:47 -0700, Phil Spector escribió: What happens if you type Sys.setlocale('LC_ALL','C') before using gsub or grep? When I do that, R hangs and don't show any message. - Phil Spector

Re: [R] User input when running R code in batch mode

2009-10-27 Thread Martin Maechler
Barry Rowlingson b.rowling...@lancaster.ac.uk on Tue, 27 Oct 2009 10:17:24 +0100 writes: On Tue, Oct 27, 2009 at 9:20 AM, Kaushik Krishnan kaushik.s.krish...@gmail.com wrote: $ r --vanilla test.r a - scan(what='character',n=1); a 1: Read 0 items character(0)

Re: [R] lasso plot using LARS

2009-10-27 Thread Kingsford Jones
Within the plot.lars code, change the type argument to matplot from 'b' to 'l': library(lars) myplot.lars - edit(plot.lars) #change type = 'b' to type = 'l' in the call to matplot data(diabetes) object - with(diabetes, lars(x,y)) myplot.lars(object, lty = 1, breaks = FALSE) hth, Kingsford

[R] Exclude rows in xyplot

2009-10-27 Thread Joel Fürstenberg-Hägg
Hi all, I'm searching for a way to exclude outliers from my dataset while making xyplots. While plotting using pairs(), I exclude specific row in my data frame and save the settings as a variable which I later include as an argument: # Discard outliers and save settings as idx

Re: [R] explalinig the output of my linear model analysis

2009-10-27 Thread Rolf Turner
On 27/10/2009, at 7:58 PM, Peter Ehlers wrote: ... if one doesn't understand the output of lm(), then one's knowledge of statistics is insufficient to warrant using lm(). I nominate this as a fortune. cheers, Rolf Turner

[R] Print several xyplots to the same page in a pdf file

2009-10-27 Thread Joel Fürstenberg-Hägg
Hello everybody, I'm using the lattice package and the xyplot to make several graphs like below. However, I can just print the three grouped plots onto one page as I'm putting them into a pdf-file, which gives me a huge amount of pages... Is it possible to put them all, or at least more than

Re: [R] Sandard deviation calculation

2009-10-27 Thread marcos carvajalino
Dear ms. The most possible reason for this: Dear all I am trying to learn R I was trying to calculate standard deviation here are the commands and outputs in R sd(Ht_cm[type=='SD'])[1] 3.283605 sd(Ht_cm[from_treeline=='above'])[1] 16.83289 sd(Ht_cm[type=='SD'][from_treeline=='above'])[1] NA

Re: [R] automatically adjusting axis limits

2009-10-27 Thread Ben Bolker
Sac-6 wrote: Dear R users, I am a newbie. Just switched from MATLAB. So thanks a lot for your patience. I have 5 spectra collected in field. Each spectra has two columns : Wavelength (56) and the actual measurement. Each measurement came in a different .txt file on disk

Re: [R] automatically adjusting axis limits

2009-10-27 Thread Ben Bolker
Sac-6 wrote: Dear R users, I am a newbie. Just switched from MATLAB. So thanks a lot for your patience. I have 5 spectra collected in field. Each spectra has two columns : Wavelength (56) and the actual measurement. Each measurement came in a different .txt file on disk

Re: [R] automatically adjusting axis limits

2009-10-27 Thread John Kane
PLEASE provide commented, minimal, self-contained, reproducible code. particularly a sample data.set. At the moment it sounds like you have one variable (Wavelength) with a length of 56 and Reflectance with a length of 5,000. What format(s) are the data in? Data.frames? Clearly this is not

Re: [R] Sandard deviation calculation

2009-10-27 Thread marcos carvajalino
Hi ms. also try this: sd(Ht_cm[type=='SD' from_treeline=='above']) cause i think you can't use [][] two apply 2 conditions to a data frame, in that cases you should use an AND () conector like above 2009/10/27 marcos carvajalino maancafe...@gmail.com: Dear ms. The most possible reason for

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Ankush Bhargava
Thank you everyone for responding. David, 3.75 in my example was equivalent to the mean of the values, which i thought was too much a coincidence... What do you think the significance of (0.03343). What is this value? Kjetil, Are you saying that mean(x) is same as the MLE for the poisson

Re: [R] Exclude rows in xyplot

2009-10-27 Thread Dieter Menne
Joel Fürstenberg-Hägg wrote: Now I'm trying to make xyplots to compare the result from three different categories: # Plot Pro against Glc for each of the three categories xyplot(Pro ~ Glc | Categories_BBCH_ID, data=fieldTrial0809, pch=°, layout=c(1, 3), aspect=1,

Re: [R] Detection Times and Poisson Distribution

2009-10-27 Thread Ben Bolker
Lorenzo Isella wrote: Dear All, Apologies if my questions are too basic for this list. I am given a set of data corresponding to list of detection times (real, non-integer numbers in general) for some events, let us say nuclear decays to fix the ideas. It is a small dataset,

Re: [R] Sandard deviation calculation

2009-10-27 Thread Duncan Murdoch
On 27/10/2009 2:58 PM, marcos carvajalino wrote: Dear ms. The most possible reason for this: Dear all I am trying to learn R I was trying to calculate standard deviation here are the commands and outputs in R sd(Ht_cm[type=='SD'])[1] 3.283605 sd(Ht_cm[from_treeline=='above'])[1] 16.83289

Re: [R] explalinig the output of my linear model analysis

2009-10-27 Thread Kingsford Jones
...which would complement the apropos fortune from John Fox: library(fortunes) fortune('dangerous') If you give people a linear model function you give them something dangerous. -- John Fox useR! 2004, Vienna (May 2004) On Tue, Oct 27, 2009 at 12:48 PM, Rolf Turner

Re: [R] Sandard deviation calculation

2009-10-27 Thread marcos carvajalino
Yup you're right, after posting i took a better look at the subsetting, found the problem and post again, my bad... 2009/10/27 Duncan Murdoch murd...@stats.uwo.ca: On 27/10/2009 2:58 PM, marcos carvajalino wrote: Dear ms. The most possible reason for this: Dear all I am trying to learn R

[R] Non-normal residuals.

2009-10-27 Thread rkevinburton
Hello, I asked a question about what the most likely process to follow if after a time-series fit is performed the residuals are found to be non-normal. One peron responded and offered to help if I supplied a sample data set. Unfortunately now that I have a sample I have lost the emai

Re: [R] sm.regression

2009-10-27 Thread Peter Ehlers
Val wrote: Hi all, I was looking for a non parametric survival analysis and I came up with the following sample from the web. However, I could not run it. Which library or function does sm.regression require? Time to install package 'sos'. It would quickly tell you that you'll find

[R] option to control the spac between columns in data frame

2009-10-27 Thread Subodh Acharya
Hello, I have a question regarding a way to control the appreance of output exported by R when I use capture.output( x, file = Directory/file.txt) , I get a text file which when I paste to a word file looks like the first table below. The following table has its clumns spaced closely so when I

[R] How to quit unwanted execution immediately?

2009-10-27 Thread Peng Yu
Occasionally, I start a command (taking long time to finish) that I did not really want to start. I type 'ctrl+C' to try to quit the execution. However, R does not quit the execution of the command immediately. I'm wondering if R could response to ctrl+C immediately.

[R] newbie: grading an histogram for several classes of the same discipline

2009-10-27 Thread Mauricio Ortiz Calvao
Hi there Sorry for the simple question, but I have been unable to find an answer. I am the coordinator of a basic discipline in a university; we have several classes taking this discipline and I would like to compare the grades, via a histogram for the distinct classes. Concretely, I have

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Peter Ehlers
Ankush Bhargava wrote: Thank you everyone for responding. David, 3.75 in my example was equivalent to the mean of the values, which i thought was too much a coincidence... What do you think the significance of (0.03343). What is this value? Kjetil, Are you saying that mean(x) is same as the

Re: [R] Cannot activate chm help in R 2.10

2009-10-27 Thread Ulrike Groemping
Uwe Ligges-3 wrote: Duncan Murdoch wrote: On 10/27/2009 10:56 AM, Ulrike Groemping wrote: Duncan Murdoch-2 wrote: Ulrike Groemping wrote: ... My system is Windows XP, German locale computer. On installation, I was asked to decide for text or html help and chose html (there was

[R] anova interaction contrasts: crossing helmert and linear contrasts

2009-10-27 Thread Katie Surrence
I am new to statistics, R, and this list, so apologies in advance for the errors etiquette I am certain to make (in spite of reading the posting guide, help on various commands, etc.).  Any help is greatly appreciated. Here is my data: fatigue =

Re: [R] Warning message when Loading package svMisc in R 2.10.0

2009-10-27 Thread Uwe Ligges
Richardson, Patrick wrote: svMisc was built and submitted to CRAN under 2.9.2. If you are trying to use it within TINN, TINN will still work. The package needs to be built under 2.10 and resubmitted. I'm sure Phillipe will get to it in due time. No, that's nonsense. It does *not* need to

Re: [R] Non-normal residuals.

2009-10-27 Thread Ben Bolker
rkevinburton wrote: Hello, I asked a question about what the most likely process to follow if after a time-series fit is performed the residuals are found to be non-normal. One peron responded and offered to help if I supplied a sample data set. Unfortunately now that I have a sample I

Re: [R] Cbind() on the right-side of a formula in xYplot()

2009-10-27 Thread Dylan Beaudette
On Monday 26 October 2009, Frank E Harrell Jr wrote: Dylan Beaudette wrote: Hi, Using the latest rms package I am able to make nice plots of model predictions +/- desired confidence intervals like this: # need this library(rms) # setup data d - data.frame(x=rnorm(100),

Re: [R] tk_choose.files without multiple selection?

2009-10-27 Thread Prof Brian Ripley
This depends on the version of Tk you have installed. The Tk 8.5.7 ChangeLog has 2008-08-25 Todd M. Helfter t...@users.sourceforge.net The patch is associated with the bug tracker id: 1936220 library/tkfbox.tcl : fix the multiple selection error for tk_getOpenFile

  1   2   >