-
6 votesanswersviews
使用dplyr计算95％-CI的长度
上次我询问如何计算每个测量时间（周）的平均分数，对于多个受访者重复测量的变量（procras） . 所以我的（简化）长格式数据集看起来像下面的例子（这里有两个学生，5个时间点，没有分组变量）： studentID week procras 1 0 1.4 1 6 1.2 1 16 1.6 1 ... -
0 votesanswersviews
使用MATLAB进行离散分布均值的置信区间
我有一个负二项分布，需要解决其平均值的95％置信区间 . 我正在使用MATLAB，这是一个 probDist 对象 . 因为这是不正常的我想我必须使用bootstrapping但我不知道如何实现它 . 我尝试了 bootci 函数，但它需要样本作为输入，而不是分发 . 编辑：这是我的代码（它不起作用）： r = 2; p = 0.5; dist_x = makedist('NegativeBino... -
1 votesanswersviews
根据使用自举方法获得的经验分布计算置信区间
我使用bootstrap方法计算了样本均值的经验分布，但现在我还需要使用我发现的经验分布计算总体均值的置信区间 . 有没有办法在我的状态下在Matlab中自动完成？如果没有，您如何找到人口平均值的95％置信区间？ -
2 votesanswersviews
Python中概率分布的置信区间
我想计算分布的置信区间（Pearson 3），我正在使用scipy.stats . 使用stats.pearson3.interval我只有平均值的置信区间，但我需要所有数据的置信区间 . 所以，我需要一个像这样的图像： 我在matlab中知道这很简单，使用cdf函数，但我没有在Python中找到这样的东西 . -
0 votesanswersviews
仅从均值和标准差绘制置信区间
我试图从数据的均值和标准差（std）绘制置信区间 . 这是我写的代码片段： MeanA=1.876; %mean of A STDA=0.018; % std of A MeanB=1.821; STDB=0.039; MeanC=1.735; STDC=0.023; MeanD=1.667; STDD=0.039; Y = [MeanA MeanB ; MeanC Mea... -
2 votesanswersviews
假设未知分布，从样本数据计算置信区间
我有样本数据，我想计算一个置信区间，假设分布不正常且未知 . 基本上，看起来分布是Pareto但我不确定 . 正态分布的答案： Compute a confidence interval from sample data Correct way to obtain confidence interval with scipy -
2 votesanswersviews
Bootstrap置信区间
我应该计算父亲和儿子之间平均高度差异的单侧90％置信度bootstrap关键置信区间（使用UsingR包中的father.son数据集） . 当我在网上搜索并实现时，我找到了smean.cl.boot函数： smean.cl.boot(father.son$sheight-father.son$fheight, conf=0.90 但我不确定这是片面的还是关键的（因为我对关键的置信区间不熟悉） .... -
2 votesanswersviews
Bonferroni同时置信区间差异的间隔
我试图在R中获得Bonferroni同时置信区间 . 我有以下数据集，我为实践做了补充： df2 <- read.table(textConnection( 'group value 1 25 2 36 3 42 4 50 1 27 2 35 3 49 4 57 1 22 2 37 3 45 4 51'), header = TRUE) ... -
0 votesanswersviews
ggplot2图中的置信区间和使用R中的预测函数获得的值不相同
这是我的数据（例子）： EXAMPLE<-data.frame( X=c(99.6, 98.02, 96.43, 94.44, 92.06, 90.08, 87.3, 84.92, 82.14, 79.76, 76.98, 74.21, 71.03, 67.86, 65.08, 62.3, 59.92, 56.35, 52.38, 45.63, 41.67, 35.71, 30.... -
1 votesanswersviews
lme4中glmer线周围的置信区间
我花了很多时间试图弄清楚如何生成一个理想的情节，并想知道是否有人可以提供帮助 . 该图用于说明二元响应变量上“时间”和“组”之间的相互作用，“组2”的“时间”和“组”之间的相互作用比“组1”更快 . 以下是使用模拟数据所需的简化版本 . 我已经正确地为每个小组绘制了一个斜率（我认为），但是现在想要在每个斜坡周围绘制一个置信度 . #simulate 41 individuals, each occ... -
0 votesanswersviews
R中线性回归模型的自变量置信区间
假设一个数据集包含单个时间步长（以秒为单位）和一个距离（以m为单位）每个时间步长管理的跑步者： time <- c(300,600,900) # time in sec dist <- c(1050,1950,3100) # dist in m 这里有一个情节来说明它： plot(x=time,y=dist,xlim=c(0,1800),ylim=c(0,5000)) 基于零截距（... -
0 votesanswersviews
手动计算多元线性回归（OLS）的置信区间
我试图了解如何手动计算多元线性回归（OLS）的置信区间 . 我的问题是我不知道如何计算所有单个系数的标准误差 . 对于只有一个自变量的回归，我遵循以下教程：http://stattrek.com/regression/slope-confidence-interval.aspx . 本教程提供以下公式： 事实证明，该公式有效 . 但是，我并没有完全理解这个公式 . 例如，为什么（-2）位于公式的... -
1 votesanswersviews
coxph中产品的置信区间
我使用coxph作为生存模型，使用三个变量变量作为独立变量 . 我知道他们之间没有互动，所以我只做了以下模型： fit <- coxph(Surv(Time, status) ~ factor(v1) + factor(v2) + factor(v3) 我的问题是我需要产品之间的置信区间 . 我可以手工估计使用该公式的两个产品中的一个 Var(v1*v2) = v1^2 Var(v2) +... -
5 votesanswersviews
多元线性回归：用置信区间绘制直线
这是我的问题： 1) 我运行了多元线性回归：假设如下： lm(attitude~quality+price+location+Income) 我主要关心态度和质量之间的关系，其他变量是控制变量 . 2) 然后我想在态度和质量之间做一个散点图 . 这很容易： Q <-ggplot(data=data, aes(x=quality, y=attitude)) Q + geom_point(s... -
1 votesanswersviews
matlab中95％置信区间的求和高斯分布
我在Matlab中的曲线拟合工具箱中使用由2个求和高斯组成的拟合拟合了一些数据 . Matlab能够非常好地为我绘制95％的上下置信区间 . 现在，我有一些新数据最适合由3个高斯组成的曲线，但在这种情况下工具箱不会绘制95％CI . 它沿曲线的所有点返回NaN，而不是95％CI的数值 . 以下是一些代码： X = [0 1 2 3 4 5 6 7 8 9 10]'; Y = [0.009 0.0... -
9 votesanswersviews
使用MATLAB进行简单的二元逻辑回归
我正在使用MATLAB进行逻辑回归，以解决一个简单的分类问题 . 我的协变量是一个介于0和1之间的连续变量，而我的分类响应是0（不正确）或1（正确）的二进制变量 . 我正在寻找运行逻辑回归来 Build 预测器，该预测器将输出某些输入观察的概率（例如，如上所述的连续变量）是正确的或不正确的 . 虽然这是一个相当简单的场景，但我在MATLAB中运行它时遇到了一些麻烦 . 我的方法如下：我有一个列向量... -
0 votesanswersviews
使用前缀在Stata中存储二项式置信区间的结果
我试图计算95％二项式威尔逊置信区间为每年完成治疗的人数（数据集列出每个人的数据集） . 我想将结果存储到矩阵中，以便我可以使用 putexcel 命令将结果导出到现有的Excel电子表格而不更改工作表的格式 . 我创建了一个二进制变量 dscomplete_binary ，如果治疗没有完成，则为0，如果治疗完成则为1 . 我尝试过以下方法： bysort year: ci dscomplete_... -
2 votesanswersviews
置信区间
我的样本为160.平均值为262，标准差为165.正态分布曲线很好地从大约-200（平均-3 *偏差）到大约800（平均3 *偏差） . 在下一步我看95％的置信区间 . 对于这个例子，我在excel中使用了函数CONFIDENCE，它说这个inteval是[236,287] . 我的问题是 - 为什么这个范围这么小？每本书陈述并显示95％置信区间为大约/ - 2 * st.dev . 在我的... -
0 votesanswersviews
R泊松置信区间
我使用泊松回归（ glm() 与 family='poisson' ）并希望以下列形式给出答案： “每个预测值的区间[低，高]的误差为0.95” . 有可能吗？如果有，怎么样？ 我看到的方法很少 . 1）计算残差并获得低值和高值 . 我应该使用哪种发行版？ 2）使用功能 predict.glm(model, type = 'link', se.fit=TRUE) . 计算正态分布的置... -
1 votesanswersviews
Azure可以计算回归的置信区间吗？
我计划尝试使用Azure ML Studio提供的不同回归方法来预测数值 . 我想知道是否有可能将预测与相应的置信区间一起得到 . 换句话说，我希望回归函数不仅可以告诉我期望值（预测），还可以告诉我它（模型）对这个值的信心 . Azure回归是否支持此功能？ ADDED 一个相关的问题 . 可以 Build “回归量”估计概率密度函数吗？例如，对于给定的情况（数据表中的一行），我希望不仅有一个数... -
0 votesanswersviews
如何使用R中使用的“facet_wrap”绘制预测概率的置信区间
我使用以下包： library(aod) library(MASS) library(ggplot2) 我正在关注以下链接中的示例R代码：http://stats.idre.ucla.edu/r/dae/logit-regression/ 这是我的GLMM的代码 str(data1) flocation <- factor(data1$location) fID <- factor(... -
0 votesanswersviews
如何对R回归模型中的Bootstrap预测和置信水平
我在R中使用一个定量预测变量（ amplitude = 40,50,60,70）在比例数据上运行β回归模型 . 我已经能够从模型中获得40,50,60,70比例的预测并绘制它 . 我从围绕这个网站和其他人的研究中了解到，从β回归模型获得预测的置信区间与其他模型不同 . 我已经读过这个以获得β回归模型的置信水平，从模型中引导预测是获得置信区间/带的一种方式（如其他帖子中所述，例如Achim Zei... -
1 votesanswersviews
使用Bebber等人2007查找估计的总物种和该数字的置信区间[已结束]
我的问题是统计数据和R相关的混合 . 我正在尝试实施Bebber等人的模型 . 2007（使用发现曲线预测未知物种数量）估算我的数据集中未发现物种的总数（通过找到预期的物种总数），以及围绕物种总数预测的置信区间 . 物种总数的估计由E（St）= k（Ntot-Nt-1）给出，其通过拟合Nt-1上的St的glm（具有拟线性误差）来计算，其中St是物种数每年，Nt-1是前一年的累积物种数，k是发现率... -
1 votesanswersviews
绘制分段回归对象时出错
我正在使用程序包'segmented'对程序R中的某些二元逻辑回归进行分段回归 . 我选择了最终模型，现在想要绘制预测的概率和置信带 . 我的最终模型有一个连续变量，我正在分析回归线以及两个二分因子变量 . 这是一些输出，因此您可以可视化模型的结构 Estimated Break-Point(s): Est. St.Err 3.799 1.117 t value for the g... -
0 votesanswersviews
使用Bootstrapping在GLMM中创建回归曲线的置信区间
创建了GLMM模型并绘制了每个因子的预测概率 . 但是，我无法理解如何使用BootMer函数创建置信区间 . 我不断收到错误消息无法模拟非整数先前权重 . 我希望有人能帮忙吗？提前致谢 . glmm1 <- glmer(cbind(Total_Dead, Total_Collected - Total_Dead) ~ Species + timeseri... -
2 votesanswersviews
使用predict.glm返回什么标准错误（...，type =“response”，se.fit = TRUE）？
在执行逻辑回归后，我将根据excellent example中提供的关于如何计算响应的95％置信区间的数据拟合模型： foo <- mtcars[,c("mpg","vs")]; names(foo) <- c("x","y") mod <- glm(y ~ x, data = foo, family... -
0 votesanswersviews
计算和显示均值的置信区间
我现在有一个图表，显示每组距离的平均值，但是我很难在图中添加置信区间并显示它们 . 我的代码是： TeamMeans<- read.csv(file = file.choose()) teamM=TeamMeans graph6<-ggplot(aes(x=Team, y=Mean), data=teamM) + geom_point() + geom_smooth(size ... -
-1 votesanswersviews
计算1维随机数据的置信区间
我想知道如何计算1维的给定置信区间（这可以简化问题）随机数据 . 情况是这样的：假设我有50个随机数据点，从1到100（不是正态分布），我想知道我们对该范围内90％数据点的最小范围 . 这与正态分布的置信区间非常相似 . 有没有人可以帮我这个？谢谢 . -
1 votesanswersviews
在最终估算中使用“Inf”的Bootstrap置信区间 - boot / dplyr包
我试图从拟合模型中引导一个值（曲线的最大斜率） . 虽然我认为代码可能更有效，但我认为它通常有效（关于如何提高效率的建议也会受到欢迎） . 我遇到的问题是，当我得到最终样本值时，其中一个单元格中有 Inf 会阻止我获得置信区间 . 我不确定这是一个引导问题还是我编写它的方式很简单 . 示例代码： library(boot) library(dplyr) df <- data.frame ... -
1 votesanswersviews
用dplyr引导置信区间
看看答案here如何估算自举时间间隔？这个问题也在ggplot2列表中被提出 . library(dplyr) mtcars %>% group_by(vs) %>% summarise(mean.mpg = mean(mpg, na.rm = TRUE), sd.mpg = sd(mpg, na.rm = TRUE), n.mpg = n()) %>% muta...