Re: [R] Help with "For" instruction

2011-01-03 Thread David Winsemius
On Jan 4, 2011, at 1:00 AM, ADias wrote: Hi, I am having a problem in doing something similar to this example: Suppose I have this vector a, and from it I wish to create 5 other vector each one with less one value than what object a has So I have "a" a<-c(1,2,3,4,5) and I want a1 that

Re: [R] Print plot to pdf, jpg or any other format when using scatter3d error

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 8:17 PM, Jurica Seva wrote: Hi, I have been trying to output my graphs to a file (jpeg, pdf, ps, it doesnt matter) but i cant seem to be able to get it to output. I tried a few things but none of them worked and am lost as what to do now. I am using the scatter3d funct

Re: [R] Help with "For" instruction

2011-01-03 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of ADias > Sent: Monday, January 03, 2011 10:00 PM > To: r-help@r-project.org > Subject: [R] Help with "For" instruction > > > Hi, > > I am having a problem in doing something simi

Re: [R] Regex to remove last character

2011-01-03 Thread Martin Morgan
On 01/03/2011 02:44 PM, David Winsemius wrote: > > On Jan 3, 2011, at 2:24 PM, rivercode wrote: > >> >> Hi, >> >> Have been having trouble trying to figure out the right regex >> parameters to >> remove the last "." in timestamp with the following format: >> >> Convert 09:30:00.377.853 to 09:30:

[R] Print plot to pdf, jpg or any other format when using scatter3d error

2011-01-03 Thread Jurica Seva
Hi, I have been trying to output my graphs to a file (jpeg, pdf, ps, it doesnt matter) but i cant seem to be able to get it to output. I tried a few things but none of them worked and am lost as what to do now. I am using the scatter3d function, and it prints out the graphs on tot he screen w

[R] Help with "For" instruction

2011-01-03 Thread ADias
Hi, I am having a problem in doing something similar to this example: Suppose I have this vector a, and from it I wish to create 5 other vector each one with less one value than what object a has So I have "a" a<-c(1,2,3,4,5) and I want a1 that shoud have (2,3,4,5) a2 that should have (1,3,4,

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread luke-tierney
On Mon, 3 Jan 2011, Frank Harrell wrote: Correct. I'm doing this because of non-exported functions in other packages, so I need ::: I'd still appreciate any insight about whether importFrom in NAMESPACE defers package loading so that if the package is not actually used (and is not installed)

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Martin Morgan
On 01/03/2011 09:06 PM, Hadley Wickham wrote: >> Correct. I'm doing this because of non-exported functions in other packages, >> so I need ::: > > But you really really shouldn't be doing that. Is there a reason that > the package authors won't export the functions? > >> I'd still appreciate an

Re: [R] matrices call a function element-wise

2011-01-03 Thread Dennis Murphy
Hi: The idea is as follows: * string your four matrices into vectors, cbinding them so that the columns correspond to what you want as the (1,1), (1, 2), (2, 1) and (2, 2) elements, respectively, of the matrix/table to be used for the Fisher test; * Operate row-wise on the constructed matr

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Bill.Venables
If you use ::: to access non-exported functions, as Frank confesses he does, then you can't complain if in the next release of the package involved the non-exported objects are missing and things are being done another way entirely. That's the deal. On the other hand, sometimes package authors

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Hadley Wickham
> Correct.  I'm doing this because of non-exported functions in other packages, > so I need ::: But you really really shouldn't be doing that. Is there a reason that the package authors won't export the functions? > I'd still appreciate any insight about whether importFrom in NAMESPACE > defers

[R] how to subset unique factor combinations from a data frame.

2011-01-03 Thread SNV Krishna
Hi All I have these questions and request members expert view on this. a) I have a dataframe (df) with five factors (identity variables) and value (measured value). The id variables are Year, Country, Commodity, Attribute, Unit. Value is a value for each combination of this. I would like to

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Frank Harrell
Correct. I'm doing this because of non-exported functions in other packages, so I need ::: I'd still appreciate any insight about whether importFrom in NAMESPACE defers package loading so that if the package is not actually used (and is not installed) there will be no problem. Thanks Frank --

Re: [R] matrices call a function element-wise

2011-01-03 Thread Jonathan Christensen
Hi, I would recommend reformatting the data as a 2x2x1000 array and using apply. Jonathan On Mon, Jan 3, 2011 at 7:57 AM, zhaoxing731 wrote: > Hello > > I have 4 1000*1000 matrix A,B,C,D. I want to use the corresponding element of > the 4 matrices. Using the "for loop" as follow: > > E<-o > f

Re: [R] unicode variable and function names?

2011-01-03 Thread Duncan Murdoch
On 11-01-03 7:44 PM, stvienna wiener wrote: Dear List, Is it possible to have function names like ∀ (unicode universal quantifier)? This question is inspired by agda source code, which supports this. http://www.cs.nott.ac.uk/~nad/listings/lib-0.4/Algebra.html It would be handy to use. My gues

Re: [R] unicode variable and function names?

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 7:44 PM, stvienna wiener wrote: Dear List, Is it possible to have function names like ∀ (unicode universal quantifier)? This question is inspired by agda source code, which supports this. http://www.cs.nott.ac.uk/~nad/listings/lib-0.4/Algebra.html It would be handy to

Re: [R] function masking and gmp questions

2011-01-03 Thread Duncan Murdoch
On 11-01-03 7:44 PM, Carl Witthoft wrote: > The problem is that I can't 'insert' the :: operator inside an existing > function. (Well, I could, but that would mean rewriting the function) > > I was hoping for a way to call the function, in this case turnpoints() > in some way that told turnpoints

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Peter Langfelder
Well, I'm pretty sure that, inside package A, calling B::functionName will not work if B has not been imported. That's why I use ::: (after spending some time trying to figure out why :: didn't work). At least that was the state of affairs as of R 2.9 or so, perhaps things have changed since then.

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Hadley Wickham
>> I think you mean packagename::functionname?  The three colon form is >> for accessing non-exported objects. > > Normally two colons suffice, but within a package you need three to > access exported but un-imported objects :) Are you sure? Note that it is typically a design mistake to use

Re: [R] function masking and gmp questions

2011-01-03 Thread Carl Witthoft
The problem is that I can't 'insert' the :: operator inside an existing function. (Well, I could, but that would mean rewriting the function) I was hoping for a way to call the function, in this case turnpoints() in some way that told turnpoints itself to look for base::apply Carl On 1/3/1

[R] unicode variable and function names?

2011-01-03 Thread stvienna wiener
Dear List, Is it possible to have function names like ∀ (unicode universal quantifier)? This question is inspired by agda source code, which supports this. http://www.cs.nott.ac.uk/~nad/listings/lib-0.4/Algebra.html It would be handy to use. My guess is, however, that it's not supportet in R.

[R] [R-pkgs] orddom package for Ordinal Dominance Statistics

2011-01-03 Thread Jens Rogmann (ABK/SK-Zentrum EPB Uni HH)
Dear R-users, a new R package has been released named orddom. This package provides ordinal estimates as alternatives to independent or paired group mean comparisons, especially for Cliff’s delta statistics. It provides basic parameters for various robust tests of stochastic equality with ord

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Peter Langfelder
On Mon, Jan 3, 2011 at 3:48 PM, Hadley Wickham wrote: > Hi Frank, > > I think you mean packagename::functionname?  The three colon form is > for accessing non-exported objects. Normally two colons suffice, but within a package you need three to access exported but un-imported objects :) Peter _

Re: [R] function masking and gmp questions

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 7:17 PM, Carl Witthoft wrote: Hi, Here's the problem I ran into: the gmp package has a method for apply() so it masks the base::apply function. With gmp installed, I tried to run the function turnpoints() from the pastecs package. It fails because it calls apply() in

Re: [R] how to invert the axes in the wireframe() plot

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 7:07 PM, Simone Gabbriellini wrote: I am trying to reproduce this graph: http://www.digitaldust.it/materiali/their.png the default axes orientation of wireframe gives me this: http://www.digitaldust.it/materiali/mine.pdf I am trying to understand if I can reproduce the ax

[R] function masking and gmp questions

2011-01-03 Thread Carl Witthoft
Hi, Here's the problem I ran into: the gmp package has a method for apply() so it masks the base::apply function. With gmp installed, I tried to run the function turnpoints() from the pastecs package. It fails because it calls apply() internally, like this: apply(mymatrix,1,max,na.rm=TRUE)

Re: [R] how to invert the axes in the wireframe() plot

2011-01-03 Thread Simone Gabbriellini
I am trying to reproduce this graph: http://www.digitaldust.it/materiali/their.png the default axes orientation of wireframe gives me this: http://www.digitaldust.it/materiali/mine.pdf I am trying to understand if I can reproduce the axes orientation of the first figure. many thanks, simone

[R] Resampling to find Confidence intervals

2011-01-03 Thread Ben Ward
Hi, I'm doing some modelling (lm) for my 3rd year dissertation and I want to do some resampling, especially as I'm working with microbes, getting them to evolve resistance to antimicrobial compounds, and after each exposure I'm measuring the minimum concentration required to kill them (which I'

Re: [R] packagename:::functionname vs. importFrom

2011-01-03 Thread Hadley Wickham
Hi Frank, I think you mean packagename::functionname? The three colon form is for accessing non-exported objects. Otherwise, I think using :: vs importFrom is functionally identical - either approach delays package loading until necessary. Hadley On Mon, Jan 3, 2011 at 9:45 PM, Frank Harrell

Re: [R] Logistic Regression Fitting with EM-Algorithm

2011-01-03 Thread Ted Harding
On 03-Jan-11 14:02:21, Robin Aly wrote: > Hi all, > is there any package which can do an EM algorithm fitting of > logistic regression coefficients given only the explanatory > variables? I tried to realize this using the Design package, > but I didn't find a way. > > Thanks a lot & Kind regards >

Re: [R] Regex to remove last character

2011-01-03 Thread Henrique Dallazuanna
Try this: format(strptime("09:30:00.377.853", "%T%OS"), "%H:%M:%OS") On Mon, Jan 3, 2011 at 5:24 PM, rivercode wrote: > > Hi, > > Have been having trouble trying to figure out the right regex parameters to > remove the last "." in timestamp with the following format: > > Convert 09:30:00.377.85

Re: [R] Inverse Gaussian Distribution

2011-01-03 Thread Dennis Murphy
Hi: "How can i manipulate the data to set the baseline of area to C? R is producing errors when I'm trying to do so." See ?relevel Dennis On Mon, Jan 3, 2011 at 12:03 PM, Louisa wrote: > > Dear, > > I want to fit an inverse gaussion distribution to a data set. > > The predictor variables are

Re: [R] randomForest speed improvements

2011-01-03 Thread apresley
I haven't tried changing the mtry or ntree at all ... though I suppose with only 6 variables, and tens-of-thousands of rows, we can probably do less than 500 tree's (the default?). Although tossing the forest does speed things up a bit, seems to be about 15 - 20% faster in some cases, I need to k

Re: [R] how to invert the axes in the wireframe() plot

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 5:06 PM, Simone Gabbriellini wrote: Dear List, I am using the wireframe function in the lattice package, and I am wondering if it is possible to invert the default axes orientation for x and y axes... what parameter should I look for? Perhaps you should define "invert"

Re: [R] Regex to remove last character

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 2:24 PM, rivercode wrote: Hi, Have been having trouble trying to figure out the right regex parameters to remove the last "." in timestamp with the following format: Convert 09:30:00.377.853 to 09:30:00.377853 gsub("()(\\.)(.{3}$)", "\\1\\3" , "09:30:00.377.853") [1]

[R] [R-pkgs] dclone 1.3-0

2011-01-03 Thread Peter Solymos
Dear R Community, I am happy to introduce the latest version 1.3-0 of the 'dclone' R package. The package provides low level functions for implementing maximum likelihood estimating procedures for complex models using data cloning and Bayesian Markov chain Monte Carlo methods with support for JAG

Re: [R] Greetings. I have a question with mixed beta regr ession model in nlme.

2011-01-03 Thread Ben Bolker
elciclopeturnio gmail.com> writes: > > *Dear R-help: > > My name is Rodrigo and I have a question with nlme package > in R to fit a mixed beta regression model. The details of the model are: > [snip] > *The question is: > How can I use nlme package in R to fit this model? > If you want to

Re: [R] iPhone 3G App For R?

2011-01-03 Thread David Duffy
Roy Mendelssohn wrote: On Jan 2, 2011, at 6:28 PM, Shige Song wrote: Does iphone even support the GNU tool chain? Check the archive for r-sig-mac (https://stat.ethz.ch/pipermail/r-sig-mac/). There has been extensive discussion about this. If memory serves (and it rarely does anymore :-) )

[R] Regex to remove last character

2011-01-03 Thread rivercode
Hi, Have been having trouble trying to figure out the right regex parameters to remove the last "." in timestamp with the following format: Convert 09:30:00.377.853 to 09:30:00.377853 Thanks, Chris -- View this message in context: http://r.789695.n4.nabble.com/Regex-to-remove-last-charact

Re: [R] Problem with uploading library

2011-01-03 Thread gaja
Woho. I did it... Thanx guys ;) -- View this message in context: http://r.789695.n4.nabble.com/Problem-with-uploading-library-tp3170455p3172572.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://

Re: [R] ARIMA simulation including a constant

2011-01-03 Thread Paolo Rossi
Sorry I am not really sure I have been taht clear. I meant ARMA which is not bound to have zero mean. More precisely, suppose I estimate y(t) = a + by(t-1) + e(t) + ce(t-1) , i.e. and ARMA(1,1). My question is how do I simulate values for yt given the values for a, b and c? My problem with arima.

[R] how to invert the axes in the wireframe() plot

2011-01-03 Thread Simone Gabbriellini
Dear List, I am using the wireframe function in the lattice package, and I am wondering if it is possible to invert the default axes orientation for x and y axes... what parameter should I look for? Best regards, Simone Gabbriellini __ R-help@r-projec

Re: [R] Saving objects inside a list

2011-01-03 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Eduardo de > Oliveira Horta > Sent: Monday, January 03, 2011 12:58 PM > To: r-help > Subject: Re: [R] Saving objects inside a list > > sapply(ls(),get) works fine. Thanks. > >

[R] packagename:::functionname vs. importFrom

2011-01-03 Thread Frank Harrell
In my rms package I use the packagename:::functionname construct in a number of places. If I instead use the importFrom declaration in the NAMESPACE file would that require the package to be available, and does it load the package when my package loads? If so I would keep using packagename::: to

Re: [R] Saving objects inside a list

2011-01-03 Thread Gabor Grothendieck
On Mon, Jan 3, 2011 at 3:58 PM, Eduardo de Oliveira Horta wrote: > sapply(ls(),get) works fine. Thanks. > > ps: the as.list and the eapply suggestions didn't work. > They work for me. Starting in a fresh session: > x <- 1; f <- function(x) x; DF <- data.frame(a = 1:3) > as.list(.GlobalEnv) $DF

Re: [R] Saving objects inside a list

2011-01-03 Thread Eduardo de Oliveira Horta
sapply(ls(),get) works fine. Thanks. ps: the as.list and the eapply suggestions didn't work. On Mon, Jan 3, 2011 at 3:56 PM, Jorge Ivan Velez wrote: > Hi Eduardo, > > Try > > r <- ls() > result <- sapply(r, get) > result > > HTH, > Jorge > > > On Mon, Jan 3, 2011 at 12:25 PM, Eduardo de Oliveira

Re: [R] Formatted output with alternating format at different rows

2011-01-03 Thread jim holtman
'sprintf' if your friend: > dummy3 = c(1.1, 2.2, 3.3) > dummy4 = c(4.4, 5.5, 6.6, 7.7) > dummy2 = c(8.8, 9.9) > > cat(sprintf("%5.1f%6.2f%7.3f\n", dummy3[1], dummy3[2], dummy3[3])) 1.1 2.20 3.300 > cat(sprintf("%5.2f%6.3f%7.4f%8.5f\n", dummy4[1], dummy4[2], dummy4[3], > dummy4[4])) 4.40 5.50

Re: [R] Inverse Gaussian Distribution

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 3:03 PM, Louisa wrote: Dear, I want to fit an inverse gaussion distribution to a data set. The predictor variables are gender, area and agecategory. For each of these variables I've defined a baseline e.g. #agecat: baseline is 3 data<-transform(data, agecat=C(factor(ageca

Re: [R] randomForest speed improvements

2011-01-03 Thread Jonathan P Daily
Have you tried adjusting: mtry - the number of parameters to try per tree ntree - the number of trees grown keep.forest - logical on whether to store tree Specifically, I found huge improvements in speed by switching keep.forest to FALSE in the past when I didn't actually need the forest post ana

[R] Inverse Gaussian Distribution

2011-01-03 Thread Louisa
Dear, I want to fit an inverse gaussion distribution to a data set. The predictor variables are gender, area and agecategory. For each of these variables I've defined a baseline e.g. #agecat: baseline is 3 data<-transform(data, agecat=C(factor(agecat,ordered=TRUE), contr.treatment(n=6,base

[R] randomForest speed improvements

2011-01-03 Thread apresley
Hi there, We're trying to use randomForest to do some predictions. The test-harness for our code is pretty straightforward: library ('randomForest'); data202 <- read.csv ("random.csv", header=TRUE); x<- data202[1:5,1:6]; y<- data202[1:5,8]; y<- y[,drop=TRUE]; x2 <- data202[

[R] Problem loading Tcl/Tk (?)

2011-01-03 Thread Peter Langfelder
Hello, a user of one of my packages has sent me the following problem that apparently appears when loading the tcltk package. I'm guessing something (Tcl/Tk? R open scripting framework?) is not installed correctly... can someone please confirm or correct my guess? The OS is Mac 10.6.5. Thanks, P

Re: [R] matrices call a function element-wise

2011-01-03 Thread Jonathan P Daily
IIRC, R is perfectly able to call matrices as vectors, so you might be able to do this: FT <- function(i) fisher.test(matrix(c(A[i],B[i],C[i],D[i]),2)) E <- sapply(1:100, FT) Though I don't know how much time you will save. -- Jonathan P. Daily Technician

[R] Formatted output with alternating format at different rows

2011-01-03 Thread Ray Laymon
Dear all, I have a simple question. I couldn't find a solution in the forums/R-user manual; I have also asked to my friends who use R, but couldn't get any answer from them either. I would appreciate any solutions. I want to write formatted text file like in Fortran. More specifically with the fo

[R] matrices call a function element-wise

2011-01-03 Thread zhaoxing731
Hello I have 4 1000*1000 matrix A,B,C,D. I want to use the corresponding element of the 4 matrices. Using the "for loop" as follow: E<-o for (i in 1:1000) {for (j in 1:1000) { E<-fisher.test(matrix(c(A[i][j],B[i][j],C[i][j],D[i][j]),2))#call fisher.test for every element } }

[R] Logistic Regression Fitting with EM-Algorithm

2011-01-03 Thread Robin Aly
Hi all, is there any package which can do an EM algorithm fitting of logistic regression coefficients given only the explanatory variables? I tried to realize this using the Design package, but I didn't find a way. Thanks a lot & Kind regards Robin Aly __

[R] Greetings. I have a question with mixed beta regression model in nlme (corrected version).

2011-01-03 Thread elciclopeturnio
*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. I'm so sorry. In the last email, I forgot to say that W is also a unknown parameter in the mixed beta regression model. In any case, here I send you the correct formulation. ** Supp

Re: [R] Please, need help with a plot

2011-01-03 Thread Victor F Seabra
although the code somehow didn't work on my Vista / R 2.8, it did work perfectly on a XP machine / R 2.10 I've been trying to fix this for days, Thank you very much for your help! __ 02/01/2011 19:30, David Winsemius < dwinsem

[R] Greetings. I have a question with mixed beta regression model in nlme.

2011-01-03 Thread elciclopeturnio
*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. The details of the model are: Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~ Beta(mu_{ij} * phi_{ij}; (1 - mu_{ij}) * phi_{ij}) y_{ij} =

Re: [R] changing method of estimation in GLM

2011-01-03 Thread ATANU
can you please explain me this with an example(preferably NEWTON RAPHSON METHOD)? -- View this message in context: http://r.789695.n4.nabble.com/changing-method-of-estimation-in-GLM-tp3170836p3172088.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Proof for computing sums of squares

2011-01-03 Thread Douglas Bates
On Mon, Jan 3, 2011 at 12:15 PM, Ethan Arenson wrote: > Hi. > > I know that R computes sums of squares based on the diagonal of > > t(Q) %*% y %*% t(y) %*% Q, > > where Q comes from the QR-decomposition of the model matrix. So, how do you know that? When I look at the code for summary.lm I see t

[R] Proof for computing sums of squares

2011-01-03 Thread Ethan Arenson
Hi. I know that R computes sums of squares based on the diagonal of t(Q) %*% y %*% t(y) %*% Q, where Q comes from the QR-decomposition of the model matrix. Does anyone know where I can find a proof for this result? All Best and Happy New Year, Ethan [[alternative HTML version deleted]

Re: [R] Saving objects inside a list

2011-01-03 Thread Jorge Ivan Velez
Hi Eduardo, Try r <- ls() result <- sapply(r, get) result HTH, Jorge On Mon, Jan 3, 2011 at 12:25 PM, Eduardo de Oliveira Horta <> wrote: > Hello there, > > any ideas on how to save all the objects on my workspace inside a list > object? > > For example, say my workspace is as follows > ls()

Re: [R] Saving objects inside a list

2011-01-03 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Eduardo de > Oliveira Horta > Sent: Monday, January 03, 2011 9:25 AM > To: r-help > Subject: [R] Saving objects inside a list > > Hello there, > > any ideas on how to save all

Re: [R] Saving objects inside a list

2011-01-03 Thread David L Lorenz
Eduardo, Try this: object.list <- sapply(objects(), function(x) get(x), simplify=F, USE.NAMES=T) Dave From: Eduardo de Oliveira Horta To: r-help Date: 01/03/2011 11:32 AM Subject: [R] Saving objects inside a list Sent by: r-help-boun...@r-project.org Hello there, any ideas on how to sa

Re: [R] Saving objects inside a list

2011-01-03 Thread Henrique Dallazuanna
Remove the pattern argument: sapply(ls(), get, simplify = F) On Mon, Jan 3, 2011 at 3:51 PM, Eduardo de Oliveira Horta < eduardo.oliveiraho...@gmail.com> wrote: > Thanks, but this is still too restrictive... I wanted something that could > actually save objects with arbitrary names inside a lis

Re: [R] Saving objects inside a list

2011-01-03 Thread Gabor Grothendieck
On Mon, Jan 3, 2011 at 12:25 PM, Eduardo de Oliveira Horta wrote: > Hello there, > > any ideas on how to save all the objects on my workspace inside a list > object? > > For example, say my workspace is as follows > ls() > [1] "x" "y" "z" > > and suppose I want to put these objects inside a list o

Re: [R] Saving objects inside a list

2011-01-03 Thread Eduardo de Oliveira Horta
Thanks, but this is still too restrictive... I wanted something that could actually save objects with arbitrary names inside a list object... On Mon, Jan 3, 2011 at 3:39 PM, Henrique Dallazuanna wrote: > Try this: > > sapply(ls(pattern = "x|y|z"), get, simplify = F) > > > On Mon, Jan 3, 2011 at 3

Re: [R] Saving objects inside a list

2011-01-03 Thread Henrique Dallazuanna
Try this: sapply(ls(pattern = "x|y|z"), get, simplify = F) On Mon, Jan 3, 2011 at 3:25 PM, Eduardo de Oliveira Horta < eduardo.oliveiraho...@gmail.com> wrote: > Hello there, > > any ideas on how to save all the objects on my workspace inside a list > object? > > For example, say my workspace is

[R] Saving objects inside a list

2011-01-03 Thread Eduardo de Oliveira Horta
Hello there, any ideas on how to save all the objects on my workspace inside a list object? For example, say my workspace is as follows ls() [1] "x" "y" "z" and suppose I want to put these objects inside a list object, say object.list <- list() without having to explicitly write down their nam

Re: [R] parLapply - Error in do.call("fun", lapply(args, enquote)) : could not find function "fun"

2011-01-03 Thread Adrienne Wootten
Alright, I found the solution on my own at least it seems that way. Here is what I discovered for future reference in the R-help archive. The issue with parLapply comes in with the docall function in parLapply, > parLapply function(cl, x, fun, ...){ docall(c, clusterApply(cl, splitList(x, lengt

[R] ARIMA simulation including a constant

2011-01-03 Thread Paolo Rossi
Hi, I have been looking at arima.sim to simulate the output from an ARMA model fed with a normal and uncorrelated input series but I cannot find a way to pass an intercept / constant into the model. In other words, the model input in the function allows only for the AR and MA components but I need

Re: [R] optimize

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 11:09 AM, Muhammad Rahiz wrote: Josh's recommendation to use predict works. At the same time, I'll work on your suggestions, David. Thanks. Muhammad Rahiz Researcher & DPhil Candidate (Climate Systems & Policy) School of Geography & the Environment University of Oxford O

Re: [R] optimize

2011-01-03 Thread Muhammad Rahiz
Josh's recommendation to use predict works. At the same time, I'll work on your suggestions, David. Thanks. Muhammad Rahiz Researcher & DPhil Candidate (Climate Systems & Policy) School of Geography & the Environment University of Oxford On Mon, 3 Jan 2011, David Winsemius wrote: On Jan 3,

Re: [R] optimize

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 9:52 AM, Muhammad Rahiz wrote: Hi all, I'm trying to get the value of y when x=203 by using the intersect of three curves. The horizontal curve does not meet with the other two. How can I rectify the code below? Thanks Muhammad ts <- 1:10 dd <- 10:1 ts <- seq(200,20

Re: [R] How to compute the density of a variable that follows a proportional error distribution

2011-01-03 Thread mahesh samtani
Hello, Please accept my apologies; the symbols from my earlier message did not make it through correctly on the mailing list. Here is my question again: I am trying to compute the density of a variable k that is either (1) Normally distributed; (2) Log-Normally distributed; or (3) follows proportio

Re: [R] optimize

2011-01-03 Thread David Winsemius
On Jan 3, 2011, at 9:52 AM, Muhammad Rahiz wrote: Hi all, I'm trying to get the value of y when x=203 by using the intersect of three curves. The horizontal curve does not meet with the other two. How can I rectify the code below? Extend the appropriate axes. xlim and ylim are the argume

Re: [R] optimize

2011-01-03 Thread Joshua Wiley
Hi Muhammad, On Mon, Jan 3, 2011 at 6:52 AM, Muhammad Rahiz wrote: > Hi all, > > I'm trying to get the value of y when x=203 by using the intersect of three > curves. The horizontal curve does not meet with the other two. How can I > rectify the code below? What is your code supposed to do? Thi

Re: [R] using " plot" with time series object - " axes = FALSE" option does not appear to work

2011-01-03 Thread Ben Bolker
Clifford Long gmail.com> writes: > > Dear R-help, > > I am told by Professor Ripley that my question (quote) wastes the time > of (and has insulted on-list) the R developers by falsely claiming > there are problems with their code. I am writing, as he instructed, > to publicly apologize to the

[R] optimize

2011-01-03 Thread Muhammad Rahiz
Hi all, I'm trying to get the value of y when x=203 by using the intersect of three curves. The horizontal curve does not meet with the other two. How can I rectify the code below? Thanks Muhammad ts <- 1:10 dd <- 10:1 ts <- seq(200,209,1) dd <- c(NA,NA,NA,NA,1.87,1.83,1.86,NA,1.95,1.9

Re: [R] using "plot" with time series object - "axes = FALSE" option does not appear to work

2011-01-03 Thread Clifford Long
Dear R-help, I am told by Professor Ripley that my question (quote) wastes the time of (and has insulted on-list) the R developers by falsely claiming there are problems with their code. I am writing, as he instructed, to publicly apologize to the R developers for any slight that my question migh

Re: [R] Probably with default library tree

2011-01-03 Thread Chiquoine, Ben
Josh and Gabor thanks for those two extremely useful responses. I played around with those environment variables (R_LIBS and R_LIBS_USER) for a bit without success but I like the rationale hat you've given for the alternate path and will probably just stick with it. Thanks, Ben -Original

Re: [R] using "plot" with time series object - "/hijack thread for quetion re R data capture

2011-01-03 Thread Mike Marchywka
> Date: Mon, 3 Jan 2011 00:03:15 -0600 > From: gnolff...@gmail.com > To: r-help@r-project.org > Subject: [R] using "plot" with time series object - "axes = FALSE" option > does not appear to work > > Dear R-help, > > I am attempting to plot data using standard R plot utilities. The > data was

Re: [R] using "plot" with time series object - "axes = FALSE" option does not appear to work

2011-01-03 Thread Victor F Seabra
I guess you should specify x and y variables: plot(x=PCECTPI$var1, y=PCECTPI$var2, axes = FALSE) __ 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-gui

Re: [R] factor names

2011-01-03 Thread Iasonas Lamprianou
Thank you Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357-22-713178 Fax: +357-22-590539 Honorary Research Fellow Department of Education The Universit

Re: [R] factor names

2011-01-03 Thread Henrique Dallazuanna
Try this: f <- factor(sample(c("Engineer", "Doctor", "Teacher"), 50, rep = TRUE)) levels(f)[levels(f) %in% c('Doctor', 'Teacher')] <- paste(abbreviate(c('Doctor', 'Teacher'), 5), collapse = " & ") On Mon, Jan 3, 2011 at 9:18 AM, Iasonas Lamprianou wrote: > Dear all > I have a factor variable wh

[R] factor names

2011-01-03 Thread Iasonas Lamprianou
Dear all I have a factor variable which holds values like "Engineer", "Doctor", "Teacher" etc. I would like to collapse those categories so that Teachers and Sociologists form one category named "Teach & Soc" etc. However, I do not know how I can do it. Recoding does not seem to work. Thank

Re: [R] personal details appearing on google such after I pose a question in the R-help forum

2011-01-03 Thread Michael Sumner
library(fortunes);fortune() John Miller: How do I prevent google search to post my questions asked here?? Martin Maechler: you don't: R-help is famous and celebrity can't be gotten rid of ;-) -- John Miller and Martin Maechler R-help (June 2004) On Mon, Jan 3, 2011 at 8:06 PM, Siri Bjon

Re: [R] personal details appearing on google such after I pose a question in the R-help forum

2011-01-03 Thread Siri Bjoner
Not use a sig when posting to the forum? It is an open forum, and as such, all details that are sent to the list will also appear in the searchable archives. I don't think it is possible to reserve oneself against appearing in the archives... Siri. Siterer "taby gathoni" : Hi all, Out o

Re: [R] error in calling source(): invalid multibyte character in parser

2011-01-03 Thread Prof Brian Ripley
On Mon, 3 Jan 2011, peter dalgaard wrote: On Jan 3, 2011, at 08:32 , Luca Meyer wrote: Being italians when writing comments/instructions we use accented letters - like à, ò, è, etc when running R scripts using such characters I get and error saying: invalid multibyte character in parse

[R] personal details appearing on google such after I pose a question in the R-help forum

2011-01-03 Thread taby gathoni
Hi all, Out of curiosity I googled my name yesterday to just see what new information  in the web there is associated with me. To  my surprise, I found in addition to the  questions i have posted on this forum, my email address and my signature details (name, address, telephone number) seem

[R] Keeping lattice plot axes in the same place when legend sizes vary

2011-01-03 Thread Wayne Rochester
Dear all, I am generating a number of lattice plots with common axes but varying legends. There is one plot in the window at a time. I would like the plot axes to always appear in the same place within the window (and saved images) so the plots do not move around when they are combined in a v

Re: [R] Errors in installing Matrix package

2011-01-03 Thread Dennis Murphy
Hi: On Sun, Jan 2, 2011 at 7:36 PM, Yingfeng Zheng wrote: > Good Morning R-Help Community: > > I cannot install the "Matrix" package on R (version 2.12.1). Here are the > errors: > > > Install or load? They're different processes... > > > > Library(Matrix) > Why is library capitalized? R *is*

Re: [R] error in calling source(): invalid multibyte character in parser

2011-01-03 Thread peter dalgaard
On Jan 3, 2011, at 08:32 , Luca Meyer wrote: > Being italians when writing comments/instructions we use accented letters - > like à, ò, è, etc when running R scripts using such characters I get and > error saying: > > invalid multibyte character in parser > > I have been looking at the he