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

news2024/9/22 15:38:47

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

点击查看原文

导师:张博士:来自中国科学院及重点高校资深专家,长期从事R语言模型、群落生态学、保护生物学、景观生态学和生态模型方面的研究和教学工作,以发表了多篇论文,拥有丰富的科研及实践经验。

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)模型及贝叶斯实现

 点击查看原文

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

<IBM DB2>《DB2 进程技术模型》

《DB2 进程技术模型》 1 概念说明2 引擎可分派单元EDU3 多线程体系结果优点4 协调代理程序5 防火墙6 客户机程序7 侦听器8 代理程序9 db2fmp10 db2vend10.1 数据库 EDU10.2 事件监视器线程的标识方式如下:10.3 备份和复原线程的标识方式如下: 11 数据库服…

苹果手机没有声音怎么回事?3分钟解决!

案例:苹果手机听不见声音怎么回事? 【朋友们,苹果手机没有声音,不知道我是不是按错了什么。】 如果你的苹果手机没有声音,可能是由于多种原因导致的。苹果手机没有声音怎么回事?看这里,下面是一…

(转载)简述马尔可夫链

赶紧记录一下,通俗易懂。 参考:https://zhuanlan.zhihu.com/p/448575579 马尔科夫链的思想:过去所有的信息都已经被保存到了现在的状态,基于现在就可以预测未来。(用数学方法就能解释自然变化的一般规律模型) 马尔科夫链为状态空…

ROS学习第二十四节——rosbag

1 rosbag使用_命令行 需求: ROS 内置的乌龟案例并操作,操作过程中使用 rosbag 录制,录制结束后,实现重放 实现: 1.准备 创建目录保存录制的文件 mkdir ./xxx cd xxx2.开始录制 -a:all,录制所有话题消息 -o:out&#xff0c…

领跑行泊一体,纵目科技剑指自动驾驶L2到L4的规模化商业落地机遇

‍数据智能产业创新服务媒体 ——聚焦数智 改变商业 2019年,通用、丰田、特斯拉等11家车企承诺自动驾驶时间表,他们大都表示在2020年底实现高级别自动驾驶。以特斯拉为例,其CEO埃隆马斯克曾承诺在2020年实现自动驾驶食言后,随后在…

【工作记录】centos7.5环境下通过源码编译方式部署mysql5.7.25

前言 本文介绍centos7.5环境下通过源码编译安装mysql5.7.25的过程及安装过程中遇到的问题解决。 一、准备工作 # 新建目录 # data 数据 config 配置 boost 引导文件 log 日志文件 mkdir -p /opt/mysql/data /opt/mysql/config /opt/mysql/boost /opt/mysql/log /opt/mysql/ …

Python小姿势 - # 基础数据结构与算法

基础数据结构与算法 Python中基础的数据结构与算法是非常重要的,它们可以帮助我们解决很多实际问题。今天我们就来学习一下Python中的基础数据结构与算法。 首先,我们先来了解一下数据结构。数据结构是一种用来存储、组织、处理数据的方法。它可以帮助我…

理解Java程序的执行

main 方法 public class Solution {public static void main(String[] args) {Person person new Person();person.hello();} }class Person {public void hello() {System.out.println("hello");} }源文件名是 Solution.java,这是因为文件名必须与 pub…

初刷leetcode题目总结 -- 数据结构

魔王的介绍:😶‍🌫️一名双非本科大一小白。魔王的目标:🤯努力赶上周围卷王的脚步。魔王的主页:🔥🔥🔥大魔王.🔥🔥🔥 ❤️‍&#x1…

Leetcode38. 外观数列

一、题目描述: 给定一个正整数 n ,输出外观数列的第 n 项。 「外观数列」是一个整数序列,从数字 1 开始,序列中的每一项都是对前一项的描述。 你可以将其视作是由递归公式定义的数字字符串序列: countAndSay(1) “…

4月19日第壹简报,星期三,农历闰二月廿九

4月19日第壹简报,星期三,农历闰二月廿九坚持阅读,静待花开1. 国家统计局:一季度GDP同比增长4.5%,当前中国经济没有通缩,下阶段也不会出现通缩,全国城镇调查失业率平均值为5.5%。2. 字节跳动发布…

海格里斯HEGERLS高速穿梭车按需定制|四向穿梭车货架和子母穿梭车货架别傻傻分不清?

随着物流行业和仓储行业的发展,越来越多的企业用户对仓储自动化程度要求越来越高。而近年来,各式各样的穿梭车AGV小车也现身各大物流展,备受各大中小型企业用户的青睐。且为了进一步提高仓库仓储的存储率,越来越多的仓储货架和仓储…

Docker之容器数据卷

Docker之容器数据卷 1. 容器数据卷的概念2. 数据卷的使用2.1 方式一2.2 测试22.3 方式二 3. 安装MySQL4. 具名挂载和匿名挂载4.1 匿名挂载4.2 具名挂载 5. 初识 Dockerfile6. 数据卷-容器之间 1. 容器数据卷的概念 将应用和环境打包成一个镜像数据?如果数据都在容器…

Java | 一分钟掌握JDK命令行工具 | 4 - 可视化分析工具

作者:Mars酱 声明:本文章由Mars酱编写,部分内容来源于网络,如有疑问请联系本人。 转载:欢迎转载,转载前先请联系我! 前言 我们其实在分析的时候,也并不是必须使用命令行工具才能可以…

八、vue-基础之列表渲染v-for、v-for中的key属性的作用

一、v-for列表渲染 在真实开发中,我们往往会从服务器拿到一组数据,并且需要对其进行渲染。 这个时候我们可以使用v-for来完成;v-for类似于JavaScript的for循环,可以用于遍历一组数据; 二、v-for基本使用 &#xff0…

Kubeadm方式搭建K8s高可用集群【1.23.0版本】

文章目录 一、集群规划及架构二、系统初始化准备(所有节点同步操作)三、安装kubeadm(所有节点同步操作)四、高可用组件安装及配置1、安装Nginx及配置2、安装keepalived及配置 五、初始化Master集群六、扩容K8S集群1、扩容master节点2、扩容node节点 七、安装网络组件Calico八、…

电脑文件剪切到U盘,为什么不见了?这4个技巧帮你找回丢失文件

文件剪切是指对计算机中的文件进行移动操作,将文件从原先的位置剪切到新的位置,然后将这些文件粘贴到新的位置,以完成文件的整个剪切操作。虽然剪切操作起来并不难,但是近日,有位小伙伴遇到了将电脑文件剪切到u盘后出现…

【C++】模板(二)

文章目录 非类型模板参数简单对容器array(c11)介绍及对比模板特化函数模板特化类模板特化全特化偏特化 模板分离模板总结 非类型模板参数 模板参数分类类型形参与非类型形参。 类型形参即:出现在模板参数列表中,跟在class或者type…

文献智能管理工具

当下载的文献非常多的时候,对文献的管理将是一个非常耗时且困难的工作。困难主要来自两个方面,一是文献的分类,当文献属于多个类别的时候,究竟应该放在哪个类别的文件夹下,还是说每个类别的文件夹下都放一个该文献的备…

传输层重点协议之【TCP协议】

1. TCP协议段格式 2. TCP原理 2.1 可靠性机制 2.1.1 确认应答 确认应答是实现可靠性的最核心机制 首先来看一个例子,下面是我和女神的对话~ 所以为了解决上述问题呢,就需要针对消息进行编号!给发送的消息分配一个”序号“,同时…