全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学
2510 6
2007-02-07
<P>我有一朋友想让我帮助设计SAS程序,可惜我也是初学者,所以特来请教各位师兄师姐:</P>
<P>配合下列模型进行最小二乘方差分析,比较微卫星标记基因型与表型性状之间的关系<br>Y<SUB>ijkI</SUB>= µ十P<SUB>i </SUB>+ S<SUB>j</SUB>+ H<SUB>I</SUB>+G<SUB>k</SUB>+ e<SUB>ijkI</SUB><br>其中,YijkI为表型性状的测定值; µ为群体平均值;Pi为i年龄的固定效应,i为1、2;S<SUB>j</SUB>为第j品系的固定效应,j为1、2;H<SUB>I</SUB> 为I性别的固定效应,I 为0,1;Gk为微卫星标记第k种基因型的固定效应,k为1,2,3,4,......;e<SUB>ijkI</SUB>为随机残差效应。用SAS(6.12或8.0版)的GLM (General I inearMode1)过程完成。</P>
105648.rar
大小:(2.86 KB)

 马上下载

本附件包括:

  • Book1.xls

<br>

[此贴子已经被作者于2007-4-5 11:38:50编辑过]

二维码

扫码加我 拉你入群

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

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

全部回复
2007-2-7 22:29:00
<P>我把你的公式换个写法,你看是不是一样的意思?</P>
<P>Y=p1+p2+s1+s2+h1+h2+k1+k2+...+kn+ U</P>
<P>U= &micro;+e<SUB>ijkI</SUB><BR>进一步变形:</P>
<P>X1=p1,X2=p2,X3=s1,X4=s2,X5=h1,X6=h2,X7=k1,X8=K2...</P>
<P>Y=A1X1+A2X2+A3X3+A4X4+A5X5+A6X6+A7X7+A8X8+...+U</P>
<P>你现在想求的就是A1-A8和U的参数估计值对吗?</P>
<P>另外,你的方程的G<SUB>k</SUB> 没有限制吗?可以无限多吗?你这样的方程回归不了的哦~~<SUB> </SUB></P>
二维码

扫码加我 拉你入群

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

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

2007-2-7 22:31:00
你还是给个数据集,不然帮不了你,你给个有齐所有变量的,样本数是200左右的数据集。这样才能帮你哦。
二维码

扫码加我 拉你入群

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

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

2007-2-17 07:23:00
<P>配合下列模型进行最小二乘方差分析, Y<SUB>ijkI</SUB>= &micro;十P<SUB>i </SUB>+ S<SUB>j</SUB>+ H<SUB>I</SUB>+G<SUB>k</SUB>+ e<SUB>ijkI,</SUB>SAS程序是:</P>
<P>proc glm data=yourdata;</P>
<P>class p s h g;</P>
<P>model y= p s h g/solution type3;</P>
<P>run;</P>
<P>Note: SAS always use the last level(highest level, last letter) as reference, if you need change reference group, you need  recode your variables or use dummy variables.</P>
<P>'solution' helps you get all the beta estimations. </P>
<P>Type III test tells you your effects of p, s, h, g are significant or not. </P>
<P>If you need help to read your data, please let me know.</P>
<P><SUB></SUB></P>
<P><SUB></SUB></P>

[此贴子已经被作者于2007-2-21 4:21:35编辑过]

二维码

扫码加我 拉你入群

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

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

2007-3-2 10:08:00
<P>由于最近回家过年去了,没有来的及上网,今天看到各位前辈的指点,真的是很开心.我会尽快将使用效果返回给您们的.谢谢各位.</P>
二维码

扫码加我 拉你入群

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

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

2007-4-5 11:33:00

我的数据

<P><STRONG><FONT face=Verdana color=#61b713>你好sakunamary:</FONT></STRONG></P>
<P><STRONG><FONT face=Verdana color=#61b713>上面是我的数据,我现在是个初学者,请你多帮忙。我还有其它几个引物的基因型,因为太多就没有一一上传。我的邮箱<a href="mailto:zengxiancun5022@163.com" target="_blank" >zengxiancun5022@163.com</A></FONT></STRONG></P>
<P><STRONG><FONT face=Verdana color=#61b713>                         曾献存</FONT></STRONG></P>
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群