similar to: beginner's problem in displaying large data

Hi, there. Is there any function in R to plot the probability-probability plot (PP plot)? Suppose I am testing some data against normal. Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Is there a simple function in R to get a moving standard deviation (i.e. for the last x samples)? My goal is to plot bollinger bands around a moving average for price data. I use kernel smoothing for the moving average. cheers and thanks! over and out -- doktora

I have a very simple problem. When using glm to fit binary logistic regression model, sometimes I receive the following warning: Warning messages: 1: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, 2: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x = X, y = Y, weights = weights, start = start,

two ARIMA models, both have several bars signicant in ACF and PACF plots of their residuals, but when run Ljung Box tests, both don't show any significant correlations... however, one model has p-value that is larger than the other model, based on the p-values, can I say the model with larger p-values should be better than the model with smaller p-values? [[alternative HTML version

# I can not draw a 45 degree line on a qqnorm() plot, jj <- sample(c(1:100), 10) qqnorm(jj) abline() don't work. Thank you.

Hi, I would like to get the coefficients of x^0, x^1, x^2, . , x^6 from expansion of (1+x+x^2)^3. The result should be 1, 3, 6, 7, 6, 3, 1; How can I calculate in R? You help will be greatly appreciated. Peter [[alternative HTML version deleted]]

Hello, has someone written by chance a function to extract the variance-covariance matrix from a lmer-object? I've noticed the VarCorr function, but it gives unhandy output. Regards, Roel de Jong

Hi, I would like to know whether there exist algorithms to compute the coefficients or, at least, the degree of the minimal polynomial of a square matrix A (over the field of complex numbers)? I don't know whether this would require symbolic computation. If not, has any of the algorithms been implemented in R? Thanks very much, Ravi. P.S. Just for the sake of completeness, a

dear R users, I need some help for multivariate analysis. I have 2 anaesthetic treatment groups (20 patients/group) where I register heart frequency and pressure for 60 min (repeated measures every 5 minutes). I would like to perform a test to check if treatments are different in controlling freq and pressures during the anaesthesia, but i would like to have also an overall measure and not

Hi Can any one inform me an r package used to fit Cubic M-splines or r functions used to have design and penalty matrices in fitting cubic  M-splines please? Thanks, Tigist [[alternative HTML version deleted]]

the following the the lower.tri matrix in a file named luxry.car and i want to read it in R as a lower.tri matrix.how can i do? i have try to use help.search("read"),but no result what i want. 1.000 0.591 1.000 0.356 0.350 1.000

hello! i don't want to test my sample data for normality, but exponential- lognormal- or gammadistribution. as i've learnt the anderson-darling-test in R is only for normality and i am not supposed to use the kolmogorov-smirnov test of R for parameter estimates from sample data, is that true? can you help me, how to do this anyway! thank you very much! nadja

I'm trying to rename the columns in a list of data.frames using the following... for(i in 1:length(filenames)) { assign(names(get(filenames[i])), c("name", "infood", "time") ) } R returns no errors, but the names are unchanged in the data.frames. The original names were things like > names(get(filenames[2])) [1] "Tc45w4.V1" "Tc45w4.V2"

How can I extract K (number of parameters) from an AIC calculation, both to report K itself and to calculate AICc? I'm aware of the conversion from AIC -> AICc, where AICc = AIC + 2K(K+1)/(n-K-1), but not sure of how K is calculated or how to extract that value from either an AIC or logLik calculation. This is probably more of a basic statistics question than an R question, but I thank

Hello, can someone tell me how to extract a partial table from a multiway table? Something else: I tried to do basic operations on cross-tables, like subtracting two cross tables with same dimensions and i was unable to do so. Is there another way? thanks! __________________________________________________________

Hello, Brand new user of R that I want to interface with Windew (http://www.pcsoft.fr), I am looking for a detailed documentation of R(D)Com. Does anybody know where I can find it, or do you have any form of experience of the use of R with Windev? Thanks for your help. PS : As I am new to the group, people interested in some details about the requestor could find his resume at address

Greetings all, Is there a way to do the equivalent of a lmList (fit separate models for each subject) in the package nlme for a glm? Thank you very much, Peter B. [[alternative HTML version deleted]]

Hi Does R have a proceedure/software for capture recapture? Thank you. Lawrence Lessner Please respond to LLessner at nycap.rr.com

Hello, I have been asked to perform a exponential correlation on a data set. Unfortunately, i did not find any simple example on how to do that with R. If the correlation is significant, one asked me to also find the thresold of the acceleration slope. Have any of you a little example or a good link to get me started? Best regards Blaise -- Mutt, c'est de la balle, quand ??a marche :-)

Hello R people, I am need some help using the try() function. Currently I am running a loop which uses arima() for some values of p and q, which sometimes crashes. When it crashes, I want the program to just ignore it and move on to the next values to loop through. I currently have this, looping through a range of values for p and q: lo = try( arima1 <- arima( y, order=c( p, 0, q ) ) if(