Skip to content

Instantly share code, notes, and snippets.

@MJacobs1985
Created September 30, 2021 18:25
Show Gist options
  • Save MJacobs1985/656ab9137f217945b52d8800489eed7a to your computer and use it in GitHub Desktop.
Save MJacobs1985/656ab9137f217945b52d8800489eed7a to your computer and use it in GitHub Desktop.
How many samples do you need?
#### DATA IMPORT & MANAGEMENT ####
ARCdata <- read_csv("ARC_crosssectional_presentation.csv")
ARCdata$gain<-ARCdata$End-ARCdata$Start
str(ARCdata)
ARCdata<-as.data.frame(ARCdata)
ARCdata$Diet<-as.factor(ARCdata$Diet)
ARCdata$Tank<-as.factor(ARCdata$Tank)
#### DATA EXPLORATION ####
## Look at data
print(tapply(ARCdata$gain, ARCdata$Tank, mean, na.rm=TRUE))
print(tapply(ARCdata$gain, ARCdata$Diet, mean, na.rm=TRUE))
histogram(~ gain | Diet, data = ARCdata)
histogram(~ gain | Tank, data = ARCdata)
# number of observations by tank
ARCdata %>%
group_by(Tank) %>%
summarise(no_rows = length(Tank))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment