search for: juvfor

...ral Transect Pole 1-2 Quadrat 1,t,Montastrea,0,HP+MEM I need to breakdown the data before I can plot it, so this is the stepwise code I use: Juv=read.csv("juvenile_density_20110506.csv") Juv$count[Juv$count<0]<- NA #Missing data is imputed as negative data, so first label it as NA JuvFor=subset(Juv,(location==" 2 Outer 10 m") | (location==" 1 Outer 10 m")) #Subset the data of interest Juv_Sum_by_quad=aggregate(count~year+location+quadrat_juvenile,data= JuvFor,sum) # Calculate sum of each quadrat Juv_Avg=aggregate(count~year+location,data=Juv_Sum_by_quad, mean) #...