[R] Robust Regression + Stepwise

2010-08-15 Thread Shubha Vishwanath Karanth
Hi R, Can I perform a stepwise method for robust regression? In other words, how do I combine the methods "step" and "rlm" together? Thanks and Regards, Shubha This e-mail may contain confidential and/or privileged i...{{dropped:13}} __ R

[R] Help for using nnet in R for NN training and testing

2010-08-15 Thread kgorahava
Hello, I want to use nnet package in R, to train and simulate a NN and get the value of MSE. I am reading in a file which has 19 input variables and one output variable and has a total of 2000 observations. The first column in the file is a column just for giving the serial numbers of the observa

[R] Using AMORE package for NN

2010-08-15 Thread kgorahava
Hi, I am running the AMORE package functions in R to compare the MSE given by R , with the results I am getting by running my Professors algorithm in Matlab. The code I am using in R is given below. I have one file which has 125 observations to be used for training the NN and one file with 12

Re: [R] How to setup a Rprofile file?

2010-08-15 Thread yehengxin
I think I got it! -- View this message in context: http://r.789695.n4.nabble.com/How-to-setup-a-Rprofile-file-tp2326328p2326339.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

[R] How to setup a Rprofile file?

2010-08-15 Thread yehengxin
http://www.youtube.com/watch?v=zml4uXQDL4U Could you take a look at the end of the video. He made a batch file to run Rcdmr. I wonder how to setup Rprofile file to realize it. Thanks!! -- View this message in context: http://r.789695.n4.nabble.com/How-to-setup-a-Rprofile-file-tp2326328p2326

Re: [R] Limited output

2010-08-15 Thread Seth Falcon
Hi, On 7/21/10 1:58 AM, confusedcius wrote: > > The details of my problem are as follows: > I have an sql that returns 2192 rows in sqlite. > In R, I typed the following: >> library("RSQLite") >> con <- dbConnect(dbDriver("SQLite"), dbname = "C:\\sqlite\\... .sqlite") >> dbListTables(con) > #[1]

[R] How can I read date format '02-Jan-02' ?

2010-08-15 Thread Dong-Hee Koh
Hello, I am trying to read a database exported from SAS. It is form of csv, and the date format reads '02-Feb-99'. I used following code to convert character to date format, db$dob<-as.Date(db$dob, format="%d-%b-%y"). but it doesn't work, only seems NA. What's wrong with this code? Thank you in

Re: [R] re-order the rows of a data frame accroding to a specified order of some column

2010-08-15 Thread Leon Yee
> > Yes, > > x[order(order(z)),] > > Two uses of order are needed, as shown. > Thank you very much, Rich. It's like a magic! Leon [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] R 64-bit Windows isn't using much memory

2010-08-15 Thread Joshua Wiley
On Sun, Aug 15, 2010 at 8:46 PM, wrote: > >  I have a 9 GB RAM Windows Vista machine. > >  I installed the 64-bit version of R 2.11.1 for Windows from here: > > http://cran.r-project.org/bin/windows64/base/ > > I am running a program now in R. However, looking at Windows Task Manager, I > see th

Re: [R] re-order the rows of a data frame accroding to a specified order of some column

2010-08-15 Thread RICHARD M. HEIBERGER
Yes, x[order(order(z)),] Two uses of order are needed, as shown. Rich On Sun, Aug 15, 2010 at 11:47 PM, Leon Yee wrote: > Dear all, > > I have a data frame with several columns, and I have the specified > order of a given column. How can I re-order the rows of my data frame > accroding

[R] re-order the rows of a data frame accroding to a specified order of some column

2010-08-15 Thread Leon Yee
Dear all, I have a data frame with several columns, and I have the specified order of a given column. How can I re-order the rows of my data frame accroding to this specified order? for example, x = data.frame( num = 1:26, alpha = letters[1:26], stringsAsFactors=FALSE) z = sample(x$alpha, re

[R] R 64-bit Windows isn't using much memory

2010-08-15 Thread ericksonsm
I have a 9 GB RAM Windows Vista machine. I installed the 64-bit version of R 2.11.1 for Windows from here: http://cran.r-project.org/bin/windows64/base/ I am running a program now in R. However, looking at Windows Task Manager, I see that Rgui.exe is only using 12% of CPU and 191,900K of mem

Re: [R] ASCI characters

2010-08-15 Thread David Winsemius
To improve your further efforts at searching you should take note of the fact that it is ASCII, not ASCI. On Aug 15, 2010, at 8:50 PM, Orvalho Augusto wrote: Thank you!! Caveman On Mon, Aug 16, 2010 at 2:41 AM, Matt Shotwell wrote: How about: rawToChar(as.raw(82)) [1] "R" -Matt

Re: [R] ASCI characters

2010-08-15 Thread Orvalho Augusto
Thank you!! Caveman On Mon, Aug 16, 2010 at 2:41 AM, Matt Shotwell wrote: > How about: > > > rawToChar(as.raw(82)) > [1] "R" > > -Matt > > On Sun, 2010-08-15 at 19:50 -0400, Orvalho Augusto wrote: > > Hello guys! > > > > Is there any function that permits me to get an ASCI character from its >

Re: [R] ASCI characters

2010-08-15 Thread Matt Shotwell
How about: > rawToChar(as.raw(82)) [1] "R" -Matt On Sun, 2010-08-15 at 19:50 -0400, Orvalho Augusto wrote: > Hello guys! > > Is there any function that permits me to get an ASCI character from its > code? Eg. ascifunction(34) would give me ' > or ascifunction(92) gives \ > > Thanks > Caveman >

Re: [R] About .RData

2010-08-15 Thread Stephen Liu
- Original Message From: Duncan Murdoch To: Stephen Liu Cc: r-help@r-project.org Sent: Mon, August 16, 2010 4:20:08 AM Subject: Re: [R] About .RData - snip - >> 2) Why on Windows R the path of the folders/directories shows the >>Administrator; >> [1] "C:/Users/Administrator/Documents

[R] ASCI characters

2010-08-15 Thread Orvalho Augusto
Hello guys! Is there any function that permits me to get an ASCI character from its code? Eg. ascifunction(34) would give me ' or ascifunction(92) gives \ Thanks Caveman [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] CRAN (and crantastic) updates this week

2010-08-15 Thread Crantastic
CRAN (and crantastic) updates this week New packages * ber (1.0) Unknown http://crantastic.org/packages/ber In this package there are functions for adjusting the data to remove batch effects * charlson (1.0) Vanessa Cox http://crantastic.org/packages/charlson Input p

Re: [R] Paired t-tests

2010-08-15 Thread Marc Schwartz
On Aug 15, 2010, at 2:48 PM, David Winsemius wrote: > > On Aug 15, 2010, at 3:31 PM, Peter Dalgaard wrote: > >> Marc Schwartz wrote: >>> On Aug 15, 2010, at 9:05 AM, R Help wrote: >>> Hello List, I'm trying to do a paired t-test, and I'm wondering if it's consistent with eq

Re: [R] how to display delta AIC of all models

2010-08-15 Thread Ben Bolker
elaine kuo gmail.com> writes: > I am making model selection of generalised linear models based on delta > AIC, > using the command step(AIC). > > Also, I would like to learn the explanatory powers of each independent > variable. > To phrase differently, it is in need to show the delta AIC of all

Re: [R] greatest common divisor

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 5:57 PM, qingchen wrote: Hello, Is there a function in R with which we could find the greatest common divisor directly? Learn to search: > rhelpSearch <- function(string, restrict = c("Rhelp10", "Rhelp08", "Rhelp02", "functions" ),

[R] how to display delta AIC of all models

2010-08-15 Thread elaine kuo
Dear all, I am making model selection of generalised linear models based on delta AIC, using the command step(AIC). Also, I would like to learn the explanatory powers of each independent variable. To phrase differently, it is in need to show the delta AIC of all models rather than the final model

[R] greatest common divisor

2010-08-15 Thread qingchen
Hello, Is there a function in R with which we could find the greatest common divisor directly? Or is there a funciton, which could help us to find the remainder of a division directly? Best wishes! Thank you! Your truly, Qing -- View this message in context: http://r.789695.n4.nabble.com

[R] Sorry, I have changed the Chinese R help mailing list name to avoid similarity problem.

2010-08-15 Thread Tengfei Yin
Sorry for my mistake, I set Chinese help mailing list as r-h...@r-china.org, the first time I sent my test mail, I use auto-complete in gmail to finish my input, so I sent chinese character to English help mailing list, anyway, that warned me the potential problem in the future. I have change the

Re: [R] fGarch: how to use garchFit() in loop?

2010-08-15 Thread Marius Hofert
Dear David, thanks for the quick response. Yes, I tried "formula", but this gives the error "Multivariate data inputs require lhs for the formula." library(fGarch) spec <- garchSpec(model = list(alpha = 0.1, beta = c(0.4, 0.4))) data <- garchSim(spec, n = 100) x <- list() for(q in 1:3){

Re: [R] 测试3

2010-08-15 Thread Tengfei Yin
Sorry! I apologize for my mistake! 2010/8/15 Tengfei Yin > ²âÊÔ3 > > -- > Tengfei Yin > MCDB PhD student > 1620 Howe Hall, 2274, > Iowa State University > Ames, IA,50011-2274 > Homepage: www.tengfei.name > > > -- Tengfei Yin MCDB PhD student 1620 Howe Hall, 2274, Iowa State University

Re: [R] Fwd: as.logical(factor) behaviour

2010-08-15 Thread Peter Dalgaard
David Winsemius wrote: > >> > > What about changing it to behave thusly: > > > as.logical.factor <- function(vec) > as.logical(as.numeric(factor(vec))-1 ) > > > as.logical(as.numeric(factor(c("TRUE", "FALSE", NA)))-1 ) > [1] TRUE FALSENA > > as.logical(as.numeric(factor(c(TRUE, FALSE

[R] 测试4测试4测试4测试4

2010-08-15 Thread Tengfei Yin
²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4²âÊÔ4 -- Tengfei Yin MCDB PhD student 1620 Howe Hall, 2274, Iowa State University Ames, IA,50011-2274 Homepage: www.tengfei.name [[alternative HTML version deleted]] __ R-help@r-project.org mai

[R] 测试3

2010-08-15 Thread Tengfei Yin
²âÊÔ3 -- Tengfei Yin MCDB PhD student 1620 Howe Hall, 2274, Iowa State University Ames, IA,50011-2274 Homepage: www.tengfei.name [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-

[R] 测试

2010-08-15 Thread Tengfei Yin
²âÊÔ -- Tengfei Yin MCDB PhD student 1620 Howe Hall, 2274, Iowa State University Ames, IA,50011-2274 Homepage: www.tengfei.name [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-h

Re: [R] fGarch: how to use garchFit() in loop?

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 4:26 PM, Marius Hofert wrote: Dear expeRts, How can I specify the order p,q of a GARCH(p,q) model within a loop? Here's a minimal example showing that an "Formula and data units do not match"-error appears: library(fGarch) spec <- garchSpec(model = list(alpha = 0.1, b

[R] fGarch: how to use garchFit() in loop?

2010-08-15 Thread Marius Hofert
Dear expeRts, How can I specify the order p,q of a GARCH(p,q) model within a loop? Here's a minimal example showing that an "Formula and data units do not match"-error appears: library(fGarch) spec <- garchSpec(model = list(alpha = 0.1, beta = c(0.4, 0.4))) data <- garchSim(spec, n = 100) x <-

Re: [R] About .RData

2010-08-15 Thread Duncan Murdoch
On 15/08/2010 10:24 AM, Stephen Liu wrote: - Original Message From: Michael Bedward To: Stephen Liu Cc: r-help@r-project.org Sent: Sun, August 15, 2010 3:43:55 PM Subject: Re: [R] About .RData Furthermore can I run following command on R? rm /home/userA/.RData .RData contains all

Re: [R] Trouble loading "saved" Rdata

2010-08-15 Thread steven mosher
During my session I write several .Rdata objects to a variety of subdirectories so replicating the exact problem wasnt very easy. In the actual program all the files get written. all the files have sizes that fit the amount of data in them. It looks like the problem was naming the files .RDat

Re: [R] Trouble loading "saved" Rdata

2010-08-15 Thread steven mosher
I think it came down to my actual program having a function that saved the objects it was passed with a .RData extension as opposed to .Rdata Rechecking the whole thing. On Sun, Aug 15, 2010 at 11:05 AM, Joshua Wiley wrote: > Steven, > > I have exited my R session and restarted and I can load th

[R] Fwd: as.logical(factor) behaviour

2010-08-15 Thread David Winsemius
Begin forwarded message: From: Peter Dalgaard Date: August 15, 2010 3:43:31 PM EDT To: phgrosj...@sciviews.org Cc: r-help@r-project.org Subject: Re: [R] as.logical(factor) behaviour Philippe Grosjean wrote: Thank you, but I already know that. I am not surprised by this behavior, but by an

Re: [R] reading a text file, one line at a time

2010-08-15 Thread jim holtman
Read in as big a chunk as you can; take a look at your memory usage and make sure you environment does not have any unnecessary large objects sitting around. On Sun, Aug 15, 2010 at 10:12 AM, Data Analytics Corp. wrote: > Hi, > > This seems like a good solution.  I was concerned about the time t

Re: [R] Paired t-tests

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 3:31 PM, Peter Dalgaard wrote: Marc Schwartz wrote: On Aug 15, 2010, at 9:05 AM, R Help wrote: Hello List, I'm trying to do a paired t-test, and I'm wondering if it's consistent with equations. I have a dataset that has a response and two treatments (here's an exampl

Re: [R] Adding colored background area to a time series plot

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 12:21 PM, James Toll wrote: Hi, I am trying to add a rectangular colored background area to a plot of a time series of relative price changes. I believe that what I'm trying to do is very similar to the question and example given here: http://www.mail-archive.com/r-h

Re: [R] as.logical(factor) behaviour

2010-08-15 Thread Peter Dalgaard
Philippe Grosjean wrote: > Thank you, but I already know that. I am not surprised by this behavior, > but by an inconsistency between that behavior and the documentation that > says "For factors, this uses the levels (labels).", which it does not. > Best, My gut feeling say that the docs need to

Re: [R] reading a text file, one line at a time

2010-08-15 Thread Data Analytics Corp.
Hi, This seems like a good solution. I was concerned about the time taken up reading one at a time. If a chuck can be read in each time, then that should be the way for me to handle the problem. Thanks, Walt Walter R. Paczkowski, Ph.D. Data Analytics Corp. 44 Hami

Re: [R] Paired t-tests

2010-08-15 Thread Peter Dalgaard
Marc Schwartz wrote: > On Aug 15, 2010, at 9:05 AM, R Help wrote: > >> Hello List, >> >> I'm trying to do a paired t-test, and I'm wondering if it's consistent >> with equations. I have a dataset that has a response and two >> treatments (here's an example): >> >> ID trt order resp >>

Re: [R] Adding colored background area to a time series plot

2010-08-15 Thread James Toll
On Aug 15, 2010, at 10:40 AM, Gabor Grothendieck wrote: > On Sun, Aug 15, 2010 at 12:21 PM, James Toll wrote: >> Hi, >> >> I am trying to add a rectangular colored background area to a plot of a time >> series of relative price changes. I believe that what I'm trying to do is >> very similar t

Re: [R] Adding colored background area to a time series plot

2010-08-15 Thread Gabor Grothendieck
On Sun, Aug 15, 2010 at 2:04 PM, James Toll wrote: > On Aug 15, 2010, at 10:40 AM, Gabor Grothendieck wrote: > >> On Sun, Aug 15, 2010 at 12:21 PM, James Toll wrote: >>> Hi, >>> >>> I am trying to add a rectangular colored background area to a plot of a >>> time series of relative price changes.

Re: [R] time of serialization

2010-08-15 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Saptarshi Guha > Sent: Sunday, August 15, 2010 9:23 AM > To: R-help@r-project.org > Subject: [R] time of serialization > > Hello, > I have question about the overhead in lapply.

Re: [R] time of serialization

2010-08-15 Thread Saptarshi Guha
Hello, I take this back, wrong implementation of the experiment. Sorry for the noise. Regards Saptarshi On Sun, Aug 15, 2010 at 12:22 PM, Saptarshi Guha wrote: > Hello, > I have question about the overhead in lapply. > x is a list of 3000 lists. Each of the i (1<=i<=3000) list elements is > pai

Re: [R] Trouble loading "saved" Rdata

2010-08-15 Thread Joshua Wiley
Steven, I have exited my R session and restarted and I can load the file without issue. I have also tried loading the saved data on some older versions of R (2.10.1 and 2.11.0) and Windows (XP). Have you tried recreating the test object, ensuring that it is not NULL itself, resaving it, and then

Re: [R] Paired t-tests

2010-08-15 Thread Marc Schwartz
On Aug 15, 2010, at 9:05 AM, R Help wrote: > Hello List, > > I'm trying to do a paired t-test, and I'm wondering if it's consistent > with equations. I have a dataset that has a response and two > treatments (here's an example): > > ID trt order resp > 17 1 0 1 0.0037513592 >

[R] unexpected behaviour with sparse.model.matrix

2010-08-15 Thread Jarrod Hadfield
Hi, I'm trying to get sparse.model.matrix to retain unused levels. I can't seem to get this working through the most obvious routes such as specifying drop.unused.levels = FALSE in the model.frame or trying to pass all levels in xlev,which is an argument to sparse.model.matrix (see code b

Re: [R] reading a text file, one line at a time

2010-08-15 Thread jim holtman
For efficiency of processing, look at reading in several hundred/thousand lines at a time. One line read/write will probably spend most of the time in the system calls to do the I/O and will take a long time. So do something like this: con <- file('yourInputFile', 'r') outfile <- file('yourOutpu

Re: [R] reading a text file, one line at a time

2010-08-15 Thread Dimitri Liakhovitski
In read.table - under parameters can't you specify the number of lines to be read in? You could create a loop and read them one by one. Dimitri On Sun, Aug 15, 2010 at 10:58 AM, Data Analytics Corp. wrote: > Hi, > > I have an upcoming project that will involve a large text file.  I want to > >  1

Re: [R] reading a text file, one line at a time

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 10:58 AM, Data Analytics Corp. wrote: Hi, I have an upcoming project that will involve a large text file. I want to 1. read the file into R one line at a time ?readLines 2. do some string manipulations on the line 3. write the line to another text file. I can h

Re: [R] reading a text file, one line at a time

2010-08-15 Thread Matt Shotwell
Walt, Something like: con <- file("your-large-file.txt", "rt") readLines(con, 1) # Read one line -Matt On Sun, 2010-08-15 at 10:58 -0400, Data Analytics Corp. wrote: > Hi, > > I have an upcoming project that will involve a large text file. I want to > >1. read the file into R one line

Re: [R] Adding colored background area to a time series plot

2010-08-15 Thread Gabor Grothendieck
On Sun, Aug 15, 2010 at 12:21 PM, James Toll wrote: > Hi, > > I am trying to add a rectangular colored background area to a plot of a time > series of relative price changes.  I believe that what I'm trying to do is > very similar to the question and example given here: > > http://www.mail-archi

[R] reading a text file, one line at a time

2010-08-15 Thread Data Analytics Corp.
Hi, I have an upcoming project that will involve a large text file. I want to 1. read the file into R one line at a time 2. do some string manipulations on the line 3. write the line to another text file. I can handle the last two parts. Scan and read.table seem to read the whole file

[R] Adding colored background area to a time series plot

2010-08-15 Thread James Toll
Hi, I am trying to add a rectangular colored background area to a plot of a time series of relative price changes. I believe that what I'm trying to do is very similar to the question and example given here: http://www.mail-archive.com/r-h...@stat.math.ethz.ch/msg73948.html http://www.mayin.or

[R] time of serialization

2010-08-15 Thread Saptarshi Guha
Hello, I have question about the overhead in lapply. x is a list of 3000 lists. Each of the i (1<=i<=3000) list elements is pair of two elements: a string vector and a data frame x is roughly 235MB. > gc() ## > z <- system.time(y <- lapply(x,function(r){ system.time(serialize(r,NULL))['elapsed

Re: [R] as.logical(factor) behaviour

2010-08-15 Thread Philippe Grosjean
Thank you, but I already know that. I am not surprised by this behavior, but by an inconsistency between that behavior and the documentation that says "For factors, this uses the levels (labels).", which it does not. Best, Philippe On 15/08/10 16:09, R Help wrote: The problem is that, underne

Re: [R] Learning ANOVA

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 10:58 AM, Stephen Liu wrote: - Original Message From: David Winsemius To: JesperHybel Cc: r-help@r-project.org Sent: Sat, August 14, 2010 5:44:00 PM Subject: Re: [R] Learning ANOVA Hi David, - snip - JH > Also how to save InsectSprays.aov? I think I can only

Re: [R] calibration plot labels

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 10:52 AM, paaventhan jeyaganth wrote: Dear all, when i do the calibration plot, Not a well defined operation i put the x label y label , there is some labels are i did not put it , like "resample optimism added ..." i want to get rid of the these label ,

Re: [R] problems with which

2010-08-15 Thread Nicola Spotorno
You took the point. Thanks a lot, Nicola David Winsemius ha scritto: On Aug 15, 2010, at 4:32 AM, Nicola Spotorno wrote: Dear all, I'm quite new in R and I have a problem with the function which. When I use it to select a subset of a dataframe it works well but somewhere R takes trace of

Re: [R] Learning ANOVA

2010-08-15 Thread Stephen Liu
- Original Message From: David Winsemius To: JesperHybel Cc: r-help@r-project.org Sent: Sat, August 14, 2010 5:44:00 PM Subject: Re: [R] Learning ANOVA Hi David, - snip - >> JH > Also how to save InsectSprays.aov? I think I can only save it as >> InsectSprays.csv. I can't find "wr

[R] calibration plot labels

2010-08-15 Thread paaventhan jeyaganth
Dear all, when i do the calibration plot, i put the x label y label , there is some labels are i did not put it , like "resample optimism added ..." i want to get rid of the these label , is any body know how can i get rid of these label. these are the following command i used cal <

Re: [R] band pass filter

2010-08-15 Thread Matt Shotwell
nuncio, If you already have a filter kernel, you can use the filter function. Of course, convolution filters can be applied directly using the discrete Fourier transform via the fft function. For an example of filtering (lowpass) with R, see http://biostatmatt.com/archives/78 , and the associated

Re: [R] problems with which

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 4:32 AM, Nicola Spotorno wrote: Dear all, I'm quite new in R and I have a problem with the function which. When I use it to select a subset of a dataframe it works well but somewhere R takes trace of the past dataframe and this creates problems with following operation

Re: [R] About .RData

2010-08-15 Thread Stephen Liu
- Original Message From: Michael Bedward To: Stephen Liu Cc: r-help@r-project.org Sent: Sun, August 15, 2010 3:43:55 PM Subject: Re: [R] About .RData >> Furthermore can I run following command on R? >> rm /home/userA/.RData > .RData contains all of the objects in a workspace (your dat

Re: [R] as.logical(factor) behaviour

2010-08-15 Thread R Help
The problem is that, underneath the factors are actually numbers (1 and 2), where as, if you extract the levels and then get the logical, it converts them to strings and then to logicals. I run into this problem ALL THE TIME with numerics in a dataset. Consider the following: > factor(c(3,6,5,2,

Re: [R] problems with which

2010-08-15 Thread Nicola Spotorno
Hi Michael, Thanks a lot! Best, Nicola Michael Bedward ha scritto: Hi Nicola, Although you have subsetted the data.frame, the factor variable will still include the empty level. For example, the following... foo <- data.frame(f = factor(c("a", "b", "c", "b", "a")), val = 1:5) foo2 <- foo[

[R] Paired t-tests

2010-08-15 Thread R Help
Hello List, I'm trying to do a paired t-test, and I'm wondering if it's consistent with equations. I have a dataset that has a response and two treatments (here's an example): ID trt order resp 17 1 0 1 0.0037513592 18 2 0 1 0.0118723051 19 4 0 1 0.0002610251

Re: [R] Trouble loading "saved" Rdata

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 3:06 AM, steven mosher wrote: Did you exit R and then return? fname<-"test.Rdata" full<-file.path("Example",fname,fsep=.Platform$file.sep) full [1] "Example/test.Rdata" load(full) test NULL I am unable to reproduce the problem (after correcting two different syntactic e

Re: [R] Rows index/colProds

2010-08-15 Thread David Winsemius
On Aug 15, 2010, at 4:47 AM, pablo.andrade wrote: Hi, Is there any function to replace colProds that finds column-wise products of a matrix? Or is there any other function that would give a better solution this problem ? I have a matrix and an array, example: a GSM

Re: [R] legend outside plot area

2010-08-15 Thread Pete B
Emily Maybe this will help # Data y1=rnorm(10) y2=runif(10) x=1:10 # Set XPD and Outer Margin par(xpd=NA,oma=c(3,0,0,0)) # Plot and Legend plot(x,y1,type="n",ylab="") lines(x,y1, col="red") lines(x,y2, col="blue") legend(par("usr")[1],par("usr")[3],c("Y1","Y2"),col=c("red","blue"),lty=1,xjust=

Re: [R] problems with which

2010-08-15 Thread Michael Bedward
Hi Nicola, Although you have subsetted the data.frame, the factor variable will still include the empty level. For example, the following... foo <- data.frame(f = factor(c("a", "b", "c", "b", "a")), val = 1:5) foo2 <- foo[ foo$f != "a", ] foo2 str(foo2) Produces this output... f val 2 b 2

Re: [R] legend outside plot area

2010-08-15 Thread Hrishi Mittal
Emily, You need to set the xpd parameter to TRUE (using par) and use a negative value for the inset argument in the legend command. For example to place a legend on the right, legend("right",inset=c(-0.3,0),legend=somelegendvector) See ?legend and play around with the values of inset to get an

Re: [R] non-linear regression for 3D data

2010-08-15 Thread szisziszilvi
Thanks a lot, this I(xx^2) ... worked. I guess, I should learn more abot the function poly itself. (so will I... :) ) Thanks again! -- View this message in context: http://r.789695.n4.nabble.com/non-linear-regression-for-3D-data-tp2320982p2325911.html Sent from the R help mailing list archive a

Re: [R] problems with which

2010-08-15 Thread Nicola Spotorno
Sorry I send again the e-mail because I receive a delivery status notification message and because I made a mistyping in the previous message: 'interaction.plot' and not 'interaction_plot' sorry... Nicola Nicola Spotorno ha scritto: Hello, thanks for your replies, but the logic operator is '!

Re: [R] problems with which

2010-08-15 Thread Nicola Spotorno
Hello, thanks for your replies, but the logic operator is '!='. If I call sentence_trial the dataframe shows no 'an' datapoints as I want but I still have problems with 'str' and 'interaction_plot'. Nicola Tim Gruene ha scritto: Isn't '!=' the operator applied here, not '='? On Sun, Aug

[R] legend outside plot area

2010-08-15 Thread Emily Forbes
Hi, please can you help me. When I add a legend to a boxplot it appears inside the plot area, how do you get it into the margins? I have already changed the parameters so there is space for it on the margin on the right hand side of the graph. Thanks, Emily

[R] Fwd: Simple problem with lm/predict

2010-08-15 Thread Michael Bedward
Hi Nick, Normally the data you provide to predict has the same var names that you used in your regression formula e.g. predict( extrapolate1, newdata=list(t3=something, t4=something.else) ) Or you could have a two col matrix containing t3 and t4 values and provide that... my.matrix <- matrix( c(n

Re: [R] problems with which

2010-08-15 Thread Tim Gruene
Isn't '!=' the operator applied here, not '='? On Sun, Aug 15, 2010 at 09:04:46PM +1000, Michael Bedward wrote: > Hello Nicola, > > > sentences_trial <- sentences[which(sentences$Cond!= "an"),] > > That should be == "an" instead of = "an" > > == is the test for equality whereas = is assignment

Re: [R] problems with which

2010-08-15 Thread Michael Bedward
Hello Nicola, > sentences_trial <- sentences[which(sentences$Cond!= "an"),] That should be == "an" instead of = "an" == is the test for equality whereas = is assignment Michael __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

Re: [R] Moving average in R

2010-08-15 Thread Tim Gruene
Dear Suman, you can download the TTR package through the "Packages"-link on http://cran.r-project.org/ and install it with R CMD INSTALL TTR_0.20-2.tar.gz (or the appropriate name if you use the MacOS X or Windows binary). Tim On Sun, Aug 15, 2010 at 02:59:18PM +0530, suman dhara wrote: >

[R] Moving average in R

2010-08-15 Thread suman dhara
Hi, I want to fit moving average trend in R. In google, I see that it is in the package 'TTR'. But, I can't install this package. I have used the following code: >install.packages("TTR") But, it says there is no package called 'TTR'. Can you help me? Regards, Suman Dhara [[alternativ

Re: [R] Rows index/colProds

2010-08-15 Thread Henrik Bengtsson
What is your sessionInfo()? /Henrik On Sun, Aug 15, 2010 at 10:47 AM, pablo.andrade wrote: > > Hi, > > Is there any function to replace colProds that finds column-wise products of > a matrix? > > Or is there any other function that would give a better solution this > problem ? > > I have a matri

Re: [R] Median abline how-to ?

2010-08-15 Thread Jim Lemon
On 08/13/2010 12:41 AM, David martin wrote: Hi, I'm newbie with R and don't really know how to add a median line to each of the groups that is not all the plot long. Here is a small working code that i have adapted for my purpose. If somebody could tell me how to draw median lines on each group

Re: [R] x-axis label print in 45 degree

2010-08-15 Thread Jim Lemon
On 08/13/2010 05:14 AM, array chip wrote: Hi how can print x-axis labels in 45 degree in boxplot() (or plot in general)? I can use las=2 to print in 90 degree, but it looks ugly. Is there a simple option to do 45 degree easily? Hi John, Have a look at staxlab in the plotrix package, particular

[R] Rows index/colProds

2010-08-15 Thread pablo.andrade
Hi, Is there any function to replace colProds that finds column-wise products of a matrix? Or is there any other function that would give a better solution this problem ? I have a matrix and an array, example: > a GSM1.CEL GSM2.CEL GSM1.CEL 1_at 1

[R] as.logical(factor) behaviour

2010-08-15 Thread Philippe Grosjean
Hello, According to ?as.logical: "as.logical attempts to coerce its argument to be of logical type. For factors, this uses the levels (labels)." However, > as.logical(factor(c("FALSE", "TRUE"))) [1] TRUE TRUE Shouldn't it be the same as: > as.logical(levels(factor(c("FALSE", "TRUE" [1]

[R] problems with which

2010-08-15 Thread Nicola Spotorno
Dear all, I'm quite new in R and I have a problem with the function which. When I use it to select a subset of a dataframe it works well but somewhere R takes trace of the past dataframe and this creates problems with following operations. For example: sentences <- read.xls("frasi.tot.march.3

Re: [R] About .RData

2010-08-15 Thread Michael Bedward
> Furthermore can I run following command on R? >> rm /home/userA/.RData .RData contains all of the objects in a workspace (your data, functions etc.). If you start R in a directory that contains a .RData file it will be loaded automatically. If you delete the file you are deleting your data. D

Re: [R] incrementing matrix elements more efficiently

2010-08-15 Thread Romain Francois
Le 15/08/10 02:43, david h shanabrook a écrit : I need to increment cells of a matrix (collusionM). The indexes to increment are in an index (matchIndex). This is sample code library(seqinr) library(Matrix) x<- "abcabcabc" mx<- s2c(x) collisionM<- Matrix(0,nrow=10, ncol=1

Re: [R] Trouble loading "saved" Rdata

2010-08-15 Thread steven mosher
Did you exit R and then return? fname<-"test.Rdata" full<-file.path("Example",fname,fsep=.Platform$file.sep) full [1] "Example/test.Rdata" load(full) test NULL > sessionInfo() R version 2.11.1 (2010-05-31) x86_64-apple-darwin9.8.0 locale: [1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8

[R] About .RData

2010-08-15 Thread Stephen Liu
Hi folks, OS Ubuntu 10.04 64 bit On starting R following line is displayed; ... [Previously saved workspace restored] On terminal; $ locate .RData /home/userA/.RData /usr/lib/R/site-library/qtl/data/badorder.RData /usr/lib/R/site-library/qtl/data/bristle3.RData /usr/lib/R/site-library/qtl/d