如何读取R 的sumary()结果
思路
step 1: sum = summary(model)
step 2: sum有好多属性,直接根据属性名称引用(\()即可, 如:
+ > sum\)call 返回 model 使用的模型语句
+ > sum\(coefficients; 返回一个列表,可以继续引用,如下
+ > sum\)coefficients[1 : 12]; 返回一个列表的一个切片,还可以继续切片
+ > sum$coefficients[1 : 12][2]; 返回一个列表的一个切片的第二个元素
下面是一些 测试代码,未整理,可以大致学习一下
用'demo()'来看一些示范程序,用'help()'来阅读在线帮助文件,或
用'help.start()'通过HTML浏览器来看帮助文件。
用'q()'退出R.
[Workspace loaded from ~/.RData]
> y=c(53,434,111,38,108,48)
> x1=c(1,2,3,1,2,3)
> x2=c(1,2,1,2,1,2)
> log.glm <-glm(y~x1+x2,family = possion(link=log))
Error in possion(link = log) : 没有"possion"这个函数
> log.glm <-glm(y~x1+x2,family = possion(link=log),data=(y,x1,x2))
错误: 意外的',' in "log.glm <-glm(y~x1+x2,family = possion(link=log),data=(y,"
> dataframe <-data.frame(y,x1,x2)
> head(dataframe)
y x1 x2
1 53 1 1
2 434 2 2
3 111 3 1
4 38 1 2
5 108 2 1
6 48 3 2
> log.glm <-glm(y~x1+x2,family = possion(link=log),data=data.frame(y,x1,x2))
Error in possion(link = log) : 没有"possion"这个函数
> log.glm <-glm(y~x1+x2,family = poisson(link=log),data=data.frame(y,x1,x2))
> summary(log.glm)
Call:
glm(formula = y ~ x1 + x2, family = poisson(link = log), data = data.frame(y,
x1, x2))
Deviance Residuals:
1 2 3 4 5 6
-3.1382 16.6806 0.8189 -11.0398 1.8210 -12.6942
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.59532 0.15792 22.767 < 2e-16 ***
x1 0.12915 0.04370 2.955 0.00312 **
x2 0.64803 0.07483 8.660 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 662.84 on 5 degrees of freedom
Residual deviance: 575.10 on 3 degrees of freedom
AIC: 619.08
Number of Fisher Scoring iterations: 5
> log.glm.x1
错误: 找不到对象'log.glm.x1'
>
>
>
>
>
>
> summary(log.glm)
Call:
glm(formula = y ~ x1 + x2, family = poisson(link = log), data = data.frame(y,
x1, x2))
Deviance Residuals:
1 2 3 4 5 6
-3.1382 16.6806 0.8189 -11.0398 1.8210 -12.6942
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.59532 0.15792 22.767 < 2e-16 ***
x1 0.12915 0.04370 2.955 0.00312 **
x2 0.64803 0.07483 8.660 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 662.84 on 5 degrees of freedom
Residual deviance: 575.10 on 3 degrees of freedom
AIC: 619.08
Number of Fisher Scoring iterations: 5
> log.glm.x1
错误: 找不到对象'log.glm.x1'
> help("glm")
> anova(log.glm)
Analysis of Deviance Table
Model: poisson, link: log
Response: y
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev
NULL 5 662.84
x1 1 8.770 4 654.07
x2 1 78.978 3 575.10
> ano= anova(log.glm)
> ano[1]
Df
NULL
x1 1
x2 1
> ano[2]
Deviance
NULL
x1 8.770
x2 78.978
> ano[3]
Resid. Df
NULL 5
x1 4
x2 3
> sum= summary(log.glm)
> sum
Call:
glm(formula = y ~ x1 + x2, family = poisson(link = log), data = data.frame(y,
x1, x2))
Deviance Residuals:
1 2 3 4 5 6
-3.1382 16.6806 0.8189 -11.0398 1.8210 -12.6942
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.59532 0.15792 22.767 < 2e-16 ***
x1 0.12915 0.04370 2.955 0.00312 **
x2 0.64803 0.07483 8.660 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 662.84 on 5 degrees of freedom
Residual deviance: 575.10 on 3 degrees of freedom
AIC: 619.08
Number of Fisher Scoring iterations: 5
> sum[1]
$call
glm(formula = y ~ x1 + x2, family = poisson(link = log), data = data.frame(y,
x1, x2))
> sum[1,1]
Error in sum[1, 1] : 量度数目不对
> sum[2]
$terms
y ~ x1 + x2
attr(,"variables")
list(y, x1, x2)
attr(,"factors")
x1 x2
y 0 0
x1 1 0
x2 0 1
attr(,"term.labels")
[1] "x1" "x2"
attr(,"order")
[1] 1 1
attr(,"intercept")
[1] 1
attr(,"response")
[1] 1
attr(,".Environment")
<environment: R_GlobalEnv>
attr(,"predvars")
list(y, x1, x2)
attr(,"dataClasses")
y x1 x2
"numeric" "numeric" "numeric"
> sum[3]
$family
Family: poisson
Link function: log
> sum[4]
$deviance
[1] 575.0954
> sum[4,1]
Error in sum[4, 1] : 量度数目不对
> sum
Call:
glm(formula = y ~ x1 + x2, family = poisson(link = log), data = data.frame(y,
x1, x2))
Deviance Residuals:
1 2 3 4 5 6
-3.1382 16.6806 0.8189 -11.0398 1.8210 -12.6942
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.59532 0.15792 22.767 < 2e-16 ***
x1 0.12915 0.04370 2.955 0.00312 **
x2 0.64803 0.07483 8.660 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 662.84 on 5 degrees of freedom
Residual deviance: 575.10 on 3 degrees of freedom
AIC: 619.08
Number of Fisher Scoring iterations: 5
> sum[5]
$aic
[1] 619.0808
> sum[6]
$contrasts
NULL
> sum[8]
$null.deviance
[1] 662.8432
> sum[9]
$df.null
[1] 5
> sum[10]
$iter
[1] 5
> sum$aic
[1] 619.0808
> sum$null.deviance
[1] 662.8432
> sum$residual.deviance
NULL
> sum$residual.devianc
NULL
> sum[11]
$deviance.resid
1 2 3 4 5 6
-3.1382350 16.6805594 0.8189003 -11.0397892 1.8209720 -12.6941833
> summary.aov()
Error in summary.aov() : 缺少参数"object",也没有缺省值
> sum
Call:
glm(formula = y ~ x1 + x2, family = poisson(link = log), data = data.frame(y,
x1, x2))
Deviance Residuals:
1 2 3 4 5 6
-3.1382 16.6806 0.8189 -11.0398 1.8210 -12.6942
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.59532 0.15792 22.767 < 2e-16 ***
x1 0.12915 0.04370 2.955 0.00312 **
x2 0.64803 0.07483 8.660 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 662.84 on 5 degrees of freedom
Residual deviance: 575.10 on 3 degrees of freedom
AIC: 619.08
Number of Fisher Scoring iterations: 5
> sum$coefficients
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.5953201 0.15791713 22.767132 9.709542e-115
x1 0.1291456 0.04370053 2.955240 3.124256e-03
x2 0.6480267 0.07482977 8.660013 4.717107e-18
> sum$coefficients[4]
[1] 0.1579171
> sum$coefficients[5]
[1] 0.04370053
> sum$coefficients[6]
[1] 0.07482977
> sum$coefficients[1]
[1] 3.59532
> sum$coefficients[1..2]
错误: unexpected numeric constant in "sum$coefficients[1..2"
> sum$coefficients[1:2]
[1] 3.5953201 0.1291456
> sum$coefficients[1:5]
[1] 3.59532005 0.12914558 0.64802675 0.15791713 0.04370053
> sum$coefficients[11:12]
[1] 3.124256e-03 4.717107e-18
> sum$coefficients[11:12][1]
[1] 0.003124256
> sum$coefficients[11:12][2]
[1] 4.717107e-18
如何读取R 的sumary()结果的更多相关文章
- R语言基础
一.扩展包的基本操作语句R安装好之后,默认自带了"stats" "graphics" "grDevices" "utils&qu ...
- R语言︱文件读入、读出一些方法罗列(批量xlsx文件、数据库、文本txt、文件夹)
笔者寄语:小规模的读取数据的方法较为简单并且多样,但是,批量读取目前看到有以下几种方法:xlsx包.RODBC包.批量转化成csv后读入. R语言中还有一些其他较为普遍的读入,比如代码包,R文件,工作 ...
- R语言手册
在R的官方教程里是这么给R下注解的:一个数据分析和图形显示的程序设计环境(A system for data analysis and visualization which is built bas ...
- 让R与Python共舞
转载:http://ices01.sinaapp.com/?p=129 R(又称R语言)是一款开源的跨平台的数值统计和数值图形化展现 工具.通俗点说,R是用来做统计和画图的.R拥有自己的脚本 ...
- Android小例子:使用反射机制来读取图片制作一个图片浏览器
效果图: 工程文件夹: 该例子可供于新手参考练习,如果有哪里不对的地方,望指正>-< <黑幕下的人> java代码(MainActivity.java): package co ...
- ArcGIS二次开发之读取遥感图像像素值的做法
作者:朱金灿 来源:http://blog.csdn.net/clever101 首先是读取遥感图像的R.G.B波段数据的做法.读取R.G.B波段数据的像素值主要通过IRaster接口的Read方法在 ...
- go语言学习笔记---读取文件io/ioutil 包
io/ioutil 包几个函数方法 名称 作用 备注 ReadAll 读取数据,返回读到的字节 slice 1 ReadDir 读取一个目录,返回目录入口数组 []os.FileInfo, 2 Re ...
- GO 文件读取常用的方法
方式1: 一行一行的方式读取 其中常用的方法就有:ReadString,ReadLine,ReadBytes ReadLine 返回单个行,不包括行尾字节,就是说,返回的内容不包括\n或者\r\n,返 ...
- [高性能MYSQL 读后随笔] 关于事务的隔离级别(一)
一.锁的种类 MySQL中锁的种类很多,有常见的表锁和行锁,也有新加入的Metadata Lock等等,表锁是对一整张表加锁,虽然可分为读锁和写锁,但毕竟是锁住整张表,会导致并发能力下降,一般是做dd ...
随机推荐
- 使用httpClient发送post请求连接restful接口
public static String httpPost(String url,String arg){ InputStream is; BufferedReader br; StringBuild ...
- Failed to start component [StandardEngine[Tomcat].StandardHost[localhost].StandardContext[]]
今天在测试项目代码时,在idea中配置tomcat7插件运行后一直报如下错误: 解决方案:看了网上大多数办法都是修改xml文件配置,感觉并不适用,最后看到比较靠谱解释如下: pom.xml中jar包发 ...
- Jmeter正则表达式提取器(转载)
转载自 http://blog.csdn.net/qq_35885203 使用jmeter来测试时,经常会碰到需要上下文传输数据的情况,如登录后生成的token,在其他页面的操作,都需传入这个toke ...
- idea搭spring boot项目
Spring Boot是由Pivotal团队提供的全新框架,设计目的是用来简化新Spring应用的初始搭建以及开发过程.它主要推崇的是'消灭配置’,实现零配置. 那么,如何在idea中创建一个spri ...
- __getattr__ __delattr__ __setattr__ __getattribute__使用(重写python提供的错误信息)
自己定义了这些attr 查找删除设置就会触发自己定义的逻辑,如果不重新,pyton会提供自己报错信息class Room: def __init__(self,name): self.name = n ...
- 模拟器 Unable to execute simctl install Error 117
ios 模拟器 Unable to execute '"/usr/bin/xcrun" simctl install "/Users/tt/PAServer/scrat ...
- 关于Bootstrap的入门知识
问:Bootstrap是什么? 答:开源的前端框架,就是一些事先写好的css.js等. 问:Bootstrap在哪儿下载? 答:官方(https://getbootstrap.com/),中文(htt ...
- 一次sendmsg的改造过程
比较蛋疼的一个改造过程,简单记录一下. 场景:用户态使用sendmsg发包,tcp报文,由于内核实现过程中存在一次kernel_read,也就是存在将pagecache中的内容拷贝一次的问题. 为了减 ...
- MySQL统计信息以及执行计划预估方式初探
数据库中的统计信息在不同(精确)程度上描述了表中数据的分布情况,执行计划通过统计信息获取符合查询条件的数据大小(行数),来指导执行计划的生成.在以Oracle和SQLServer为代表的商业数据库,和 ...
- hello1以及hello2的部分代码分析
(一)1.GreetingServlet.java源码文件: @WebServlet("/greeting") //以@WebServlet注释开头,注释指定相对于上下文根的URL ...