plot.window(...) 中的错误：使用 descdist 函数时需要有限的“xlim”值

Question

我正在处理一个包含 33 个变量和 1084 个观察值的数据集。现在我想建立一个模型。为此，我想使用 fitdistrplus 包中的 descist 函数来查看哪个分布最适合我的模型。 我的因变量如下所示：

1   15.0
2   5.0
3   5.0
4   NA
5   NA
6   100.0
7   5.0
8   10.0
9   50.0
10  240.0
11  20.5
12  NA
13  0
14  NA
15  60.0
16  NA
17  5.0
18  5.0
19  25.0
20  20.0
21  50.0
22  75.0
23  NA
24  10.0

我正在运行以下命令：

library(fitdistrplus)
library(MASS)
library(logspline)
descdist(mydata$var1, discrete = FALSE)

但我收到此错误消息并且没有情节：

Error in plot.window(...) : need finite 'xlim' values

我确实有几个缺失值 (NA) 和一些零 (0)，但还有很多其他观察结果不是 0 或缺失。我没有负值。 有人可以帮助我理解错误消息并可能帮助我解决问题吗？

那太好了！

最好的 安卡

Answer 1

您需要使用上面的表格删除缺失的数据：

mydata = structure(list(var1 = c(15, 5, 5, NA, NA, 100, 5, 10, 50, 240, 
20.5, NA, 0, NA, 60, NA, 5, 5, 25, 20, 50, 75, NA, 10)), class = "data.frame", row.names = c(NA, 
24L))

library(fitdistrplus)
library(MASS)
library(logspline)
descdist(mydata$var1, discrete = FALSE)

我得到错误：

Error in plot.window(...) : need finite 'xlim' values

试试这个：

descdist(mydata[complete.cases(mydata),"var1"], discrete = FALSE)