[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
入門(寫過其他程式，只是對語法不熟悉)
[問題敘述]:
想請教各位大大
目前我有一筆降雨資料(ncdf檔)，資料的維度為經度41格，緯度77格，總天數為9862天，
現在我想要將所有網格所有天數(41x77x9862)的降雨值出來畫累積分布圖(CDF)，我所用
的指令為ecdf()，但圖始終畫不出來，Console右上角的紅點出現了快半小時後消失且沒
有出現錯誤訊息，但畫圖區還是一片空白，不確定是不是要畫的值太多所導致還是另有
其他原因，所以想請教各位大大有什麼方法可以解決這個問題，謝謝。
降雨資料放置於此:http://0rz.tw/JI056
[程式範例]:
library(ncdf4)
library(data.table)
library(tidyverse)
memory.limit(size=50000)
#讀取資料
HIRAM_WRF_file <- "C:\\Users\\TOM\\Desktop\\197901-200512_pr_axis_time_domain.nc"
HIRAM_WRF_data <- nc_open(HIRAM_WRF_file)
print(HIRAM_WRF_data)
hiram_wrf_lon <- ncvar_get(HIRAM_WRF_data,"lon")
hiram_wrf_lat <- ncvar_get(HIRAM_WRF_data,"lat")
hiram_wrf_time <- ncvar_get(HIRAM_WRF_data,"time")
hiram_wrf_pr <- ncvar_get(HIRAM_WRF_data,"pr")
hiram_wrf_pr
dim(hiram_wrf_pr)
#畫CDF曲線
hiram_wrf_pr_cdf <- ecdf(hiram_wrf_pr[hiram_wrf_pr>=0])
summary(hiram_wrf_pr_cdf)
plot(hiram_wrf_pr_cdf,col="black",lwd=2
,xlim=c(0,2000),ylim=c(0,1),xlab="pr(mm)",ylab="CDF")
[環境敘述]:
[關鍵字]:
累積分布函數 CDF ncdf檔