搅屎棍茅于轼在5月26日发表了《50-60年饿死人数的估算方法》，随后这篇博文就给和谐掉了。文中这老头不知道在哪里找到的国家估计局的统计数据，弄了一个简单的二次回归就得出大饥荒时期我国饿死的人数大约是3600万。
疑惑一：估计局的数据和模拟的数据出奇的一致，这不得不让人怀疑数据真是估计出来的。
疑惑二：3600万的计算方法。
【根据回归方程
第N年的人口数（千万） = 54.074 + 1.0684N + 0.02933N2
1950年的N为1； 1951年的N为2，如此等等。
推断1960年底（1961年初）的人口数应为71118万人。这是正常人口变化在1960年底应有的人口数。
而1960年底的人口数为65859万人，这是非正常人口变化的结果。二者相差 5259 万人。
差数是由于灾荒期间（1959和1960两年内）非正常人口出生和非正常人口死亡造成的。
将此差数5259万减去1959和1960 两年内因灾荒而减少的出生人口数（非正常出生），就能得出这两年中因灾 荒而死亡的人口数（非正常死亡）。也就是饿死的人数（万人）。
5259 – 65994（3.5% - 2.478%）- 67207（3.5% - 2.086%）
= 5259 – 675 – 950 = 3634
得出1959和1960两年内饿死的人为3634万。
计算中65994和67207是1959,1960年底的实际人口数（万人）。此处假定正常的出生率是3.5%，
是1950-1958九年出生率的平均数。2.478%和2.086%是1959,1969年公布的实际出生率。】
用50-58的平均出生率3.5%减去59和60的实际出生率做为非正常出生率给人的感觉有点像为了数据而数据的味道。
<br />
mao<-read.table('c:/temp/mao.txt')<br />
names(mao)<-c('year','pop','birthrate','mortality','fitpop')<br />
attach(mao)<br />
mao58<-mao[1:9,1:2]<br />
mao58<-transform(mao58,nb=1:9)<br />
fm<-lm(pop~nb+I(nb^2),data=mao58)<br />
plot(year,pop,pch=16,col=4,xlab='年份',ylab='人口')<br />
points(year,fitpop,col='red',pch=2)<br />
lines(year,fitpop,col='red')<br />
</p>
[attachment=219663,1151]
