tmpdata<-read.csv("HadCRUT.5.0.1.0.analysis.summary_series.global.monthly.csv",header=FALSE,skip=1) tmp_median<-ts(tmpdata$V2,start = 1850,frequency = 12) tmp_l95<-ts(tmpdata$V3,start = 1850,frequency = 12) tmp_u95<-ts(tmpdata$V4,start = 1850,frequency = 12) source('~/SA/iceemdanP.R', echo=TRUE) imfs_median<-iceemdanP(tmp_median) imfs_l95<-iceemdanP(tmp_l95) imfs_u95<-iceemdanP(tmp_u95) # source('~/Temperature/extract_spectra_temps.R', echo=TRUE) # plot(spectrum_median$instantfreq[,6]) # plot(spectrum_l95$instantfreq[,6]) # plot(spectrum_u95$instantfreq[,6]) # plot(abs(spectrum_u95$instantfreq[,6])) source('~/Temperature/extract_final_imfs_temp.R', echo=TRUE) plot(final_imfs_median) plot(final_imfs_l95) plot(final_imfs_u95) trend_median<-final_imfs_median[,9] trend_l95<-final_imfs_l95[,9] trend_u95<-final_imfs_u95[,9] ts.plot(trend_l95,trend_median,trend_u95,gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) ts.plot(final_imfs_l95[,6],final_imfs_median[,6],final_imfs_u95[,6],gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) ts.plot(final_imfs_l95[,7],final_imfs_median[,7],imfs_u95[,7],gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) ts.plot(final_imfs_l95[,8],final_imfs_median[,8],final_imfs_u95[,8],gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) # plot(dtmp_median_dt<-ts(ddt(trend_median,1/12),start=1850,frequency=12),ylab="Rate of Change of Temperature Anomaly, ºC per Year") # plot(dtmp_l95_dt<-ts(ddt(trend_l95,1/12),start=1850,frequency=12),ylab="Rate of Change of Temperature Anomaly, ºC per Year") # plot(dtmp_u95_dt<-ts(ddt(trend_u95,1/12),start=1850,frequency=12),ylab="Rate of Change of Temperature Anomaly, ºC per Year") # plot(ts(rowSums(final_imfs_median[,5:8]),start = 1850,frequency = 12),ylab="Sum of IMFs 4-7 and Residual") source('~/Temperature/extract_final_spectra_temps.R', echo=TRUE) median_periods=1/final_spectrum_median$instantfreq[,4:8] median_periods[median_periods<0]=NA median_periods[median_periods[,1]>10,1]=NA median_periods[median_periods[,2]>20,2]=NA Measured=tmp_median_nh decomposition_global=ts.union(Measured,final_imfs_median_nh) colnames(decomposition_global)=c("Measured",paste0("IMF ",1:8),"Residual")