Re: [R] Save multiple plots in a single pdf file when the plots are generated by a single plot command

2012-06-15 Thread Peter Ehlers
On 2012-06-15 16:14, Debs Majumdar wrote: I thought that might be the case and did a couple of dev.off() even though I started a new R session. Each time I try to use the plot, it comes up with "pdf 2". I think that the code for plot.lordif may have a bug; it contains this line: if (Sys.in

Re: [R] Replication of linear model/autoregressive model

2012-06-15 Thread R. Michael Weylandt
On Fri, Jun 15, 2012 at 6:56 AM, Al Ehan wrote: > Hi, > > I would like to make a replication of 10 of a linear, first order > Autoregressive function, with respect to the replication of its innovation, > e. for example: > > #where e is a random variables of innovation (from GEV distribution-that >

Re: [R] time zones and the chron to POSIXct conversion

2012-06-15 Thread R. Michael Weylandt
It does seem that as.POSIXct.date doesn't respect the tz= argument the generic suggests it would. I'd think this is a bug that could be changed without breaking back-compatibility, but I don't have the power to make such things happen. R-Core ruling? Michael On Fri, Jun 15, 2012 at 2:25 AM, Jann

Re: [R] Loop Help

2012-06-15 Thread Özgür Asar
Hi, Do you mean printing the results of degree_w(net.101031) and apply(net.101031, 2, mean) in a loop? see ?print to print the results and ?cat as a suporting function. For instance see the toy example below for (i in 1:5){ print('') cat("model",i,"\

Re: [R] Cholesky decomposition error

2012-06-15 Thread nataraj
Dear Mr.Kjetil, Thanks for your comment. You have already pointed me the article in reply to one of my earlier post to this list and I am following the paper. Now I am checking for condition for positive definiteness for original matrix using a simple script (got from earlier posting in the lis

Re: [R] time zones and the chron to POSIXct conversion

2012-06-15 Thread David Winsemius
On Jun 15, 2012, at 11:18 PM, R. Michael Weylandt wrote: It does seem that as.POSIXct.date doesn't respect the tz= argument the generic suggests it would. Why would a date function have a TZ? `as.Date` wouldn't. I'd think this is a bug that could be changed without breaking back-compatibili

Re: [R] time zones and the chron to POSIXct conversion

2012-06-15 Thread Jeff Newmiller
> Why would a date function have > a TZ? `as.Date` wouldn't. This statement seems nonsensical to me. POSIXt objects have tzone attributes. Date and chron objects do not. Since they do not include tzone, logically the user should be able to supply it during a conversion from Date or chron to PO

Re: [R] time zones and the chron to POSIXct conversion

2012-06-15 Thread R. Michael Weylandt
On Fri, Jun 15, 2012 at 11:26 PM, David Winsemius wrote: > > On Jun 15, 2012, at 11:18 PM, R. Michael Weylandt wrote: > >> It does seem that as.POSIXct.date doesn't respect the tz= argument the >> generic suggests it would. > > > Why would a date function have a TZ? `as.Date` wouldn't. > It seems

Re: [R] wgcna

2012-06-15 Thread Fabrice Tourre
Peter, Thank you. In fact, I am also very interesting to WGCNA. On Sat, Jun 16, 2012 at 3:29 AM, Peter Langfelder wrote: > On Fri, Jun 15, 2012 at 8:04 AM, Ingezz wrote: >> Dear Peter, >> >> I am trying to apply the WGCNA meta-analysis for two (or more) microarray >> datasets-tutorial to my own

[R] MCMCfactanal inquiries

2012-06-15 Thread Haillie Lee
Dear list, Hi, I am emailing to see if it would be possible to get some help on running MCMCfactanal on R. I am trying to derive a single score measuring a likelihood to recieve IMF loans from 8 different variables related to IMF representation but for reason, my codes do not seem to take me to wh

[R] Flexmix package

2012-06-15 Thread kimk
Hi, I am using the package "flexmix" and would like get some assistance. I am trying to run two equations jointly Y1=X1B+E1 Y2=X2G+E2 So that I have X and Y in a matrix format and would like to run the latent class model using flexmix. Though, my problem here is that Flexmix automatically genera

Re: [R] Wrong computation of time differenze in POSIXct - additional digits

2012-06-15 Thread arun
Hi, I checked the same with strptime.  It is the rounding issue. Try this: dat1<-data.frame(datetime=c("2012-06-15 16:32:39.0025 CEST","2012-06-15 16:32:39.0086 CEST"))  op<-options(digits.secs=4)  dat1$datetime<- strptime(dat1$datetime, "%Y-%m-%d %H:%M:%OS") formatC(as.numeric(difftime(dat1[

<    1   2