SAS t检验_sas 单样本t检验_Tiaaaaa的博客-程序员秘密

/*----------------- 一.单样本t检验 -----------------*/
/*有原始数据的t检验*/
DATA TEMP;
	INPUT X@@;
	Y=X-20.7;
	CARDS;
	20.99 20.41 20.10 20.00 20.91 22.60
	20.99 20.42 20.90 22.99 23.12 20.89
	;
TITLE 'T检验 - PROC MEANS过程';
PROC MEANS 
     mean            /*求平均数*/
     stderr          /*求平均数的标准误*/
     t              /*t检验，检验样本的平均数等于0的假设是否成立，H0：均值为0*/
     probt;         /*或prt，t检验的显著性*/
	 var y;        /*参与分析的变量名称*/
RUN;
TITLE 'T检验 - PROC Univariate过程';
PROC UNIVARIATE DATA=TEMP
                 normal;  /*要求univariate程序检定输入资料是否呈现常态分配，并且输出检验结果*/
				 var y;
RUN;



/*无原始数据的t检验*/
DATA TEMP1;
	INPUT u0 u1 s n;        /*依次输入总体均值、样本均值、样本标准差、样本数*/
	t=(u1-u0)/(s/sqrt(n));  /*样本t检验的计算公式*/
	v=n-1;                  
	p=probt(t,v);           /*调用标准正态分布函数查看t值对应的p值*/
	cards; 
72 74.2 6.5 25
;
proc print;   /*输出t检验的结果*/
run;

/*----------------- 二.配对样本t检验 -----------------*/
/*有原始数据的配对t检验*/
/*ttest过程*/
data temp;
     input id a b;
	 cards;
1 2.41 2.80
2 2.90 3.04
3 2.75 1.88
4 2.23 3.43 
5 3.67 3.81
6 4.49 4.00
7 5.16 4.44
8 5.45 5.41
9 2.06 1.24
10 1.64 1.83
11 1.06 1.45
12 0.77 0.92
;
proc ttest;
     paired a*b;
run;

/*means过程*/
data temp;
	set temp;
	y=a-b;
	keep id a b y;
run;
proc means mean stderr t probt;
	var y;
run;

/*univariate过程*/
proc univariate normal;
	var y;
run;

/*无原始数据的配对t检验*/
data temp1;
	input d s n;  /*依次输入配对样本的差值、标准差、配对数*/
	t=d/(s/sqrt(n));
	v=n-1;
	p=probt(t,v);
	cards;
0.625 0.78 8
;
proc print;
run;


/*----------------- 三.独立样本t检验 -----------------*/
/*有原始数据的两独立t检验*/
data temp;
      input group x@@;
	  cards;
1 134 1 146 1 104 1 119 1 124 1 161 1 107 1 83 1 113 1 129 1 97 1 123
2 70 2 118 2 101 2 85 2 107 2 132 2 94

proc ttest;
	class group;
	var x;
run;

/*无原始数据的两独立t检验*/
data temp;
     input x1 x2 n1 n2 s1 s2;
	 sc=sqrt((1/n1+1/n2)*((n1-1)*s1**2+(n2-1)*s2**2)/(n1+n2-2));   /*计算两独立样本共同的标准差*/
	 t=(x2-x1)/sc;
	 v=n1+n2-2;
	 p=probt(t,v);
	 cards;
20.95 21.79 20 20 5.89 3.43
;
proc print;
run;

本文链接：https://blog.csdn.net/Tiaaaaa/article/details/58612860

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

细谈Verilog中$signed()和$unsigned()用法_不吃人的坤坤坤坤坤的博客-程序员秘密

之前在做Verilog，有符号数和无符号数运算时，对这俩个一直没弄的很透彻。本文就以比较简短的篇幅去捋一下其用法。

opencv的机器学习的bug_使用opencv中的机器学习算法,进行预测predict时内存错误_hnsdgxylh的博客-程序员秘密

在VS跟opencv做图像处理时，win32平台使用Opencv的ML库加载数据时，放在当前目录下不会有任何错误但是假如使用MFC做界面时，采用绝对路径，比如C:\data.xml读取数据不会出错，但是放在当前路径下会有读取不到数据的情况，在使用SVM时就发生了采用的方法是采用这样的当前路径“../data.xml”可以避免错误

TTGO可编程手表入门_库 esp insights 中的分类 '' 不可用。设置为 'uncategorized_流鱼a_ec的博客-程序员秘密

一、简介前几天淘了一款可编程手表，于是准备玩一下这个小玩意。二、环境搭建在开始之前，首先就是搭建相关的环境。按照说明书上的来，不可避免的，我又遇到了麻烦。于是开始在百度上寻找答案。终于找到一篇博文：(16条消息) ESP32环境搭建(arduino)_风度青年的博客-程序员秘密_arduino esp32...

Spark面试相关_aeluwl2038的博客-程序员秘密

Spark Core面试篇01一、简答题1.Spark master使用zookeeper进行HA的，有哪些元数据保存在Zookeeper？答：spark通过这个参数spark.deploy.zookeeper.dir指定master元数据在zookeeper中保存的位置，包括Worker，Driver和Application以及Executors。standby节点要从zk中获...

gaussdb 数据库用户和安全管理【数据加密】【03】_高斯数据库关闭ssl加密_Zhao.Mr的博客-程序员秘密

1. 数据加密数据加密的密钥由系统对钱夹密码的复杂变化逻辑生成。钱夹密码由用户指定输入，必须符合密码复杂度要求。注意：GaussDB不提供修改钱夹密码的接口，一旦设定，请牢记。01.设定用于数据加密的钱夹密码。通过以下SQL命令设定钱夹密码，以设定钱夹密码为gaussdb_123为例：ALTER SYSTEM SET ENCRYPTION KEY IDENTIFIED BY 'gaussdb_123'; 说明：钱夹密码设定后，GaussDB通过对钱夹密码的一系列复杂变化，生成用于数

ATM模拟（查余额，存款，取款）_一键生成atm机余额视频_是小光a~的博客-程序员秘密

ATM模拟该代码只能执行简单的查余额，存款，取款操作（卡号信息由用户输入，默认密码为"000000"，默认余额为0）该代码缺点:功能较少,注释信息，提示信息较多,虽然部分体现了类的思想,但是主函数中仍然是大的模拟,可以优化为操作类或外部函数来使代码具有可读性(这个代码读起来还是挺烦的)。。。#include<cstdio>#include<iostream>#in...