模拟从 R 中截断的逻辑分布生成随机数答案

【问题标题】：Simulation to generate random numbers from a truncated logistic distribution in R模拟从 R 中截断的逻辑分布生成随机数
【发布时间】：2014-03-27 10:58:53
【问题描述】：

我编写了一个代码来在 R 中生成逻辑随机变量，但截断值介于 0.2 和 0.5 之间。那里的任何人都可以帮助正确编写它，尤其是在 if 条件下。谢谢。

<##The Logistic distribution##
ns=10
for(i in 1:100){
rtruncl=function(u,a,b,alpha1,alpha2)
#error: {
qalpha1 = log(alpha1/(1-alpha1))
qalpha2 = log((1-alpha2)/alpha2)
X=a+b*log(u/(1-u))
u=runif(ns)
if((X >= qalpha1) && (X <= qalpha2)){
# error: }
a=0; b=2 ##a is a location parameter and b is a scale parameter##
t1=rtruncl(u,a,b, 0.2,0.5)
t1
}
}  #closed the for loop

【问题讨论】：

您能修正一下格式吗？
@krlmlr 希望现在没事。
非常感谢任何其他方式。
如何将rlogis() 与拒绝抽样一起使用？（假设“截断”意味着您只需删除这些值。）
如这里所写，有不平衡的括号，语句 if((X >= qalpha1) && (X

标签： r simulation

【解决方案1】：

您可以使用 R 中的内置逻辑分布函数（例如，rlogis(...)，以及 truncdist 包。

library(truncdist)
location <- 0
scale    <- 2
X <- rlogis(1000, location=location, scale=scale)
Y <- rtrunc(1000, spec="logis", a=-5,b=5, location=location, scale=scale)
par(mfrow=c(1,2))
hist(X, xlim=c(-10,10), breaks=30, freq=F, col="lightgreen")
hist(Y, xlim=c(-10,10), breaks=30, freq=F, col="lightblue")

如果这是某种家庭作业，以展示如何计算截断分布，那么我建议您查看 rtrunc(...) 函数的代码。

【讨论】：