[问题类型]:
程式咨询(我想用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")
[环境叙述]:
https://imgur.com/psT3fMx
[关键字]:
累积分布函数 CDF ncdf档