R语言混合效应(多水平/层次/嵌套)模型及贝叶斯实现技术

回归分析是科学研究中十分重要的数据分析工具。随着现代统计技术发展,回归分析方法得到了极大改进。混合效应模型(Mixed effect model),即多水平模型(Multilevel model)/分层模型(Hierarchical Model)/嵌套模型(Nested Model),无疑是现代回归分析中应用最为广泛的统计模型,代表了现代回归分析主流发展方向。混合效应模型形式灵活可以应对现代科学研究中各种数据情况,与传统回归模型相比具有更为强大数据分析能力,且结果更为可信。将分为复杂数据的回归及混合效应模型概述及数据探索;回归与混合效应模型,包括一般线性回归(lm)、广义线性回归(glm);线性混合效应模型(lmm)及广义线性混合效应模型(glmm);贝叶斯(brms)回归与混合效应模型;相关数据回归与混合效应模型及贝叶斯实现,包括嵌套数据、时间自相关数据,空间自相数据及系统发育数据分析;非线性数据回归分析及贝叶斯实现,包括广义可加(混合)模型和非线性(混合)模型等。将通过大量实例讲解,使大家能应对科研工作中复杂数据局面,选择合适模型,提高数据分析能力。        

查看原文

内容详情:

课前预习:复杂数据的回归及混合效应模型概述及数据探索由Ai尚研修提供教材】

1复杂数据回归模型的选择策略

1)科学研究中数据及其复杂性

2)回归分析历史、理论基础

3)回归分析基本假设和常见问题

4)复杂数据回归模型选择策略

2如何通过数据探索避免常见统计问题

1)数据缺失(missing value)

2)零值(zero trouble)

3)奇异值/离群值(outliers)

4)异质性(heterogeneity)

5)数据分布正态性(normality)

6)响应变量与预测变量间关系(relationships)

7)交互作用项(interaction)

8)共线性(collinearity)

9)样本独立性(independence)

专题一:回归与混合效应(多水平/层次/嵌套)模型

1.1一般线性模型(lm)

1)基本形式、基本假设、估计方法、参数检验、模型检验

2)一般线性回归、方差分析及协方差分析

3)一般线性回归模型验证

4)一般线性回归模型选择-逐步回归

案例1:鱼类游速与水温关系的回归及协方差分析;

案例2:施肥和种植密度对作物产量的影响

案例3:决定海洋植食性鱼类多样性的决定因子-模型验证

案例4:淡水鱼丰度的环境因子的筛选-逐步回归

1.2广义线性模型(glm)

1) 基本形式、基本假设、估计方法、参数检验、模型检验

2) 0,1数据分析:伯努利分布、二项分布及其过度离散问题

3)计数数据各种情况及模型选择:泊松、伪泊松、负二项、零膨胀泊松、零膨胀负二项、零截断泊松及零截断负二项模型

4) 广义线性模型的模型比较和选择-似然比LR和AIC

案例1:动物身体特征与患病与否(0,1)的关系的逻辑斯蒂回归

案例2:海豹年龄与攻击行为的关系-0,1数据转化为比率数据分析

案例3:不同实验处理下蚜虫多度的差异分析-计数数据泊松回归

其他案例:零膨胀、零截断数据分析。。。。。。。。。。。。。。。。。。。。

1.3线性混合效应模型(lmm)

1) 线性混合效应模型基本原理

2) 线性混合效应模型建模步骤及实现

3) 线性混合效应模型的预测和模型诊断

4) 线性混合效应模型的多重比较

案例1:睡眠时间与反应速度关系

案例2:多因素实验(分层数据)的多重比较

1.4广义线性混合效应模型(glmm)

1)广义线性混合效应模型基本原理

2)广义线性混合效应模型建模步骤及流程

3)广义线性混合效应模型分析0,1数据

4)广义线性混合效应模型分析计数数据及模型选择:泊松、伪泊松、负二项、零膨胀泊松、零膨胀负二项、零截断泊松及零截断负二项模型

案例1:蝌蚪“变态”与否(0,1)的多因素分析-逻辑斯蒂混合效应模型

案例2:虫食种子多度影响因素的多变量分析-泊松混合效应模型

案例3:模拟计数数据-零膨胀、零截断、过度离散等广义混合效应模型

专题二:贝叶斯(brms)回归与混合效应(多水平/层次/嵌套)模型

2.1贝叶斯回归及混合效应模型上

1)贝叶斯回归分析简介

2)利用brms实现贝叶斯回归分析简介

3)贝叶斯回归分析的模型诊断、交叉验证、预测和作图

4)贝叶斯广义线性模型实现:gamma分布、伯努利分布、二项分布等

案例1:鱼游速与温度关系的贝叶斯回归-结果解读、模型验证、模型诊断

案例2:森林生物量与林龄关系贝叶斯回归-gamma分布、brms参数调整

案例3:动物身体特征与患病与否(0,1)的关系的贝叶斯回归-伯努利分布

案例4:海豹年龄与攻击行为的关系-0,1数据转化为比率数据分析-二项分布

其他案例:贝叶斯分析计数数据过度离散、零膨胀等问题

2.2贝叶斯回归及混合效应模型下

1)贝叶斯线性混合效应模型:实现步骤、模型验证、多重比较

2)贝叶斯广义混合效应模型-计数数据分析:泊松、负二项、零膨胀泊松、零膨胀负二项等

案例1:睡眠时间与反应速度关系的贝叶斯线性混合效应模型

案例2:教师受欢迎程度的多变量预测-贝叶斯线性混合效应模型

案例3:虫食种子多度(计数数据)影响因素的多变量分析-贝叶斯广义混合效应模型

其他案例:贝叶斯分析计数数据过度离散、零膨胀等问题

专题三:相关数据回归分析:嵌套、时间、空间、系统发育相关数据分析

3.1嵌套型随机效应混合效应模型分析及贝叶斯实现

1)数据分层问题及嵌套型随机效应混合效应模型介绍

2)嵌套型随机效应混合效应模型分析步骤及流程及模型选择(MuMIn)

3)嵌套型随机效应混合效应模型的方差分解:ICC、varcomp及贝叶斯法

4)经典方差分解案例讲解

案例1:不同种类海豚年龄多因素预测模型及模型选择(MuMIn)- 嵌套结构

案例2:纲/科/属/种型嵌套随机效应的方差分解及贝叶斯方法

案例3:物种属性可塑性和基因多样性对物种丰富度影响的相对贡献-全模型变差分解

3.2时间相关数据分析及贝叶斯实现

1)回归模型的方差异质性问题及解决途径

2)时间自相关分析:线性及混合效应模型及贝叶斯方法

3)时间自相关+方差异质性分析及贝叶斯实现

案例1:模拟数据方差异质性问题-gls,lmm及brms方法比较

案例2:鸟类多度变化的时间自相关分析-gls vs brms

案例3:资源脉冲与食谱关系分析:方差异质性+时间相关-lmm vs brms

3.3空间相关数据分析及贝叶斯实现

1)空间自相关概述

2)空间自相关问题解决方式:自相关修正参数、空间距离权重法、空间邻接权重法

3)空间自相关问题修正基本流程-gls和lme

4)空间自相关贝叶斯修正-空间距离权重 VS 空间邻接权重

案例1:北方林物种多样性与气候关系-一般线性回归模型空间自相关问题修正

案例2:全球水鸟巢穴捕食率影响因素分析-混合效应模型空间自相关问题修正

3.4系统发育相关数据分析及贝叶斯实现

1、系统发育简介:系统发育假说、系统发育信号及系统发育树

2、系统发育树及系统发育距离矩阵构建

3、系统发育信息纳入回归模型-广义最小二乘(gls)

4、系统发育信息纳入混合效应模型(lmm/glmm)及贝叶斯方法实现案例

案例1:模拟数据-系统发育相关对物种属性影响-gls vs brms

案例2:全球水鸟巢穴捕食率影响因素分析-系统发育混合效应模型:lmm vs brms

专题四:非线性关系数据分析:广义可加(混合)模型(GAM/GAMM)和非线性(混合)(NLM/NLMM)模型

4.1“线性”回归的含义及非线性关系的判定

4.2广义可加(混合效应)(GAM/GAMM)模型及贝叶斯实现

4.3非线性(混合效应)(NLM/NLMM)模型及贝叶斯实现

 

推荐:

A类

 基于R语言的现代贝叶斯统计学方法(贝叶斯参数估计、贝叶斯回归、贝叶斯计算)

B类

基于R语言BIOMOD2模型的物种分布模拟实践技术应用

C类

基于R语言的Copula变量相关性分析及应用

D类

基于R语言的Meta分析【全流程、不确定性分析】方法与Meta机器学习

E类

基于R语言:INLA下的贝叶斯回归\生存分析\随机游走、广义可加模型\极端数据的贝叶斯分析等实践

点击查看原文 

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

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

相关文章

如何快速运用R语言实现生物群落(生态)数据统计分析与绘图

R 语言作的开源、自由、免费等特点使其广泛应用于生物群落数据统计分析。生物群落数据多样而复杂,涉及众多统计分析方法。本次以生物群落数据分析中的最常用的统计方法回归和混合效应模型、多元统计分析技术及结构方程等数量分析方法为主线,通过多个来自…

Linux:查看进程。

Linux:查看进程。 windows linux TTY如果是?说明是不是终端(控制台)启动的,而是系统内部自己启动的。 TIME是启动Linux后,这个进程一共占用了cpu多少时间00…

QT 设计ROS GUI界面订阅和发布话题

QT 设计ROS GUI界面订阅和发布话题 主要参考下面的博客 ROS项目开发实战(三)——使用QT进行ROS的GUI界面设计(详细教程附代码!!!) Qt ROS 相关配置请看上一篇博客 首先建立工作空间和功能包&a…

【探索】机器指令翻译成 JavaScript

前言 前些时候研究脚本混淆时,打算先学一些「程序流程」相关的概念。为了不因太枯燥而放弃,决定想一个有趣的案例,可以边探索边学。 于是想了一个话题:尝试将机器指令 1:1 翻译 成 JavaScript,这样就能在浏览器中&am…

Java程序设计入门教程-- if 条件语句

目录 单分支选择语句(if) 双分支选择语句(if…else) 嵌套if语句 单分支选择语句(if) 情形 当判断条件满足时,执行语句体S,而不满足则什么都不做。 格式 if (条件判断表…

【计算机视觉 | 目标检测】术语理解6:ViT 变种( ViT-H、ViT-L ViT-B)、bbox(边界框)、边界框的绘制(含源代码)

文章目录 一、ViT & ViT变种1.1 ViT的介绍1.2 ViT 的变种 二、bbox(边界框)三、边界框的绘制 一、ViT & ViT变种 1.1 ViT的介绍 ViT,全称为Vision Transformer,是一种基于Transformer架构的视觉处理模型。传统的计算机视…

java企业工程项目管理系统平台源码(三控:进度组织、质量安全、预算资金成本、二平台:招采、设计管理)

工程项目管理软件(工程项目管理系统)对建设工程项目管理组织建设、项目策划决策、规划设计、施工建设到竣工交付、总结评估、运维运营,全过程、全方位的对项目进行综合管理 工程项目各模块及其功能点清单 一、系统管理 1、数据字典&#…

回调函数与钩子函数的区别,另QT中connect函数的实现,lambda的使用

1、钩子函数是回调函数的一种 广泛来说两者都是一样的 严格来说 钩子函数的函数名早已被定义好,只是函数内部需要用户在应用层来定义, 1)可以完全通过宏来实现系统是否调用该函数(底层不封闭,修改宏的参数实现是否编…

【2023 · CANN训练营第一季】MindSpore模型快速调优攻略 第二章——MindSpore调试调优

1.生态迁移 生态迁移工具使用示例 生态迁移工具技术方案 不同框架间模型定义前端表达差别巨大(相同算子的API技术难点 、 算子功能、模型构建方式差别较大); 对于同一框架,不管前端表达差异如何,最终对应的计算 图是相似的。因此提出&#x…

Kubernetes部署+kubesphere管理平台安装

Kubernetes官网;kubesphere官网 不论是Kubernetes官网还是找的其它部署步骤,基本都是推荐搭建集群的方式,是为了实现高可用.....等等,这样一来至少需要两台或三台的服务器来搭建,这样对我们的成本也是非常大的&#xf…

Axure教程——直方图(中继器)

本文将教大家如何用AXURE用中继器制作直方图 一、效果介绍 如图: 预览地址:https://yjkepz.axshare.com 下载地址:https://download.csdn.net/download/weixin_43516258/87842701 二、制作方法 (1)制作刻度表 设计5个刻…

CSDN上海城市开发者社区线下活动纪实

引言 5月27号中午,很高兴能和现CSDN副总裁、前微软 Azure 工程团队首席研发经理、技术畅销书《编程之美》及《构建之法》的作者邹欣邹老师,以及CSDN的 “上海城市开发者社区” 的部分成员齐聚一堂,参加CSDN上海城市开发者社区自5月初成立以来…

Ctfshow基础二刷(1)

前言: 前两天的信安给我整emo了,头一回打正经比赛,结果发现基础太差,代码审计烂得一踏糊涂。 寻思寻思,从头整一遍基础。又买了安恒出的新书。争取7号去吉林打省队选拔不给导儿丢脸吧呜呜 文件包含 web78: 这题一…

curl 命令-接口测试

curl 命令-接口测试 JUST DO IT 温暖春日 在linux/Unix 为代表的os上, 对后端进行测试, 模拟连接请求都会书写脚本 场景: 在Linux 上接口测试工具有ab, restClient, postman等, 最常用的方法是curl进行简单测试 curl是非常方便的Rest 客户端, 可以很方便的完成 Rest API测…

AcWing算法提高课-1.3.11二维费用的背包问题

宣传一下算法提高课整理 <— CSDN个人主页&#xff1a;更好的阅读体验 <— 本题链接&#xff08;AcWing&#xff09; 点这里 题目描述 有 N N N 件物品和一个容量是 V V V 的背包&#xff0c;背包能承受的最大重量是 M M M。 每件物品只能用一次。体积是 v i v_…

vcruntime140.dll丢失怎么办?怎么解决vcruntime140.dll丢失的问题

当您运行一个需要此文件的程序时&#xff0c;如果您的系统中不存在这个文件&#xff0c;会提示出错信息“找不到vcruntime140.dll”或“vcruntime140.dll丢失”。这种情况下&#xff0c;您需要解决这个问题&#xff0c;才能继续运行此应用程序。我们将介绍vcruntime140.dll丢失…

5.27下周黄金行情走势预测及开盘操作策略

近期有哪些消息面影响黄金走势&#xff1f;下周黄金多空该如何研判&#xff1f; ​黄金消息面解析&#xff1a;周五(5月26日)黄金大幅下跌&#xff0c;主要受到美国数据影响&#xff0c;美国公布的4月PCE和耐用品订单数据向好&#xff0c;再次强化市场对美联储的鹰派押注。现货…

(四)ArcGIS空间数据的转换与处理——数据结构转换

ArcGIS空间数据的转换与处理——数据转换 空间数据的来源很多&#xff0c;如地图、工程图、规划图、航空与遥感影像等&#xff0c;因此空间数据也有多种格式。根据应用需要&#xff0c;需对数据进行格式转换&#xff0c;不同数据结构间的转换主要包括矢量数据到栅格数据的转换…

戏曲APP软件开发需具备哪些功能呢?

戏曲是我国的国粹&#xff0c;传统戏曲文化源远流长&#xff0c;博大精深&#xff0c;数千年以来一直都是深受大众喜欢的文化生活的重要环节。随着时代的推进&#xff0c;娱乐形式更加多样化&#xff0c;传统的剧场演出形式的戏曲传播方式已经跟不上时代发展以及人们的需求了。…

Office project 2016安装

哈喽&#xff0c;大家好。今天一起学习的是project 2016的安装&#xff0c;Microsoft Office project项目管理工具软件&#xff0c;凝集了许多成熟的项目管理现代理论和方法&#xff0c;可以帮助项目管理者实现时间、资源、成本计划、控制。有兴趣的小伙伴也可以来一起试试手。…