Re: [R] "na.strings" and the like; suspending interpretation of "NA"

2009-08-03 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Jan > Theodore Galkowski > Sent: Monday, August 03, 2009 10:23 PM > To: R Project > Subject: [R] "na.strings" and the like; suspending > interpretation of "NA" > > Can someone

Re: [R] RServe - How to use 'createReference' method?

2009-08-03 Thread Romain Francois
Hi, The Rserve implementation of REngine does not support references (yet?). Anyway, I don't think references are what you need here. You probably want the methods createFile and openFile that create RFileInputStream and RFileOutputStream which you can use to transfer files through the R ser

Re: [R] Add columns in a dataframe and fill them from another table according to a criteria

2009-08-03 Thread Meenu Sahi
Thank you very much !!! It works !!! Sorry it took me some time to understand your solution and apply it to my extended dataset. Just got it working !!! Yes as.numeric is the trick.. Many thanks again for your help and time. Regards Meenu On Sun, Aug 2, 2009 at 6:06 AM, David Winsemius wrote: >

[R] "na.strings" and the like; suspending interpretation of "NA"

2009-08-03 Thread Jan Theodore Galkowski
Can someone point me to the proper place in the documentation or on the Wiki where I can learn how to get R to stop interpreting the string "NA" as something special? I have a table in a database which contains (among other things) country codes and continent codes. The standard set of two-letter

Re: [R] Save model and predictions from svm

2009-08-03 Thread Steve Lianoglou
2) I'm coming from RapidMinder, so am still getting used to how data is handled in R. (In RM, everything is like the R data.frame and predictions are automatically appended as new columns to the data.) What I'd like is this: Starting with data frame of: label, v1, v2, v3 After svm prediction

Re: [R] Save model and predictions from svm

2009-08-03 Thread Noah Silverman
Thanks Steve, 1) That helps. Exactly what I need. 2) I'm coming from RapidMinder, so am still getting used to how data is handled in R. (In RM, everything is like the R data.frame and predictions are automatically appended as new columns to the data.) What I'd like is this: Starting with da

Re: [R] Save model and predictions from svm

2009-08-03 Thread Steve Lianoglou
Hi, On Aug 3, 2009, at 10:55 PM, Noah Silverman wrote: Hello, I'm using the e1071 package for training an SVM. It seems to be working well. This question has two parts: 1) Once I've trained an SVM model, I want to USE it within R at a later date to predict various new data. I see the

[R] status of boolean package

2009-08-03 Thread Ronggui Huang
FYI -- Forwarded message -- From: Bear Braumoeller Date: 2009/8/4 Subject: Re: Package boolean on CRAN To: "c...@datanalytics.com" Cc: Ronggui Huang , Ben Goodrich Gentlemen, Thanks for your enquiries.  Boolean is under active development, and we should have a new release on

[R] Save model and predictions from svm

2009-08-03 Thread Noah Silverman
Hello, I'm using the e1071 package for training an SVM. It seems to be working well. This question has two parts: 1) Once I've trained an SVM model, I want to USE it within R at a later date to predict various new data. I see the write.svm command, but don't know how to LOAD the model back

[R] graph data mining

2009-08-03 Thread Weiwei Shi
Hi I am wondering if R or some others provide free-to-use graph mining tools, like mining some frequent structure in a x-y plot data? thanks, Weiwei -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. "Did you always know?" "No, I did not. But I believed..." ---Matrix III [[alternativ

Re: [R] lme funcion in R

2009-08-03 Thread Hongwei Dong
Thanks. I tried to set a higher tolerance for the convergence, such as changing tolerance from 1e-6 to 1e-5, msTol from 1e-7 to 1e-6, but R still does not converge. Any more suggestions? Harry On Mon, Aug 3, 2009 at 10:36 AM, Jason Morgan wrote: > On 2009.08.03 10:15:46, Hongwei Dong wrote: >

Re: [R] Character from Symbol Font on rgl Plot

2009-08-03 Thread Duncan Murdoch
On 02/08/2009 4:49 PM, Tom La Bone wrote: I want to print characters from the symbol font (or perhaps even Wingdings) in an rgl 3d plot, but I am having no luck. So, what do I have to do in order to get this snippet to print out a character from the symbol font? library(rgl) open3d() text3d(1,1,

Re: [R] Help with reshaping data.frame

2009-08-03 Thread Gabor Grothendieck
Try this: > tst2 <- tst[-(2:3)] > tst2$timevar <- with(tst, interaction(K2, K3)) > reshape(tst2, dir = "wide", idvar = "K1", timevar = "timevar") K1 V1.D.a V2.D.a V3.D.a V1.E.a V2.E.a V3.E.a V1.D.b V2.D.b V3.D.b V1.E.b V2.E.b V3.E.b 1 10 0.08 99105 NA NA NA NA 76

Re: [R] Merge partially duplicated rows

2009-08-03 Thread David Winsemius
On Aug 3, 2009, at 7:12 PM, David Winsemius wrote: On Aug 3, 2009, at 9:24 AM, Rnewbie wrote: Dear all, I have a dataset, and I wanted to merge the rows with duplicated IDs by calculating the means or medians from the duplicate rows. I tried using the command duplicated(x), but it only

Re: [R] Help with data type

2009-08-03 Thread Don MacQueen
It's because Delta isn't a number. Or, more accurately, Delta[1] and Delta[i-1] are not numbers. The are objects that have a class, and therefore an internal structure: class(Delta[1]) [1] "xts" "zoo" To understand what's going on, you need to understand the structure of these classes, an

Re: [R] Installing package for i386

2009-08-03 Thread Sarah Hawley
Hi Mark, Thanks for the suggestion and for setting me straight about my operating system. I will contact my collaborator. Best, Sarah On 8/3/09 1:59 PM, "Marc Schwartz" wrote: > On Aug 3, 2009, at 3:01 PM, Sarah Hawley wrote: > >> I am trying to install a package from a collaborator that was

Re: [R] Help with Ecdf function

2009-08-03 Thread Mateus Teixeira
David, Actually, I have used the example given in Ecdf help page, see below a part of the script: comparacao<-c(obs,simul[,,1],simul[,,2],simul[,,3],simul[,,4]) grupos<-c( rep('obs',length(obs)), rep('exp1',length(simul[,,1])), rep('exp2',length(simul[,,2])), rep('exp3',length(simul[,,3])), rep('

Re: [R] Merge partially duplicated rows

2009-08-03 Thread David Winsemius
On Aug 3, 2009, at 9:24 AM, Rnewbie wrote: Dear all, I have a dataset, and I wanted to merge the rows with duplicated IDs by calculating the means or medians from the duplicate rows. I tried using the command duplicated(x), but it only tells where the duplicated rows are. You might wa

Re: [R] Help with Ecdf function

2009-08-03 Thread David Winsemius
On Aug 3, 2009, at 4:19 PM, Mateus Teixeira wrote: Dear R users, I'm using Ecdf (Hmisc library) to plot four cdf in a same graphic. In this graphic I also plot the 0.99 quantile for these cdf. How did you do this? Experimenting with the example on hte Ecdf help page it appears that addi

Re: [R] Help with data type

2009-08-03 Thread David Winsemius
On Aug 3, 2009, at 10:48 AM, ehxpieterse wrote: Using a quantmod function, I calculate the daily change between two points in a time series. However, I don't think I am using the data correctly. Code: getSymbols("^GSPC", src="yahoo") CloseData <- Cl(GSPC) Delta <- diff(CloseData, lag=1) for

Re: [R] scatterplot3d bug??

2009-08-03 Thread David Winsemius
Not sure you can call it a bug when the help page says for angles in range (180, 360) that some support functions may "not work properly". On Aug 3, 2009, at 2:12 PM, Vivek Ayer wrote: Hey guys, Not sure if I encountered a bug with the scatterplot3d function. Here's the calls I made: s3d1

Re: [R] What does this error message mean?

2009-08-03 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Hongwei Dong > Sent: Monday, August 03, 2009 2:26 PM > To: r-help@R-Project.org > Subject: [R] What does this error message mean? > > Hi, I used R to run a linear regression and

Re: [R] Truncating based on attribute range and serial no

2009-08-03 Thread Jorge Ivan Velez
Hi JR, There are few ways to do what you asked for, so here is one suggestion using tapply: Age <- c(48, 57, 56, 76, 76, 66, 70, 14, 7, 3, 62, 62, 30, 10, 7, 53, 44, 29, 46, 47, 15, 13, 84, 77, 26) SerialNo <- c(001147, 005979, 005979, 006128, 006128, 007004, 007004, 007004, 007004, 007004, 00

Re: [R] Truncating based on attribute range and serial no

2009-08-03 Thread Marc Schwartz
On Aug 3, 2009, at 2:04 PM, PDXRugger wrote: COnsider the following: Age<-c(48, 57, 56, 76, 76, 66, 70, 14, 7, 3, 62, 62, 30, 10, 7, 53, 44, 29, 46, 47, 15, 13, 84, 77, 26) SerialNo<-c(001147, 005979, 005979, 006128, 006128, 007004, 007004, 007004, 007004, 007004, 007438, 007438,0094

Re: [R] What does this error message mean?

2009-08-03 Thread David Winsemius
On Aug 3, 2009, at 5:26 PM, Hongwei Dong wrote: Hi, I used R to run a linear regression and keep getting the following error message. I do not understand it very well. Anyone can help out? Thanks. Error in storage.mode(y) <- "double" : invalid to change the storage mode of a factor In add

Re: [R] Help with reshaping data.frame

2009-08-03 Thread hadley wickham
On Mon, Aug 3, 2009 at 5:23 PM, wrote: > I'm having trouble reshaping a data.frame from long to wide. > (I think that's the right terminology; feel free to educate me.) > I've looked at the reshape function and package and plyr package, > but I can't quite figure out how to do this after a dozen v

Re: [R] Help with reshaping data.frame

2009-08-03 Thread rmailbox
How about: tst.long <- melt ( tst, id.vars = c("K1", "K2", "K3"), variable_name = "V" ) tst2 <- transform ( tst.long, K23Vx = paste ( K2, K3, V, sep = ".") ) tst.wide <- cast ( tst2 , K1 ~ K23Vx ) tst.wide > tst.wide K1 D.a.V1 D.a.V2 D.a.V3 D.b.V1 D.b.V2 D.b.V3 E.a.V1 E.a.V2 E.a.V3 E.b.V1

Re: [R] penalized logistic regression

2009-08-03 Thread Marc Schwartz
On Aug 3, 2009, at 10:50 AM, annie Zhang wrote: Hi, R users, Is there any package for penalized logistic regression with more than two response classes? I read the manual for stepPlr, but it seems it's only for binary case. Thank you, Annie Take a look at the lrm() function in Frank's

[R] fitting a truncated power law

2009-08-03 Thread John Sanders
Dear all, How can I fit a truncated power law to a vector? I can't find a function to do that. If the function provides an AIC, even better. John __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the po

[R] Truncating based on attribute range and serial no

2009-08-03 Thread PDXRugger
COnsider the following: Age<-c(48, 57, 56, 76, 76, 66, 70, 14, 7, 3, 62, 62, 30, 10, 7, 53, 44, 29, 46, 47, 15, 13, 84, 77, 26) SerialNo<-c(001147, 005979, 005979, 006128, 006128, 007004, 007004, 007004, 007004, 007004, 007438, 007438,009402,009402, 009402, 012693, 012693, 012693, 014063,0140

Re: [R] boxplot( ) headers with Greek letters, values, and text

2009-08-03 Thread Johnson, Roger W.
Thank you! (This works.) Roger From: Gabor Grothendieck [ggrothendi...@gmail.com] Sent: Monday, August 03, 2009 12:29 PM To: Johnson, Roger W. Cc: r-help@R-Project.org Subject: Re: [R] boxplot( ) headers with Greek letters, values, and text Try this: box

Re: [R] how do i retain decimal values

2009-08-03 Thread Manisha Brahmachary
Thanks Patrick,the options(digits=) works. Manisha Here is the code: rand.val<- list() rand.val<-data.frame(a=runif(11,0.1,0.0001),b=runif(11,0. 1,0.0001), c=runif(11,0.1,0.0001),d=runif(11,0.1,0.0001),e=runif (11,0.1,0.0001),f=runif(11

[R] RODBC fails to read ArcGIS mdb shape table at 2.9.1; works at

2009-08-03 Thread Phil Hurvitz
I've found what might be a bug in RODBC 1.3-0. When reading (with `sqlFetch') an ESRI format personal geodatabase table containing binary shape data (MS Access mdb format), R crashes. With RODBC 1.2-6, R does not crash. Details, code, sample data, and errors are posted here: http://gis.washi

Re: [R] how do i retain decimal values

2009-08-03 Thread Manisha Brahmachary
Thanks Peter, the options(digits=) works. The matrix I send as example is the x values in the correlation I will do. The Y values are different and are not tied. Thanks again. manisha -Original Message- From: P Ehlers [mailto:ehl...@ucalgary.ca] Sent: Sunday, August 02, 2009 10:32 AM To

[R] Merge partially duplicated rows

2009-08-03 Thread Rnewbie
Dear all, I have a dataset, and I wanted to merge the rows with duplicated IDs by calculating the means or medians from the duplicate rows. I tried using the command duplicated(x), but it only tells where the duplicated rows are. Any suggestions will be appreciated. -- View this message in cont

[R] RServe - How to use 'createReference' method?

2009-08-03 Thread joaodaniel
I need to input a txt, or xls, file from a client to R, using RServe. >From what I've been reading, the best way to do this, is using the 'createReference' method, from REngine package. But I couldn't find any documents exemplifying it's use. I got to upload a file from java? And then? How do i

[R] Some SQL Challenges

2009-08-03 Thread JoK LoQ
Hi, Im trying to complete a list of jobs using SQL Querries and some "if else" commands but im stucked in some steps. Could any of you give me some help? -in COLUMN1 change the format 20JAN2000:00:00:00 to 20JAN2000 and exclude every row that date is different from 20. -extract the first charac

[R] Help with data type

2009-08-03 Thread ehxpieterse
Hi there, Using a quantmod function, I calculate the daily change between two points in a time series. However, I don't think I am using the data correctly. Code: getSymbols("^GSPC", src="yahoo") CloseData <- Cl(GSPC) Delta <- diff(CloseData, lag=1) for (i in 3:length(Delta)) { if (Delta[i]>Del

Re: [R] Variable names as inputs...

2009-08-03 Thread Horacio Samaniego
That's a great guide!!! Thanks I guess, I wasn't 2 clear about my question. I have succeed naming my vars, but I can't seem to be able to create a c() of objects named after the string. So far I have MyList$lm1 MyList$lm2, ... but can't use them as an input to AIC() that is, AIC(Mylist) won

[R] What does this error message mean?

2009-08-03 Thread Hongwei Dong
Hi, I used R to run a linear regression and keep getting the following error message. I do not understand it very well. Anyone can help out? Thanks. Error in storage.mode(y) <- "double" : invalid to change the storage mode of a factor In addition: Warning message: In model.response(mf, "numeric"

[R] Help with reshaping data.frame

2009-08-03 Thread davidr
I'm having trouble reshaping a data.frame from long to wide. (I think that's the right terminology; feel free to educate me.) I've looked at the reshape function and package and plyr package, but I can't quite figure out how to do this after a dozen variations. I have a data.frame with more level

Re: [R] session logging

2009-08-03 Thread Marc Schwartz
Jacob, As far as I am aware, the OSX GUI has no logging/transcript feature per se, short of saving the console output to an external file. The history file will, as you note, only contain the commands entered at the R prompt, not the additional output. One of the solutions would be to use

Re: [R] Open data-file in R - from Rserve

2009-08-03 Thread joaodaniel
Do I have to use to openFile method? How can I assign this file to a dataframe? joaodaniel wrote: > > How can I open a data file from the client on the RServe server using > Java? > -- View this message in context: http://www.nabble.com/Open-data-file-in-R---from-Rserve-tp24795557p24798244.

Re: [R] Installing package for i386

2009-08-03 Thread David Winsemius
On Aug 3, 2009, at 4:01 PM, Sarah Hawley wrote: I am trying to install a package from a collaborator that was given to me as a zip file. I'm running R 2.9.0 using Mac OSX 10.5.7 [R.app GUI 1.28 (5395) i386-apple-darwin8.11.1] foo.zip contains all the usual package directories (man, R, R-ex,

Re: [R] Installing package for i386

2009-08-03 Thread Marc Schwartz
On Aug 3, 2009, at 3:01 PM, Sarah Hawley wrote: I am trying to install a package from a collaborator that was given to me as a zip file. I'm running R 2.9.0 using Mac OSX 10.5.7 [R.app GUI 1.28 (5395) i386-apple-darwin8.11.1] foo.zip contains all the usual package directories (man, R, R-ex, e

Re: [R] hist2d

2009-08-03 Thread Steve Lianoglou
Hi, On Aug 3, 2009, at 4:44 PM, simon.b.firest...@frb.gov wrote: I'd like to use the hist2d function. What library contains it? Thanks! Try: R> RSiteSearch('hist2d') Several signs point to the gplots package ... try that one. -steve -- Steve Lianoglou Graduate Student: Computational Sy

Re: [R] if confusion

2009-08-03 Thread Rolf Turner
On 4/08/2009, at 8:40 AM, rkevinbur...@charter.net wrote: Simple question: Why doesn't the following work? Or what 'R' rule am I missing? tclass <- "Testing 1 2 3" if(tclass == "Testing 1 2 3") { cat("Testing", tclass, "\n") } else { cat(tclass, "\n") } I get an error 'else' is unex

Re: [R] if confusion

2009-08-03 Thread Erik Iverson
This is because you could imagine a program like the following: tclass <- "Testing 1 2 3" if(tclass == "Testing 1 2 3") { cat("Testing", tclass, "\n") } cat("more stuff") However, in your program, the next line is "else", which is not a command you can just put in the middle of your progra

Re: [R] if confusion

2009-08-03 Thread David Huffer
Because the line if... Is syntactically complete. Move the else statement to the line above it like tclass <- "Testing 1 2 3" if ( tclass == "Testing 1 2 3" ) { cat ( "Testing" , tclass , "\n" ) } else { cat ( tclass , "\n" ) } david -- David   ---

[R] hist2d

2009-08-03 Thread Simon . B . Firestone
I'd like to use the hist2d function. What library contains it? Thanks! Simon Firestone Economist, Fair Lending Enforcement Division of Consumer and Community Affairs Federal Reserve Board Washington, DC 20551 simon.b.firest...@frb.gov (202) 785-6056 [[alternative HTML version deleted]]

[R] if confusion

2009-08-03 Thread rkevinburton
Simple question: Why doesn't the following work? Or what 'R' rule am I missing? tclass <- "Testing 1 2 3" if(tclass == "Testing 1 2 3") { cat("Testing", tclass, "\n") } else { cat(tclass, "\n") } I get an error 'else' is unexpected. Thank you. Kevin __

Re: [R] session logging

2009-08-03 Thread Jacob Wegelin
Maybe I'm missing something in the history help file. I read it and subsequently tried > savehistory("trythis.txt") but the resulting file does not contain the text that scrolled past the console. It contains only a few uninteresting timestamps. Is there some more useful application of history()

[R] Help with Ecdf function

2009-08-03 Thread Mateus Teixeira
Dear R users, I'm using Ecdf (Hmisc library) to plot four cdf in a same graphic. In this graphic I also plot the 0.99 quantile for these cdf. I successfully plot cdfs using different types of line to distinguish them, but I can't determine the type of lines showing 0.99 quantile. Is there a way t

[R] Installing package for i386

2009-08-03 Thread Sarah Hawley
I am trying to install a package from a collaborator that was given to me as a zip file. I'm running R 2.9.0 using Mac OSX 10.5.7 [R.app GUI 1.28 (5395) i386-apple-darwin8.11.1] foo.zip contains all the usual package directories (man, R, R-ex, etc) After unzipping foo, I used the following comman

[R] session logging

2009-08-03 Thread Jacob Wegelin
Consider all the text that one sees on the console during an R session. Is there a way, within R, that all this text--i.e., both the "output" and the "messages"--can be automatically copied to a single file, in addition to its being visible on the console? If I remember to save the console to a f

Re: [R] status of boolean?

2009-08-03 Thread Uwe Ligges
Ronggui Huang wrote: package of boolean has been removed from CRAN. Anyone know why? Thanks. Yes, because the maintainer has not fixed the issues the package had with R CMD check and got ORPHANED. You can try to contact the maontainer in order to take over maintainer ship. Unfortunately, th

[R] MSM package and qmatrix

2009-08-03 Thread Ross Culloch
Hi R-helpers, I am having a problem understanding how to construct the qmatrix, i have read the help menu for msm but I am still a bit lost. Below is an example of one of my transition matrices: statetable.msm(BEH, ID, data = A1) to from 1 2 3 4 5 23 200 201 203 999 1 86 11 2 20

[R] session logging

2009-08-03 Thread Jacob Wegelin
Consider all the text that one sees on the console during an R session. Is there a way, within R, to make all this text--both the "output" and the "messages"--automatically get copied to a single text file, in addition to seeing it on the console? If I remember to save the console to a file a

Re: [R] boxplot( ) headers with Greek letters, values, and text

2009-08-03 Thread Gabor Grothendieck
Try this: boxplot(1:5, main = bquote(mu == .(mu) * "," ~ lambda == .(lambda) * "," ~ trunction == .(truncation) * "," ~ N[T] == .(n))) On Sun, Aug 2, 2009 at 8:37 PM, Johnson, Roger W. wrote: > Hi - > > I've been using the option > > main=bquote(paste(mu==.(mu),", ",lambda==.(lambda),", > ",tru

Re: [R] Lattice xyplot: same scales within one factor

2009-08-03 Thread OB
Thank you! I ended up hunting down a few previous posts you made and getting it to to work, although yes, it was not pretty. -Orion On 8/3/09 10:30 am, "Deepayan Sarkar" wrote: > On Mon, Jul 13, 2009 at 8:33 AM, OB wrote: >> I am using R 2.8.1 and lattice to produce xyplots conditioned on >> tw

Re: [R] wireframe color

2009-08-03 Thread Deepayan Sarkar
On Mon, Aug 3, 2009 at 9:36 AM, Marlin Keith Cox wrote: > Hello, hopefully simple question, but I cannot find the answer.  I need to > change the color from the standard default.  Still want the scaled colors, > but need different colors for different graphs. > > Code is: > > wireframe(z ~ y*x, mat

[R] Open data-file in R - from Rserve

2009-08-03 Thread joaodaniel
How can I open a data file from the client on the RServe server using Java? -- View this message in context: http://www.nabble.com/Open-data-file-in-R---from-Rserve-tp24795557p24795557.html Sent from the R help mailing list archive at Nabble.com. __ R

Re: [R] Defining different objects within a loop

2009-08-03 Thread Deepayan Sarkar
On Thu, Jul 30, 2009 at 3:09 PM, jimdare wrote: > > Dear R Users, > > I have a dataset that I have split by group.  For each group I plot a figure > using: > > for (i in splitdf){ > plot<-xyplot() > print(plot) > savePlot(filename=paste(i$Group[1]),type="pdf") > } > > This gives me

Re: [R] lattice shingle plot axis annotation

2009-08-03 Thread Deepayan Sarkar
On Thu, Jul 30, 2009 at 12:43 AM, Armin Goralczyk wrote: > Hello (R-)Experts > > I hope someone can help with this problem concerning axis annotation > of a lattice shingle plot. I want a plot with three shingles to > display some laboratory value over time. In the first panel over the > first few

Re: [R] boxplot( ) headers with Greek letters, values, and text

2009-08-03 Thread cameron.bracken
Johnson, Roger W. wrote: > > Any suggestions on how I can produce the same title using the boxplot > command? > If you use LaTeX you could try the tikzDevice package ( https://r-forge.r-project.org/projects/tikzdevice/ ). In that case you would use: main = paste('\\mu =',mu,', \\lambda =',la

Re: [R] dot plot with several points for 2 categories

2009-08-03 Thread Deepayan Sarkar
On Wed, Jul 29, 2009 at 9:11 AM, jaregi wrote: > > Hi Michael, Steve, and 1Rnwb, > > I'm very impressed by the quick replies on the mailer. Thanks a lot for your > suggestions. They worked very well. > > In general, I have to say that I'm a bit disappointed that in R, like in > Excel, one basically

Re: [R] principal component analysis for class variables:somenone has experience with made 4 and plotgenes?

2009-08-03 Thread andreiabb
Alain Guillet-2 wrote: > > The PCA doesn't work with class variables so the error is normal. You > should try to work with a discriminant factorial analysis (see > discrimin.coa in ade4). > > Alain thanks for the tip but it did not work, well then I manaege to put > it to work but the res

Re: [R] lme funcion in R

2009-08-03 Thread Hongwei Dong
Thanks for the replies above. Here are my script and data structure: library(nlme) tlevel<-lme(fixed = LN_unitlandval ~ MH_D+APT_D+ResOth_D+NonRes_D+Vacant_D+access_emp1+pct_vacant+transit_D+park_dum,data=lusdrdata,random = ~MH_D+APT_D+ResOth_D+NonRes_D+Vacant_D | TAZ) str: $ TAZ : int 100 100 10

Re: [R] lme funcion in R

2009-08-03 Thread Jason Morgan
On 2009.08.03 10:15:46, Hongwei Dong wrote: > Hi, R users, > I'm using the "lme" function in R to estimate a 2 level mixed effects > model, in which the size of the subject groups are different. It turned out > that It takes forever for R to converge. I also tried the same thing in SPSS > and SPS

Re: [R] Lattice xyplot: same scales within one factor

2009-08-03 Thread Deepayan Sarkar
On Mon, Jul 13, 2009 at 8:33 AM, OB wrote: > I am using R 2.8.1 and lattice to produce xyplots conditioned on > two factors. What I would like is to have the scales be free between values > of one factor, but some within. Thus, in this example, > > xyplot(mpg ~ disp | factor(gear) + factor(cyl), mt

Re: [R] lme funcion in R

2009-08-03 Thread Chuck Cleland
On 8/3/2009 1:15 PM, Hongwei Dong wrote: > Hi, R users, > I'm using the "lme" function in R to estimate a 2 level mixed effects > model, in which the size of the subject groups are different. It turned out > that It takes forever for R to converge. I also tried the same thing in SPSS > and SPSS c

[R] lme funcion in R

2009-08-03 Thread Hongwei Dong
Hi, R users, I'm using the "lme" function in R to estimate a 2 level mixed effects model, in which the size of the subject groups are different. It turned out that It takes forever for R to converge. I also tried the same thing in SPSS and SPSS can give the results out within 20 minutes. Anyone c

Re: [R] principal component analysis for class variables

2009-08-03 Thread Hongwei Dong
Hi, R users, I'm using the "lme" function in R to estimate a 2 level mixed effects model, in which the size of the groups are different. It turned out that It takes forever for R to converge. I also tried the same thing in SPSS and SPSS can give the results out within 20 minutes. Anyone can give

Re: [R] principal component analysis for class variables

2009-08-03 Thread Alain Guillet
The PCA doesn't work with class variables so the error is normal. You should try to work with a discriminant factorial analysis (see discrimin.coa in ade4). Alain andreiabb wrote: Dear Forum, I have a class variable 1 (populations A-E), and two other class variables, variable 2 and variabl

[R] wireframe color

2009-08-03 Thread Marlin Keith Cox
Hello, hopefully simple question, but I cannot find the answer. I need to change the color from the standard default. Still want the scaled colors, but need different colors for different graphs. Code is: wireframe(z ~ y*x, mat.df, drape = TRUE, zlab = list("Water mass error

[R] principal component analysis for class variables

2009-08-03 Thread andreiabb
Dear Forum, I have a class variable 1 (populations A-E), and two other class variables, variable 2 and variable 3. What I want is to see if the combination of var 2 and var 3, will give me a pattern that allows to distinguish populations. I found several packages like ade4, with pcaiv function a

Re: [R] Cox ridge regression

2009-08-03 Thread ljubomir
Thank you Terry, that answered all questions. As a suggestion, help page for ridge() might indicate that the ridge term simultaneously introduces predictors and penalizes them. Ljubomir From: Terry Therneau To: Ljubomir Buturovic Cc: r-help@r-project.org Subject: Re: Cox ridge regression Dat

Re: [R] For Loop Iteration: Variable Counter

2009-08-03 Thread Vivek Ayer
Got it..Thanks Vivek On Mon, Aug 3, 2009 at 9:13 AM, Steve Lianoglou wrote: > Hi, > > On Aug 3, 2009, at 12:09 PM, Vivek Ayer wrote: > >> Hey guys, >> >> How do I iterate such that I add 100 to the counter every time? >> >> Suppose: for (i in c(1:100)) >> >> I want i to be 1, 10, 20, 30, ... inst

Re: [R] For Loop Iteration: Variable Counter

2009-08-03 Thread Steve Lianoglou
Hi, On Aug 3, 2009, at 12:09 PM, Vivek Ayer wrote: Hey guys, How do I iterate such that I add 100 to the counter every time? Suppose: for (i in c(1:100)) I want i to be 1, 10, 20, 30, ... instead of 1,2,3,4,5 ... How can this be done? ?seq -steve -- Steve Lianoglou Graduate Student: Com

[R] For Loop Iteration: Variable Counter

2009-08-03 Thread Vivek Ayer
Hey guys, How do I iterate such that I add 100 to the counter every time? Suppose: for (i in c(1:100)) I want i to be 1, 10, 20, 30, ... instead of 1,2,3,4,5 ... How can this be done? Thanks, Vivek __ R-help@r-project.org mailing list https://stat.e

Re: [R] How to catch an error using try

2009-08-03 Thread jim holtman
gotError <- try(gene.seq <- getSequence (id=gene.map[,"ensembl_transcript_id"], type="ensembl_transcript_id", seqType="3utr", mart=hmart)) if (inherits(gotError, 'try-error')){.error occurred, do recovery} On Mon, Aug 3, 2009 at 11:53 AM, wrote: > Sometimes the following function call

Re: [R] How to catch an error using try

2009-08-03 Thread Steve Lianoglou
Hi, On Aug 3, 2009, at 11:53 AM, wrote: Sometimes the following function call causes a database exception: gene.seq <- getSequence (id=gene.map[,"ensembl_transcript_id"], type="ensembl_transcript_id", + seqType="3utr", mart=hmart) I understand the above functi

Re: [R] selectively altering variable value

2009-08-03 Thread jim holtman
This should set the first two ia==4 entries to 1 in saccade column: > x <- read.table('/tempxx.txt.r', header=TRUE) > x vp trial ia saccade 1 1 1 5 0 2 1 1 5 0 3 1 1 5 0 4 1 1 5 1 5 1 1 5 1 6 1 1 9 1 7 1 1 9

Re: [R] plm summary error

2009-08-03 Thread Achim Zeileis
Ulrich, as usual with such specific questions, it is a good idea to have the package authors at least Cc (done now). In April I had already forwarded your message privately to the plm authors with the following comment: I've quickly checked what is going on. If you do something like foo <

[R] survdiff for left-truncated data?

2009-08-03 Thread Rajen Shah
Hi Does anyone know if there is a function like survdiff which can also handle left-truncated and right-censored data? When I use it on left-truncated and right-censored data I get an error message saying Right censored data only. Many thanks Rajen [[alternative HTML version deleted]]

[R] How to catch an error using try

2009-08-03 Thread mauede
Sometimes the following function call causes a database exception: > gene.seq <- getSequence (id=gene.map[,"ensembl_transcript_id"], > type="ensembl_transcript_id", + seqType="3utr", mart=hmart) I understand the above function must be called by try to capture the ev

[R] penalized logistic regression

2009-08-03 Thread annie Zhang
Hi, R users, Is there any package for penalized logistic regression with more than two response classes? I read the manual for stepPlr, but it seems it's only for binary case. Thank you, Annie [[alternative HTML version deleted]] __ R-help@r-

[R] Variable names as inputs...

2009-08-03 Thread Horacio Samaniego
Hi, I have been using R for a bit, but never came across this simple issue... So, I wonder if anyone could give me a hint on it. I have calculated a bunch o models and now need to compare them using AIC(). The issue seem to be that I have placed each output in a list object. So, the question is h

[R] plm summary error

2009-08-03 Thread Ulrich Morawetz
Dear "plm"-Package insiders, [I posted the following observation is April already but unfortunately I am not aware of any answers. With the hope that someone found an answer in the mean time, I ask again:] I realized the following difficulty with the summary.plm function (demonstrated wit

[R] problem with predict(...,type="terms")

2009-08-03 Thread Antonio.Gasparrini
Hello R users, I'm trying to understand the meaning of the results returned by the command 'predict' with type="terms". I thought it was the matrix with the fitted values for each term (value of predictor * related coef), but the calculation is different. I include an example below: y <- rnor

[R] selectively altering variable value

2009-08-03 Thread Jens Bölte
Hello, I have run an eye-tracking experiment for which I now like to analyse the saccades. Participants looked from a fixation cross (ia = 5) to the target area (ia = 4) in following example of a data frame. ia = 9 stands for everything else. A saccade is indicated by saccade = 1. Sometimes the s

Re: [R] min frequencies of categorical predictor variables in GLM

2009-08-03 Thread Marc Schwartz
On Aug 3, 2009, at 12:06 AM, Thomas Mang wrote: Hi, Suppose a binomial GLM with both continuous as well as categorical predictors (sometimes referred to as GLM-ANCOVA, if I remember correctly). For the categorical predictors = indicator variables, is then there a suggested minimum frequen

Re: [R] Table to Data-frame

2009-08-03 Thread Jorge Ivan Velez
Hi Olivier, How about this? # some data x <- rpois(100,5) # the table table(x) # converting the table as.data.frame.table(table(x)) HTH, Jorge On Mon, Aug 3, 2009 at 10:35 AM, Olivier Peron wrote: > Hello > > How to pass from table to data-frame ? > > Thanks > >

Re: [R] Table to Data-frame

2009-08-03 Thread Steve Lianoglou
Hi, On Aug 3, 2009, at 10:35 AM, Olivier Peron wrote: Hello How to pass from table to data-frame ? R> a <- sample(LETTERS[1:10], 20, replace=T) R> b <- table(a) R> c <- as.data.frame(b) R> c a Freq 1 A4 2 B2 3 C1 4 D4 5 E1 6 F1 7 H1 8 J6 -steve -- Steve Li

Re: [R] Table to Data-frame

2009-08-03 Thread Paulo Eduardo Cardoso
as.data.frame(table())? 2009/8/3 Olivier Peron > Hello > > How to pass from table to data-frame ? > > Thanks > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-pro

Re: [R] apply function to named numeric object

2009-08-03 Thread Paulo Eduardo Cardoso
Yes, mean of a one element object is meanless, but I'm trying to avoid a ifelse in the sentence. The result of mean values of x and y is added at each loop to a data.frame (or any other object that could speed up the loop. 2009/8/3 ONKELINX, Thierry > Dear Paulo, > > Your object is a named vec

Re: [R] Table to Data-frame

2009-08-03 Thread David Winsemius
On Aug 3, 2009, at 10:35 AM, Olivier Peron wrote: Hello How to pass from table to data-frame ? Copied from the help page for table: "The as.data.frame method for objects inheriting from class "table" can be used to convert the array-based representation of a contingency table to a data

[R] Table to Data-frame

2009-08-03 Thread Olivier Peron
Hello How to pass from table to data-frame ? Thanks __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-con

Re: [R] apply function to named numeric object

2009-08-03 Thread ONKELINX, Thierry
Dear Paulo, Your object is a named vector. Apply() and colMeans() only work on matrices and dataframes. So you need to convert the vector to a matrix or data.frame. See the example below. Please not that calculating the mean of one element is a waiste of time. # a named numerical vector object <-

Re: [R] apply function to named numeric object

2009-08-03 Thread Paulo Eduardo Cardoso
The objects are > centroid x y 669.5 313.5 > str(centroid) Named num [1:2] 670 314 - attr(*, "names")= chr [1:2] "x" "y" or > centroid x y [1,] 659.5 472.5 [2,] 659.5 473.5 > str(centroid) num [1:2, 1:2] 660 660 472 474 - attr(*, "dimnames")=List of 2 ..$ : NULL ..$ :

  1   2   >