search for: saji

...ass(mydata)="numeric" but after using "fitdistr", I got the result below >fitdistr(mydata,"normal") mean sd NA NA (NA) (NA) the help doc of "fitdistr" does not mention anything about that, thus I need your help. Thank you in advanced, Saji from Shanghai -- View this message in context: http://n4.nabble.com/Help-with-function-fitdistr-in-MASS-tp997609p997609.html Sent from the R help mailing list archive at Nabble.com.

...etric and the mean is close to 0. It just that the spread is large (there are outliers). The right is the transformed data, and the distribution is obviously no normal. Can anyone explain that to me? Thank you in advanced. ----- ------------------------------------------------------------------ Saji Ren from Shanghai China GoldenHeart Investment Group ------------------------------------------------------------------ -- View this message in context: http://n4.nabble.com/Problem-about-Box-Cox-transformation-topic-in-html-form-tp1011015p1011015.html Sent from the R help mailing list archive at...

...ave a data series of 500 data, and I want to limit the value of it to be less than 1. Below is my code: >for (i in 1:500) +if( x[i] > 1) +x[i] = 1 but the system told me it's wrong. Can anyone told me the reason? ----- ------------------------------------------------------------------ Saji Ren from Shanghai China GoldenHeart Investment Group ------------------------------------------------------------------ -- View this message in context: http://r.789695.n4.nabble.com/How-to-perform-a-substitution-in-a-loop-tp2325048p2325048.html Sent from the R help mailing list archive at Nabble....

...ere the 'x' is the one generated above in 1). And the R told me: "mistakes in match.fun(FUN) : 'quantile(x, 0.8)' is not a function, character or symbol" Can anyone help? Thank you in advanced. ----- ------------------------------------------------------------------ Saji Ren from Shanghai China GoldenHeart Investment Group ------------------------------------------------------------------ -- View this message in context: http://n4.nabble.com/Confusion-in-quantile-and-getting-rolling-estimation-of-sample-quantiles-tp1015893p1015893.html Sent from the R help mailing...

Hello: I want to get a rolling estimation of the stdev of my data. Searching the document, I found the function "rollapply" in the zoo package. For example, my series is "c", and i want get a period of 10 days, so i write the command below: roll.sd = rollapply( c, 10, sd, na.pad = TRUE, align = 'right' ) but there is an error in it ,and the computing cannot be

...ma1*sma1*a(t-46) But when I use this forcast formula in Excel, it gives a totally different predict from R. And I don't know why? I guess the expression of the forcast formula of this SARMA(1,1)*(0,1)45 is wrong, but I don't know the right form. Can anybody help me with this?Thank, again! saji from Shanghai _________________________________________________________________ ÔõÑùÂò³µÆ±¸ü·½±ã£¿Î¢Èí´º½ÚËÑË÷£¬µÚһʱ¼äΪÄúÌṩ»ð³µÆ±ÐÅÏ¢£¡ [[alternative HTML version deleted]]

Guys: Is it possible to simulate a seasonal ARIMA model in R? Which package can do this job? saji from Shanghai

...the forecasting. Now I just want to know that How does R perform the predict procedure(the predict formula, the initial setting of errors,etc.)? I run the following commands and get the original code of the "predict" command, but I can't read it. Can anybody explain it to me? Thanks! saji from Shanghai the code: > getS3method("predict","Arima") function (object, n.ahead = 1, newxreg = NULL, se.fit = TRUE, ...) { myNCOL <- function(x) if (is.null(x)) 0 else NCOL(x) rsd <- object$residuals xr <- object$call$xreg xreg...