《数学建模第一次论文研读报告》由会员分享，可在线阅读，更多相关《数学建模第一次论文研读报告（55页珍藏版）》请在金锄头文库上搜索。
1、第一篇第一问：1、（1）CD4浓度未做变换时图形为：Regression analysis - curve estimation (回归分析曲线估计)为了使拟合曲线相对平滑，进行变换由于测得的的浓度有0值，而log的真数必须大于0，故将其加1。具体操作：transformcompute（2）HIV浓度未做变换时图形为了使拟合曲线相对平滑，进行变换轻度患者二次曲线回归CD4变化情况：时达到最大。时达到最大。综上：取重度患者当时达到最大。当时达到最大综合可知第二问：1、 多因素方差分析由于做主成分分析单位不统一，需要进行归一化，因此考虑进行多因素方差分析，具体步骤如下：(1)本例数据保存为dier
2、wen.txt，将dierwen.txt导入SPSS，General Linear Model - Univariate的顺序打开Univariate主对话框，将nongdu选入Dependent Variable框；将liaofa shijian nianling 选入Fixed Factor框;(2)单击Model按钮， 选择Custom选项，在Buil Term下拉菜单中选择Main effects，然后在Factor and Covariates中将liaofa shijian nianling liaofa*shijian 引入Model列表框， Sum of squares下拉菜单
3、中选择Type III， 勾选Include Intercept in mode框，再单击Continue按钮返回主对话框;(3)Contrasts对话框选择默认Plots对话框选择默认Save对话框选择默认Options对话框选择默认(4)单击Post Hoc按钮，在Univariate: Post Hoc Multiple Comparisions Observed Means对话框的Factor列表框中选择变量liaofa shijian nianling liaofa*shijian ，移动到Post Hoc Tests for列表框，并且在Equal Variances Assume
4、d中选择S-N-K、LSD和Scheffe选项，再单击Continue按钮返回主对话框单击OK按钮，得到分析结果.得到表格如下：2、 单因素方差分析单独对疗法进行单因素方差分析如下：(1)本例数据保存为dierwen.txt，Analyze-Compare Means-One-Way ANOVA的顺序打开One-Way ANOVA主对话框，打开dierwen.txt并将其导入SPSS，并将nongdu选入Dependent List框；将liaofa选入Factor框;(2)单击Contrasts按钮， 勾选Polynomial复选框，在Degree框中选择5次方曲线，再单击Continue按
5、钮返回主对话框;(3)单击Post Hoc按钮，在Equal Variance Assumed复选框勾选LSD和Scheffe等，再单击Continue按钮返回主对话框;(4)单击Options 按钮，在statistics复选框勾选所有选项，再选中means plots框，最后单击Continue按钮返回主对话框;(5)单击OK按钮，得到分析结果.3、 确定最佳终止治疗时间第三问1、 增量成本效用模型2、 预测第二篇第一问1、考虑首次测试的特殊性较好。从第三篇可以看出年龄对CD4和HIV也具有一定影响，我认为忽略其影响不是特别合理。2、 MATLAB程序clearA=load(3.txt);
6、CD=A(:,3); %第三列为CD4浓度HIV=A(:,5); %第五列为HIV浓度t=A(:,2); %第二列为CD4的测量时间tA=A;23991 0 0 0 0; %在A矩阵中再加入一行i=1;ID=1;C0=;H0=; %从编号为1的开始，while ID=23990 %编号从1到23990时 i0=i; %此时i0=1 c0=A(i0,3); %c0表示第三列编号为i0的病人CD4浓度的第一个数值 C0=C0;c0; %把每个病人CD4的初始浓度放在一个矩阵C0里 h0=A(i0,5); %h0表示第三列编号为i0的病人HIV浓度的第一个数值 H0=H0;h0; %把每个病人HIV
7、的初始浓度放在一个矩阵H0里 while A(i,1)=A(i+1,1) i=i+1; c0=A(i0,3); C0=C0;c0; h0=A(i0,5); H0=H0;h0; end %选出每个病人的初始浓度 ID=A(i+1,1);i=i+1;endX=C0,H0,t,C0.*H0,C0.*t,C0.2,H0.*t,H0.2,t.2;stepwise(X,CD)a=ones(1569,1) C0 H0 t t.2b,bint,r,rint,stats=regress(CD,a) X=C0,H0,t,C0.*H0,C0.*t,C0.2,H0.*t,H0.2,t.2;stepwise(X,HIV
8、)aa=ones(1569,1) H0 t t.2;b,bint,r,rint,stats=regress(HIV,aa)3、 结果：b = -125.9014 1.0582 26.7680 6.1818 -0.0973b = 0.4322 0.7884 -0.16190.00334、 预测：利用Excel将t=0时的CD4和HIV浓度筛选出来，求平均值C0=85.898，H0=5.027 y=-125.9014+89.6533+116.1745+6.1818.*t-0.0973.*t.2plot(t,y,o)title(CD4随时间的变化曲线)y=0.4322+3.9633-0.1619.*
9、t+0.0032.*t.2plot(t,y,o)title(HIV随时间的变化曲线)y=(79.9264+6.1818.*t-0.0973.*t.2)./(4.3955-0.1619.*t+0.0032.*t.2);plot(t,y,o)title(CH随时间的变化曲线)第三篇第一问1、 利用单因素方差分析求平均浓度和方差23725-255-24.8为异常数据，应剔除。利用SPSS里的analyseMean命令可以得到以下结果： 注：画蓝线处与答案差距较大2、 U检验法（1）原理：检验假设：;：;选取统计量在为真的条件下： 对给定的，选取，, ,故 ,这表明在为真的条件，是一小概率事件，由此可
10、以得出如下判定方法：计算的试验值；若，则拒绝原假设,接受；否则接受原假设。其中：检验假设：;： (事先算出样本观察值,才提这样检验问题)选取统计量在为真的条件下 对给定的，选取，,故 ,这表明在为真的条件下，是一小概率事件，由此可以得出如下判定方法：计算的试验值；若，则拒绝原假设,接受；否则接受原假设。（2） 具体计算：注意：计算时，和指的是第段时间的平均浓度和人数，而和指的是第段时间内的平均浓度和标准差。1-5 中的CD4同理可计算得：时间段（周）CD4显著性HIV显著性1512.4344显著提高-56.4170显著提高695.0669显著提高-5.988显著提高10253.1832显著提高
11、-0.9735显著提高26401.3897无显著性0.6743显著提高40以上-3.2979无显著性1.7731无显著性第二问1、 双因子方差分析Analysecompare meanmeans排序：4132432142313421总：4321年龄12CD4显著性CD4显著性31-402.5574显著提高7.7937显著提高41-50-1.8190无显著性1.0612无显著性50以上2.6073显著提高4.1988显著提高年龄34CD4显著性CD4显著性31-406.6576显著提高3.5821显著提高41-50-1.8592无显著性0.8223无显著性50以上4.0046显著提高1.8371
12、显著提高时间1234CD4显著性CD4显著性CD4显著性CD4显著性8-1.4200无0.8136无2.2960有4.2571有16-2.2089无-2.5289无-1.4527无0.3698无24-2.4363无-3.0108无-3.1955无-2.7784无32-0.3975无0.9124无0.8863无0.2426无40-1.6445无-1.8817无-0.2126无-0.4267无第四篇第一问：1、 lagrange线性插值法并求总体浓度变化量：思路：根据插值后的数据，由第四周所有病人的平均浓度减去第0周所有病人的平均浓度即为浓度平均变化量。（1）CD4clearclca=xlsread(1.xls);%修改编号，23777，23780，23990.把2370改为23778，23990改为23779;d=a(any(isnan(a(:,2),2),:);x,y=size(d);b=;B=cell(356,1);%创建单元数组，使23424编号5行5列为B的第一个元素t=1;k=0;for
