今天刚刚碰到一个SAS群里面的人问起求函数最优解的问题以前刚刚好写过这方面的代码,现在重新写发现生疏了很多干脆写一篇博客记下来,省得到时候又到处找资料。。。。在SAS里面optmodel过程都能解决关于求最优的所有问题,当然还包括一些其他的过程例如proc nlp,proc lp等等但是这些过程步骤全部包含茬proc optmodel下面,proc optmodel求最优过程数据来源有两种一种是自己创建,一种是去读取sas的表自己创造数据很简单,proc optmodel过程如下:
另外一种是去读取表变量和数据都在sas表里面,proc optmodel过程具体例子如下:然后在这个数据集里 id是主键 factor1跟2都是每个id不同的两个因子值
然后给每个id分配权重w(id) (一共多少个id我不知道)
内容提示:用t检验比较多个均数犯Ⅰ类错误概率的SAS电脑实验
文档格式:PDF| 浏览次数:50| 上传日期: 22:53:06| 文档星级:?????
全文阅读已结束如果下载本文需要使用
SAS 医学统计学实习上机笔记含t检驗,正态性检验方差分析,秩和检验卡方检验,直线、回归相关分析等格式详见附件 (a.class表示建在永久库a中,不加a.表示建在临时库中) /*缺失的数据用.表示但SAS系统在储存及处理数据时文本型的缺失值用空格表示数值型的缺失值用“.”表示,*/ 直接从外部数据文件中读取 直接导入导出文件-导入- 数据集的修改,见sas基础教程实用版 一、定量资料的统计推断—t检验 备注:正态性检验,样本量小(n≤50)看S-W值,樣本大时看K-S值 仅用于平衡设计资料(即按所分析的变量水平将研究对象分组每个组内的观察对象数目应相等) 功能较ANOVA过程强大,除可进行方差分析外还可进行回归等分析 用于平衡或非平衡设计资料均可,当不能肯定资料类型是否符合ANOVA的适用条件时最好改用GLM过程进行方差分析 1. 完全随机设计的方差分析 /*输出为文本文件*/ 1、独立样本设计2 ×2列联表资料的卡方检验 /*P142 例11.1血塞通和银杏达莫治疗急性脑梗死的疗效比较*/ /*例 11.2 两種方法治疗肱骨外上髁炎的效果比较*/ 2、独立样本设计R×C表资料的卡方检验 3、配对2 ×2列联表的卡方检验 /*例11.6 两种方法诊断低血钾的结果(阳性率)比较*/ 四、直线、回归相关分析 /*i=rl表示出现直线回归线,width=1表示线的宽度height表示散点的大小*/ model y=x/clm cli clb; /*以y为应变量,x为自变量建立线性回归模型并输絀应变量条件均数的95%可信区间,输出个体y值的95%预测区间回归系数的95%可信区间*/ label x="过去2周内抗生素使用情况" y="多重耐药肠杆菌分离结果"; ? 成组设計资料的秩和检验 ? 随机区组设计的秩和检验 1、配对设计的两样本比较 2、单一样本与总体中位数比较 title "单一样本与总体中位数比较"; 3、成组设計资料的秩和检验 3.1成组设计定量资料两样本比较 分别使用含25%蔗糖的奶嘴和普通奶嘴,测得新生儿足底采血时疼痛评分问两组疼痛评分有無差别? title "成组设计两样本比较秩和检验"; if为分支(条件)语句结合语句"g=1",这条分支语句的意思是:g的初始值等于1如果"_n_>10"不成立,则不执行then後面的程序g还是等于1; 如果"_n_>10"成立,则执行then后面的程序g就等于2。程序执行的结果是:前10条数据为第1组第10条以后的数据为第2组就是分组 */ 3.2等级资料两样本比较 4、成组设计资料的秩和检验 4.1成组设计定量资料多样本比较 /*以下4条语句对变量x作秩变换:求x的秩次,并存入新变量rx中*/ 4.1成組设计等级资料多样本比较 5、随机区组设计的秩和检验 title "随机区组设计的秩和检验"; |
版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。