| Line 1... |
Line 1... |
| 1 |
|
1 |
|
| 2 |
R : Copyright 2000, The R Development Core Team
|
2 |
R : Copyright 2001, The R Development Core Team
|
| 3 |
Version 1.2.0 Under development (unstable) (2000-10-09)
|
3 |
Version 1.3.0 Under development (unstable) (2001-03-06)
|
| 4 |
|
4 |
|
| 5 |
R is free software and comes with ABSOLUTELY NO WARRANTY.
|
5 |
R is free software and comes with ABSOLUTELY NO WARRANTY.
|
| 6 |
You are welcome to redistribute it under certain conditions.
|
6 |
You are welcome to redistribute it under certain conditions.
|
| 7 |
Type "?license" or "?licence" for distribution details.
|
7 |
Type `license()' or `licence()' for distribution details.
|
| 8 |
|
8 |
|
| 9 |
R is a collaborative project with many contributors.
|
9 |
R is a collaborative project with many contributors.
|
| 10 |
Type "?contributors" for a list.
|
10 |
Type `contributors()' for more information.
|
| 11 |
|
11 |
|
| 12 |
Type "demo()" for some demos, "help()" for on-line help, or
|
12 |
Type `demo()' for some demos, `help()' for on-line help, or
|
| 13 |
"help.start()" for a HTML browser interface to help.
|
13 |
`help.start()' for a HTML browser interface to help.
|
| 14 |
Type "q()" to quit R.
|
14 |
Type `q()' to quit R.
|
| 15 |
|
15 |
|
| 16 |
> options(digits=7)
|
16 |
> options(digits=7)
|
| 17 |
>
|
17 |
>
|
| 18 |
> ## powers, including complex ones
|
18 |
> ## powers, including complex ones
|
| 19 |
> outer(a <- -4:12, -2:7, "^")
|
19 |
> outer(a <- -4:12, -2:7, "^")
|
| Line 132... |
Line 132... |
| 132 |
|
132 |
|
| 133 |
n= 29 : 435+0i -14.5+133.3253i -14.5+65.87416i -14.5+43.03447i -14.5+31.34124i -14.5+24.09919i -14.5+19.07442i -14.5+15.30746i -14.5+12.31641i -14.5+9.831244i -14.5+7.687413i -14.5+5.777328i -14.5+4.025905i -14.5+2.377154i -14.5+0.7861672i -14.5-0.7861672i -14.5-2.377154i -14.5-4.025905i -14.5-5.777328i -14.5-7.687413i -14.5-9.831244i -14.5-12.31641i -14.5-15.30746i -14.5-19.07442i -14.5-24.09919i -14.5-31.34124i -14.5-43.03447i -14.5-65.87416i -14.5-133.3253i
|
133 |
n= 29 : 435+0i -14.5+133.3253i -14.5+65.87416i -14.5+43.03447i -14.5+31.34124i -14.5+24.09919i -14.5+19.07442i -14.5+15.30746i -14.5+12.31641i -14.5+9.831244i -14.5+7.687413i -14.5+5.777328i -14.5+4.025905i -14.5+2.377154i -14.5+0.7861672i -14.5-0.7861672i -14.5-2.377154i -14.5-4.025905i -14.5-5.777328i -14.5-7.687413i -14.5-9.831244i -14.5-12.31641i -14.5-15.30746i -14.5-19.07442i -14.5-24.09919i -14.5-31.34124i -14.5-43.03447i -14.5-65.87416i -14.5-133.3253i
|
| 134 |
|
134 |
|
| 135 |
n= 30 : 465+0i -15+142.7155i -15+70.56945i -15+46.16525i -15+33.69055i -15+25.98076i -15+20.64573i -15+16.65919i -15+13.50606i -15+10.89814i -15+8.660254i -15+6.67843i -15+4.873795i -15+3.188348i -15+1.576564i -15+0i -15-1.576564i -15-3.188348i -15-4.873795i -15-6.67843i -15-8.660254i -15-10.89814i -15-13.50606i -15-16.65919i -15-20.64573i -15-25.98076i -15-33.69055i -15-46.16525i -15-70.56945i -15-142.7155i
|
135 |
n= 30 : 465+0i -15+142.7155i -15+70.56945i -15+46.16525i -15+33.69055i -15+25.98076i -15+20.64573i -15+16.65919i -15+13.50606i -15+10.89814i -15+8.660254i -15+6.67843i -15+4.873795i -15+3.188348i -15+1.576564i -15+0i -15-1.576564i -15-3.188348i -15-4.873795i -15-6.67843i -15-8.660254i -15-10.89814i -15-13.50606i -15-16.65919i -15-20.64573i -15-25.98076i -15-33.69055i -15-46.16525i -15-70.56945i -15-142.7155i
|
| 136 |
>
|
136 |
>
|
| - |
|
137 |
> ## lowess() {incl. sort, etc}:
|
| - |
|
138 |
> options(digits = 5)
|
| - |
|
139 |
>
|
| - |
|
140 |
> lowess(c(3,2,6,3,8,4))$y # this used to differ on Linux
|
| - |
|
141 |
[1] 2.4259 3.4293 3.9896 5.2832 5.3386 4.9822
|
| - |
|
142 |
>
|
| - |
|
143 |
> y1 <- c(3,1:2,5:2,4,1:3,3)
|
| - |
|
144 |
> lowess(y1)$y
|
| - |
|
145 |
[1] 2.0723 2.4362 2.7508 3.0038 3.2036 3.2435 2.8485 2.4832 2.4197 2.5231
|
| - |
|
146 |
[11] 2.6286 2.7813
|
| - |
|
147 |
> lowess(y1, f = .4)$y
|
| - |
|
148 |
[1] 2.7030 1.9593 2.4351 3.4933 4.0000 3.0000 2.4651 1.8579 1.7357 2.0000
|
| - |
|
149 |
[11] 2.7103 3.1184
|
| - |
|
150 |
>
|
| - |
|
151 |
> lowess(c(y1,100), f = .4)$y
|
| - |
|
152 |
[1] 2.2135 2.2468 2.5459 3.7984 4.0000 3.0000 2.8427 2.4999 2.1512 2.0000
|
| - |
|
153 |
[11] 2.7129 3.1195 3.4172
|
| - |
|
154 |
>
|
| - |
|
155 |
> ## this is the test sample from Cleveland's original lowess.doc:
|
| - |
|
156 |
> x <- c(1:5, rep(6,10),8,10,12,14,50)
|
| - |
|
157 |
> y <- c(18,2,15,6,10,4,16,11,7,3,14,17,20,12,9,13,1,8,5,19)
|
| - |
|
158 |
> lowess(x,y, f = .25, iter = 0, delta = 0)$y
|
| - |
|
159 |
[1] 13.6588 11.1446 8.7012 9.7220 10.0000 11.3000 11.3000 11.3000 11.3000
|
| - |
|
160 |
[10] 11.3000 11.3000 11.3000 11.3000 11.3000 11.3000 13.0000 6.4399 5.5959
|
| - |
|
161 |
[19] 5.4557 18.9982
|
| - |
|
162 |
> lowess(x,y, f = .25, iter = 0, delta = 3)$y
|
| - |
|
163 |
[1] 13.6588 12.3466 11.0343 9.7220 10.5110 11.3000 11.3000 11.3000 11.3000
|
| - |
|
164 |
[10] 11.3000 11.3000 11.3000 11.3000 11.3000 11.3000 13.0000 6.4399 5.5959
|
| - |
|
165 |
[19] 5.4557 18.9982
|
| - |
|
166 |
> lowess(x,y, f = .25, iter = 2, delta = 0)$y
|
| - |
|
167 |
[1] 14.8112 12.1154 8.9838 9.6757 10.0000 11.3460 11.3460 11.3460 11.3460
|
| - |
|
168 |
[10] 11.3460 11.3460 11.3460 11.3460 11.3460 11.3460 13.0000 6.7345 5.7437
|
| - |
|
169 |
[19] 5.4147 18.9981
|
| - |
|
170 |
>
|
| - |
|
171 |
>
|
The R Project SVN
R