search for: fomula

Dear all, I found really difficult with the time series questions, please help me with this monthly airline series! I have run the following r code, and there is an error appeared at the end. The data files was enclosed in the email. I'm sorry the errors message appeared in chinese, but it says "plot.xy(xy.coords(x, y), type = type, ...) : errors in argument has more than 3

...an app that gets about 5000 unique visitors per day. I figured I''d give it plenty of Mongrels -- twenty to be specific. After running out of memory and hitting the swap periodically, I scaled it back to five and it still seems to serve up visitors fine. So, is there some super-secret math fomula I can use to determine the ideal number of Mongrels I should start in a cluster for an app? Joe ____________________________________________________________________________________ Do you Yahoo!? Everyone is raving about the all-new Yahoo! Mail beta. http://new.mail.yahoo.com

Hi, don' t understand why the function fomula have this error, i enclose the parameter "a" with the function I() Thank Ruben x<-1:5 y<-c( 2 ,4 , 6 , 8 ,11) formu<-y~I(a*x) form<-formula(formu) dummy<-data.frame(x=x,y=y) fm<-lm(form,data=dummy) Error in unique(c("AsIs", oldClass(x))) : Object "a&quo...

...rep(x = 0:1, times = 4)) Call: lm(formula = (1:8) ~ rep(x = 0:1, each = 4) + rep(x = 0:1, times = 4)) Coefficients: (Intercept) rep(x = 0:1, each = 4) 2.5 4.0 So while the model is identified because the two regressors are not the same, terms.fomula does not recognize this and drops the second regressor. What I would have wanted can be obtained by switching the arguments: R> lm((1:8) ~ rep(each = 4, x = 0:1) + rep(x = 0:1, times = 4)) Call: lm(formula = (1:8) ~ rep(each = 4, x = 0:1) + rep(x = 0:1, times = 4)) Coefficients:...

Dear R-help, This is an example in the {Hmisc} manual under rcorr.cens function: > set.seed(1) > x <- round(rnorm(200)) > y <- rnorm(200) > round(rcorr.cens(x, y, outx=F),4) C Index Dxy S.D. n missing uncensored Relevant Pairs Concordant Uncertain 0.4831 -0.0338 0.0462 200.0000

...alanced Data confidence interval calculation from Pg 62. The data are fill weights from bottles sampled randomly from a sample of four filling machines. There are 12 values, and the confidence intervals are for 1-2a = 95%. I have got the same results as the book but using slightly different fomulae (see variables for H1, G1 and H12 and G12). I'd appreciate any help, and any comments on whether their is a better way to do this. Thanks Paul. > BGBottles Machine weight 1 1 14.23 2 1 14.96 3 1 14.85 4 2 16.46 5 2 16.74 6 2 15.94 7...

...h = 4) + rep(x = 0:1, times = 4)) > Call: lm(formula = (1:8) ~ rep(x = 0:1, each = 4) + rep(x > = 0:1, times = 4)) > Coefficients: (Intercept) rep(x = 0:1, each = 4) 2.5 4.0 > So while the model is identified because the two > regressors are not the same, terms.fomula does not > recognize this and drops the second regressor. What I > would have wanted can be obtained by switching the > arguments: R> lm((1:8) ~ rep(each = 4, x = 0:1) + rep(x = 0:1, times =4)) > Call: lm(formula = (1:8) ~ rep(each = 4, x = 0:1) + rep(x &gt...

hi ,everyone, I have a script shown as below: > bj2=bjerrdata$tyerr[bjyearnum[2]+1:bjyearnum[3]] > length(bj2) [1] 448 while > b=bjyearnum[3] > a=bjyearnum[2]+1 > bj1=bjerrdata$tyerr[a:b] > length(bj1) [1] 169 it is different with bj2 and bj1 . and the array bjyearnum is [1] 0 279 448 633 1021 1365 1813 2237 2839 3314 3798 4157 12 why bj1 is 169 while bj2 is 448