Skip to content

Instantly share code, notes, and snippets.

@johnstantongeddes

johnstantongeddes/df_fit_lm.R

Created June 13, 2014 16:22
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save johnstantongeddes/70eb26acd8353792166b to your computer and use it in GitHub Desktop.
Save johnstantongeddes/70eb26acd8353792166b to your computer and use it in GitHub Desktop.
Download ZIP
Challenge to fit a lm to every "gene" in df
Raw
df_fit_lm.R
df <- data.frame(gene = rep(c(paste("gene", 1:10, sep="")), 5),
treatment = rep(paste("trt", 1:5, sep=""), each=10),
temp = rep(1:5, each=10),
exp = rnorm(50))
head(df)
# lm for a single gene
lmout <- lm(exp ~ temp, data = df[df$gene == "gene1", ])
summary(lmout)
# note "p-value"
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment