search for: means2

Assume a C program invoked by .Call, that returns a list. Near the top of the program we allocate space for all the list elements. (It is my habit to use "xyz2" for the name of the R object and "xyz" for the pointer to its contents.) PROTECT(means2 = allocVector(REALSXP, nvar)); means = REAL(means2); PROTECT(u2 = allocVector(REALSXP, nvar)); u = REAL(u2); PROTECT(loglik2 = allocVector(REALSXP, 2)); loglik = REAL(loglik2); PROTECT(rlist = mknamed(VECSXP, outnames)); Can I assign the individual elements into rlis...

...e intervals for the means And then again for a negative binomial GLM require(MASS) mydata2 <- data.frame('treatments'=gl(4, 4, 64, LETTERS[1:4]), 'species'=gl(4, 1, 64, letters[1:4]), 'value'=rnegbin(64,1,1)) model2 <- glm.nb(value ~ treatments*species, data=mydata2) means2 <- with(cbind('x'=gl(16,4),mydata2), tapply(value, x, mean)) Can I get the confidence intervals for each mean in means1 and means2, according to the family used? Thanks in advance, all the best! Rafael. ___________________________________________________________________________...

...ave nothing to do with each other. I also tried calling plotCI with argument add=TRUE, which didn't seem to work (that is actually what I wanted I think). (It should look the same as if I called plotCI twice with same labels/xlim/ylim/etc.) plotCI(x = cbind(x1,x2), y = cbind(means1,means2), # means1 == ci1["Estimate",] xlim = c(0,100), #ylim = c(0.2,0.5), ylab = "System welfare", pch = 7, col = c("red","blue"), type = "b", uiw = cbind(uiw1,uiw2)) Thanks in advance, JeeBee.

Hello, Here is a general question that I think many of you will have insight to. I have 40 independent dataframes. Each has a dimension of 1359 x 15. For example >dim(meanAnnualWaterDepth) > 1539 15 I need to calculate the average value for each dataframe but only for columns 7: 15. I can do this with apply as follows: meanAnnualWaterDepth$Average <-