排除自身的平均命令怎么写

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

最近看一篇文章，要计算排除了本家庭的信任度之外的其他样本家庭的平均信任度，求教这个命令是怎么写的？
数据简介： COUN省县代码  VILL  村庄代码 HOUS    住户代码   trust信任度
一个家庭所在村庄层面的平均信任度：by  COUN VILL: egen trustmean=mean(trust)，可是如何计算排除此家庭之外其他样本家庭的平均信任度呢？先谢谢了。

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：trust tmean mean egen 信任度 命令 排除

1# 天涯在何方
*看是否合要求
bys COUN VILL :egen t=total(trust)
bys COUN VILL :gen tmean=(t-trust)/(_N-1)

jzhyue 发表于 2010-9-14 11:22
1# 天涯在何方
*看是否合要求
bys COUN VILL :egen t=total(trust)
bys COUN VILL :gen tmean=(t-trust)/(_N-1)

知道这个思想，就是不知道命令怎么写，呵呵，多谢兄台指点！

levelsof HOUS, local(household)

foreach x of local household {
        gen trust2= trust
        replace trust2 = . if HOUS == `x'
        by COUN VILL: egen trustmean`x' = mean(trust2)
        drop trust2
}

ctx5518 发表于 2010-9-14 11:41
levelsof HOUS, local(household)

foreach x of local household {
        gen trust2= trust
        replace trust2 = . if HOUS == `x'
        by COUN VILL: egen trustmean`x' = mean(trust2)
        drop trust2
}

佩服佩服。也感谢jzhyue的解答，我本来和您是一样的想法，但这位兄台的命令也很有道理。特奉上统计描述，以供大家比较。
   Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
  trustmean1 |      9195    3.513549    .8065215          1          5
       tmean |      9194    3.518055    .8044878          1          5
       trust |      9194     3.51849    1.323353          1          5

天涯在何方 发表于 2010-9-14 12:06

ctx5518 发表于 2010-9-14 11:41
levelsof HOUS, local(household)

foreach x of local household {
        gen trust2= trust
        replace trust2 = . if HOUS == `x'
        by COUN VILL: egen trustmean`x' = mean(trust2)
        drop trust2
}

佩服佩服。也感谢jzhyue的解答，我本来和您是一样的想法，但这位兄台的命令也很有道理。特奉上统计描述，以供大家比较。
   Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
  trustmean1 |      9195    3.513549    .8065215          1          5
       tmean |      9194    3.518055    .8044878          1          5
       trust |      9194     3.51849    1.323353          1          5

*按照ctx5518的思路
levelsof HOUS, local(household)
gen ttmean=.
foreach x of local household {
        gen trust2= trust
        replace trust2 = . if HOUS == `x'
        by COUN VILL: egen trustmean`x' = mean(trust2)
        drop trust2
  replace ttmean=trustmean`x' if HOUS == `x'
}
drop trustmean*

jzhyue 发表于 2010-9-14 15:48

天涯在何方 发表于 2010-9-14 12:06

ctx5518 发表于 2010-9-14 11:41
levelsof HOUS, local(household)

foreach x of local household {
        gen trust2= trust
        replace trust2 = . if HOUS == `x'
        by COUN VILL: egen trustmean`x' = mean(trust2)
        drop trust2
}

佩服佩服。也感谢jzhyue的解答，我本来和您是一样的想法，但这位兄台的命令也很有道理。特奉上统计描述，以供大家比较，从样本数目来看，9194是和原样本一致的。
   Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
  trustmean1 |      9195    3.513549    .8065215          1          5
       tmean |      9194    3.518055    .8044878          1          5
       trust |      9194     3.51849    1.323353          1          5

*按照ctx5518的思路
levelsof HOUS, local(household)
gen ttmean=.
foreach x of local household {
        gen trust2= trust
        replace trust2 = . if HOUS == `x'
        by COUN VILL: egen trustmean`x' = mean(trust2)
        drop trust2
  replace ttmean=trustmean`x' if HOUS == `x'
}
drop trustmean*

谢谢，这回我把这个命令输入，得到第一种思路和第二种思路的描述统计，他们虽然不一致，但是很接近了。
    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
      ttmean |      9195    3.518543    .8046014          1          5
       tmean |      9194    3.518055    .8044878          1          5

*问题的关键是trust是否有缺失值，而你的数据就有一个：下面排除缺失值（同6楼）
bys COUN VILL :egen t=total(trust)
bys COUN VILL :egen tn=count(trust)
gen tmean=(t-trust)/(tn-1)
sum if trust~=.
*如果将缺失值作为观察值之一，下面比较合理（因为没有删除缺失值，而上面的方法与6楼的方法都不考虑缺失值）
bys COUN VILL :egen t=total(trust)
bys COUN VILL :gen tmean=(t-trust)/(_N-1)

天涯在何方 发表于 2010-9-14 16:14

jzhyue 发表于 2010-9-14 15:48

天涯在何方 发表于 2010-9-14 12:06

ctx5518 发表于 2010-9-14 11:41

谢谢，这回我把这个命令输入，得到第一种思路和第二种思路的描述统计，他们虽然不一致，但是很接近了。
    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
      ttmean |      9195    3.518543    .8046014          1          5
       tmean |      9194    3.518055    .8044878          1          5

我想应这样统计：
sum if trust~=.

Variable	Obs	Mean	Std. Dev.	Min	Max
trust	9194	3.51849	1.323353	1	5
ttmean	9194	3.51849	.8046295	1	5
tmean	9194	3.51849	.8046295	1	5

jzhyue 发表于 2010-9-14 16:50
*问题的关键是trust是否有缺失值，而你的数据就有一个：下面排除缺失值（同6楼）
bys COUN VILL :egen t=total(trust)
bys COUN VILL :egen tn=count(trust)
gen tmean=(t-trust)/(tn-1)
sum if trust~=.
*如果将缺失值作为观察值之一，下面比较合理（因为没有删除缺失值，而上面的方法与6楼的方法都不考虑缺失值）
bys COUN VILL :egen t=total(trust)
bys COUN VILL :gen tmean=(t-trust)/(_N-1)

这回这个问题看的透彻多了，多谢兄台指点！