Skip to content

Instantly share code, notes, and snippets.

@jcheong0428
Created March 20, 2021 05:33
Show Gist options
  • Save jcheong0428/22bfe8767c5a76b4af1f428c68fa53b5 to your computer and use it in GitHub Desktop.
Save jcheong0428/22bfe8767c5a76b4af1f428c68fa53b5 to your computer and use it in GitHub Desktop.
Manipulating data in R from rpy2
%%R -i data
# Specify Evit as a factor and set a polynomial contrast.
data$Evit <- as.factor(data$Evit)
contrasts(data$Evit) <- contr.poly
m<-lmer('Weight ~ Time * Evit + (1|Pig)', data=data)
print(summary(m))
# Relevel categorical data so base level is Evit100
m<-lmer('Weight ~ Time * relevel(Evit, ref="Evit100") + (1|Pig)', data=data)
print(summary(m))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment