，xp，用向量表示X＝（x1，x2，…，βp 为待估参数．Logistic回归模型中的参数估计常用极大似然估计法得到．设y 是0－1变量，x1，x2，…，n），取P＝πi，P＝1－πi，则yi 的联合概率函数为，yi＝0，1；i＝1，2，…，yn 的似然函数为对数似然函数为对于Logistic回归，将代入，得令＝0，可以用数值计算求待估参数β0，β1，β2，…，βp 的极大似然估计．用R 软件可以解决Logistic回归模型中的参数估计、检验等问题．......
2023-11-18
应用R实验：淋巴细胞白血病人生存数据Logistic模型
【摘要】：50名急性淋巴细胞白血病病人，在入院治疗时取得了外猿血中的细胞数x1（千个/mm3），淋巴结浸润等级x2（分为0，1，2，3级），出院后有无巩固治疗x3（1表示有巩固治疗，0 表示无巩固治疗）．通过随访取得病人的生存时间，并用变量y＝0表示生存1年以内，y＝1表示生存1年或1年以上．关于x1，x2，x3 和y 的观测数据，见表5-1．试用Logistic回归模型分析病人生存时间长短的概率与x1，x
50名急性淋巴细胞白血病病人，在入院治疗时取得了外猿血中的细胞数x1（千个/mm3），淋巴结浸润等级x2（分为0，1，2，3级），出院后有无巩固治疗x3（1表示有巩固治疗，0 表示无巩固治疗）．通过随访取得病人的生存时间，并用变量y＝0表示生存1年以内，y＝1表示生存1年或1年以上．关于x1，x2，x3 和y 的观测数据，见表5-1．试用Logistic回归模型分析病人生存时间长短的概率与x1，x2，x3 的关系．
表5-1 50名急性淋巴细胞白血病人生存数据
（续表）
以下用Logistic回归模型分析病人生存时间长短的概率与x1，x2，x3 的关系，其代码如下：
结果如下：
从以上结果可以看到，用三个变量x1，x2，x3 进行回归是无法通过检验的．
以下用“step（ ）”函数对变量进行筛选．代码如下：
＞glm．new＜－step(glm．sol)
结果如下：(www.chuimin.cn)
用“summary（ ）”函数显示细节，代码如下：
＞summary(glm．new)
结果如下：
从以上计算结果可以看出，去掉了x1 后，所有参数通过了检验（显著性水平为0．1时）．于是回归方程为
进行预测：
以上结果说明，巩固治疗比没有巩固治疗提高了0．4325522/0．04496518＝9．619715倍．
有关应用多元统计分析：基于R的实验的文章
- 详细阅读
-
广义线性模型概述：应用多元统计分析基于R实验详细阅读
，βp的线性函数．因此，对μ 作变换，则可得到下面几种分布的连接函数的形式：正态分布m（μ）＝μ＝∑βi xi．二项分布．Poisson分布m（μ）＝lnμ＝∑βi xi．上述推广体现在以下两个方面：通过一个连接函数，将响应变量的期望与解释变量建立线性关系m[E]＝β0＋β1 x1＋β2 x2＋…......
2023-11-18
-
应用R实验：曲线回归分析方法详细阅读
曲线回归分析的基本任务是通过两个变量x 和y 的实际观测数据建立曲线回归方程，以揭示x 和y 间的曲线关系的形式．常用的一种方法是：通过变量替换，把一元非线性回归问题转化为一元线性回归问题．曲线回归分析首要的工作是确定因变量y 与自变量x 之间曲线关系的类型．通常通过两个途径来确定：（1）利用有关专业知识，根据已知的理论规律和实践经验；（2）如果没有已知的理论规律和实践经验可以利用，可在直角坐标系......
2023-11-18
-
基于R的实验：iris数据集的聚类分析结果详细阅读
在实验2．3．2中曾对iris数据集进行描述和展示，在实验2．3．4中曾对iris数据集进行可视化．以下将对iris数据集进行进行聚类分析．根据实验2．3．4，我们只知道数据集内有三个品种鸢尾花而不知道每朵花的真正分类，只能凭借花萼及花瓣的长度和宽度分类．以下对鸢尾花iris数据集进行聚类分析，代码如下：＞data(iris);attach(iris)＞iris．hc1＜-hclust(dist(......
2023-11-18
-
数据库开发教程：E-R模型设计方法详细阅读
由图1.4.2可知，E-R图表达的是一个实体内部的联系和两个实体之间多个数量的联系。图1.4.2E-R图例1.4.2E-R模型的设计方法要设计一个E-R模型，应先确定实体集、实体属性和联系集。根据这些信息即可绘制E-R图。......
2023-12-04
-
R实验：对smoke数据集进行对应分析详细阅读
以下对smoke数据集进行对应分析．（1）首先查看smoke数据集的信息＞library(ca)＞data(＂smoke＂)＞smoke结果如下：这个数据集来自Greenacre（1984），被应用于多个统计软件作为对应分析的说明案例数据．它的内容是一个5行（阶层：SM，JM，SE，JE和SC）4列（吸烟习惯：none，light，medium 和heavy）的列联表，给出了一个虚构的公司内各阶层......
2023-11-18
-
多元统计分析：R实验中的数据变换方法详细阅读
，n）．记则其中，1p＝（1，1，…，pn·）表示对角线元素为p1·，p2·，…，pn· 的对角矩阵．因此，经过变换后所得到的新数据矩阵B，可以看成是由对应矩阵P 经过中心化和标准化后得到的矩阵．设用于检验行与列是否不相关的χ2 统计量为其中，表示第（i，j）单元在检验行与列两个属性变量否不相关时对总χ2 统计量的贡献，有其中，χ2＝＝T[tr]＝T[tr]＝T[tr]，tr表示方阵SQ 的迹．......
2023-11-18
-
使用R进行实验：分析ability.cov数详细阅读
ability．cov数据集是R 软件自带的数据集，以下对该数据集进行因子分析．（1）查看ability．cov数据集中的信息ability．cov数据集提供了Ability and Intelligence Tests（能力和智力测试）中，112个人参加的六个测试指标general（普通），picture（画图），blocks（积木），maze（迷津），reading（阅读），vocab（词汇）......
2023-11-18
相关推荐