Re: [R] speed issue in simulating a stochastic process

2014-11-06 Thread Rolf Turner
On Thu, Nov 6, 2014 at 2:05 PM, Matteo Richiardi wrote: Final question: in your code you have mean(M[t-1L,]): what is the 'L' for? I removed it at apparently the code produces the same output... The constant "1L" is stored as an integer; the constant "1" is stored as double precisio

Re: [R] R version 2.14.2 installation on mac OS X Yosemite

2014-11-06 Thread Uwe Ligges
On 06.11.2014 20:20, Dalela, Deepansh wrote: Hi, I noticed your email on the R support page. I am currently using SPSS 21 for mac, and need to install a macro plug-in which requires me o install R 2.14 (I already have R 3.0). When i download the R 2.14 package from the CRAN website, and try

Re: [R] Lme4 Package Help!

2014-11-06 Thread Ben Bolker
Daniel Mello ucmerced.edu> writes: > > Hello, all! > > So, as stated in the title, the Lme4 package used > to output p-values for the > fixed effects. What happened?! > > Literally 2 weeks ago, I ran code, got output with no errors, and had > p-values listed for my fixed effects. > > Now, ru

Re: [R] R 2.14.2 installation on mac OS X yosemite

2014-11-06 Thread Uwe Ligges
On 06.11.2014 20:21, Deepansh Dalela wrote: Hi, I noticed your email on the R support page. I am currently using SPSS 21 for mac, and need to install a macro plug-in which requires me o install R 2.14 (I already have R 3.0). When i download the R 2.14 package from the CRAN website, and try

[R] R version 2.14.2 installation on mac OS X Yosemite

2014-11-06 Thread Dalela, Deepansh
Hi, I noticed your email on the R support page. I am currently using SPSS 21 for mac, and need to install a macro plug-in which requires me o install R 2.14 (I already have R 3.0). When i download the R 2.14 package from the CRAN website, and try to install it on my Macintosh HD, it gives me an

[R] R 2.14.2 installation on mac OS X yosemite

2014-11-06 Thread Deepansh Dalela
Hi, I noticed your email on the R support page. I am currently using SPSS 21 for mac, and need to install a macro plug-in which requires me o install R 2.14 (I already have R 3.0). When i download the R 2.14 package from the CRAN website, and try to install it on my Macintosh HD, it gives me a

[R] Lme4 Package Help!

2014-11-06 Thread Daniel Mello
Hello, all! So, as stated in the title, the Lme4 package used to output p-values for the fixed effects. What happened?! Literally 2 weeks ago, I ran code, got output with no errors, and had p-values listed for my fixed effects. Now, running THE SAME CODE with THE SAME DATASET (nothing at all ha

Re: [R] VGAM package : Frechet distribution - 2 parameter estimation

2014-11-06 Thread Thomas Yee
Hello, a day or two ago I submitted VGAM 0.9-5 to CRAN, which has myriads of changes to family functions---their names, their arguments, and their order thereof. Especially regarding family functions for discrete and continuous distributions. In a nutshell, I found lots of inconsistencies while

[R] Differences between MTDFReml and kinship2

2014-11-06 Thread David Duffy
From: "silvano" I fitted a genetic model using kinship2 and I compared it with MTDFReml program output. The residual variance of both are very close but the genetic variance are very differents. The output are: MTDFReml: genetic variance = 1.24015 residual variance = 5.93424 R (Kinship2): g

Re: [R] speed issue in simulating a stochastic process

2014-11-06 Thread William Dunlap
Loops are not slow, but your code did a lot of unneeded operations in each loop. E.g, you computed D$id==i & D$t==t for each row of D. That involves 2*nrow(D) equality tests for each of the nrow(D) rows, i.e., it is quadratic in N*T. Then you did a data.frame replacement operation D[k,]$y

Re: [R] Daylight Saving Time

2014-11-06 Thread Jeff Newmiller
?timezones You probably need to recompile R. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#..

Re: [R] Daylight Saving Time

2014-11-06 Thread Vasantha Kumar Kesavan
In LINUX, I don't want R to use the operating system zoneinfo(Olsen database) instead of that I like to point different path which has the latest zoneinfo(latest Olsen database). On Fri, Nov 7, 2014 at 2:34 AM, Jeff Newmiller wrote: > Working now... after what action? > > AFAIK on *NIX systems R

Re: [R] Daylight Saving Time

2014-11-06 Thread Jeff Newmiller
Working now... after what action? AFAIK on *NIX systems R uses the OS installation of the Olsen database, so on a fresh login R should pick up any OS update you have installed. --- Jeff NewmillerThe

Re: [R] Daylight Saving Time

2014-11-06 Thread Vasantha Kumar Kesavan
Thanks Jeff Newmiller, it is working now. I would like to know, the same kind of configuration can be done in Linux and Solaris platform. Instead of R is mapping to operating system(/usr/share/. /usr/share/lib/) zoneinfo directory. Thanks Vasanth On Fri, Nov 7, 2014 at 1:18 AM, Jeff Newmiller

Re: [R] Daylight Saving Time

2014-11-06 Thread Jeff Newmiller
R on Windows uses the Olsen timezone database, a copy of which is stored with R in the Program Files directory (e.g. R/R-3.1.1/share/zoneinfo). You could update the file yourself if you can find a corrected version, or download an updated version of R. ---

Re: [R] Timezone Upgrade

2014-11-06 Thread Vasantha Kumar Kesavan
I am log out and login many times. but no luck. Thanks Vasanth On Fri, Nov 7, 2014 at 1:04 AM, Jeff Newmiller wrote: > Log out and log back in again. For many people this may happen in the > normal course of daily use. > --

Re: [R] Timezone Upgrade

2014-11-06 Thread Jeff Newmiller
Log out and log back in again. For many people this may happen in the normal course of daily use. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go..

[R] Timezone Upgrade

2014-11-06 Thread Vasantha Kumar Kesavan
Hi, I would like to know, once the operating system timezone information is updated, what step should be carried out at R for reflecting the operating system(LINUX X64 SOLARIS X64 and WINDOWS X64) update. Thanks Vasanth [[alternative HTML version deleted]] _

[R] Daylight Saving Time

2014-11-06 Thread Vasantha Kumar Kesavan
Hi, I am working R on windows 2012 R2 platform, I have updated the latest hotfixes for time zone information(Microsoft KB 2981580. ). But still R is not populating correct date time values(Standard and Daylight saving). Could you please advise me, how to

Re: [R] speed issue in simulating a stochastic process

2014-11-06 Thread William Dunlap
I find that representing the simulated data as a T row by N column matrix allows for a clearer and faster simulation function. E.g., compare the output of the following two functions, the first of which uses your code and the second a matrix representation (which I convert to a data.frame at the e

Re: [R] speed issue in simulating a stochastic process

2014-11-06 Thread Thomas Adams
Matteo, Ah — OK, N=20, I did not catch that. You have nested for loops, which R is known to be exceedingly slow at handling — if you can reorganize the code to eliminate the loops, your performance will increase significantly. Tom On Thu, Nov 6, 2014 at 7:47 AM, Matteo Richiardi wrote: > I wis

Re: [R] cluster + tt terms in coxph

2014-11-06 Thread Henric Winell
On 2014-11-05 14:50, Therneau, Terry M., Ph.D. wrote: This is fixed in version 2.37-8 of the survival package, which has been in my "send to CRAN real-soon-now" queue for 6 months. Your note is a prod to get it done. I've been updating and adding vignettes. Is your fixed code publicly availa

Re: [R] speed issue in simulating a stochastic process

2014-11-06 Thread Thomas Adams
Matteo, I tried your example code using R 3.1.1 on an iMac (24-inch, Early 2009), 3.06 GHz Intel Core 2 Duo, 8 GB 1333 MHz DDR3, NVIDIA GeForce GT 130 512 MB running Mac OS X 10.10 (Yosemite). After entering your code, the elapsed time from the time I hit return to when the graphics appeared was

[R] speed issue in simulating a stochastic process

2014-11-06 Thread Matteo Richiardi
I wish to simulate the following stochastic process, for i = 1...N individuals and t=1...T periods: y_{i,t} = y_0 + lambda Ey_{t-1} + epsilon_{i,t} where Ey_{t-1} is the average of y over the N individuals computed at time t-1. My solution (below) works but is incredibly slow. Is there a faster

Re: [R] limit of cmdscale function

2014-11-06 Thread David L Carlson
You avoid the call to cmdscale() by supplying your own starting configuration (see the manual page for the y= argument). You could still hit other barriers within isoMDS() or insufficient memory on your computer. - David L Carlson Department of Anthropology Te

Re: [R] VGAM package : Frechet distribution - 2 parameter estimation

2014-11-06 Thread Rolf Turner
Dear Ms. Gobin, There would appear to be a typo in the package manual that appears on CRAN. Doing ?Frechet (it would be nice to have this aliased also to "frechet") points you immediately to frechet. With "frechet2" substituted for "frechet" your code, everything works. cheers, Rolf Turn

Re: [R] VGAM package : Frechet distribution - 2 parameter estimation

2014-11-06 Thread Katherine Gobin
Dear Mr Michael, Thanks a lot for your guidance. The pdf file describing VGAM package has mentioned 'frechet' in the example, so I got the error. Regards Katherine On Thursday, 6 November 2014 2:54 PM, Michael Dewey wrote: On 06/11/2014 06:04, Katherine Gobin wrote: > Dear R forum, > >

Re: [R] VGAM package : Frechet distribution - 2 parameter estimation

2014-11-06 Thread Michael Dewey
On 06/11/2014 06:04, Katherine Gobin wrote: Dear R forum, I am trying to execute following code (Page no 259 - VGAM.pdf) # . library(VGAM) set.seed(123) fdata <- data.fra

[R] Time dependant covariate

2014-11-06 Thread Paula Iso-Markku
Hello, Can I get help? I have the attached kind of material and I should create a time-dependant covariate from the variate "Optime"? I have tried using survsplit function but I don't know exactly how to do it. Or I can't get it to work. How can I split every id's follow-up time into two rows acco

[R] limit of cmdscale function

2014-11-06 Thread Kawashima, Masayuki
Hi We have a few questions regarding the use of the "isoMDS" function. When we run "isoMDS" function using 60,000 x 60,000 data matrix, we get the following error message: cmdscale(d, k) : invalid value of 'n' Calls: isoMDS -> cmdscale --

[R] VGAM package : Frechet distribution - 2 parameter estimation

2014-11-06 Thread Katherine Gobin
Dear R forum, I am trying to execute following code (Page no 259 - VGAM.pdf) # . library(VGAM) set.seed(123) fdata <- data.frame(y1 = rfrechet(nn <- 1000, shape = 2 + exp(1))