r语言求置信区间的函数_r语言区间估计

r语言求置信区间的函数_r语言区间估计转载链接:https://www.codesd.com/item/confidence-interval-of-coefficients-using-the-generalized-estimation-equation-ege.html.链接:https://stackoverflow.com/questions/21221280/confidence-interval-of-coefficients-using-generalized-estimating-equation-gee.需要明确,

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

转载
链接: https://www.codesd.com/item/confidence-interval-of-coefficients-using-the-generalized-estimation-equation-ege.html.
链接: https://stackoverflow.com/questions/21221280/confidence-interval-of-coefficients-using-generalized-estimating-equation-gee.

需要明确,在置信度(检验水准为0.05)确定的情况下置信区间的计算需要Estimate(系数样本估计值)和Std.err(标准误)
即upr=Estimate+1.96Std.err
lwr=Estimate-1.96
Std.err

方法一
这个方法不能得到哑变量的系数,只能得到这个变量的系数;比如说Time这个变量是3个水平,两个哑变量,按照方法一只能得到Time这个变量的,无法得到Time2和TIme3参照于TIme1的。

library(geepack)
data(dietox)
dietox$Cu     <- as.factor(dietox$Cu)
mf1 <- formula(Weight~Cu*poly(Time,3))
gee1 <- geeglm(mf1, data=dietox, id=Pig,
               family=poisson("identity"),corstr="ar1")
cc <- coef(summary(gee1))
citab <- with(as.data.frame(cc),
     cbind(lwr=Estimate-1.96*Std.err,
           upr=Estimate+1.96*Std.err))
rownames(citab) <- rownames(cc)

Jetbrains全家桶1年46,售后保障稳定

方法二
自己定义confint函数,因为confint函数在geeglm函数中无法使用

confint.geeglm <- function(object, parm, level = 0.95, ...) {
    cc <- coef(summary(object))
    mult <- qnorm((1+level)/2)
    citab <- with(as.data.frame(cc),
                  cbind(lwr=Estimate-mult*Std.err,
                        upr=Estimate+mult*Std.err))
    rownames(citab) <- rownames(cc)
    citab[parm,]
}

confint(gee1)
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/234256.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • Stopwatch用法

    Stopwatch用法获取系统时间计算System.currentTimeMillis()Stopwatch对程序部分代码进行计时(ms级别),适用于同步单线程代码块。StopWatch实例一次只能开启一个task,不能同时start多个task,并且在该task未stop之前不能start一个新的task,必须在该taskstop之后才能开启新的task,若要一次开启多个,需要new不同的StopWatch实例//只输出运行多少秒Stopwatchstopwatch=Stopwatch.creat.

  • HTML的注释与JavaScript注释及快捷键「建议收藏」

    HTML的注释与JavaScript注释及快捷键「建议收藏」1、单行注释:在代码前输入两个反斜杠//,或用快捷键ctrl+/2、多行注释:使用/**/把要注释的代码块包起来,或用快捷键shift+alt+a3、若觉得2、的快捷键用着不方便,可以按自己喜欢可以在vscode中修改注释快捷键。…

  • jar运行发生ClassNotFound的完整解决方案

    jar运行发生ClassNotFound的完整解决方案原文地址:http://www.java2000.net/viewthread.jsp?tid=6053转载请注明上述链接或者CSDN的链接1今天彻底测试了jar程序TestJar.javapackagenet.java2000.test.jar;importjavax.swing.JOptionPane;importorg.springframework.beans.fa

  • java写入文件

    java写入文件

  • 播放.avi后缀视频报出0xc00d5212,编码格式不支持

    播放.avi后缀视频报出0xc00d5212,编码格式不支持以avi后缀的格式视频文件,在win10系统上播放可能会报如下如下错误:最普遍的现象就是高版本Windows媒体播放器播放不了采用早期编码编辑的AVI格式视频,而低版本Windows媒体播放器又播放不了采用最新编码编辑的AVI格式视频解决方案:这里我总结了两种方案:第一种:安装一个插件名字叫格式工厂,这款插件可以很好的支持大批量的文件格式转换,它会把avi视频转换成mp4格式视频,…

  • 使用Fastjson生成Json字符串少字段属性(数据丢失)「建议收藏」

    在开发中经常要使用到fastJson来转换对象为json 串,但是最近发现在一个对象转换的时候,总是丢失了一个字段的值,(数据丢失).就很纳闷。到网上找了很多方法总是不行,最后总算是在一篇博文中看到问题的关键!现在整理如下,希望可以帮助到更多的人。package per.eblink.pojo; public class Node { private String id

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号