similar to: variance components in lme

Hello, I have two data frames, NameAndConc.df and WhichOnes.df. NameAndConc.df has two columns, the first column is the names of some material, and the second column is the concentration of active ingredient. > NameAndConc.df name conc 1 material1 0.8 2 material2 0.5 3 material3 0.4 WhichOnes.df has two columns, each of which specifies which material. > WhichOnes.df

Hello, I can't seem to make qqmath plot text. v <- rnorm(8) qqmath(~v) # this is ok, but qqmath(~v, panel=function(x, y) { points(x, y) text(0, 0, 'print me', cex=2)}) isn't ok. The axis frame is printed properly (ie, as in qqmath(~v)), but neither 'print me' nor the points show up on the graph. I can, however, see the upper half of the word 'me' in the

Hello, I'm somewhat new to R. I've searched the archive for the last year and tried to consult the manual pages for the following problem, but did not find an answer. I want to sort an array by the index values. Here is the array >acc.gp.bl.wtmn Gp17 Gp4 1 0.5703125 0.6406250 10 0.7812500 0.7109375 11 0.8046875 0.7343750 12 0.8359375 0.7890625 13 0.8515625 0.7109375

below is my data frame. I would like to compute summary statistics for mgl for each river mile (mean, median, mode). My apologies in advance- I would like to get something like the SAS print out of PROC Univariate. I have performed an ANOVA and a tukey LSD and I would just like the summary statistics. thanks stephen RM mgl 1 215 0.9285714 2 215 0.7352941 3 215 1.6455696 4 215

Buenas, una pregunta. Si yo estoy calculando un modelo lineal, el caso más simple, 1 variable respuesta y una variable explicativa y creo un modelo, me da un R2 del 80% y quiero ver como es esa relacion entre las variables, para calcular el error de predicción del modelo, basta con ver el intervalo de confianza del modelo e irme a los extremos? Por si no me he expresado bien, un ejemplo tonto:

Enun ejemplo real estoy viendo como el intervalo de confianza usando lo que me comentas me ha salido mucho más pequeño de lo que la realidad luego refleja. ¿Cómo es esto posible?? Es decir, veo que para valores de 2,70 obtengo una respuesta de entre 2,69 y 2,90 y sin embargo luego en la realidad tengo valores entre 2,20 y 3 Gracias Jesús From: jorgeivanvelez en gmail.com Date: Thu, 21 Apr

How would one go about putting titles in each of several plots that are generated from within a call to tapply? For example I'd like the following two barplots to have titles 'Group 1' and 'Group 2', where '1' and '2' come from the levels of 'group'. group <- gl(2, 10) result <- sample(c('A', 'B'), size=length(group),

Hi R fans, I was trying to write some code in R when I realized that I was just writing a FORTRAN 77 style program in the R syntax, and hence coming nowhere near tapping the potential of R. I'm wondering what is a good reference that would help me get up to speed in programming in R (i.e., using these classes, methods, objects, accessor functions, and many other things I don't

Dear colleagues, I'm conducting a meta-analysis of studies evaluating adherence of HIV-positive drug users into AIDS treatment, therefore I'm looking for some advice and syntax suggestion for running the meta-regression using proportions, not the usual OR/RR frequently used on RCT studies. Have already searched already several handbooks, R-manuals, mailing lists, professors, but... not

Hello, model.frame changes the factor labels when na.action = na.omit. > f <- gl(3, 2, 6, paste('m', 1:3, sep = '')) > r <- c(NA, NA, 3:6) > mf <- model.frame(~ r + f, na.action = na.omit) > mf r f 3 3 m1 4 4 m1 5 5 m2 6 6 m2 But it seems like it should be this > mf r f 3 3 m2 4 4 m2 5 5 m3 6 6 m3 intead. How do I get the second mf? Thanks a lot,

Hello all. I have two data frames. Group Start End G1 200 700 G2 500 1000 G3 2000 3000 G4 4000 6000 G5 7000 8000 and Pos C0 C1 200 0.9 0.6 500 0.8 0.8 800 0.9 0.7 1000 0.7 0.6 2000 0.6 0.4 2500 1.2 0.8