2

我正在尝试拟合一条曲线来模拟直邮活动随时间推移的响应。使用 R，我能够使用 fitdistr() 函数获得形状和比例因子。然后我在 weibull() 函数中使用形状和比例作为参数。但是，我们的广告系列通常会持续 63 天（8 周），并且拟合的威布尔曲线的“长度”会过早地被截断。有没有办法设置“长度”？

...或者有没有更好的方法来模拟 Direct Mail 营销活动的响应？？？

谢谢！

set.seed(5)
install.packages("MASS")
library("MASS")
responses <-c(4,5,1,12,24,16,16,15,5,18,7,12,5,13,6,2,9,2,5,1,4,4,5,3,3,4,7,3,9,2,2,4,3,2,5,4,3,2,2,2,2,2,2,2,2,2,1,1,1,1,1,2,3,2,1,1,1,1,1,1,1,1,1)                                                                                        

f <- fitdistr(responses,'weibull')

f #check the shape and scale

#plug in the shape and scale. 284 is the number of total responders that we're trying to fit the curve to. 
weibulldraws <- as.data.frame(table(round(.5 +   rweibull(284,1.0753863,4.6579543))))
weibulldraws
4

1 回答 1

0

当您用 对向量进行制表时table()，只有向量中出现的值才会出现在结果中；如果要包含其他值，则需要创建一个因子，其中包含这些值作为级别。因此，如果您想在最后包含所有零，请在制表之前将您的结果转换为具有适当水平的因子：

set.seed(5)
wshape <- 1.0753863
wscale <- 4.6579543
n <- 284
rvals <- round(rweibull(n,wshape,wscale)+0.5)
frvals <- factor(rvals,levels=0:63)
weibulldraws <- as.data.frame(table(frvals))

如果您想绘制与 Weibull 相关的理论曲线（即密度曲线而不是对应于从分布中随机抽取的表格），则使用curve()以下dweibull()函数：

curve(n*dweibull(x,wshape,wscale),from=0,to=63)

（用于add=TRUE将曲线添加到现有绘图）。

于 2015-03-30T17:08:44.770 回答