从R中截断Logistic分布生成随机数的模拟 [英] Simulation to generate random numbers from a truncated logistic distribution in R

问题描述

我已经编写了一段代码来生成R中的Logistic随机变量，但截断在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

推荐答案

您可以使用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(...)函数的代码。

这篇关于从R中截断Logistic分布生成随机数的模拟的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！