泊松抽样问题代码咨询

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

N = nrow(data)
ID = 1:N
l = data.frame(ID,h,data)#h为另一种方法算出的抽样概率
Tot=l$h
n =nrow(data)*0.1
pik=inclusionprobabilities(Tot,n)
m=UPpoisson(pik)  
请问，m运行出的结果中“1”的个数（1表示被抽中）为什么不等于n

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Inclusion Abilities poisson Frame Data

是不是少了几个就是这么神奇用泊松抽样就是会这样，别管他为什么了

是多了3行，不会是因为代码不合适的问题吗

llb_321 发表于 2021-1-27 11:29
是不是少了几个就是这么神奇用泊松抽样就是会这样，别管他为什么了

多了3行，与代码无关吗，我一直以为是代码哪里不合适

不是代码的问题，是UPpoisson这个函数抽样时是按照pik进行的，而pik是inclusionprobabilities给出的每个观察值纳入取样的概率，虽然sum(pik)=n，但最后table(m)确实会出现1的结果不等于n，而且同样的pik，多进行几次UPpoisson，你会发现每次都不一样。
正常。貌似只有UPpoisson会这样，其他的抽样就不会。

llb_321 发表于 2021-1-27 19:35
不是代码的问题，是UPpoisson这个函数抽样时是按照pik进行的，而pik是inclusionprobabilities给出的每个观察 ...

谢谢，学到了

llb_321 发表于 2021-1-27 19:35
不是代码的问题，是UPpoisson这个函数抽样时是按照pik进行的，而pik是inclusionprobabilities给出的每个观察 ...

您好，我再请教一个问题可以不，想要对矩阵的行进行均匀抽样，在R中可以用哪个函数包呢，找了好久了，

sampling::srswor() 可以试试

llb_321 发表于 2021-1-28 10:55
sampling::srswor() 可以试试

[em23]谢谢！

llb_321 发表于 2021-1-28 10:55
sampling::srswor() 可以试试

跟您再请教一下，在R中怎样做可以提取矩阵中n行与对应的n列相交位置的数据？