R软件操作实例

R version 3.3.2 (2016-10-31) -- "Sincere Pumpkin Patch"

Copyright (C) 2016 The R Foundation for Statistical Computing

Platform: i386-w64-mingw32/i386 (32-bit)

R是自由软件，不带任何担保。

在某些条件下你可以将其自由散布。

用'license()'或'licence()'来看散布的详细条件。

R是个合作计划，有许多人为之做出了贡献.

用'contributors()'来看合作者的详细情况

用'citation()'会告诉你如何在出版物中正确地引用R或R程序包。

用'demo()'来看一些示范程序，用'help()'来阅读在线帮助文件，或

用'help.start()'通过HTML浏览器来看帮助文件。

用'q()'退出R.

[原来保存的工作空间已还原]

> getwd()

[1] "C:/Users/A/Documents"

> getwd()

[1] "D:/360极速浏览器下载"

> library(base)

> library(mgcv)

载入需要的程辑包：nlme

This is mgcv 1.8-15. For overview type 'help("mgcv-package")'.

>

> y <- c(rep(1,10),rep(0,20),rep(1,40),rep(0,10),rep(1,40),rep(0,40))

> y

[1] 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1

[37] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0

[73] 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

[109] 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

[145] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

> y <- c(rep(1,10),rep(0,20),rep(1,40),rep(0,10),rep(1,40),rep(0,40))

> x <- 1:160

> x

[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 [19] 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 [37] 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 [55] 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 [73] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 [91] 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 [109] 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126

[127] 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 [145] 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160

> b <- glm(y~x,family=binomial)

> b

Call: glm(formula = y ~ x, family = binomial)

Coefficients:

(Intercept) x

1.4548 -0.0146

Degrees of Freedom: 159 Total (i.e. Null); 158 Residual

Null Deviance: 219.3

Residual Deviance: 203 AIC: 207

> par(mfrow=c(2,2))

> qq.gam(b)

> qq.ga (mb,rep=50,level=1)

错误: 没有"qq.ga"这个函数

> qq.gam(b,rep=50,level=1

+ )

> libraray(meta)

错误: 没有"libraray"这个函数

> library(meta)

Loading 'meta' package (version 4.6-0).

Type 'help("meta-package")' for a brief overview.

> ?metabin

starting httpd help server ... done

> search()

[1] ".GlobalEnv" "package:meta" "package:mgcv"

[4] "package:nlme" "package:stats" "package:graphics"

[7] "package:grDevices" "package:utils" "package:datasets"

[10] "package:methods" "Autoloads" "package:base"

> qq.gam()

Error in qq.gam() : 缺少参数"object",也没有缺省值

> data(amlodipine)

> amlodipine

study n.amlo mean.amlo var.amlo n.plac mean.plac var.plac

1 Protocol 154 46 0.2316 0.2254 48 -0.0027 0.0007

2 Protocol 156 30 0.2811 0.1441 26 0.0270 0.1139

3 Protocol 157 75 0.189

4 0.1981 72 0.0443 0.4972

4 Protocol 162A 12 0.0930 0.1389 12 0.2277 0.0488

5 Protocol 163 32 0.1622 0.0961 34 0.005

6 0.0955

6 Protocol 166 31 0.183

7 0.1246 31 0.0943 0.1734

7 Protocol 303A 27 0.6612 0.7060 27 -0.0057 0.9891