tmpdata<-read.csv("HadCRUT.5.0.1.0.analysis.summary_series.180W_30S_180E_30N.monthly.csv",skip=1,header = FALSE) tmp_median_tropics<-ts(tmpdata$V2,start = 1850,frequency = 12) tmp_l95_tropics<-ts(tmpdata$V3,start = 1850,frequency = 12) tmp_u95_tropics<-ts(tmpdata$V4,start = 1850,frequency = 12) source('~/Temperature/iceemdan.R', echo=TRUE) imfs_median_tropics<-iceemdan(tmp_median_tropics) imfs_l95_tropics<-iceemdan(tmp_l95_tropics) imfs_u95_tropics<-iceemdan(tmp_u95_tropics) # source('~/Temperature/extract_spectra_temps_tropics.R', echo=TRUE) # plot(spectrum_median_tropics$instantfreq[,6]) # plot(spectrum_l95_tropics$instantfreq[,6]) # plot(spectrum_u95_tropics$instantfreq[,6]) # plot(abs(spectrum_u95_tropics$instantfreq[,6])) source('~/Temperature/extract_final_imfs_temp_tropics.R', echo=TRUE) plot(final_imfs_median_tropics) plot(final_imfs_l95_tropics) plot(final_imfs_u95_tropics) trend_median_tropics<-final_imfs_median_tropics[,8] trend_l95_tropics<-final_imfs_l95_tropics[,8] trend_u95_tropics<-final_imfs_u95_tropics[,8] ts.plot(trend_l95_tropics,trend_median_tropics,trend_u95_tropics,gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) ts.plot(final_imfs_l95_tropics[,6],final_imfs_median_tropics[,6],final_imfs_u95_tropics[,6],gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) ts.plot(final_imfs_l95_tropics[,7],final_imfs_median_tropics[,7],final_imfs_u95_tropics[,7],gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) ts.plot(final_imfs_l95_tropics[,8],final_imfs_median_tropics[,8],final_imfs_u95_tropics[,8],gpars=list(xlab="Time",ylab="Monthly Anomaly (ºC)",lty=c("dashed","solid","dashed"))) # plot(dtmp_median_tropics_dt<-ts(ddt(trend_median_tropics,1/12),start=1850,frequency=12),ylab="Rate of Change of Temperature Anomaly, ºC per Year") # plot(dtmp_l95_tropics_dt<-ts(ddt(trend_l95_tropics,1/12),start=1850,frequency=12),ylab="Rate of Change of Temperature Anomaly, ºC per Year") # plot(dtmp_u95_tropics_dt<-ts(ddt(trend_u95_tropics,1/12),start=1850,frequency=12),ylab="Rate of Change of Temperature Anomaly, ºC per Year") # plot(ts(rowSums(final_imfs_median_tropics[,5:8]),start = 1850,frequency = 12),ylab="Sum of IMFs 4-7 and Residual") source('~/Temperature/extract_final_spectra_temps_tropics.R', echo=TRUE) Measured=tmp_median_tropics decomposition_tropics=ts.union(Measured,final_imfs_median_tropics) colnames(decomposition_tropics)=c("Measured",paste0("IMF ",1:8),"Residual")