Fitdistrplus不支持的分发错误
Fitdistrplus是一个R语言包，用于拟合和评估各种概率分布模型。它提供了一种方便的方式来拟合数据并估计分布参数。然而，有时在使用Fitdistrplus时可能会遇到不支持的分布错误。
不支持的分布错误意味着Fitdistrplus无法拟合或评估所选择的特定概率分布模型。这可能是由于数据不适合所选择的分布，或者Fitdistrplus包本身不支持该特定分布。
在遇到不支持的分布错误时，可以考虑以下解决方案：
- 检查数据的适应性：首先，确保所选择的分布模型适用于数据。可以通过绘制数据的直方图、密度图或Q-Q图来检查数据的分布特征。如果数据不符合所选择的分布模型的假设，那么Fitdistrplus可能无法成功拟合。
- 尝试其他分布模型：如果所选择的分布模型不适合数据，可以尝试使用其他分布模型进行拟合。Fitdistrplus支持多种常见的概率分布模型，如正态分布、指数分布、伽玛分布等。可以根据数据的特征选择适当的分布模型进行尝试。
- 使用其他R包：如果Fitdistrplus无法支持所需的分布模型，可以尝试使用其他R语言包来进行拟合和评估。R语言生态系统中有许多其他包可用于概率分布拟合，如fitdistr、actuar、VGAM等。
- 自定义拟合算法：如果没有现成的R包支持所需的分布模型，还可以考虑自定义拟合算法。这需要对概率分布的数学特性和参数估计方法有一定的了解，并使用R语言编写自己的拟合函数。
总之，当遇到Fitdistrplus不支持的分发错误时，需要检查数据的适应性，尝试其他分布模型，使用其他R包或自定义拟合算法来解决该问题。
相关·内容
我正在尝试将一个发行版与我的数据(ssdist；ssi是值)相匹配，但是我一直被抛出下面的错误代码，我不知道为什么。值得注意的是，错误只抛给pois和nbinom，它们是我认为最适合数据的发行版。
错误：：fitdist(ss$ssiss$name == s，distr = "nbinom")：函数mle无法估计参数，错误代码为100
这是我的代码：
ssdist <- read_excel(here("Data/Raw/ssdist.xlsx")) %>%
verify(has_all_names("name", "s
浏览 5提问于2022-05-19得票数 0
回答已采纳
3回答
我正在尝试从复制代码，但是我在这样做时遇到了问题。我正在使用来自package VGAM和fitdistrplus的gumbel发行版。当这样做时，问题就出现了：
fit = fitdist(data1, 'gumbel', start = list(location = 0, scale = 1))
Error in mledist(data, distname, start, fix.arg, ...) :
'start' must specify names which are arguments to 'distr'.
好像lo
浏览 2提问于2016-03-18得票数 4
回答已采纳
通过运行以下代码： set.seed(12)
eta <- 10
nu <- 2
n <- 100
t <- eta*(((-1)*log(1-runif(n)))^(1/nu))# generates T_i
t
library("fitdistrplus")
fitll <- fitdist(t, "llogis") 在R中，我收到了这个错误： Error in fitdist(t, "llogis") : The dllogis function must be defined 我之前安装过包"f
浏览 31提问于2020-10-07得票数 2
回答已采纳
1回答
我试图找到的发行版。我试过用fitdistrplus包
data <- data.matrix(Book1)
descdist(data, discrete = FALSE)
但是得到这个错误：
下降错误(数据，离散=假)：数据必须是数字向量
浏览 3提问于2016-01-26得票数 1
我想在观测数据的直方图上绘制一组观测得到的伽马密度函数。我能够为伽马拟合生成直方图和参数估计。正在对来自主数据集的多个数据子集执行此操作。如何在这个循环中创建的每个直方图上绘制伽马密度函数？
我目前有：
library(MASS)
species <- c('acesac', 'acesac', 'acesac', 'acesac', 'acesac', 'acesac',
'polbif', 'polbif', 'polbif', '
浏览 3提问于2015-10-09得票数 1
回答已采纳
我正在尝试拟合我的数据的多个分布，以找到最适合的，从库( fit Fitdistrplus)中找到最佳匹配。我使用库中的dgamma3 (FAdist)。如何为每个分发设置start参数？ library(FAdist)
library(fitdistrplus)
x <- rnorm(30,5,1)
fitdist(x, distr = 'gamma3',start = list(meanlog=0,sdlog=1)) 下面是错误： $start.arg
$start.arg$meanlog
[1] 0
$start.arg$sdlog
[1] 1
$fix.a
浏览 27提问于2021-07-14得票数 0
如何使用fitdistrplus来估计t分布的位置和尺度参数？我知道我需要提供初始值(总体上它工作得很好)，但在这个包中只允许使用df。你有什么解决方案吗？
非常感谢。
浏览 44提问于2016-07-30得票数 2
我试图通过使用来自fitdist包的fitdistrplus函数将我的数据与瑞利分布相匹配。
x <- c(19.000000,23.000000,26.000000,45.000000,8.050000,46.900000,1.268333,30.000000,
1.466667,3.733333,1.683333,4.000000,3.950000,1.850000,42.000000,1.333333,
1.550000,1.000000,2.066667,1.566667,1.216667,1.850000,1.400000,8.366667,
浏览 6提问于2021-01-27得票数 0
回答已采纳
我试图找出Beta分布的参数。它是一个更大的循环的一部分(在Stackoverflow中找到的)，它可以很好地适应其他几个发行版，但是，这个特别的循环给我带来了困难。问题在于适配函数。
我没有零值。在stats包的文档中，它指出使用参数ncp=0可以帮助解决太接近于零的值问题。我试过了但也没用。
我做错了什么？
Tks，
P.D .我用其他方法估计了参数。问题是，这个功能更方便、更灵活，因为它是功能的一部分。
...
else if(distrib[[i]] == "beta"){
gf_shape = "beta"
fd_b <
浏览 1提问于2015-06-17得票数 2
回答已采纳
1回答
以下代码每次运行时和在两个不同的R版本上都会崩溃。我敢打赌我做错了什么(也许我无意中把优化工作推到了它不能去的地方？)但我看不出是什么。我在这里的想法是通过提供几乎好的a、mean和sd来最大限度地帮助MLE。如您所见，我生成的数据样本非常接近拟合参数(故意)。
library(truncnorm)
library(fitdistrplus)
fitdist(rtruncnorm(100, a=8, mean=10, sd=1),
"truncnorm", fix.arg=list(a=8),
start = list(mean = 10+0.5
浏览 19提问于2022-04-27得票数 1
回答已采纳
我试图覆盖一个威布尔概率密度函数(PDF)的右删失数据的直方图上的数据使用包fitdistrplus，但一直无法这样做。
df是一个包含虚拟数据集的data.frame，其格式为我的更大的数据集。
> library(fitdistrplus)
> left <- c(rep(1,12),rep(2,5),rep(3,3),rep(4,2),rep(5,1))
> right <- c(rep(1,12),rep(2,5),rep(3,3),rep(4,2),NA)
> (df <- data.frame(left,right)
浏览 4提问于2022-08-07得票数 0
尝试在设备上调试应用程序时，出现错误“获取进程12567的任务失败”。只有当我的应用程序尝试在相同的设备和系统上运行任何其他应用程序时，才会发生这种情况，它运行良好。不知道会出什么问题。在我的应用程序中，我使用了zbar sdk和tapku库框架。请指导，如果需要任何其他信息，请随时询问。谢谢。
浏览 0提问于2013-04-26得票数 0
我想对我拥有的一些数据点进行对数皮尔逊III拟合。然而，每次我尝试这样做时，我都会收到错误消息，我真的不知道该怎么办。也许我应该补充说，我只是在几天前才开始使用R，所以，我不是它的专家。
重要的代码部分，没有导入内容的部分等等如下：
pIIIpars<-list(shape=1, location=1, scale=1)
dPIII<-function(x, shape, location, scale) PearsonDS::dpearsonIII(x, shape=1, location=1, scale=1, params=pIIIpars, log=FALSE)
pP
浏览 1提问于2013-02-21得票数 2
当我输入代码plot(x.logis)时，输出是四张图片一起显示，但我只是在寻找其中之一。我怎么把它们分别画出来？
这是我的代码：
x<-c(2800,3260,66.5,195,420,840,1380,469,260,50,209.8,370,27,420,157)
y<-log10(x)
# load fitdistrplus package for using fitdist function
library(fitdistrplus)
# fit logistic distribution using MLE method
x.logis <- fitdist(y,
浏览 4提问于2016-06-11得票数 2
回答已采纳
这是我的代码：
x<-c(1300,541,441,35,278,167,276,159,126,60.8,160,5000,718,9740,135,3480,264.6,379,170,251.3,155.84,187.01,850)
y<-log10(x)
# load fitdistrplus package for using fitdist function
library(fitdistrplus)
# fit logistic distribution using MLE method
x.logis <- fitdist(y, "logis"
浏览 5提问于2016-10-12得票数 0
回答已采纳
我想执行一个符合我在R中的一些数据点的日志Pearson，我遵循了这个准则，。但是当我的偏斜度(g)是负的(当然参数“比例”也是负的，因为在计算标度时“符号(G)”)，我会遇到一个问题。来自"FAdist“的发行版不适用于负刻度参数，我需要这个值作为"fitdist”(从fitdistrplus)的起始值。有些页面说，参数的形状和规模只有在pearson III (或概括伽玛分布)中是正的，而其他没有，我的想法用完了。守则是：
library(fitdistrplus)
library(FAdist)
library(e1071)
#data
df <-(92.8, 5
浏览 0提问于2020-03-06得票数 0
一旦我厌倦了构建我的项目，.I在下面的代码中得到了构建错误。
创建过程dlm.uspScriptTableInsertionForADW
(
---一些代码
创建表#selectColListTable
WITH
(
DISTRIBUTION = ROUND_ROBIN
)
错误：
SQL80001:WITH‘.Excepting '(’，或AS. )附近的语法不正确。
浏览 0提问于2017-04-17得票数 1
我正在尝试使用来自fitdist包的fitdistrplus函数来使我的数据适合Burr发行版。
library(fitdistrplus)
library(actuar)
fit_bur <- fitdist(data=c(6,8,6,3,3,6,6,6,7,8,4,5,1,4,4,0,0,0,0,7,5,6,6,5,3,4,8,6,5,6,3,6,5,4,3,9,4,2,4), method = "mle",distr="burr",
start = list(shape1 = 0.3, shape2 = 1
浏览 0提问于2018-06-21得票数 0
回答已采纳
我已经在我的Windows 7计算机上安装了cassandra的最新版本，目录如下：
/cygdrive/c/Development/Cassandra/apache-cassandra-1.1.4
在/var/lib/cassandra下为data_file_directories、saved_caches_direcory和commitlog_directory创建了所有必需的目录，但当我尝试使用以下方法启动它时：
bin/cassandra -f
我得到以下错误：
$ bin/cassandra -f xss = -ea -javaagent:bin/../lib/jamm-0.2.5
浏览 2提问于2012-09-10得票数 4
回答已采纳
我现在使用包fitdistrplus来构造Gamma分布，我的问题是如何提取拟合的值来计算根均方误差？谢谢你的帮助。
library(fitdistrplus)
Sev = c(1.42,5.15,2.5,2.29,12.36,2.82,1.4,3.53,1.17,1.0,4.03,5.26,1.65,1.41,3.75,1.09,
3.44,1.36,1.19,4.76,5.58,1.23,2.29,7.71,1.12,1.26,2.78,1.13,3.87,15.43,1.19,
4.95,7.69,1.17,3.27,1.44,1.05,3.94,1.58,2.29,2.
浏览 0提问于2018-12-22得票数 1
回答已采纳
1回答
在R中绘制对数似然函数
、、、
嗨，我在我的代码中收到了x和y长度不同的错误。datasim是一个大小为1000的模拟样本。请帮我弄一下这个。
x <- datasim
loglik <- function(theta){
k<- theta[1]
lambda<- theta[2]
out <- sum(dweibull(x,shape = k, scale=lambda, log = TRUE) )
return(out)
}
theta<- c(0.5,1.5)
plot(theta, loglik(theta), type="l", lwd=3,
浏览 50提问于2020-05-23得票数 1
回答已采纳
这是为我崩溃R的代码，没有任何有用的错误消息。有什么想法吗？我在这快疯了。谢谢,
library(fitdistrplus)
library(truncnorm)
set.seed(0)
x = rtruncnorm(n=30, a=0, b=Inf, mean=1, sd=0.45)
my_fit = fitdist(x, "truncnorm", method='mle', fix.arg=list(a=0, b=Inf),
start=list(mean=mean(x), sd=sd(x)))
浏览 9提问于2020-08-16得票数 2
回答已采纳
2回答
library(VGAM)
library(fitdistrplus)
fitdist(u_NI$k_u, 'truncpareto',
start = list(lower=1,
upper=42016,
shape=1)) -> fit.k_u
length(u_NI$k_u) = 637594
我发现了一个错误：
<simpleError in optim(par = vstart, fn = fnobj, fix.arg = fix.arg, obs
浏览 10提问于2022-02-05得票数 0
回答已采纳
我正在使用来自fitdistrplus包的某些函数作为我正在创建的包的一部分。
当函数运行时，我试图阻止任何错误消息在控制台中显示，而是希望将错误消息记录在我正在创建的错误日志中。
在大多数情况下，使用tryCatch()使我能够做到这一点。
但是对于fitdist()函数，即使消息被写入错误日志(意味着tryCatch()表达式正在运行)，我也不能阻止在控制台中打印的错误消息。
我在下面的代码中复制了我的问题。
library(fitdistrplus)
file.create("error_log.txt")
func_desc<-function(x){
t
浏览 2提问于2018-07-02得票数 2
回答已采纳
我试着在metRology包的fitdist函数中使用t-scaled函数。然而，即使我正在运行这两个包，我仍然收到错误‘必须定义dt缩放函数’。
我不知道现在该怎么办。
任何帮助都将不胜感激！
library(metRology)
library(fitdistrplus)
z = rnorm(1000, mean = 0, sd = 1)
x = fitdist('z', "t-scaled", start = list(df = 3, mean = mean(z), sd = sd(z)) )
浏览 0提问于2018-05-28得票数 0
3回答
我从R中的fitdist()包中知道了fitdistrplus函数，但是我不能用它来预测高斯分布。我可以预测正常，逻辑，威布尔等。我如何使用它的高斯？还有其他方法可以预测吗？
浏览 0提问于2016-11-08得票数 -1
回答已采纳
1回答
请问如何解决这个错误？这个版本还有别的包吗？
install.packages("fitdistrplus"); library(fitdistrplus)
Installing package into ‘/home/linux/R/x86_64-pc-linux-gnu-library/3.4’
(as ‘lib’ is unspecified)
Warning in install.packages :
package ‘fitdistrplus’ is not available (for R version 3.4.4)
Error in library(fitd
浏览 1提问于2021-10-14得票数 0
我想定义我自己的分布函数，以便在R中使用fitdist或fitdistr函数，例如在fitdistrplus包中使用fitdist。我定义了一个名为sgamma的定制发行版，如下所示：
dsgamma<-function(x,shape){return(dgamma(x,shape,scale=1));}
qsgamma<-function(p,shape){return(qgamma(p,shape,scale=1));}
psgamma<-function(q,shape){return(pgamma(q,shape,scale=1));}
rsgamma<-fun
浏览 2提问于2014-07-24得票数 5
回答已采纳
我使用fitdistrplus软件包中的fitdist函数来获得最优分布，并绘制ppcomp图，我使用?fitdistrplus的示例代码来替换我的数据和代码。
data(groundbeef)
serving <- groundbeef$serving
fitW <- fitdist(serving, "weibull")
fitg <- fitdist(serving, "gamma")
fitln <- fitdist(serving, "lnorm")
ppcomp(list(fitW, fitg, fitln),
浏览 2提问于2016-07-04得票数 3
我有一个第三方sdk作为.dll和tlb提供。它适用于WPF，但是当我将它添加到UWP项目时，一旦我到达使用sdk的屏幕，它就会抛出。
使用CLSID {DBA84AA0-3342-4AF0-96ED-3D8D08D544F3}的COM组件实例由于以下错误而失败: 80040154类未注册(HRESULT异常: 0x80040154 (REGDB_E_CLASSNOTREG))。请确保COM对象在允许的CoCreateInstanceFromApp列表中。
是否存在.dll或tbl，与我正在开发的UWP应用程序一起工作。我有x64和x86版本的.dll。
浏览 1提问于2019-03-22得票数 0
1回答
为了为我的数据集拟合几种概率分布，我尝试了使用car和MASS包的几种分布，例如对于泊松分布：
require(car); require (MASS)
#Data set
X1 <- rnorm(100, 175, 7)
Y1 <- rbinom(100,1, 0.35)
##Fitting Poisson:
poisson <- fitdistr(Y1, "Poisson")
qqp(Y1, "pois", poisson$estimate)
但是，如果我的选项是二项式分发，MASS包不支持这种分发，那么我使用fitdistrplus包：
浏览 3提问于2017-01-27得票数 1
回答已采纳
我试图在R中编写一些代码，从Excel中获取样本数据，并确定最适合于数据的分布，然后是所述分布的参数。在Googling上搜索了一段时间后，我决定尝试fitdistrplus来拟合发行版，并发现gofstat是一个可以用来检查是否合适的函数。我想比较GOF的统计数据，找出循环中最合适的分布。
我的代码的最初部分实际上只是从Excel导入我的示例数据(我在Excel中创建了一个1000个值，通常分布在Excel中，并将其保存为CSV格式的单个列)，并试图将其与一个发行版相匹配并绘制结果。
library(fitdistrplus)
testData = read.table("C:\\U
浏览 2提问于2015-03-18得票数 2
回答已采纳
1回答
所以，我有一个离散的数据集，my_dat，我试图拟合一条曲线，以便能够生成基于my_dat的随机变量。我在连续数据上使用fitdistrplus取得了很大的成功，但是尝试将它用于离散数据时有很多错误。
表设置：
library(fitdistrplus)
my_dat <- c(2,5,3,3,3,1,1,2,4,6,
3,2,2,8,3,4,3,3,4,4,
2,1,5,3,1,2,2,4,3,4,
2,4,1,6,2,3,2,1,2,4,
5,1,2,3,2)
我首先看一下数据的直方图：
浏览 8提问于2017-10-19得票数 2
2回答
假设我有一个变量x，它是使用以下方法生成的：
x <- rgamma(100,2,11) + rnorm(100,0,.01) #gamma distr + some gaussian noise
head(x,20)
[1] 0.35135058 0.12784251 0.23770365 0.13095612 0.18796901 0.18251968
[7] 0.20506117 0.25298286 0.11888596 0.07953969 0.09763770 0.28698417
[13] 0.07647302 0.17489578 0.02594517 0.
浏览 27提问于2017-08-07得票数 7
回答已采纳
1回答
为什么我得到线而不是标准气泡在我的q-q图？
我的代码：
data <- read.csv("C:\\Users\\anton\\SanFrancisco.csv")
x <- data$ï..San.Francisco
head(x)
library("fitdistrplus")
fitnor <- fitdist(x, "norm")
fitlogis <- fitdist(x, "logis")
qqcomp(list(fitnor, fitlogis), legendtext=c(&#
浏览 4提问于2020-10-13得票数 1
回答已采纳
用户问我如何才能批量使用InAppPurchase。 我的应用程序是免费的，所以他们可以通过以下链接获得大量应用程序：AppStore Volume 但他们不能获得所有设备的InAppPurchase ...有没有办法批量购买IAP？ 做这件事最好的方法是什么？
浏览 32提问于2020-12-22得票数 0
3回答
我正试图通过Jenkyns和Crashlytics自动化我的应用程序构建发行版。我按照文档中的指示进行了配置，但crashlyticsUploadDistributionRelease失败了，出现了以下错误：“分发上传失败”和“警告- Crashlytics停止编译，因为它无法分发未签名的APK:”。
我怎么才能修好它？
浏览 11提问于2016-06-22得票数 10
回答已采纳
1回答
我试图在MC55上运行MS 2003项目。我发现这个错误在下面-
有部署错误。继续吗？
MS 2003支持MC55吗？操作如何？请和我分享一下。
浏览 1提问于2011-08-15得票数 1
1回答
android原生有没有API可以使用路线图API来获取平滑的方向路线。我在谷歌地图中找到了JavaScript ，但对谷歌安卓地图库不太确定
浏览 17提问于2015-05-08得票数 3
广播时，星火可能会因错误org.apache.spark.sql.errors.QueryExecutionErrors#notEnoughMemoryToBuildAndBroadcastTableError而失败(Spark3.2.1)：
为什么BroadcastExchange需要更多的驱动程序内存？广播不是向所有驱动程序发送数据吗？为什么驱动程序内存是瓶颈？
谢谢。
浏览 7提问于2022-10-20得票数 1
回答已采纳
我想在单个图中比较不同分布与我的数据的拟合程度。fitdistrplus包中的qqcomp函数几乎完成了我想要做的事情。然而，我唯一的问题是，它主要是使用基本R图编写的，而我所有的其他图都是用ggplot2编写的。我基本上只想自定义qqcomp图，使其看起来像是在ggplot2中创建的。 从文档(https://www.rdocumentation.org/packages/fitdistrplus/versions/1.0-14/topics/graphcomp)中我了解到，通过设置plotstyle="ggplot"，这是完全可能的。然而，如果我这样做，绘图上没有显示任何
浏览 47提问于2019-12-30得票数 3
iOS企业应用程序能够使用清单文件在web上分发。tvOS可以以类似的方式分发吗？我知道tvOS没有网络浏览器，但我想知道是否有办法将文件空投到苹果电视上，以安装内部tvOS应用程序？
到目前为止，唯一的方法是使用MDM服务器或使用Apple Configurator通过Mac拖放IPA。
浏览 6提问于2018-08-23得票数 5
我已经在voro++上下载了最新版本，然后我在“包含目录”中添加了voro++源路径。但是，发生了以下错误：
LNK2019未解析的外部符号"public：__thiscall voro::voronoicell_base::voronoicell_base(void)“(??0voronoicell_base@voro@@QAE@XZ)在函数"public：__thiscall voro::voronoicell::voronoicell(void)”(？？0VOLONOICELL@voro@@QAE@XZ)中引用
我该如何解决这个问题？
浏览 1提问于2017-05-03得票数 1
提前为我糟糕的英语道歉:p
我试着用一个我称之为“短信”的新频道来传播竞选活动，
我在通道类型选择列表中添加了SMS选项，
并更新了boutton (在垃圾上)的分发操作，并使它调用一个新的javascript函数(ActionToDistributeSMS)，如果频道是短信，
我开发了一个插件(C#)来开发CampaignActivity
ActionToDistributeSMS()函数将显示插件
..。(DistributeCampaignActivityResponse)service.Execute(distrubuteCampaign)；//插件的代码.
我
浏览 0提问于2014-06-05得票数 0
在安装时，您可以获得Visual模板，用于创建“使用代码修复(NuGet +VSIX)的分析器”和“代码重构(VSIX)”。
为什么没有用NuGet传送代码重构的模板(比如“代码重构(NuGet + VSIX)")？
是否可以从NuGet包中添加代码重构提供程序(通过手动创建NuGet包)？
在每个项目的基础上安装代码重构提供程序不是很有意义吗？
编辑:我的用例是一种代码重构，它从非类型化数据(如JSON)中生成C#类型。这基本上就是F#类型提供者的特性(当然，它们的功能要强大得多)。这意味着用户创建一个类，如下所示
[JsonProvider("http://example.
浏览 4提问于2015-10-14得票数 3
回答已采纳
我正处于设置W3总缓存的最后阶段。在配置W3TC CDN时，当我单击Test S3 Upload和Cloudfront分发按钮时，会得到此错误：
Error: Distribution for origin "bucket-name.s3.amazonaws.com" not found.
我只有一个发行版，这是我的政策：
"Effect": "Allow",
"Action": "cloudfront:ListDistributions",
"Resource": "*"
我认为造
浏览 3提问于2022-03-25得票数 0
有没有一种方法可以在Google Cloud中为API端点设置适当的监控和警报？我可以在这里看到中位数和百分位数延迟：，但我想为它设置警报。例如，如果延迟中值在1小时内超过100ms，我将收到一封电子邮件。我也一直在尝试Stackdriver，甚至是高级版本，但没有成功。我已经能够绘制延迟，但不能设置警报。
有什么建议吗？
编辑:我刚刚找到了Apigee选项来部署包含监控的API。你有使用它的经验吗？
浏览 1提问于2018-01-26得票数 3
我想要定义我自己的发行版，用它来拟合我的每月降水数据，从现在起称为“月”。我正在使用“lmomco”函数来帮助我定义发行版，但无法使它工作。例如，我定义了广义极值(gev)分布，如下所示：
dgev<-pdfgev #functions which are included in lmomco
pgev<-cdfgev
qgev<-quagev
由于"fitdistrplus“需要参数"start”(由所需分布的初始参数值组成)，因此我估计这些初始值如下：
lmom=lmoms(month,nmom=5) #from lmomco package
浏览 8提问于2015-04-27得票数 7
回答已采纳
2回答
以下是我的代码：
library(fitdistrplus)
s <- c(11, 4, 2, 9, 3, 1, 2, 2, 3, 2, 2, 5, 8,3, 15, 3, 9, 22, 0, 4, 10, 1, 9, 10, 11,
2, 8, 2, 6, 0, 15, 0 , 2, 11, 0, 6, 3, 5, 0, 7, 6, 0, 7, 1, 0, 6, 4, 1, 3, 5,
2, 6, 0, 10, 6, 4, 1, 17, 0, 1, 0, 6, 6, 1, 5, 4, 8, 0, 1, 1, 5, 15, 14, 8, 1,
3, 2, 9, 4, 4, 1,
浏览 1提问于2018-08-25得票数 6
1回答
我有一个计数数据的向量，它是强烈过度分散和零膨胀的。
向量看起来像这样：
i.vec=c(0,63,1,4,1,44,2,2,1,0,1,0,0,0,0,1,0,0,3,0,0,2,0,0,0,0,0,2,0,0,0,0,
0,0,0,0,0,0,0,0,6,1,11,1,1,0,0,0,2)
m=mean(i.vec)
# 3.040816
sig=sd(i.vec)
# 10.86078
我想要拟合一个分布，我强烈怀疑它将是一个零膨胀泊松(ZIP)。但我需要执行一个显着性测试，以证明ZIP分布符合数据。
如果我有一个正态分布，我可以使用vcd包中的goodfit()函数进行卡方拟合优度测试
浏览 6提问于2011-08-23得票数 12
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
