#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)