[R] R problem

2016-12-08 Thread Grace King
program or batch file.” Here is how R is called to run the macro in the script: *R --save < IPFPUMS.R* R now resides on the same server location as EMME. Our IT person created a symbolic link for R but still the same error occurred. Any thoughts will be much appreciated. Gr

Re: [R] Rolling window linear regression

2014-10-04 Thread Grace Shi
I have to do roll regression based on the Daily data. I use the past three weeks of daily returns as the estimation window and the regression is estimated rolling forward one week at a time generating time series estimates of beta. I know I should use the rollapply in zoo package. but I am

Re: [R] Lattice z colours when using formula class

2013-10-08 Thread Grace Garner
I've found one solution, which is to transform each matrix into a 1d array. This gives me the output I wanted. Grace On 08/10/2013 09:55, Grace Garner wrote: Dear all, I'm using Windows 7, R version 2.15.0 and lattice_0.20-6. I'm plotting observed data in lattice using wirefr

[R] Lattice z colours when using formula class

2013-10-08 Thread Grace Garner
rame(Temp~Dist*Time, shade=F, drape=T, colorkey = T, col.regions=mycols(200), col="transparent") Thank you, Grace __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://

[R] using alternative models in glmulti

2012-09-11 Thread Midgley, Meghan Grace
All, I am working on a multiple-regression meta-analysis and have too many alternative models to fit by hand. I am using the "metafor" package in R, which generates AIC scores among other metrics. I'm using a simple formula to define these models. For example, rma(Effect_size,variance, mods

Re: [R] mgcv package, problems with NAs in gam

2012-08-25 Thread grace
Deviance:0 > Residual Deviance: 0 AIC: -Inf > > What model are you trying to fit? > > best, > Simon > > > On 08/24/2012 10:36 PM, grace wrote: > > > Hi there, > > > > I'm using presence-absence data in a gam (i.e. 0 or 1 as

Re: [R] mgcv package, problems with NAs in gam

2012-08-25 Thread grace
I've attached the data to this message. http://r.789695.n4.nabble.com/file/n4641291/2005pip.csv 2005pip.csv -- View this message in context: http://r.789695.n4.nabble.com/mgcv-package-problems-with-NAs-in-gam-tp4641253p4641291.html Sent from the R help mailing list archive at Nabble.com.

[R] mgcv package, problems with NAs in gam

2012-08-24 Thread grace
Hi there, I'm using presence-absence data in a gam (i.e. 0 or 1 as values) I am trying to run a gam with 'dummy covariates' i.e. 1~1 unfortunately my model: * model<-gam(1~1, data=bats, family=negbin)* keeps putting out: * Error in gam(1 ~ 1, data = bats, family = negbin) : Not enough (non-NA

[R] Error; contrasts can be applied only to factors with 2 or more levels

2012-08-14 Thread grace
~x1, umf)) Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : contrasts can be applied only to factors with 2 or more levels > (occuRN3<-occuRN(~x3+x1 ~x1, umf)) Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : contrasts can be applied only to fa

[R] bivariate localpcf in Spatstat

2012-07-31 Thread Grace Hensen
f there is way to do it either in Spatstat or other R-packages. Thanks for your suggestions. Regards, Grace Hensen [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE d

[R] 2D image correlation/Packages for 2D CCA?

2012-07-05 Thread Grace, Miriam
Hello, I was just wondering if anyone could recommend a package for R that calculates 2D image correlation, e.g. by 2D canonical correspondence analysis? Thanks a lot! MG [[alternative HTML version deleted]] __ R-help@r-project.org mailing l

[R] survival analysis simulation question

2012-05-10 Thread Grace Ma
Hi, I am trying to simulate a regression on survival data under a few conditions: 1. Under different error distributions 2. Have the error term be dependent on the covariates But I'm not sure how to specify either conditions. I am using the Design package to perform the survival analysis using the

[R] Calculating the correlation between two matrices?

2012-04-14 Thread Grace, Miriam
relation coefficient using the rcorr function from the Hmisc package, but I'm not at all confident that this is suitable. Any advice on how to approach this problem would be really helpful. Thanks a lot, M. Grace [[alternative HTML version deleted]] __

[R] R-Censoring

2012-02-06 Thread Grace Kambey
Hi there, can somebody give me a guide as to how to generate data from weibull distribution with censoring for example, the code below generates only failure data, what do i add to get the censored data, either right or interval censoring q<-rweibull(100,2,10). Thank you Grace Kam stud

[R] R-Censoring

2012-02-05 Thread Grace Kambey
Hi there, can somebody give me a guide as to how to generate data from weibull distribution with censoring for example, the code below generates only failure data, what do i add to get the censored data, either right or interval censoring q<-rweibull(100,2,10). Thank you Grace Kam stud

[R] R-Censoring

2012-02-05 Thread Grace Kambey
Hi there, can somebody give me a guide as to how to generate data from weibull distribution with censoring for example, the code below generates only failure data, what do i add to get the censored data, either right or interval censoring q<-rweibull(100,2,10). Thank you Grace Kam stud

[R] r-help, survival

2012-02-04 Thread Grace Kambey
Hi, i have used the code below to come out with censored data with real values but if i want to generate data from q<-rweibull(25,2,4) that involves percent eg 5% or 10% censored i do not know how to do it. Can somebody show me how to do it. Thank you in advance d <- Surv(time=c(27,100,100,53,78

[R] Counting number of datasets and appending them

2010-12-29 Thread Li, Grace
you help me on this? I hope this is not some simplest R question. I am a beginner. Thanks a ton! Grace Confidentiality Notice: This e-mail message including at...{{dropped:12}} __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo