similar to: means comparison with seasonal time series?

Dear R Help, I've just installed 1.7.1 and have found that the console now exhibits odd behavior when I paste commands into it. I'm on a PC running Windows XP Professional. I use Word as my script editor, and then paste commands into the R console using the standard Windows "copy" within Word then "paste" within Rgui. This worked without a hitch in 1.7.0, but now in

Dear R list, What's the best way in R to fill a matrix who's entries depend on some function of the indices? I'm currently doing: Q <- matrix(0, k, k) for (A in 1:k) { for (B in 1:k) { Q[A,B] <- my.function(A,B) } } but I wonder if there is a more terse way. Regards, Douglas Scofield Department of Biology d.scofield at

Hello! I''ve a ubuntu server running over a HVM DomU. After a hang, the domU doesn''t start again. The domU was not able to mount the virtual machine disk so I started with a live cd and I tried to mount the disk root part / with no success so I ran fsck.ext3 and it gave me input/output error in the domU and in dom0 the kernel says: Feb 14 18:22:28 scofield last message

I have a zooreg object and I want to be able to generate a value for seasons and 11-day composites paste it onto my zoo data frame, along with year, month and days. Right now I have the following to work from: eg. dat.zoo.mdy <- with(month.day.year(time(dat.zoo)), cbind(dat.zoo, year, month, day, quarter = (month - 1) %/% 3 + 1, dow = as.numeric(format(time(dat.zoo), "%w")))) For

Dear list, I am relatively new to ordinal models and have been working through the example given by Frank Harrell in the predict.lrm {Design} help All of this makes sense to me, except for the responses, i,e how do i interpret them? i would be extremely grateful if someone could explain the results? First i establish the date and model - > y <- factor(sample(1:3, 400, TRUE), 1:3,

Hi, Does anyone know if the D'agostino test is available with R ? Alex

Dear All, I am trying to fit a structural time series model using the StructTS function (package stats) with only 2 seasons (summer and winter). More than 2 seasons work fine but with 2 seasons I get this error: > fit <- StructTS(y.ts, type="BSM") Error in T[cbind(ind + 1L, ind)] <- 1 : subscript out of bounds I have looked at Prof. Ripley's 2002 RNews article but cannot

I am an R newbie, and I am facing what is for me a challenging problem. I am trying to write functions to solve the following problems. It seems as though R offers so many options I am not sure how to approach the problem. I have a data frame with dates in the first column and returns in the remaining columns. I would like to give a user a variety of options. I designate Series1 as the

Hi Irucka, No problem. I guess this could be done using: Sample data: Dailydo<- structure(list(Date = structure(c(11231, 11232, 11233, 11234, 11235, 11236, 11207, 11208, 11209, 11179, 11180, 11181, 11151, 11152, 11123, 11093, 11064, 11065, 11035, 11036, 11008, 11009, 10979, 10980, 11347, 11348, 11380, 11381, 11410, 11406, 11438, 11439, 11470, 11502, 11528, 11560, 11561, 11593, 11594,

Hello R People: When using the arima function with the seasonal option, are the seasonal options only good for monthly and quarterly data, please? Also, I believe that weekly and daily data are not appropriate for seasonal parm estimation via arima. Is that correct, please? Thanks, Sincerely, Laura Holt mailto: lauraholt_983 at hotmail.com download!

I'm setting up a cluster with multiple ups's using nut 2.2. The ups's are Belkin 1500 (ID 050d:0751) which don't appear to report a serial number. Is there any other way I can setup ups.conf with multiple of these ups over usb on Linux? If I just setup the vendor and product id's, every definition in ups.conf appears to match the first one, and therefore all of them but

Dear John, Now, I want to do dynamic analysis by dynamic slicing based on LLVM. I just see that you have implemented dynamic slicing for LLVM. I also check the Giri project, and I don't find the code. I hope I can learn something from your code. I am very appreciated if you can tell me how to get the code or send me one copy. Thank you very much! Pengfei -- View this message in context:

Dear R users, I am looking for some advice on the proper construction of a mixed model in R, using the difference in means as the response and treating within-means residuals as a random effect. I have a dataframe (my own, a snippet of which is given below) that is composed of observations of pollen viability in flowers along tree branches. Flowers (1 to 3 per position) were collected from

Hi there I'm working with ncdf data. I have different dataset for 4 runs, 4 seasons and 3 timeslices (48 datasets in total). The datasets have the following dimensions: 96 longitudes, 48 latitudes and 30 time steps. To read all of them in, I wrote the following loop: runs <- c("03","04","05","06") years <- c(1851,1961,2061) seasons <-

All, I've coded a function and it works manually if I copy it line by line into R. However, when I try to "load" (copy and paste) the entire function into R, I get the following error after the listed line of code: + N.j.list = lapply(rej.hyp, length) Error: symbol print-name too long Does anyone you know what this error means? Strangely, when I copy the same line verbatim

Hi, I was wondering if someone could please help me. I am doing a logistic regression to compare size at maturity between 3 seasons. My model is: fit <- glm(Mature ~ Season * Size - 1, family = binomial, data=dat) where Mature is a binary response, 0 for immature, 1 for mature. There are 3 Seasons. The Season * Size interaction is significant. I would like to compare the size at 50%

Sorry for the last email, sent too early. I have a small data set that has a hierarchical structure. It has both temporal (year, months) and spatial (treatment code and zone code). The following explains the data: WSZ_Code the water supply zone code (1 to 8) Treatment_Code the treatment plant which supplies each water supply zone (1 to 4)

All, I have a repeated measures data frame and was wondering if the covariance matrix can be calculated via some created indexing or built-in R function. Specifically, say there are 3 variables, where potassium concentration is measured 6 times on each patient. Patient number (discrete) Time (1 to 6, discrete) Potassium (continuous variable) I want the covariance/correlation matrix for the

Dear List Members, I am in the process of designing a study examining pay and performance in Major League Baseball across several seasons, and before I get too deeply into it, I'd like to ask whether the group members think that performance across seasons is independent, or if it needs to be treated like a time-series variable so that lack of independence can be controlled. Any ideas or

Hi, I was hoping to get some advice regarding the testing of interactions, when one factor is modelled as a random effect... I have a model with binomial error structure where the response variable is the proportion of time spent at the main sett (animals were tracked for 28 consecutive days in each season, and were recorded either at the main sett or an outlier sett, so the response variable is