【SAS求解多元回归方程】REG多元回归分析-多元二次回归

多元二次回归是线性回归的扩展形式,通过引入自变量的二次项和交互项来捕捉非线性关系。适用于因变量与多个自变量之间存在曲线或交互效应的情况。其模型形式比简单线性回归更复杂,能够拟合更灵活的数据模式。

 

目录

【示例】

数据预处理

基本语法 

SAS代码

逐步选择

​编辑

参数估计

方差分析

回归统计量

y的拟合诊断

y的回归变量值-残差值


【示例】

设Y与x_1,x_2,x_3有相关关系,考虑二次回归模型y=\beta _0+\beta _1x_1+\beta _2x_2+\beta _3x_3+\beta _4x_1^{2}+\beta _5x_2^{2}+\beta _6x_3^{3}+\beta _7x_1x_2+\beta _8x_1x_3+\beta _9x_2x_3+\varepsilon

,8组观测数据如下:

 

 数据预处理

根据二次模型要求处理数据集得到新数据集

 

基本语法 

 PROC REG data = 数据集;
MODEL 因变量 = 自变量列表 </可选项>;
< restrict 自变量的等式约束;>

 

 SAS代码

针对处理后的数据

data d1;input x1-x9 y ;cards;
38 47.5 23 1444 2256.25 529 1805.00 1805.00 1092.50 66.0 
41 21.3 17 1681  453.69 289	 873.30  873.30  362.10 43.0 
34 36.5 21 1156 1332.25 441 1241.00 1241.00  766.50 36.0 
35 18.0 14 1225	 324.00 196	 630.00  630.00  252.00 23.0 
31 29.5 11  961  870.25 121	 914.50  914.50  324.50 27.0 
34 14.2  9 1156	 201.64 81   482.80  482.80  127.80 14.0 
29 21.0  4  841  441.00 16   609.00  609.00  84.00 	12.0 
32 10.0  8 1024	 100.00 64   320.00  320.00  80.00 	7.6 
;
proc print;
run;
proc reg data=d1;model y=x1-x9 /selection=stepwisesle=0.05 sls=0.05;
run;
proc reg data=d1;model y=x4 x7;
run;
quit;

 

逐步选择

第一步

第二步

直到留在模型中的所有变量的显著性水平都为 0.0500,同时没有其他变量满足 0.0500 显著性水平。

得到x4、x7,对应于x_1^{2},x_1x_2

 

参数估计

所以二次回归方程为:

\hat{Y}=-30.0098+0.02672x_1^{2}+0.03130x_1 x_2

参数估计表不仅给出回归方程的系数,还给出检验H_0的结果(显著性概率p值)

比如给定\alpha =0.05, 若常数项和自变量的p值均\geqa,意味着与回归方程高度显著产生矛盾,为了得到最优回归方程,应从方程中删除最不重要的自变量,重新建立Y与其余自变量的回归方程后再检验,这就是变量筛选的意义所在。

方差分析

回归平方和:U=2597.3352

残差平方和:Q=27.1798

均方误差:MSE=27.1798/5=5.4360

均方误差是模型中误差方差\sigma ^{2}的估计

检验统计量F=238.90,显著性概率p值小于0.0001,这表示拟合的模型是高度显著的,该模型解释了这组数据总变差中的主要部分。

回归统计量

决定系数:R^{2}=0.9896

复相关系数:R=\sqrt{0.9896}

标准差\sigma的估计量:ROOTMSE=2.3315

y的拟合诊断

y的回归变量值-残差值

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mfbz.cn/a/439.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

多参表达式Hive UDF

支持的操作符 &#xff1a;跳过&#xff0c;即无条件筛选&#xff1a;等于!&#xff1a;不等于range&#xff1a;区间内&#xff0c;range[n,m]表示 between n and mnrange&#xff1a;区间外&#xff0c;即not between andin&#xff1a;集合内&#xff0c;in(n,m,j,k)表示 in…

Flink读取Kafka写入Paimon

Flink SQL -- 1&#xff09;注册 Paimon 源 CREATE CATALOG paimon_hive WITH(type paimon,warehouse hdfs://xxxxx/paimon,metastore hive,hive-conf-dir /xxxxx/conf,uri thrift://域名1:9083,thrift://域名2:9083);-- 2&#xff09;声明 Kafka 源 create table kafkaS…

【开源工具】:基于PyQt5的智能网络驱动器映射工具开发全流程(附源码)

&#x1f517; 【开源工具】&#xff1a;基于PyQt5的智能网络驱动器映射工具开发全流程 &#x1f308; 个人主页&#xff1a;创客白泽 - CSDN博客 &#x1f525; 系列专栏&#xff1a;&#x1f40d;《Python开源项目实战》 &#x1f4a1; 热爱不止于代码&#xff0c;热情源自每…

MySQL 索引学习笔记

1.二叉树&#xff0c;红黑树&#xff0c;B 树&#xff0c;B树 二叉树&#xff1a;就是每个节点最多只能有两个子节点的树&#xff1b; 红黑树&#xff1a;就是自平衡二叉搜索树&#xff0c;红黑树通过一下五个规则构建&#xff1a; 1.节点只能是红色或黑色&#xff1b; 2.根…

嵌入式通信模块实战新范式:基于虚拟仿真平台的NB-IoT核心技能训练——零硬件损耗的全栈式实验方案,重构物联网通信教学逻辑

在万物智联时代&#xff0c;NB-IoT通信模块已成为低功耗广域网的基石。BC260Y作为行业主流模组&#xff0c;其AT指令控制与网络诊断能力是嵌入式开发者的必备技能。传统教学受限于硬件采购成本、设备管理难度及实验风险&#xff0c;难以开展规模化训练。嵌入式仿真实验教学平台…

docker compose的变量使用说明

澄清一下 x-shared-env 和 &shared-api-worker-env 的作用范围&#xff1a; 核心概念&#xff1a;Docker Compose 配置 vs 容器环境 x-shared-env: &shared-api-worker-env 是 Docker Compose 配置的一部分 这些定义仅在 Docker Compose 解析 YAML 文件时 有效它们定义…

美团完整面经

面试岗位 面试的岗位 - 2025春季校招 【转正实习】软件服务工程师-后端方向&#xff08;成都 - 软硬件服务-SaaS事业部&#xff09; 一面&#xff08;业务初试 - 30min&#xff09; 问题 自我介绍 Java基础 HashMap底层用的数据结构是什么&#xff1f;是线程安全的吗&…

JAVA毕业设计227—基于SpringBoot+hadoop+spark+Vue的大数据房屋维修系统(源代码+数据库)

毕设所有选题&#xff1a; https://blog.csdn.net/2303_76227485/article/details/131104075 基于SpringBoothadoopsparkVue的大数据房屋维修系统(源代码数据库)227 一、系统介绍 本项目前后端分离&#xff0c;分为业主、维修人员、管理员三种角色 1、业主&#xff1a; 登…

uniapp 页面栈一定深度后,回首页导航到新页面的解决方案

uniapp 页面栈一定深度后&#xff0c;回首页导航到新页面的解决方案 uniapp 页面导航解决方案 在 uniapp 中&#xff0c;要实现先弹出页面栈回到首页&#xff0c;然后再跳转到指定页面。 /*** description 后台选择链接专用跳转*/ interface Link {path: string;name?: stri…

java实现Google邮箱SMTP协议

一、开通Google的SMTP协议 在谷歌邮箱中开启IMAP访问 到google的设置中开启两步验证功能 在到 创建和管理应用专用密码 二、java中实现 引入maven <!--邮件--><dependency><groupId>com.sun.mail</groupId><artifactId>javax.mail</artif…

【2025最新】Adobe Illustrator下载保姆级安装教程(附官方下载链接)

文章目录 Adobe Illustrator 2024新功能介绍如何提高Adobe Illustrator的运行效率 Adobe Illustrator 这款神器相信不用我多介绍了吧&#xff0c;设计师们的得力助手&#xff01;最新的2025版据说功能和体验都提升了不少。这篇呢&#xff0c;算是我个人整理的一个超详细adobe i…

2025.06.11【Ribo-seq】|根据注释文件获取外显子及ORF序列

文章目录 一、准备材料二、提取外显子区间为BED文件1. 提取GTF中exon为BED 三、用bedtools提取外显子fasta四、后续拼接外显子为ORF序列五、流程总结 一、准备材料 基因组fasta&#xff08;如&#xff1a;genome.fa&#xff09;RiboCode生成的GTF文件&#xff08;如&#xff1…