r语言中,spei能直接计算的年值吗

Python024

r语言中,spei能直接计算的年值吗,第1张

如果是统计整个矩阵,比如y

那么直接table(y)

>y

[,1] [,2] [,3]

[1,]123

[2,]222

[3,]321

>table(y)

y

1 2 3

2 5 2

所以如果统计所有列,那就自己写个循环吧

for(i in 1:ncol(y)){

table(y[,i])

}

如果只是画图,用curve()函数就好了

画正态密度:curve(dnorm,xlim=c(-3,3),col=2)

xlim是控制x轴显示从哪儿到哪儿,col是控制曲线颜色

画指数密度:curve(dexp(x,rate=1),xlim=c(0,5))

画指数分布:curve(pexp(x,rate=1),xlim=c(0,5))

你的方法是生成很多点x=seq(-6,6,0.1)

逐一算出函数值

t1[[i]]=dnorm(x,u[i],sigma)

t2[[i]]=pnorm(x,u[i],sigma)

最后在plot出来,用type="l"和lty=2的虚线弄出来。

curve这些功能都可以做到。

curve(dexp(x,rate=1),xlim=c(0,5),lty=2,add=T)就有虚线,

add=T可以一图多线