数学模型9.1代码记录

clear; clc;  % 不吸烟与吸烟孕妇的新生儿体重和怀孕期的参数估计和假设检验
M1=dlmread('data0901.txt');  % 读取有分隔符的ASCII数据文件
M=M1(:,2:end);%将纯数据读出来
k0=0;k1=0; %k0记录不吸烟，k1记录吸烟
for i=1:1236if M(i,7)==0       % 吸烟状况,0~不吸烟k0=k0+1;y0(k0)=M(i,1); % 不吸烟孕妇的新生儿体重x0(k0)=M(i,2); % 不吸烟孕妇的怀孕期else  if M(i,7)==1 % 吸烟状况,1~吸烟k1=k1+1;y1(k1)=M(i,1); % 吸烟孕妇的新生儿体重x1(k1)=M(i,2); % 吸烟孕妇的怀孕期endend
end
% k0+k1=1226,有10名孕妇吸烟状况=9，缺失
% 新生儿体重无缺失
[y0m,s0,y0c]=normfit(y0); % 估计不吸烟孕妇新生儿体重y0的均值y0m，标准差s0和均值的置信区间y0c% 显著水平α（输入参数2）,缺省时默认为0.05,即置信度为95%
[y1m,s1,y1c]=normfit(y1); % 估计吸烟孕妇新生儿体重元有的均值y1m，标准差s1和均值的置信区间y1c
[yh,pyh]=ttest2(y0,y1,'tail','right'); % 检验 H0:y0m<=y1m, H1:y0m>y1m% T检验,主要用于样本含量较小（例如n<30），总体标准差σ未知的正态分布资料% ttest 两个样本数量一致，ttest2 不一致% tail 备择假设右侧检验% yh=0 接受原假设H0，yh=1，拒绝H0% pyh p值（H0真时，得到观察值的概率）， p>Alpha 时，接受原假设 H0 ； p<=Alpha 时，拒绝H0
% 估计、检验不吸烟与吸烟孕妇的新生儿体重低（小于2500克=88.2盎司）的比例
n01=0;
for i=1:k0if y0(i)<88.2  % n01 不吸烟孕妇的新生儿体重低n01=n01+1;end
end
q0=n01/k0; % 估计不吸烟孕妇的新生儿体重低的比例
n11=0;
for i=1:k1if y1(i)<88.2 % 吸烟孕妇的新生儿体重低n11=n11+1;end
end
q1=n11/k1; % 估计吸烟孕妇的新生儿体重低的比例
% 孕妇吸烟可能导致新生儿体重低
qs=[q0*(1-q0)*(k0-1)+q1*(1-q1)*(k1-1)]/(k1+k0-2)*(1/k1+1/k0);
qt=(q1-q0)/sqrt(qs); % 两个0-1分布总体均值的检验 H0:q0=q1
%% 怀孕期均值的假设检验（H0:不吸烟早孕，H1：吸烟早孕）% 怀孕期均值皆有缺失值
j0=0;
for i=1:k0if x0(i)<999; % 不吸烟孕妇的怀孕期数据不缺失j0=j0+1;xx0(j0)=x0(i);end
end
j1=0;
for i=1:k1if x1(i)<999; % 吸烟孕妇的怀孕期数据不缺失j1=j1+1;xx1(j1)=x1(i);end
end
[x0m,sx0,x0c]=normfit(xx0), % 估计不吸烟孕妇怀孕期的均值x0m,标准差sx0和均值的置信区间x0c
[x1m,sx1,x1c]=normfit(xx1), % 估计吸烟孕妇怀孕期的均值x1m,标准差sx1和均值的置信区间x1c
[xh,pxh]=ttest2(xx0,xx1), % 检验 H0:x0m=x1m
[xhh,pxhh]=ttest2(xx0,xx1,0.05,'right'),% 检验 H0:x0m>x1m
%% 早产率的假设检验
% 估计、检验不吸烟与吸烟孕妇早产（怀孕期小于37周）的比例
m01=0;
for i=1:j0if xx0(i)<259 % 不吸烟孕妇怀孕期小于37周m01=m01+1;end
end
r0=m01/j0, % 估计不吸烟孕妇怀孕期小于37周的比例
m11=0;
for i=1:j1if xx1(i)<259 % 吸烟孕妇怀孕期小于37周m11=m11+1;end
end
r1=m11/j1,  % 估计吸烟孕妇怀孕期小于37周的比例
rs=[r0*(1-r0)*(j0-1)+r1*(1-r1)*(j1-1)]/(j1+j0-2)*(1/j1+1/j0);  % 0-1分布 % 两样本联合估计的方差 % 两独立样本均值的T检验
rt=(r1-r0)/sqrt(rs), % 两个0-1分布总体均值的检验 H0:r0=r1

数学模型9.1代码记录相关推荐

uiautomator 代码记录：BT接收测试
uiautomator 代码记录 :BT接收测试 package rom;import java.lang.*; import java.util.Random; import java.io.Fil ...
自定义起始时间的时间戳计算（个人代码记录）
个人代码记录,这个时间戳计算,由于项目数据特殊需求,自己设置的时间戳起始时间为2019年1月1日开始, 然后根据时间戳数据计算出现在时间,具体使用如下,这里setFullYear的月份是0到11,所以 ...
今日代码(20210313)--美赛代码记录
代码记录第1及第6题(PageRank+Lasso) my_pagerank <- function(M, r, n, b) {N <- dim(M)[2]r <- r/sum(a ...
insightface人脸识别代码记录(一)（数据前期准备）
一.前言这部分主要讲训练数据的制作.一是我们直接采用作者提供好的数据集,二就是制作我们自己所需要的数据集. 目录地址:insightface人脸识别代码记录(总)(基于MXNet) 二.主要内容 1 ...
git 代码记录单条合并的方法
第一种方法 **步骤一:** 1,在开发分支 git log 找提交记录找到要合并的记录一串编码**步骤二**: 2,切换到要合并的分支例如:git checkout master **步骤三:* ...
TM1650代码记录（基于51单片机）
文章目录一.TM1650 0.介绍 1.典型应用电路 2.原理图 3.PCB 4.实物图 5.键盘扫描码二.代码记录 1.C文件: 2.H文件三.结束语一.TM1650 0.介绍功能特点: ...
DSP学习（4）—— 看代码记录
DSP学习(4)-- 看代码记录文章目录 DSP学习(4)-- 看代码记录前言一.一些缩写二.C语言忘记的知识条件编译中使用的预编译指令前言本文用于记录看代码中间遇到的不认识的东西. 一 ...
详解鲸鱼优化算法原理、数学模型和实例代码
鲸鱼优化算法 (whale optimization algorithm,WOA)是 2016 年由澳大利亚格里菲斯大学的Mirjalili 等提出的一种新的群体智能优化算法,其优点在于操作简单,调整 ...
【我的笔记】VScode遇到的小问题以及个人配置代码记录
笔记记录bug,arduino选不了板个人配置代码记录环境搭建及编译时中文乱码问题下载页面的使用说明库的配置记录bug,arduino选不了板开发arduino我用VScode,但是最近 ...

数学模型9.1代码记录

数学模型9.1代码记录相关推荐

最新文章

热门文章