Re: [R] Can't get sample function from "An Introduction to R" to work

Thank you-all very much for your help, your responses and help has been very encouraging. The following doesn't close the case but it tables it... First I copied Ken's code into my R Console and...it worked great! That was baffling as it looked identical to mine. I did not explicitly say ear

Re: [R] question from environmental statistics

?qqplot Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Friday, July 15, 2005 9:03 PM To: r-help@stat.math.ethz

[R] question from environmental statistics

thanks Fran. that was useful but Im still in a fix. its a real life data which looks like this: 0.9 10.9 24.0 6.7 0.6 1.0 2.4 12.4 7.9 15.8 1.4 7.9 11000.0 (benzene conc. taken after WTC attacks)..its just a small chunk of data i pasted for you to look at. its neither normal nor lognormal. someone

Re: [R] 3d scatter plot

Hey Sundar I guess it works indeed. I had made a stupid mistake. Thanks. Eric On Jul 15, 2005, at 11:24 AM, Sundar Dorai-Raj wrote: > Hi, Eric, > > I plotted a surface using the first two columns as x and y and the > third column for z: > > library(lattice) > v <- read.table("clipboard") > name

[R] R on HP-UX

Have any successfully compiled R on HP-UX 11iv2? I was able to compile R 2.0.1 successfully but kept getting memory fault error when I ran the gmake check: gmake check gmake[1]: Entering directory `/opt/R-2.0.1/tests' gmake[2]: Entering directory `/opt/R-2.0.1/tests' gmake[3]: Entering directory `

Re: [R] Can't get sample function from "An Introduction to R" to work

"Bret Collier" <[EMAIL PROTECTED]> writes: > David, > If below is exactly what you typed, check your code again, I think you > are missing a '}' after the last 2 parentheses. That's not supposed to cause a syntax error, just another '+'. I can copy and paste the code as written and not get an err

[R] rmpi in windows

Hi Folks, Has anyone been able to get rmpi to work under windows ? Thanks, Tolga Please follow the attached hyperlink to an important disclaimer == Please access

Re: [R] Can't get sample function from "An Introduction to R" to work

Also, are you doing this with a version of R that will allow you to keep your scripts in separate files and run them one line at a time? Rgui allows this, and many people use different editors, e.g, XEmacs with ESS ("EMacs Speaks Statistics"). Some editors (like ESS) help with parent

Re: [R] Newbie guide for plot & graphics functions?

Significant goof in my first message. Here is the corrected version: x <- 1:10 y <- 1:15 z <- outer(x,y) image(x,y,z) xb <- c(.5,x+.5) yb <- c(.5,y+.5) image(xb,yb,z) ## xb, ybNOTx, y Fortunately, corrected one also succeeds. -Don At 8:06 AM -0700 7/15/05, Don MacQueen wro

Re: [R] Can't get sample function from "An Introduction to R" to work

On Fri, 2005-07-15 at 15:40 -0500, David Groos wrote: > I'm trying to figure out R, a piece at a time, hours at a time... I > was trying to copy the sample function in, "An Introduction to R" (for > version 2.1.0) by W. N. Venables, D. M. Smith, page 42. Section 10.1 > "Simple examples" provi

Re: [R] Can't get sample function from "An Introduction to R" to work

David, If below is exactly what you typed, check your code again, I think you are missing a '}' after the last 2 parentheses. HTH, Bret >>> David Groos <[EMAIL PROTECTED]> 7/15/2005 3:40:01 PM >>> I'm trying to figure out R, a piece at a time, hours at a time... I was trying to copy the sample

[R] Can't get sample function from "An Introduction to R" to work

I'm trying to figure out R, a piece at a time, hours at a time... I was trying to copy the sample function in, "An Introduction to R" (for version 2.1.0) by W. N. Venables, D. M. Smith, page 42. Section 10.1 "Simple examples" provides a sample function which I tried to duplicate (I'm using M

Re: [R] pdf() black&white; panel titles in lattice

On 7/15/05, Chuck Cleland <[EMAIL PROTECTED]> wrote: > [EMAIL PROTECTED] wrote: > > Thank you so much for your help. Now I think I can use R for > > everything I need so far. However I think that I haven't found > > a solution for my first question. trellis.device(color=F) only > > affects the tre

Re: [R] pdf() black&white; panel titles in lattice

[EMAIL PROTECTED] wrote: > Thank you so much for your help. Now I think I can use R for > everything I need so far. However I think that I haven't found > a solution for my first question. trellis.device(color=F) only > affects the trellis device, but not the pdf device, doesn´t > it? At least that

[R] Cannot update some packages after upgrade to 2.1.1

I just upgraded to version 2.1.1 (from 2.0.1) today. > R.version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major2 minor1.1 year 2005 month06 day 20 language R I am using SuSE 9.2 and did the upgrade using rpm -U w

Re: [R] 3d scatter plot

Eric Hu wrote: > Hi, I ran into a dilemma trying to plot the following data in a 3d > scatter fashion. My data are not always increasing as persp() expects. > For now I use scatterplot3d to get points in a 3d scatter plot. I > wonder if I have any way to plot the surfaces. Thanks! > > Eric >

Re: [R] pdf() black&white; panel titles in lattice

Thank you so much for your help. Now I think I can use R for everything I need so far. However I think that I haven't found a solution for my first question. trellis.device(color=F) only affects the trellis device, but not the pdf device, doesn´t it? At least that is my impression. I manage to get

[R] 3d scatter plot

Hi, I ran into a dilemma trying to plot the following data in a 3d scatter fashion. My data are not always increasing as persp() expects. For now I use scatterplot3d to get points in a 3d scatter plot. I wonder if I have any way to plot the surfaces. Thanks! Eric sample data -0.50 0.40 1

[R] 2D contour predictions

Hi All I have been fitting regression models and would now like to produce some contour & image plots from the predictors. Is there an easy way to do this? My current (newbie) experience with R would suggest there is but that it's not always easy to find it! f3 <- lm( fc ~ poly( speed, 2 ) +

Re: [R] pdf() black&white; panel titles in lattice

On 7/15/05, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > Hello, > > I have two questions: > > 1. How can I get black and white graphics using the pdf() > device? I have tried pdf(Name.pdf, bg="white", fg"black"), but > I still get full-color graphics. See ?trellis.device, specifically the 'co

Re: [R] Passing character strings from C code to R

Why invent your own function to do the same thing as the builtin one mkString? Please discuss C programming questions on the R-devel list *as the posting guide asks* On Fri, 15 Jul 2005, mkondrin wrote: > Tuszynski, Jaroslaw W. wrote: > >> Hi, >> >> I have a C code which produces array of inte

[R] [R-pkgs] new package: elliptic

Hi I recently uploaded package "elliptic" to CRAN. From the DESCRIPTION file, it is: A suite of elliptic and related functions including Weierstrass and Jacobi forms. Also includes various tools for manipulating and visualizing complex function

[R] New RVTK package is available for download and testing.

Hello! I have just completed R-language binding for VTK library from Kitware. The library is used for 3D modelling and data visualization. RVTK allows you to do that sort of things in R-session. Because RVTK is in quite an alpha stage it is not yet on CRAN and can be downloaded from my homepage: ww

[R] nlme and spatially correlated errors

Dear R users, I am using lme and nlme to account for spatially correlated errors as random effects. My basic question is about being able to correct F, p, R2 and parameters of models that do not take into account the nature of such errors using gls, glm or nlm and replace them for new F, p, R2

Re: [R] glm(family=binomial(link=logit))

Robin Hankin <[EMAIL PROTECTED]> writes: > Hi > > I am trying to make glm() work to analyze a toy logit system. > > I have a dataframe with x and y independent variables. I have > > L=1+x-y (ie coefficients 1,1,-1) > > then if I have a logit relation with L=log(p/(1-p)), > p=1/(1+exp

Re: [R] Passing character strings from C code to R

If you want to install comment attribute on return value set it this way setAttrib(ans, R_CommentSymbol, charToSexp("Some comment")); __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Passing character strings from C code to R

Tuszynski, Jaroslaw W. wrote: >Hi, > >I have a C code which produces array of integers and potentially a string, >and I have problems passing the string out. Here is the relevant part of the >code: > > 1 PROTECT(Ret = allocVector(INTSXP, n)); > 2 ret = (int*) INTEGER(Ret); /* get pointer

[R] New RVTK package is available for download and testing.

Hello! I have just completed R language binding for VTK library from Kitware. The library is used for 3D modelling and data visualization. RVTK allows you to do that sort of things in R-session. Because RVTK is in quite an alpha stage it is not yet on CRAN and can be downloaded from my homepag

Re: [R] glm(family=binomial(link=logit))

On Fri, 15 Jul 2005, Robin Hankin wrote: > I am trying to make glm() work to analyze a toy logit system. > > I have a dataframe with x and y independent variables. I have > > L=1+x-y (ie coefficients 1,1,-1) > > then if I have a logit relation with L=log(p/(1-p)), > p=1/(1+exp(L)). Not

[R] pdf() black&white; panel titles in lattice

Hello, I have two questions: 1. How can I get black and white graphics using the pdf() device? I have tried pdf(Name.pdf, bg="white", fg"black"), but I still get full-color graphics. 2. When using lattice, how do I change the title for each panel? By default it uses the categories of the grou

Re: [R] Newbie guide for plot & graphics functions?

I think you have it backwards (the error message not withstanding). Read the documentation for image, note the beginning of the Details section, and then compare these two examples, both of which succeed on my system (R 2.1.1). x <- 1:10 y <- 1:15 z <- outer(x,y) image(x,y,z) xb <- c(.5,x+.5)

[R] glm(family=binomial(link=logit))

Hi I am trying to make glm() work to analyze a toy logit system. I have a dataframe with x and y independent variables. I have L=1+x-y (ie coefficients 1,1,-1) then if I have a logit relation with L=log(p/(1-p)), p=1/(1+exp(L)). If I interpret "p" as the probability of success in a

Re: [R] Padding in lattice plots

On 7/15/05, Federico Gherardini <[EMAIL PROTECTED]> wrote: > On Friday 15 July 2005 14:42, you wrote: > > Hi all, > > I've used the split argument to print four lattice plots on a single page. > > The problem now is that I need to reduce the amount of white space between > > the plots. I've read ot

Re: [R] Newbie guide for plot & graphics functions?

If you have a little bit of time, check out Ross Ihaka's Lectures on "Information Visualisation". The specific examples he gives are built with R: http://www.stat.auckland.ac.nz/~ihaka/120/ Time could be required because this stuff is so great that you want to read all of it... Christophe. P

Re: [R] Coxph with factors

On Fri, 15 Jul 2005, Kylie-Anne Richards wrote: > > FIRST Q: The default uses the mean of 'vo' and mean of 'po', but what is it > using for the factors?? Is it the sum of the coef of the factors divided by > the number of factors?? It uses the mean of each factor variable. The $means component of

Re: [R] Newbie guide for plot & graphics functions?

Pete Cap wrote: > Hello all, > > Can anyone point me to a decent introduction to using the plotting and > assorted graphics functions in R? > > I keep getting simple errrors and I can't figure out why, for example: > > >>image(x,y,z) > > Error in image.default(x, y, z) : dimensions of z ar

[R] Newbie guide for plot & graphics functions?

Hello all, Can anyone point me to a decent introduction to using the plotting and assorted graphics functions in R? I keep getting simple errrors and I can't figure out why, for example: > image(x,y,z) Error in image.default(x, y, z) : dimensions of z are not length(x)(+1) times length(y)(+1

Re: [R] Dividing a vector into ntiles

Here is an example that shows how to divide a vector into "quartiles" and create a categorical (factor) variable with 4 levels ("A", "B", "C", "D") from it: x <- rnorm(100) xcat <- factor(cut(x, quantile(x), include.lowest = TRUE), labels = LETTERS[1:4]) table(xcat) Ravi. > --

[R] Coxph with factors

My apologies for lack of clarity. Sample data: time.sec po vo done f.pom 55 0 4402 1 -0.25 378 -0.5 5 1 -0.25 51 0 1000 1 -0.25 43 0 71581 1 -0.25 etc >DATA<-data.frame(time.sec,done,f.pom=factor(f.pom),po,vo) >final<-coxph(Surv(time.sec,done)~f

Re: [R] crossed random fx nlme lme4

On 7/13/05, Simon Blomberg <[EMAIL PROTECTED]> wrote: > At 09:35 AM 14/07/2005, Emilio A. Laca wrote: > >I need to specify a model similar to this > > > >lme.formula(fixed = sqrt(lbPerAc) ~ y + season + y:season, data = cy, > > random = ~y | observer/set, correlation = corARMA(q = 6)) > > > >e

Re: [R] Partek has Dunn-Sidak Multiple Test Correction. Is this the same/similar to any of R's p.adjust.methods?

> "Earl" == Earl F Glynn <[EMAIL PROTECTED]> > on Thu, 14 Jul 2005 12:22:49 -0500 writes: Earl> The Partek package (www.partek.com) allows only two selections for Multiple Earl> Test Correction: Bonferroni and Dunn-Sidak. Can anyone suggest why Partek Earl> implemented

[R] Ordinal data - Regression Trees & Proportional Odds

Dear Dr. Fieberg, you used a regression tree approach to explore ordinal data set in addition to the proportinal odds model. I find this very interesting. I would like to know, how good the results of the regression tree approach turned out in comparison to the proportional odds model. Since p

Re: [R] unexpected par('pin') behaviour

> "MM" == Martin Maechler <[EMAIL PROTECTED]> > on Fri, 15 Jul 2005 14:36:16 +0200 writes: MM> which I can simplify to MM> graphics.off() MM> layout(matrix(1:42,7,6)) MM> par("pin") MM> ## [1] 0.2918620 -0.2974408 --- when using x11() "default"

Re: [R] Dividing a vector into ntiles

You could try using "cut": > cut(1:9, breaks = 3) [1] (0.992,3.66] (0.992,3.66] (0.992,3.66] (3.66,6.34] (3.66,6.34] [6] (3.66,6.34] (6.34,9.01] (6.34,9.01] (6.34,9.01] Levels: (0.992,3.66] (3.66,6.34] (6.34,9.01] > Quantile is another possible solution to get break points: > quantile(1:9, p

[R] Passing character strings from C code to R

Hi, I have a C code which produces array of integers and potentially a string, and I have problems passing the string out. Here is the relevant part of the code: 1 PROTECT(Ret = allocVector(INTSXP, n)); 2 ret = (int*) INTEGER(Ret); /* get pointer to R's Ret */ 3 for(i=0; i0) {//

Re: [R] Dividing a vector into ntiles

?quantile quantile(runif(1000), probs=seq(0,1,1/3)) John Sorkin wrote: > R 2.1.1 > Win 2k > > Would someone suggest a method (or methods) that can be used to > determine ntile cutpoints of a vector, i.e. to determine values that can > be used to divide a vector into thirds such as 0-33 centile,

Re: [R] Fwd: Re: East Asian language

I'm not sure what other applications use environment variable like "LANGUAGE". I'm sorry I was so assertive. However, what I said is right in the sense that environment variables are supposed to be used NOT by one specific application. I happen to use lots of GnuWin32 and Cygwin applications such a

[R] Dividing a vector into ntiles

R 2.1.1 Win 2k Would someone suggest a method (or methods) that can be used to determine ntile cutpoints of a vector, i.e. to determine values that can be used to divide a vector into thirds such as 0-33 centile, 34-66 centile, 67-100 centile. If for example I had a vector: 1,2,3,4,5,6,7,8,9 and

Re: [R] unexpected par('pin') behaviour

Thank you, Joerg, for the reproducable example > "joerg" == joerg van den hoff <[EMAIL PROTECTED]> > on Thu, 14 Jul 2005 17:19:51 +0200 writes: .. .. joerg> ===cut joerg> graphics.off() joerg> f <-

Re: [R] Fwd: Re: East Asian language

On Fri, 15 Jul 2005, Uwe Ligges wrote: > Pingping Zheng wrote: > >> I think what Nan Lin wanted is "How to change the language used in R >> console by re-setting something within R?" >> >> For a non-english user, sometime we prefer to set the system default >> language to a native one and use engl

Re: [R] Fwd: Re: East Asian language

On Fri, 15 Jul 2005, Peter Dalgaard wrote: > Gabor Grothendieck <[EMAIL PROTECTED]> writes: > >> I am not sure how generally this works but I was able to get it to switch >> back and forth between two languages within a single R >> session like this: The posting guide does ask people to read the

Re: [R] Another simple q - removing negative values

To leave x with only the non-negative elements, you can use x[x >= 0]. Also see the function "subset". Hope this helps, Matt Wiener -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Michael Hopkins Sent: Friday, July 15, 2005 8:00 AM To: r-help@stat.math.

Re: [R] Another simple q - removing negative values

Michael Hopkins <[EMAIL PROTECTED]> writes: > Sorry folks but again I have failed in my understanding of how to do a very > simple thing. I've read the various texts and searched the help archives > but no positive result so far. > > I want to remove all the rows in a data frame where one of the

Re: [R] Padding in lattice plots

On Friday 15 July 2005 14:42, you wrote: > Hi all, > I've used the split argument to print four lattice plots on a single page. > The problem now is that I need to reduce the amount of white space between > the plots. I've read other mails in this list about the new trellis > parameters layout.heig

[R] Another simple q - removing negative values

Sorry folks but again I have failed in my understanding of how to do a very simple thing. I've read the various texts and searched the help archives but no positive result so far. I want to remove all the rows in a data frame where one of the variables has negative values. In approx Stata pseud

Re: [R] PLS: problem transforming scores to variable space

rainer grohmann writes: > However, when I try to map the scores back to variable space, I ran into > problems: [...] > cbind(t$scores[,1],(t$scores%*%(t$loadings)%*%t$projection)[,1]) You need to transpose the loadings: > all.equal(unclass(t$scores), + t$scores %*% t(t$loadings) %*% t$

Re: [R] Fwd: Re: East Asian language

Gabor Grothendieck <[EMAIL PROTECTED]> writes: > I am not sure how generally this works but I was able to get it to switch > back and forth between two languages within a single R > session like this: > > > Sys.getenv("LANGUAGE") > LANGUAGE > "" > > q() # note English prompt; c means conti

Re: [R] Fwd: Re: East Asian language

I am not sure how generally this works but I was able to get it to switch back and forth between two languages within a single R session like this: > Sys.getenv("LANGUAGE") LANGUAGE "" > q() # note English prompt; c means continue rather than exiting Save workspace image? [y/n/c]: c > Sys.p

[R] Padding in lattice plots

Hi all, I've used the split argument to print four lattice plots on a single page. The problem now is that I need to reduce the amount of white space between the plots. I've read other mails in this list about the new trellis parameters layout.heights and layout.widhts but I haven't been able to

Re: [R] Calculation of group summaries

Perhaps I lost track of what the original question was, but on my homepage http://genetics.agrsci.dk/~sorenh/misc/ there is a package called doBy in which there is a function called summaryBy (which mimics proc summary from sas). For example summaryBy(cbind(Weight,Feed)~Evit+Cu, data=di

Re: [R] Fwd: Re: East Asian language

That's pretty convenient. Set LANGUAGE in .Renviron under current directory is convenient for both Linux and windows users. Thanks. Pingping Zheng Department of Mathematics and Statistics Fylde College Lancaster University Lancaster LA1 4YF UK Uwe Ligges wrote: > Pingping Zheng wrote: > > >>D

Re: [R] exact values for p-values

> there is no such thing as an "R lib". > Do you mean "R package"? Yes. I usually refer to them R packages as as libraries, hence the use of lib. Regards. Stephen. > In any case, not it is not (yet). > As David clearly says (a few lines below), it is a 'f90' > library; which means it's "fo

Re: [R] Fwd: Re: East Asian language

Pingping Zheng wrote: > Dear all, > > I think what Nan Lin wanted is "How to change the language used in R > console by re-setting something within R?" > > For a non-english user, sometime we prefer to set the system default > language to a native one and use english in individual programs when

Re: [R] Fwd: Re: East Asian language

Dear all, I think what Nan Lin wanted is "How to change the language used in R console by re-setting something within R?" For a non-english user, sometime we prefer to set the system default language to a native one and use english in individual programs when we think necessary. a) setting the e

Re: [R] exact values for p-values

> "S" == S O Nyangoma <[EMAIL PROTECTED]> > on Wed, 13 Jul 2005 11:06:05 +0200 writes: S> Hi David, Since I am looking at very extreme values, it appears I will S> need FMLIB. Is it an R lib? there is no such thing as an "R lib". Do you mean "R package"? In any case, not

Re: [R] Coxph with factors

Yes, and please show us a reproducible example or small section of the data as well as the error output. On Thu, 2005-07-14 at 15:36 +, Dieter Menne wrote: > Kylie-Anne Richards itga.com.au> writes: > > > > > I am fitting a coxph model with factors. I am running into problems when > > usi

Re: [R] Variance components from lm?

In any case, the question was about lm, not lme. The short answer is that this is easiest from aov, a wrapper for lm, and for balanced data will be equivalent to as using lme. Steve can find worked examples for interlaboratory data in MASS (the book, see the FAQ). On Thu, 14 Jul 2005, Spencer