作者andrew43 (Myrmarachne)
看板Statistics
标题Re: 常态分配的部分平均值
时间Mon Jul 6 12:05:34 2009
以下是以 R 来模拟。
我也还处於练习 R 的新手，写得笨或丑请多多指点。
取得左侧 5% 的关键值
th <- qnorm(0.05,mean=40,sd=15)
取得模拟之常态分布
sample <- rnorm(100000,mean=40,sd=15)
以回圈及判断式将小於关键值者设定为 NA 值
for (i in 1:length(sample)) {
if (sample[i] < th) sample[i] <- NA
}
取得新的平均值
mean(sample, na.rm=TRUE)
--
http://apansharing.blogspot.com/
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.128.117.116
1F：推 mili:感谢，我在excel上也是用相同方式计算，但我想微积分应该能 07/06 12:12
2F：→ mili:求出一确定值(根据pdf)。此外，用R呈现mean的移动容易吗? 07/06 12:14
3F：→ mili:我想绘制的是类似这个图: 07/06 12:14
5F：→ andrew43:看了一些R的demo，肯定画得出来，但可能不是太直接 07/06 12:24
6F：→ andrew43:你很急的话，建议用R画曲线，再用别的软体加工可能比较快 07/06 12:27
7F：推 mili:感谢 我会研究一下 07/06 22:01
8F：→ andrew43:练习了一下，可以用积分求啦哈哈哈。基本原理不可忘啊我 07/06 23:13