正在生成范围内的rnorm（）数据

如何在10-15的范围内生成随机分布的数据？ 我试过了，但大部分数据都是10秒

library(MCMCglmm)

x=matrix(rtnorm(n = 100, lower = 10, upper = 15), nrow=10, ncol=10)

这是矩阵的一部分：

        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]     [,9]    [,10]
 [1,] 10.01323 10.09077 10.12173 10.21363 10.14710 10.07159 10.19811 10.16515 10.14877 10.13405
 [2,] 10.10269 10.08853 10.07679 10.21932 10.15413 10.00907 10.00788 10.06476 10.06058 10.09998
 [3,] 10.00339 10.02658 10.05871 10.03135 10.07163 10.03250 10.08346 10.17975 10.10160 10.12360
 [4,] 10.00334 10.07949 10.00456 10.04519 10.11282 10.06135 10.05492 10.10164 10.13717 10.00696
 [5,] 10.12996 10.07651 10.25793 10.10667 10.14016 10.02615 10.02722 10.16351 10.28760 10.06861

---

如图所示，大多数数据是10秒到15秒

＞P>你可能想要改变平均值在范围的中间；默认为零

library(MCMCglmm)
set.seed(101)
par(mfrow=c(1,2))
hist(rtnorm(1e5,lower=10,upper=15),col="gray",main="mean=0")
hist(rtnorm(1e5,mean=12.5,lower=10,upper=15),col="gray",main="mean=12.5")

也许是为了回答你修改过的问题

rtnorm(n, mean = x%*%beta+8, sd = 2, lower=10, upper=15)

---

会起作用。

< P>你可能想把平均值改变到范围的中间；默认为零

library(MCMCglmm)
set.seed(101)
par(mfrow=c(1,2))
hist(rtnorm(1e5,lower=10,upper=15),col="gray",main="mean=0")
hist(rtnorm(1e5,mean=12.5,lower=10,upper=15),col="gray",main="mean=12.5")

也许是为了回答你修改过的问题

rtnorm(n, mean = x%*%beta+8, sd = 2, lower=10, upper=15)

---

将起作用。

您使用哪种rtnorm？有许多具有截断正态分布的包。@我正在使用此包：libraryMCMCglmm@RuiBarradas我会更好地编辑我的问题。你能解释一下为什么你完全改变了你的问题吗？我回答了你先前的问题；如果那不是你真正的问题，你能问一个新问题吗？我很想把你的编辑内容放回去……我可以看出它们之间的联系，但在已经发布了答案之后，对你的问题做如此多的修改，感觉有点粗鲁。我已经给出了我认为是对你新问题的回答，但是在将来发布一个新问题会更好，即使它是密切相关的-你可以链接到上一个问题…你使用的是哪种rtnorm？有许多具有截断正态分布的包。@我正在使用此包：libraryMCMCglmm@RuiBarradas我会更好地编辑我的问题。你能解释一下为什么你完全改变了你的问题吗？我回答了你先前的问题；如果那不是你真正的问题，你能问一个新问题吗？我很想把你的编辑内容放回去……我可以看出它们之间的联系，但在已经发布了答案之后，对你的问题做如此多的修改，感觉有点粗鲁。我已经给出了我认为是对你新问题的回答，但是在将来发布一个新问题会更好，即使它是密切相关的-你可以链接到上一个问题…我编辑了我的问题，因为它不够精确，我仍然有一个噪声项，它改变了分布。我的问题是因为它不够精确，我仍然有一个噪声项，它改变了分布