求助：perl程序

joson12345

求助一perl程序，不知哪位朋友能帮帮忙。多谢了！（实验需要，有点着急）
如有不清楚的请及时发帖。条件如下：

n1＝100 n2＝200 n3＝100 m1＝25 m2＝50 m3＝50 m4＝100 m5＝25
m6＝50
Step1:smok个数n1, nonsmok个数n2, unknow个数n3。体重。
产生n1个smok个体的表型值（浓度）使其分布满足 Y(n1)— N(110,50),
产生n2个nonsmok个体的表型值（浓度）使其分布满足 Y(n2)－ N(100,50),
产生n3个unknow个体的表型值（浓度）使其分布满足 Y(n3)— N(90,50),
产生所有的个体(n1＋n2＋n3)表型分布满足 Y(n1＋n2＋n3)—  N(102, 700),
格式如下（例如）：
分类1             浓度
smok            102
smok            100
smok            125
.....                ....   
n1个          ...

nonsmok      150
nonsmok      180
nonsmok      110
.....
n2个          ...
....

unknow        100
unknow        150
unknow         98
.....                 ...
n3个          ....


Step2：满足step1条件后，
从smok组里面随机拿出m1个放到A组里面, m2个放到B组里面, 剩下的[n1－(m1+m2)]个放到C组里面。
从nonsmok组里面随机拿出m3个放到A组里面, m4个放到B组里面, 剩下的[n2－(m3+m4)]个放到C组里面。
从unknow组里面随机拿出m5个放到A组里面, m6个放到B组里面, 剩下的[n3－(m5+m6)]个放到C组里面。

Step3：按照A、B、C重新分组，此时样本数
A组有（m1＋m3＋m5）
B组有（m2＋m4＋m6）
C组有（n1＋n2＋n3－m1－m2－m3－m4－m5－m6）

格式如下（例如）：
分类2                          浓度
A                               150
A                               100
A                               125
.....                             ....
(m1＋m3＋m5)个     ....


B                                102
B                                110
B                                 98
.....                               ....
（m2＋m4＋m6）个          ....

C                                180
C                                100
C                                150
....
(n1+n2+n3-m1-m2-m3-m4-m5-m6)个        ....
....

按照A、B、C三组个体的浓度进行方差分析，求出F值（F=MSA/MSE）。

Step4：输出F值放入excel表格里面， 重复Step2、Step3步骤10次。

兰花仙子

致lz:
这里不是script writer，这里只是help你解决你的script problem，首先你应提出自己的思路，并给出自己的code，然后再提出疑问，然后我们才能help你resolve你的problem。

joson12345

兰花仙子朋友，我刚刚学perl也是出于实验的需要才学的。以前专门做基因操作上的研究，现在没有办法才学习的perl的。暂时还编写不了像样的程序。还情大家能帮帮忙。着急呀！

笨狗

原帖由 "兰花仙子" 发表：
致lz:
这里不是script writer，这里只是help你解决你的script problem，首先你应提出自己的思路，并给出自己的code，然后再提出疑问，然后我们才能help你resolve你的problem。

支持

duanwenhui

没有看懂啊！详细一点啊

joson12345

duanwenhui朋友，给个联系方式，我发email给你。

duanwenhui

duanwenhui@gmail.com
偶明天正好有空！

joson12345

N(110,50), 这是要求数据服从正态分布，均值为110， 方差为50。
是不是这里不懂呀。我的要求就是产生很多个体的浓度值，这些个体要求放在指定的分组内。然后重新分组，把个体放在指定的组内（重新分组）。

joson12345

我马上发给你

joson12345

已经发给你了！注意查收，多谢朋友，如果后面的方差分析加不进去的话，只写前面的程序也可以。