I need to figure out how to get rid of the NA for mean and standard error
summary(Copy_of_US_Vcm_daily_aflx$season)
winter spring monsoon
1804 1830 1845
season.summary<-data.frame(Copy_of_US_Vcm_daily_aflx %>%
-
group_by(season) %>%
-
summarize(GPP_int.mean=mean(GPP_int), se=sd(GPP_int)/(n()^0.5),n=n()))
look at your summary data
season.summary
season GPP_int.mean se n
1 winter NaN NA 1804
2 spring NaN NA 1830
3 monsoon NaN NA 1845