首页 > 学术百科

lsd两两比较(均数标准差)

data a;

input gp$ n X S;

cards;

A 5 2.46 0.37

B 5 2.41 0.18

C 5 3.02 0.04

D 5 4.03 0.90

;

/*采用sql过程进行方差分析计算*/

proc sql;

吉尔吉斯族create table b as

select sum(n) as N, /*计算总例数N*/

count(n) as k, /*计算组数k*/

sum(n*X**2)-sum(n*X)**2/sum(n) as SSZJ, /*计算组间变异*/

sum(S**2*(n-1)) as SSZN /*计算组内变异*/

from a;

quit;

data c;

佛山金马影剧院set b;

MSZJ=SSZJ/(k-1); /*计算组间均方*/

MSZN=SSZN/(N-k); /*计算组内均方*/

F=MSZJ/MSZN; /*计算F值*/

pubmed数据库P=1-probf(F,k-1,N-k); /*计算P值*/

proc print data=c;

var F P; /*显示F值，P值*/

run;

/*以下进行两两比较计算，采用LSD-t检验*/

data d;

set c;猎德小学

keep N k MSZN;

proc sql;

create table e as

select a1.gp as gp1,a2.gp as gp2, /* 交叉匹配生成两两比较数据集*/

a1.n as n1,a2.n AS n2,

a1.X AS X1,a2.X as X2

from a as a1, a as a2

where a1.gp<a2.gp;

create table f as

select e.*,d.* from e ,d; /*生成两两比较计算表 */

quit;

data g;

set f;

羟基磷酸钙dmean=x1-x2;

地球物理学报

sd=sqrt(MSZN*(1/n1+1/n2));

t=dmean/sd;

p=(1-probt(abs(t),n-k))*2;

proc print data=g;

var gp1 gp2 n1 n2 x1 x2 dmean t p; /*显示两两比较结果*/

run;

本文发布于:2024-09-24 14:23:08，感谢您对本站的认可！

本文链接：https://www.17tex.com/xueshu/93939.html

上一篇：【临床研究】左眶额皮质低频重复经颅磁刺激强迫症的初步研究

下一篇：医学统计学常用方法小结5篇

标签：计算生成采用进行

留言与评论（共有 0 条评论）