BpHom100Eves <- read.delim("D:/ELTETATK/Idosor_gyak/BpHom100Eves.txt") attach(BpHom100Eves) ido=seq(1,36500) evnap=rep( seq(1,365), 100) negyedev=trunc(evnap/92)+1 plot(koz[1:(365*5)], type="l") acf(koz, lag=1000) #install.packages("GeneCycle") library(GeneCycle) plot(periodogram(koz)$freq, periodogram(koz)$spec, type="l") plot(periodogram(koz)$freq, periodogram(koz)$spec, type="l", xlim=c(0.002,0.003)) szezdummy<-lm(koz~ido+factor(evnap)) summary(szezdummy) plot(szezdummy$coefficients[3:366], type="l") #szezhettrend<-lm(koz~interaction(factor(negyedev)napok*factor(negyedev)+factor(evnap)) #summary(szezhettrend) szeztrig<-lm(koz~ido+cos(2*pi*ido/365)+sin(2*pi*ido/365)) summary(szeztrig) plot(periodogram(szeztrig$res)$freq, periodogram(szeztrig$res)$spec, type="l") plot(periodogram(szeztrig$res)$freq, periodogram(szeztrig$res)$spec, type="l", xlim=c(0.004, 0.006)) szeztrig3<-lm(koz~ido+cos(2*pi*ido/365)+sin(2*pi*ido/365)+cos(2*2*pi*ido/365)+sin(2*2*pi*ido/365)+cos(3*2*pi*ido/365)+sin(3*2*pi*ido/365)) summary(szeztrig3) plot(periodogram(szeztrig3$res)$freq, periodogram(szeztrig3$res)$spec, type="l") plot(ido[1:365], szezdummy$fitted.values[1:365], type="l") lines(szeztrig3$fitted.values[1:365], type="l", col=2) acf(szeztrig3$res) acf(koz) ar1fit<-arima(szeztrig3$res, order=c(1,0,0)) acf(ar1fit$res) ar2fit<-arima(szeztrig3$res, order=c(2,0,0)) qqline(ar2fit$res) shapiro.test(ar2fit$res[seq(10,36500,10)]) t.test(koz[1:18250], koz[18251:36500])