0

我正在尝试使用 EnvStat 包和循环来一次计算多个列来拟合 R 中的概率分布。列的长度不同，并且发生了一些代码错误。数据框不会保留为数字格式。错误消息：“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")
4

1 回答 1

0

您的代码有两个主要变化：首先，使用列表而不是数据框（这样您就可以容纳不相等的向量长度）：

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
于 2022-01-07T23:28:46.517 回答