r code help

  1. U

    Efficient R coding

    Consider the following data frame: dat <- data.frame(x=c(0,128.96,0,397.37,56.80,0,801.79,340.79), y=c(40,97,172,186,210,246,327,432), d=c(1,0,0,0,0,1,0,1), A=c(0,1,0,1,1,0,1,1), B=c(0,1,0,1,1,0,1,1)) For each value of the column y which satisfies...
  2. C

    How is to define the lme model syntax .

    The following model assumes that the covariance between the random intercepts and random slopes across Subjects is Zero . lmer(Reaction ~ Days + (Days||Subject), sleepstudy) How can I write the same model by "lme" function under nlme pacakge assuming zero covariance between the random...
  3. C

    Specifying the variance component model "varComp" in R .

    I am trying to fit a random slope model by "varComp" in R . For the following example , that is in "lmer" syntax , how can I write it in "varComp" syntax : library(lme4) library(varComp) fm1 <- lmer(Reaction ~ Days + (Days||Subject), sleepstudy) I am not understanding...
  4. O

    Indices with data bound between -1 and +1

    Hi all, I would really appreciate some advice on how to deal with my data. I would like to compare data calculated from an index which is bounded by -1 and +1. I was surprised to see that previously published use of this index has been analysed using standard ANOVA, because I was led to...
  5. C

    How can I calculate Standard Errors of Variance Estimates .

    I want to extract the "standard error" of variance component from the output of "lmer" . In Chapter 12 , Experiments with Random Factors , of the book Design and Analysis of Experiments, written by Douglas C. Montgomery , at the end of the chapter , Example 12-2 is done by SAS . In Example...
  6. C

    Specifying formula of lmer() function in R .

    Am I specifying my lmer model correctly ? Say I have one independent variable (X) at individual level and one independent variable (Z) at group level. Both are continuous variable . If the model is...
  7. C

    How to compute the residual standard deviation from `glmer()` function in R?

    I want to extract standard deviation of residual from `glmer()` function in R . So I wrote : lmer_obj = glmer(Y ~ X1 + X2 + (1|Subj), data = D, family = binomial) sigma(lmer_obj) I noticed that the last command `sigma(lmer_obj)` returns always "1" irrespective of data ...
  8. B

    Simple probability in R

    Hello, newbie here. I have just started with R and am having trouble with how to enter a seemingly very simple probability question in R. The question is: There is a big court case, and 20% of the adult population believe the accused is innocent prior to jury selection. Assume that the 12...
  9. R

    Fitting confidence interval lines in R using the predict function

    Hi everyone, I am currently attempting to fit confidence interval lines to a regression model in R, but every time I try to for example add a +2 confidence interval line, I end up with two of them! I have been over my code many times and can't find any obvious mistakes. This is the code I...