R中不同长度EnvStat循环数据帧概率分布的拟合

Question

我试图在R中拟合概率分布，使用EnvStat包和循环一次计算多个列。列有不同的长度，并且正在发生一些代码错误。数据帧不保持数字格式。错误信息：'x‘必须是一个数字向量，我无法识别错误。有人能帮忙吗？非常感谢

遵循代码：

x = runif(n = 50, min = 1, max = 12)
y = runif(n = 70, min = 5, max = 15)
z = runif(n = 35, min = 1, max = 10)
m = runif(n = 80, min = 6, max = 18)

length(x) = length(m)
length(y) = length(m)
length(z) = length(m)

df = data.frame(x=x,y=y,z=z,m=m)
df

library(EnvStats)

nproc = 4
cont = 1
dfr = data.frame(variavel = character(nproc), 
                 locationevd= (nproc), scaleevd= (nproc),
                 stringsAsFactors = F)


# i = 2

for (i in 1:4) {
  
  print(i) 
  
  nome.var=colnames(df)
  
  df = df[,c(i)]
  df = na.omit(df)
  
  variavela = nome.var[i]
  dfr$variavel[cont] = variavela
  
  evd = eevd(df);evd
  locationevd = evd$parameters[[1]]
  dfr$locationevd[cont] = locationevd
  scaleevd = evd$parameters[[2]]
  dfr$scaleevd[cont] = scaleevd

  cont = cont + 1
  
}

writexl::write_xlsx(dfr, path = "Results.xls")

Answer 1

代码的两个主要更改:第一，使用一个列表而不是一个dataframe (因此您可以容纳不等的向量长度)：

x = runif(n = 50, min = 1, max = 12)
y = runif(n = 70, min = 5, max = 15)
z = runif(n = 35, min = 1, max = 10)
m = runif(n = 80, min = 6, max = 18)


vl = list(x=x,y=y,z=z,m=m)
vl
if (!require(EnvStats){ install.packages('EnvStats'); library(EnvStats)}
   
nproc = 4
# cont = 1   Not used
dfr = data.frame(variavel = character(nproc), 
                 locationevd= (nproc), scaleevd= (nproc),
                 stringsAsFactors = F)

第二:使用一个循环索引，而不是使用"cont“索引。

for ( i in 1:length(vl) ) {
    # print(i)    Not needed
    nome.var=names(vl)  # probably should have been done before loop
    var = vl[[i]]
    variavela = nome.var[i]
    dfr$variavel[i] = variavela  # all those could have been one step
    
    evd = eevd( vl[[i]] )   # ;evd
    locationevd = evd$parameters[[1]]
    dfr$locationevd[i] = locationevd
    scaleevd = evd$parameters[[2]]
    dfr$scaleevd[i] = scaleevd
}

这将为您提供所需的结构：

 dfr
  variavel locationevd scaleevd
1        x    5.469831 2.861025
2        y    7.931819 2.506236
3        z    3.519528 2.040744
4        m   10.591660 3.223352