拓端小助手的博客_原创(37)

R语言实现有限混合模型建模分析

原文链接：http://tecdat.cn/?p=6129介绍有限混合模型在应用于数据时非常有用，其中观察来自不同的群体，并且群体隶属关系未知。模拟数据首先，我们将模拟一些数据。让我们模拟两个正态分布 - 一个平均值为0，另一个平均值为50，两者的标准差为5。m1 <- 0m2 <- 50sd1 <- sd2 <- 5N1 <- ...

R语言

编程开发

原创 2021-05-12 14:13:04 281 阅读

ggplot2如何在R语言中绘制表格

原文连接：http://tecdat.cn/?p=6154偶尔我想在R中的图表旁边绘制一个表格，例如，以显示图表本身的摘要统计数据。这非常简单。该函数tableGrob创建像一个数据帧的曲线图的表，安排ggplot2图形对象的网页上。这是一个小例子：会话信息R version 3.2.1 (2015-06-18)Platform: x86_64-apple-...

R语言

编程开发

原创 2021-05-12 14:12:52 491 阅读

R语言回归中的Hosmer-Lemeshow拟合优度检验

原文链接：http://tecdat.cn/?p=6166在依赖模型得出结论或预测未来结果之前，我们应尽可能检查我们假设的模型是否正确指定。也就是说，数据不会与模型所做的假设冲突。对于二元结果，逻辑回归是最流行的建模方法。在这篇文章中，我们将看一下 Hosmer-Lemeshow逻辑回归的拟合优度检验。Hosmer-Lemeshow拟合优度检验Hosmer-Lemeshow...

R语言

编程开发

原创 2021-05-12 14:12:51 4759 阅读

R语言stan进行基于贝叶斯推断的回归模型

原文连接：http://tecdat.cn/?p=6252R的Stan可以从许多统计软件包中运行Stan。到目前为止，我一直在从R运行Stan，首先按照快速入门指南中的说明安装并运行所有内容。简单线性回归第一步是为Stan模型编写文件。这包含一个文件linreg.stan：data { intN; [N] x; vector[N] y; } parameter...

R语言

编程开发

原创 2021-05-12 14:12:50 865 阅读

R语言t检验和非正态性的鲁棒性

原文连接：http://tecdat.cn/?p=6261t检验是统计学中最常用的检验之一。双样本t检验允许我们基于来自两组中的每一组的样本来测试两组的总体平均值相等的零假设。这在实践中意味着什么？如果我们的样本量不是太小，如果我们的数据看起来违反了正常假设，我们就不应过分担心。此外，出于同样的原因，即使X不正常（同样，当样本量足够大时），组均值差异的95％置信区间也将具有正确的覆盖...

R语言

编程开发

原创 2021-05-12 14:12:49 914 阅读

R语言Wald检验 vs 似然比检验

原文链接：http://tecdat.cn/?p=6895在开展基于概率推理的课程时，关键主题之一是基于似然函数的检验和置信区间构建。通常包括Wald，似然比和分数检验。在这篇文章中，我将修改Wald和似然比检验的优缺点。我将重点关注置信区间而不是检验。示例我们将X表示观察到的成功次数的随机变量，x表示其实现的值。似然函数只是二项式概率函数，但参数是模型参数。所以ML...

R语言

编程开发

原创 2021-05-12 14:12:48 2789 阅读

R语言用于线性回归的稳健方差估计

原文链接：http://tecdat.cn/?p=6274在这篇文章中，我们将看看如何在实践中使用R 。为了说明，我们首先从线性回归模型中模拟一些简单数据，其中残差方差随着协变量的增加而急剧增加： n < - 100 x < - rnorm（n）residual_sd < - exp（x）y < - 2 * x + residual_sd ...

R语言

编程开发

原创 2021-05-12 14:12:47 801 阅读

R语言在逻辑回归中求R square R方

原文链接：http://tecdat.cn/?p=6295并非所有结果/因变量都可以使用线性回归进行合理建模。也许第二种最常见的回归模型是逻辑回归，它适用于二元结果数据。如何计算逻辑回归模型的R平方？麦克法登R平方在R中，glm（广义线性模型）命令是用于拟合逻辑回归的标准命令。据我所知，拟合的glm对象并没有直接给你任何伪R平方值，但可以很容易地计算出McFadden的...

R语言

编程开发

原创 2021-05-12 14:12:46 3698 阅读 1收藏

R语言Poisson回归的拟合优度检验

原文链接：http://tecdat.cn/?p=6304在这篇文章中，我们将看一下Poisson回归的拟合优度测试与个体计数数据。许多软件包在拟合Poisson回归模型时在输出中提供此测试，或者在拟合此类模型（例如Stata）之后执行此测试，这可能导致研究人员和分析人员依赖它。在这篇文章中，我们将看到测试通常不会按预期执行，因此，我认为，应该谨慎使用。偏差拟合度检验...

R语言

编程开发

原创 2021-05-12 14:12:45 626 阅读

R语言ROC曲线下的面积 - 评估逻辑回归中的歧视

原文链接：http://tecdat.cn/?p=6310在讨论ROC曲线之前，首先让我们在逻辑回归的背景下考虑校准和区分之间的区别。良好的校准是不够的对于模型协变量的给定值，我们可以获得预测的概率。如果观察到的风险与预测的风险（概率）相匹配，则称该模型已被很好地校准。也就是说，如果我们要分配一组值的大量观察结果，这些观察结果的比例应该接近20％。如果观察到的...

R语言

编程开发

原创 2021-05-12 14:12:43 862 阅读

R语言使用 LOWESS技术图分析逻辑回归中的函数形式

原文链接：http://tecdat.cn/?p=6322当我们在回归模型中包含连续变量作为协变量时，重要的是我们使用正确的（或近似正确的）函数形式。例如，对于连续结果Y和连续协变量X，可能是Y的期望值是X和X ^ 2的线性函数，而不是X的线性函数。一种简单但通常有效的方法是简单地查看Y对X的散点图，以直观地评估。对于我们通常使用逻辑回归建模的二元结果，事情并不那么容易（至少在...

R语言

编程开发

原创 2021-05-12 14:12:42 340 阅读

R语言中固定与随机效应Meta分析 - 效率和置信区间覆盖

原文链接：http://tecdat.cn/?p=6330Meta分析是综合现有证据的重要工具。它通常用于医学和临床环境中，以评估关于治疗或暴露对感兴趣结果的影响的现有证据。基本思想是将先前研究中感兴趣的影响的估计汇总在一起。在进行Meta分析时必须做出的选择是固定效应和随机效应之间的选择。在这篇文章中，我们将看看这种选择的一些后果，实际上研究正在测量不同的效果。设置我们...

R语言

编程开发

原创 2021-05-12 14:12:41 354 阅读

R语言对苏格兰独立民意调查的Meta分析

原文链接：http://tecdat.cn/?p=63352014年9月18日，苏格兰居民将投票进行公民投票决定是否独立于英国。虽然运动之前已经保持了相当健康的领先优势，但是根据投票意向的民意调查，竞选已经大大收紧。特别是，两个民意调查现在已经显示出更大的比例，表示他们将投票赞成投票比例。由于上周进行了一系列民意调查，结果略有不同，我决定对其进行简单的Meta分析。根据现有证据，调查结果...

R语言

编程开发

原创 2021-05-12 14:12:40 183 阅读

在Amazon Web Services中使用R语言运行模拟

原文链接：http://tecdat.cn/?p=6339我最近一直致力于R中的一些模拟研究，涉及计算机密集型MCMC采样。通常我会使用我所在机构的计算集群来做这些，利用大量的计算机核心，但暂时缺乏这种可用性导致我调查使用亚马逊的Web服务（AWS）系统。在这篇文章中，我将描述我在R中进行模拟的步骤。使用Amazon Web Services进行初始设置首先注册https...

R语言

编程开发

原创 2021-05-12 14:12:39 193 阅读

R语言分析协变量之间的非线性关系

原文链接：http：//tecdat.cn/？p = 6366最近我被问到我的 - [R和Stata的软件包是否能够适应协变量之间的非线性关系。答案是肯定的，在这篇文章中，我将说明如何做到这一点。为了说明，我们将模拟具有两个协变量X1和X2以及连续结果ý的非常大的数据集。 set.seed（123）n < - 10000 x1 < - rnorm（...

R语言

编程开发

原创 2021-05-12 14:12:38 180 阅读

R语言用多重插补法估算相对风险

原文链接：http://tecdat.cn/?p=6379在这里，我将用R中的一个小模拟示例进行说明。首先，我们使用X1和X2双变量法线和Y模拟大型数据集，其中Y遵循给定X1和X2的逻辑模型。首先，我们模拟一个非常大的完整数据集：#simulate完整数据expit < - function（x）{ EXP（X）/（1 + EXP（X））}n < -...

R语言

编程开发

原创 2021-05-12 14:12:36 357 阅读

R语言调整随机对照试验中的基线协变量

原文链接：http://tecdat.cn/?p=6386随机对照试验构成通常被认为是用于评估某些干预或感兴趣治疗效果的金标准设计。参与者被随机分配到两个（有时更多）的群体这一事实确保了，至少在期望中，两个治疗组在测量的，重要的是可能影响结果的未测量因素方面是平衡的。因此，两组之间结果的差异可归因于随机化治疗而不是对照（通常是另一种治疗）的效果。如果随机化没有受到影响，即使不调...

R语言

编程开发

原创 2021-05-12 14:12:35 412 阅读

R语言利用基线协变量提高随机对照试验的效率

原文链接：http://tecdat.cn/?p=6391在这篇文章中，我将介绍基于半参数理论的最近开发的改进边际治疗效果估计的方法之一。边际处理与条件处理效应未经调整和调整后的治疗效果不同的一个重要例子是使用逻辑回归来模拟二元结果。也就是说，治疗效果的边际或未调整比值比不同于以一个或多个基线协变量为条件的治疗效果。这意味着如果调整基线测量，真实治疗效果估计实际上与边际未调...

R语言

编程开发

原创 2021-05-12 14:12:34 243 阅读

R语言在RCT中调整基线时对错误指定的稳健性

原文链接：http://tecdat.cn/?p=6400众所周知，调整一个或多个基线协变量可以增加随机对照试验中的统计功效。调整分析未被更广泛使用的一个原因可能是因为研究人员可能担心如果基线协变量的影响在结果的回归模型中没有正确建模，结果可能会有偏差。建立我们假设我们有关于受试者的双臂试验的数据。对于第i个主题，我们记录基线协变量和结果。我们让表示受试者是否被随机分配到...

R语言

编程开发

原创 2021-05-12 14:12:33 157 阅读

R语言使用倾向评分提高RCT（随机对照试验）的效率

原文链接：http://tecdat.cn/?p=6404倾向评分已成为观察性研究中混杂因素调整的常用方法。基本思想是模拟接受治疗或暴露的概率如何取决于混杂因素，即要治疗的“倾向”。首先要注意的是，人们不会认为倾向评分在RCT中起作用。如上所述，倾向评分用于调整观察性研究中的混淆。在RCT中，随机化确保治疗和其他基线变量在统计学上是独立的，即没有混淆。那么倾向得分有什么用呢？...

R语言

编程开发

原创 2021-05-12 14:12:32 214 阅读

R语言多臂试验 - 我们应该考虑多重性吗？

原文链接：http://tecdat.cn/?p=6408我们应该担心多重性吗？事实证明（或许不足为奇）在多臂试验的设定中，某种多样性调整是否合适的问题一直备受争议。一种观点认为，如果不同的假设代表不同的研究问题，那么不允许进行多重比较是合理的。模拟研究为了进行实证研究，我根据Parmar等人的假设进行了一项非常小的模拟研究（在R中）。运行模拟mu &...

R语言

编程开发

原创 2021-05-12 14:12:31 154 阅读

用R语言制作交互式图表和地图

原文链接：http://tecdat.cn/?p=8032可以直接从R / RStudio制作在线交互式图表和地图。配置启动RStudio，创建一个新的RScript，然后将工作目录设置为下载的数据文件夹。使ggplot2图表成为交互式Plotly图表制作交互式点线图以下代码将安装并加载程序包（该程序还将自动加载ggplot2），加载readr和dplyr，...

R语言

编程开发

原创 2021-05-12 14:12:05 206 阅读

用R语言模拟随机服务排队系统

原文链接：http://tecdat.cn/?p=8159M / M / c / k系统用肯德尔的表示法，M / M / c / k系统具有指数到达（M/ M / c / k），cC具有指数服务时间（M /M/ c / k）和k−c的服务器（M / M /c/ k）ķ-C队列中的位置（M / M / c /k）。这是M / M / 2/3系统的模拟。...

R语言

编程开发

原创 2021-05-12 14:12:04 220 阅读

用R语言模拟M / M / 1随机服务排队系统

原文链接：http://tecdat.cn/?p=8199本文中我在R中构造一个简单的M / M / 1队列的离散事件模拟。模拟变量像往常一样，我们从模拟及其检测所需的变量开始。 t.end <- 10^5 # duration of simt.clock <- 0 # sim timeTa <- 1.3333 # interarri...

R语言

编程开发

原创 2021-05-12 14:12:03 387 阅读

用R语言模拟混合制排队随机服务排队系统

原文链接：http://tecdat.cn/?p=8214M / M / 1系统该系统的基本参数::使用M / M / 1系统进行仿真非常简单。 lambda <- 2mu <- 4rho <- lambda/mu # = 2/4.. 例如，可以快速可视化随时间变化的资源使用情况。在下面，我们可以看到仿真如何收敛到系统中理论上的平...

R语言

编程开发

原创 2021-05-12 14:12:02 556 阅读

在R语言中进行缺失值填充：估算缺失值

原文链接：http://tecdat.cn/?p=8287介绍缺失值被认为是预测建模的首要障碍。因此，掌握克服这些问题的方法很重要。估算缺失值的方法的选择在很大程度上影响了模型的预测能力。在大多数统计分析方法中，按列表删除是用于估算缺失值的默认方法。但是，它不那么好，因为它会导致信息丢失。在本文中，我列出了5个R语言方法。链式方程进行的多元插补通过链...

R语言

编程开发

原创 2021-05-12 14:12:01 1008 阅读

主成分分析(PCA)原理及R语言实现及分析实例

原文链接：http://tecdat.cn/?p=8491主成分分析（PCA）是一种数据降维技巧，它能将大量相关变量转化为一组很少的不相关变量，这些无关变量称为主成分。例如，使用PCA可将30个相关（很可能冗余）的环境变量转化为5个无关的成分变量，并且尽可能地保留原始数据集的信息。主成分分析模型，变量（X1到X5）映射为主成分（PC1,PC2）　　PCA分析的一般步骤...

R语言

编程开发

原创 2021-05-12 14:11:58 684 阅读

R语言OutliersO3软件包异常值检测方法比较

原文链接：http://tecdat.cn/?p=8502识别异常值的方法有很多种，R中有很多不同的方法。关于异常值方法的文章使用了理论和实践的混合体。理论一切都很好，但异常值是异常值，因为它们不遵循理论。实践涉及数据的测试方法，有时用基于理论模拟的数据，更好地使用“真实”数据集。如果一种方法发现我们都同意的异常值，那么它可以被认为是成功的，但是我们是否都同意哪些个案是异常值？...

R语言

编程开发

原创 2021-05-12 14:11:57 199 阅读

R语言NLP案例：LDA主题文本挖掘优惠券推荐网站数据

原文链接：http://blog.sina.com.cn/s/blog_14154cb430102x1wb.html介绍Groupon是一个优惠券推荐服务，您可以免费注册Groupon，并且Groupon每天都会向您发送包含该地区当天交易的电子邮件。如果您喜欢这笔交易，那么您可以立即从Groupon购买，并在餐馆/商店兑换。数据这些数据是从Groupon网站的纽约市区域获得的...

R语言

编程开发

原创 2021-05-12 14:11:56 122 阅读

【大数据部落】R语言犯罪率回归模型报告Regression model on crimerate report

原文链接：http://blog.sina.com.cn/s/blog_14154cb430102wue4.htmlObjection：We attempts to explore the relationship between different demographic factors to crime rate, find out the important factor...

R语言

编程开发

原创 2021-05-12 14:11:55 343 阅读

官方博客	全部文章	热门标签	班级博客
了解我们	网站地图	意见反馈

鸿蒙开发者社区	51CTO学堂
51CTO	软考资讯

51CTO博客

拓端小助手