#Odds: runs<-100 game<-function(){ results<-sample(c(-1,0,1,2),10,replace=T) return(sum(results)>5) } sum(replicate(runs,game()))/runs #Time Series Analysis require("astsa") plot(gtemp,xlab="Year",ylab="Deviation",main="Global Temperature Deviations,1880-2015") acf2(gtemp) plot(diff(gtemp)) acf2(diff(gtemp)) sarima(gtemp,0,1,1) sarima.for(gtemp,5,0,1,1)