search for: nilaya

Hi R community, I am new bird to R and moved recently from SAS. I am no means expert on either but very curious learner. So your help crucial for me to learn R. I have already got positive expression. I was trying to fit a mixed model in animal experiment but stuck at simple point. The following similar example is from SAS mixed model pp 212. # data genetic_evaluation <-

...# mydf[x[1]] = mydf[x[2]] and value of mydf[x[1]] < median(mydf[x[1]]), the value of a1cf is -1 Cf <- d1cf + a1cf return(Cf) } Cim <- apply(nmat, 1, coffin) I tried hard could get a appropriate way to do the job. Can you help me? NIL ---------- Forwarded message ---------- From: Nilaya Sharma <nilaya.sharma@gmail.com> Date: Wed, May 4, 2011 at 11:32 AM Subject: simple question To: r-help@r-project.org Dear R experts I have simple question, please execuse me: #example data, the real data consists of 20000 pairs of variables K1 <- c(1,2,1, 1, 1,1); K2 <- c(1, 1,2,2,...

Dear All Sorry for this simple question, I could not solve it by spending days. My data looks like this: # data set.seed(1234) clvar <- c( rep(1, 10), rep(2, 10), rep(3, 10), rep(4, 10)) # I have 100 level for this factor var; yvar <- rnorm(40, 10,6); var1 <- rnorm(40, 10,4); var2 <- rnorm(40, 10,4); var3 <- rnorm(40, 5, 2); var4 <- rnorm(40, 10, 3); var5 <- rnorm(40, 15,

Dear all Sorry for simple question: I want to put the following option into look as number of X is large 1000 variables X1 <- sample(c(1,2, 3, 4),10, replace = T, prob = c(0.4, 0.2, 0.2, 0.2)) cv1 <- round(runif(2, 1, 10)) # X2 is copy of X1 X2 <- X1 # now X2 is different in cv1 random positions X2[cv1] <- 5 cv2 <- round(runif(2, 1, 10)) # X3 is copy of X2 X3 <- X2