帮助设计程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

<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 为0，1；Gk为微卫星标记第k种基因型的固定效应，k为1,2,3,4,......；e<SUB>ijkI</SUB>为随机残差效应。用SAS(6.12或８.０版)的GLM (General I inearMode1)过程完成。</P> 105648.rar (2.86 KB) 本附件包括：

• Book1.xls

2007-4-5 11:39:00 上传

帮助设计程序

<br>

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

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：设计程序 General Micro sas程序 固定效应 设计程序 General 初学者 基因型 平均值

<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>

你还是给个数据集，不然帮不了你，你给个有齐所有变量的，样本数是200左右的数据集。这样才能帮你哦。

<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编辑过]

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

<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>

105645.rar (2.86 KB) 本附件包括：

• Book1.xls

2007-4-5 11:35:00 上传

帮助设计程序

<BR>