Dear R users

I am new R user, execuse me I bother you, but I worked hard to find a
solution:

# data
ID <- c(1:100)
set.seed(21)
y <- rnorm(100, 10,2)
x1 <- rnorm(100, 10,2)
x2 <- rnorm(100, 10,2)
x3 <- rnorm(100, 10,2)
x4 <- rnorm(100, 10,2)
x5 <- rnorm(100, 10,2)
x6 <- rnorm(100, 10,2)
mydf <- data.frame(ID,y, x1,x2, x3, x4, x5, x6)
# just seperate analyis
require(nlme)
mod1<- lme(fixed= y ~ x1 + x2, random = ~ 1 | ID) # i want to put subject /
ID as random is it correct??
m1 <- anova(mod1)
m1 [1,4] # I am not getting exact value below <.0001???, how to get one?
m1 [2,4]
m1 [3,4]
# putting in a loop
for(i in length(mydf)){
  mylme <- NULL
  print(m1[i+1,]<- lme(fixed= mydf$y ~ mydf$x[,i] + mydf$x[,i+1], random = ~
1 | ID))}

could not help myself to work ! However I have the following output in my
mind
# The output in my mind a data fram with the following
model   p-intercept      variable1   p-value1   variable2      p-value2
1            <.0001           x1            0.9452     x2             0.5455
2           <.0001            x3            0.3301      x 4           0.9905
3          <.0001              x5            0.9971     x6            0.0487

I need a solution as I have tons of variables to work on. I am sowing these
six just for an example.

Thank you in advance;

ram sharma

	[[alternative HTML version deleted]]

