Description Usage Arguments Details Value Author(s) Examples
Interfaces to MASS
functions that can be used
in a pipeline implemented by magrittr
.
1 2 3 4 5 6 7 8 9 10 | ntbt_corresp(data, ...)
ntbt_glm.nb(data, ...)
ntbt_lda(data, ...)
ntbt_lm.gls(data, ...)
ntbt_lm.ridge(data, ...)
ntbt_loglm(data, ...)
ntbt_logtrans(data, ...)
ntbt_polr(data, ...)
ntbt_qda(data, ...)
ntbt_rlm(data, ...)
|
data |
data frame, tibble, list, ... |
... |
Other arguments passed to the corresponding interfaced function. |
Interfaces call their corresponding interfaced function.
Object returned by interfaced function.
Roberto Bertolusso
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 | ## Not run:
library(intubate)
library(magrittr)
library(MASS)
## corresp
## Original function to interface
corresp(~ Age + Eth, data = quine)
## The interface reverses the order of data and formula
ntbt_corresp(data = quine, ~ Age + Eth)
## so it can be used easily in a pipeline.
quine %>%
ntbt_corresp(~ Age + Eth)
## glm.nb
## Original function to interface
glm.nb(Days ~ Sex/(Age + Eth*Lrn), data = quine)
## The interface reverses the order of data and formula
ntbt_glm.nb(data = quine, Days ~ Sex/(Age + Eth*Lrn))
## so it can be used easily in a pipeline.
quine %>%
ntbt_glm.nb(Days ~ Sex/(Age + Eth*Lrn))
## lda
Iris <- data.frame(rbind(iris3[,,1], iris3[,,2], iris3[,,3]),
Sp = rep(c("s","c","v"), rep(50,3)))
## Original function to interface
lda(Sp ~ ., Iris)
## The interface reverses the order of data and formula
ntbt_lda(Iris, Sp ~ .)
## so it can be used easily in a pipeline.
Iris %>%
ntbt_lda(Sp ~ .)
stackloss %>%
ntbt_lda(stack.loss ~ .) %>%
summary()
## lm.gls
## Original function to interface
lm.gls(conc ~ uptake, CO2, W = diag(nrow(CO2)))
## The interface reverses the order of data and formula
ntbt_lm.gls(CO2, conc ~ uptake, W = diag(nrow(CO2)))
## so it can be used easily in a pipeline.
CO2 %>%
ntbt_lm.gls(conc ~ uptake, W = diag(nrow(CO2)))
## lm.ridge
## Original function to interface
lm.ridge(GNP.deflator ~ ., longley)
## The interface reverses the order of data and formula
ntbt_lm.ridge(longley, GNP.deflator ~ .)
## so it can be used easily in a pipeline.
longley %>%
ntbt_lm.ridge(GNP.deflator ~ .)
## loglm
## Original function to interface
xtCars93 <- xtabs(~ Type + Origin, Cars93)
loglm(~ Type + Origin, xtCars93)
## The interface reverses the order of data and formula
xtCars93 <- ntbt_xtabs(Cars93, ~ Type + Origin)
ntbt_loglm(xtCars93, ~ Type + Origin)
## so it can be used easily in a pipeline.
Cars93 %>%
ntbt_xtabs(~ Type + Origin) %>%
ntbt_loglm(~ Type + Origin)
## logtrans
## Original function to interface
logtrans(Days ~ Age*Sex*Eth*Lrn, data = quine,
alpha = seq(0.75, 6.5, len=20))
## The interface reverses the order of data and formula
ntbt_logtrans(data = quine, Days ~ Age*Sex*Eth*Lrn,
alpha = seq(0.75, 6.5, len=20))
## so it can be used easily in a pipeline.
quine %>%
ntbt_logtrans(Days ~ Age*Sex*Eth*Lrn,
alpha = seq(0.75, 6.5, len=20))
## polr
op <- options(contrasts = c("contr.treatment", "contr.poly"))
## Original function to interface
polr(Sat ~ Infl + Type + Cont, housing)
## The interface reverses the order of data and formula
ntbt_polr(housing, Sat ~ Infl + Type + Cont)
## so it can be used easily in a pipeline.
housing %>%
ntbt_polr(Sat ~ Infl + Type + Cont)
options(op)
## qda
set.seed(123) ## make reproducible
tr <- sample(1:50, 25)
iris3df <- data.frame(cl = factor(c(rep("s",25), rep("c",25), rep("v",25))),
train = rbind(iris3[tr,,1], iris3[tr,,2], iris3[tr,,3]))
## Original function to interface
qda(cl ~ ., iris3df)
## The interface reverses the order of data and formula
ntbt_qda(iris3df, cl ~ .)
## so it can be used easily in a pipeline.
iris3df %>%
ntbt_qda(cl ~ .)
## rlm
## Original function to interface
rlm(stack.loss ~ ., stackloss)
## The interface reverses the order of data and formula
ntbt_rlm(stackloss, stack.loss ~ .)
## so it can be used easily in a pipeline.
stackloss %>%
ntbt_rlm(stack.loss ~ .) %>%
summary()
stackloss %>%
ntbt_rlm(stack.loss ~ ., psi = psi.hampel, init = "lts") %>%
summary()
stackloss %>%
ntbt_rlm(stack.loss ~ ., psi = psi.bisquare) %>%
summary()
## End(Not run)
|
First canonical correlation(s): 0.05317534
Age scores:
F0 F1 F2 F3
-0.3344445 1.4246090 -1.0320002 -0.4612728
Eth scores:
A N
-1.0563816 0.9466276
First canonical correlation(s): 0.05317534
Age scores:
F0 F1 F2 F3
-0.3344445 1.4246090 -1.0320002 -0.4612728
Eth scores:
A N
-1.0563816 0.9466276
First canonical correlation(s): 0.05317534
Age scores:
F0 F1 F2 F3
-0.3344445 1.4246090 -1.0320002 -0.4612728
Eth scores:
A N
-1.0563816 0.9466276
Call: glm.nb(formula = Days ~ Sex/(Age + Eth * Lrn), data = quine,
init.theta = 1.597990733, link = log)
Coefficients:
(Intercept) SexM SexF:AgeF1 SexM:AgeF1
3.01919 -0.47541 -0.70887 -0.72373
SexF:AgeF2 SexM:AgeF2 SexF:AgeF3 SexM:AgeF3
-0.61486 0.62842 -0.34235 1.15084
SexF:EthN SexM:EthN SexF:LrnSL SexM:LrnSL
-0.07312 -0.67899 0.94358 0.23891
SexF:EthN:LrnSL SexM:EthN:LrnSL
-1.35849 0.76142
Degrees of Freedom: 145 Total (i.e. Null); 132 Residual
Null Deviance: 234.6
Residual Deviance: 167.6 AIC: 1093
Call: glm.nb(formula = Days ~ Sex/(Age + Eth * Lrn), init.theta = 1.597990733,
link = log)
Coefficients:
(Intercept) SexM SexF:AgeF1 SexM:AgeF1
3.01919 -0.47541 -0.70887 -0.72373
SexF:AgeF2 SexM:AgeF2 SexF:AgeF3 SexM:AgeF3
-0.61486 0.62842 -0.34235 1.15084
SexF:EthN SexM:EthN SexF:LrnSL SexM:LrnSL
-0.07312 -0.67899 0.94358 0.23891
SexF:EthN:LrnSL SexM:EthN:LrnSL
-1.35849 0.76142
Degrees of Freedom: 145 Total (i.e. Null); 132 Residual
Null Deviance: 234.6
Residual Deviance: 167.6 AIC: 1093
Call: glm.nb(formula = Days ~ Sex/(Age + Eth * Lrn), init.theta = 1.597990733,
link = log)
Coefficients:
(Intercept) SexM SexF:AgeF1 SexM:AgeF1
3.01919 -0.47541 -0.70887 -0.72373
SexF:AgeF2 SexM:AgeF2 SexF:AgeF3 SexM:AgeF3
-0.61486 0.62842 -0.34235 1.15084
SexF:EthN SexM:EthN SexF:LrnSL SexM:LrnSL
-0.07312 -0.67899 0.94358 0.23891
SexF:EthN:LrnSL SexM:EthN:LrnSL
-1.35849 0.76142
Degrees of Freedom: 145 Total (i.e. Null); 132 Residual
Null Deviance: 234.6
Residual Deviance: 167.6 AIC: 1093
Call:
lda(Sp ~ ., data = Iris)
Prior probabilities of groups:
c s v
0.3333333 0.3333333 0.3333333
Group means:
Sepal.L. Sepal.W. Petal.L. Petal.W.
c 5.936 2.770 4.260 1.326
s 5.006 3.428 1.462 0.246
v 6.588 2.974 5.552 2.026
Coefficients of linear discriminants:
LD1 LD2
Sepal.L. -0.8293776 0.02410215
Sepal.W. -1.5344731 2.16452123
Petal.L. 2.2012117 -0.93192121
Petal.W. 2.8104603 2.83918785
Proportion of trace:
LD1 LD2
0.9912 0.0088
Call:
lda(Sp ~ ., data = Iris)
Prior probabilities of groups:
c s v
0.3333333 0.3333333 0.3333333
Group means:
Sepal.L. Sepal.W. Petal.L. Petal.W.
c 5.936 2.770 4.260 1.326
s 5.006 3.428 1.462 0.246
v 6.588 2.974 5.552 2.026
Coefficients of linear discriminants:
LD1 LD2
Sepal.L. -0.8293776 0.02410215
Sepal.W. -1.5344731 2.16452123
Petal.L. 2.2012117 -0.93192121
Petal.W. 2.8104603 2.83918785
Proportion of trace:
LD1 LD2
0.9912 0.0088
Call:
lda(Sp ~ ., data = .)
Prior probabilities of groups:
c s v
0.3333333 0.3333333 0.3333333
Group means:
Sepal.L. Sepal.W. Petal.L. Petal.W.
c 5.936 2.770 4.260 1.326
s 5.006 3.428 1.462 0.246
v 6.588 2.974 5.552 2.026
Coefficients of linear discriminants:
LD1 LD2
Sepal.L. -0.8293776 0.02410215
Sepal.W. -1.5344731 2.16452123
Petal.L. 2.2012117 -0.93192121
Petal.W. 2.8104603 2.83918785
Proportion of trace:
LD1 LD2
0.9912 0.0088
Length Class Mode
prior 14 -none- numeric
counts 14 -none- numeric
means 42 -none- numeric
scaling 9 -none- numeric
lev 14 -none- character
svd 3 -none- numeric
N 1 -none- numeric
call 3 -none- call
terms 3 terms call
xlevels 0 -none- list
$coefficients
(Intercept) uptake
73.71000 13.27633
$residuals
[1] -191.131260 -302.310396 -285.726242 -217.589432 -42.364407 80.857911
[7] 399.219746 -159.268071 -261.153776 -316.261799 -278.660544 -112.728950
[13] 51.649987 338.148634 -193.786526 -328.863053 -358.746051 -282.643443
[19] -143.264506 18.459165 322.217039 -167.233868 -218.669524 -225.982763
[25] -183.070977 -5.190686 131.307960 412.496075 -102.179857 -261.153776
[31] -288.381508 -238.831558 -86.176292 103.427670 363.373658 -179.182564
[37] -177.512905 -329.538128 -175.105179 -90.159191 75.547379 376.649987
[43] -119.439085 -153.615513 -171.549814 -121.999864 16.051440 171.136947
[49] 454.980327 -138.025945 -190.789233 -229.965661 -145.897256 -3.863053
[55] 188.396174 508.085643 -128.732515 -156.270778 -166.239283 -94.119573
[61] 47.914629 228.225161 557.208059 -118.111452 -96.527299 -64.011551
[67] 25.367386 167.401589 306.555501 635.538400 -80.937731 -50.060148
[73] 12.991156 103.697726 260.335891 419.404296 735.110866 -119.439085
[79] -137.683918 -61.356285 38.643715 188.643715 350.367386 662.091057
$effects
(Intercept) uptake
-3986.8408546 1308.0368214 100.4113518 -49.5886482 -149.5886482
-225.6474908 -227.2931404 637.4708418 319.8827398 157.5888507
2.2946378 -90.2934642 -155.7638810 -196.5867059 558.0576489
229.8811212 83.4698707 -60.1770739 -151.7063333 -192.8233710
-226.2342979 495.5859372 167.4094095 -11.8259608 -155.4729054
-233.2372114 -240.4712868 -234.7050384 456.4087621 135.6441323
-53.1208211 -196.7677657 -284.0616548 -268.0011936 -242.1169364
414.0550592 121.8791789 -37.2381825 -177.7085994 -237.4725816
-238.3536016 -227.2931404 351.5833475 45.6425137 -121.9455883
-220.0623022 -363.4748477 -257.4127679 -274.9410562 340.9949217
67.8782077 -118.7690605 -270.8867457 -330.6507280 -324.1198499
-213.5281870 380.1720969 90.1139018 -54.1796637 -226.4153577
-280.8851271 -290.2368876 -265.4114730 308.1708020 0.1122832
-164.2992911 -305.8285506 -386.7693843 -378.1208211 -286.5883245
320.8769129 26.5833475 -139.9459120 -302.6520228 -352.8864220
-324.1198499 -259.0584176 369.5836712 49.8778840 -83.8272557
-253.9452645 -328.5330428 -356.9439697 -284.4706393
$rank
[1] 2
$fitted.values
[1] 286.1313 477.3104 535.7262 567.5894 542.3644 594.1421 600.7803 254.2681
[9] 436.1538 566.2618 628.6605 612.7289 623.3500 661.8514 288.7865 503.8631
[17] 608.7461 632.6434 643.2645 656.5408 677.7830 262.2339 393.6695 475.9828
[25] 533.0710 505.1907 543.6920 587.5039 197.1799 436.1538 538.3815 588.8316
[33] 586.1763 571.5723 636.6263 274.1826 352.5129 579.5381 525.1052 590.1592
[41] 599.4526 623.3500 214.4391 328.6155 421.5498 471.9999 483.9486 503.8631
[49] 545.0197 233.0259 365.7892 479.9657 495.8973 503.8631 486.6038 491.9144
[57] 223.7325 331.2708 416.2393 444.1196 452.0854 446.7748 442.7919 213.1115
[65] 271.5273 314.0116 324.6326 332.5984 368.4445 364.4616 175.9377 225.0601
[73] 237.0088 246.3023 239.6641 255.5957 264.8891 214.4391 312.6839 311.3563
[81] 311.3563 311.3563 324.6326 337.9089
$assign
NULL
$qr
$qr
(Intercept) uptake
[1,] -9.1651514 -2.494121e+02
[2,] 0.1091089 9.852398e+01
[3,] 0.1091089 8.722411e-02
[4,] 0.1091089 8.722411e-02
[5,] 0.1091089 8.722411e-02
[6,] 0.1091089 8.620913e-02
[7,] 0.1091089 1.613177e-01
[8,] 0.1091089 1.227485e-01
[9,] 0.1091089 1.298533e-01
[10,] 0.1091089 1.420331e-01
[11,] 0.1091089 1.369582e-01
[12,] 0.1091089 1.440631e-01
[13,] 0.1091089 1.531979e-01
[14,] 0.1091089 1.907522e-01
[15,] 0.1091089 4.662486e-02
[16,] 0.1091089 4.357992e-02
[17,] 0.1091089 7.098441e-02
[18,] 0.1091089 7.707430e-02
[19,] 0.1091089 8.519415e-02
[20,] 0.1091089 1.176736e-01
[21,] 0.1091089 1.623327e-01
[22,] 0.1091089 -1.325903e-02
[23,] 0.1091089 -1.630398e-02
[24,] 0.1091089 -2.036390e-02
[25,] 0.1091089 -1.427401e-02
[26,] 0.1091089 7.040592e-03
[27,] 0.1091089 7.199939e-02
[28,] 0.1091089 1.542129e-01
[29,] 0.1091089 -5.081334e-02
[30,] 0.1091089 -4.675342e-02
[31,] 0.1091089 -5.994817e-02
[32,] 0.1091089 -5.385828e-02
[33,] 0.1091089 -4.167851e-02
[34,] 0.1091089 4.560988e-02
[35,] 0.1091089 1.471080e-01
[36,] 0.1091089 -9.141259e-02
[37,] 0.1091089 -5.994817e-02
[38,] 0.1091089 -4.472345e-02
[39,] 0.1091089 -3.558862e-02
[40,] 0.1091089 2.980667e-03
[41,] 0.1091089 7.402936e-02
[42,] 0.1091089 1.613177e-01
[43,] 0.1091089 -1.512965e-01
[44,] 0.1091089 -1.330268e-01
[45,] 0.1091089 -1.259220e-01
[46,] 0.1091089 -7.618787e-02
[47,] 0.1091089 -1.178021e-01
[48,] 0.1091089 5.575969e-02
[49,] 0.1091089 1.156436e-01
[50,] 0.1091089 -1.614463e-01
[51,] 0.1091089 -1.117122e-01
[52,] 0.1091089 -1.228770e-01
[53,] 0.1091089 -1.249070e-01
[54,] 0.1091089 -8.633769e-02
[55,] 0.1091089 -8.184127e-03
[56,] 0.1091089 1.745125e-01
[57,] 0.1091089 -1.238920e-01
[58,] 0.1091089 -9.039761e-02
[59,] 0.1091089 -6.096315e-02
[60,] 0.1091089 -8.227776e-02
[61,] 0.1091089 -3.863357e-02
[62,] 0.1091089 2.429527e-02
[63,] 0.1091089 1.247784e-01
[64,] 0.1091089 -1.929107e-01
[65,] 0.1091089 -1.766710e-01
[66,] 0.1091089 -1.665212e-01
[67,] 0.1091089 -1.584014e-01
[68,] 0.1091089 -1.401317e-01
[69,] 0.1091089 -5.994817e-02
[70,] 0.1091089 1.044788e-01
[71,] 0.1091089 -1.807309e-01
[72,] 0.1091089 -1.512965e-01
[73,] 0.1091089 -1.431766e-01
[74,] 0.1091089 -1.553564e-01
[75,] 0.1091089 -1.076523e-01
[76,] 0.1091089 -8.184127e-03
[77,] 0.1091089 1.308683e-01
[78,] 0.1091089 -1.340418e-01
[79,] 0.1091089 -1.289669e-01
[80,] 0.1091089 -8.938263e-02
[81,] 0.1091089 -1.086673e-01
[82,] 0.1091089 -8.430772e-02
[83,] 0.1091089 -3.964855e-02
[84,] 0.1091089 1.065088e-01
$qraux
[1] 1.109109 1.077074
$pivot
[1] 1 2
$tol
[1] 1e-07
$rank
[1] 2
attr(,"class")
[1] "qr"
$df.residual
[1] 82
$terms
conc ~ uptake
attr(,"variables")
list(conc, uptake)
attr(,"factors")
uptake
conc 0
uptake 1
attr(,"term.labels")
[1] "uptake"
attr(,"order")
[1] 1
attr(,"intercept")
[1] 1
attr(,"response")
[1] 1
attr(,".Environment")
<environment: R_GlobalEnv>
attr(,"predvars")
list(conc, uptake)
attr(,"dataClasses")
conc uptake
"numeric" "numeric"
$call
lm.gls(formula = conc ~ uptake, data = CO2, W = diag(nrow(CO2)))
$xlevels
named list()
attr(,"class")
[1] "lm.gls"
$coefficients
(Intercept) uptake
73.71000 13.27633
$residuals
[1] -191.131260 -302.310396 -285.726242 -217.589432 -42.364407 80.857911
[7] 399.219746 -159.268071 -261.153776 -316.261799 -278.660544 -112.728950
[13] 51.649987 338.148634 -193.786526 -328.863053 -358.746051 -282.643443
[19] -143.264506 18.459165 322.217039 -167.233868 -218.669524 -225.982763
[25] -183.070977 -5.190686 131.307960 412.496075 -102.179857 -261.153776
[31] -288.381508 -238.831558 -86.176292 103.427670 363.373658 -179.182564
[37] -177.512905 -329.538128 -175.105179 -90.159191 75.547379 376.649987
[43] -119.439085 -153.615513 -171.549814 -121.999864 16.051440 171.136947
[49] 454.980327 -138.025945 -190.789233 -229.965661 -145.897256 -3.863053
[55] 188.396174 508.085643 -128.732515 -156.270778 -166.239283 -94.119573
[61] 47.914629 228.225161 557.208059 -118.111452 -96.527299 -64.011551
[67] 25.367386 167.401589 306.555501 635.538400 -80.937731 -50.060148
[73] 12.991156 103.697726 260.335891 419.404296 735.110866 -119.439085
[79] -137.683918 -61.356285 38.643715 188.643715 350.367386 662.091057
$effects
(Intercept) uptake
-3986.8408546 1308.0368214 100.4113518 -49.5886482 -149.5886482
-225.6474908 -227.2931404 637.4708418 319.8827398 157.5888507
2.2946378 -90.2934642 -155.7638810 -196.5867059 558.0576489
229.8811212 83.4698707 -60.1770739 -151.7063333 -192.8233710
-226.2342979 495.5859372 167.4094095 -11.8259608 -155.4729054
-233.2372114 -240.4712868 -234.7050384 456.4087621 135.6441323
-53.1208211 -196.7677657 -284.0616548 -268.0011936 -242.1169364
414.0550592 121.8791789 -37.2381825 -177.7085994 -237.4725816
-238.3536016 -227.2931404 351.5833475 45.6425137 -121.9455883
-220.0623022 -363.4748477 -257.4127679 -274.9410562 340.9949217
67.8782077 -118.7690605 -270.8867457 -330.6507280 -324.1198499
-213.5281870 380.1720969 90.1139018 -54.1796637 -226.4153577
-280.8851271 -290.2368876 -265.4114730 308.1708020 0.1122832
-164.2992911 -305.8285506 -386.7693843 -378.1208211 -286.5883245
320.8769129 26.5833475 -139.9459120 -302.6520228 -352.8864220
-324.1198499 -259.0584176 369.5836712 49.8778840 -83.8272557
-253.9452645 -328.5330428 -356.9439697 -284.4706393
$rank
[1] 2
$fitted.values
[1] 286.1313 477.3104 535.7262 567.5894 542.3644 594.1421 600.7803 254.2681
[9] 436.1538 566.2618 628.6605 612.7289 623.3500 661.8514 288.7865 503.8631
[17] 608.7461 632.6434 643.2645 656.5408 677.7830 262.2339 393.6695 475.9828
[25] 533.0710 505.1907 543.6920 587.5039 197.1799 436.1538 538.3815 588.8316
[33] 586.1763 571.5723 636.6263 274.1826 352.5129 579.5381 525.1052 590.1592
[41] 599.4526 623.3500 214.4391 328.6155 421.5498 471.9999 483.9486 503.8631
[49] 545.0197 233.0259 365.7892 479.9657 495.8973 503.8631 486.6038 491.9144
[57] 223.7325 331.2708 416.2393 444.1196 452.0854 446.7748 442.7919 213.1115
[65] 271.5273 314.0116 324.6326 332.5984 368.4445 364.4616 175.9377 225.0601
[73] 237.0088 246.3023 239.6641 255.5957 264.8891 214.4391 312.6839 311.3563
[81] 311.3563 311.3563 324.6326 337.9089
$assign
NULL
$qr
$qr
(Intercept) uptake
[1,] -9.1651514 -2.494121e+02
[2,] 0.1091089 9.852398e+01
[3,] 0.1091089 8.722411e-02
[4,] 0.1091089 8.722411e-02
[5,] 0.1091089 8.722411e-02
[6,] 0.1091089 8.620913e-02
[7,] 0.1091089 1.613177e-01
[8,] 0.1091089 1.227485e-01
[9,] 0.1091089 1.298533e-01
[10,] 0.1091089 1.420331e-01
[11,] 0.1091089 1.369582e-01
[12,] 0.1091089 1.440631e-01
[13,] 0.1091089 1.531979e-01
[14,] 0.1091089 1.907522e-01
[15,] 0.1091089 4.662486e-02
[16,] 0.1091089 4.357992e-02
[17,] 0.1091089 7.098441e-02
[18,] 0.1091089 7.707430e-02
[19,] 0.1091089 8.519415e-02
[20,] 0.1091089 1.176736e-01
[21,] 0.1091089 1.623327e-01
[22,] 0.1091089 -1.325903e-02
[23,] 0.1091089 -1.630398e-02
[24,] 0.1091089 -2.036390e-02
[25,] 0.1091089 -1.427401e-02
[26,] 0.1091089 7.040592e-03
[27,] 0.1091089 7.199939e-02
[28,] 0.1091089 1.542129e-01
[29,] 0.1091089 -5.081334e-02
[30,] 0.1091089 -4.675342e-02
[31,] 0.1091089 -5.994817e-02
[32,] 0.1091089 -5.385828e-02
[33,] 0.1091089 -4.167851e-02
[34,] 0.1091089 4.560988e-02
[35,] 0.1091089 1.471080e-01
[36,] 0.1091089 -9.141259e-02
[37,] 0.1091089 -5.994817e-02
[38,] 0.1091089 -4.472345e-02
[39,] 0.1091089 -3.558862e-02
[40,] 0.1091089 2.980667e-03
[41,] 0.1091089 7.402936e-02
[42,] 0.1091089 1.613177e-01
[43,] 0.1091089 -1.512965e-01
[44,] 0.1091089 -1.330268e-01
[45,] 0.1091089 -1.259220e-01
[46,] 0.1091089 -7.618787e-02
[47,] 0.1091089 -1.178021e-01
[48,] 0.1091089 5.575969e-02
[49,] 0.1091089 1.156436e-01
[50,] 0.1091089 -1.614463e-01
[51,] 0.1091089 -1.117122e-01
[52,] 0.1091089 -1.228770e-01
[53,] 0.1091089 -1.249070e-01
[54,] 0.1091089 -8.633769e-02
[55,] 0.1091089 -8.184127e-03
[56,] 0.1091089 1.745125e-01
[57,] 0.1091089 -1.238920e-01
[58,] 0.1091089 -9.039761e-02
[59,] 0.1091089 -6.096315e-02
[60,] 0.1091089 -8.227776e-02
[61,] 0.1091089 -3.863357e-02
[62,] 0.1091089 2.429527e-02
[63,] 0.1091089 1.247784e-01
[64,] 0.1091089 -1.929107e-01
[65,] 0.1091089 -1.766710e-01
[66,] 0.1091089 -1.665212e-01
[67,] 0.1091089 -1.584014e-01
[68,] 0.1091089 -1.401317e-01
[69,] 0.1091089 -5.994817e-02
[70,] 0.1091089 1.044788e-01
[71,] 0.1091089 -1.807309e-01
[72,] 0.1091089 -1.512965e-01
[73,] 0.1091089 -1.431766e-01
[74,] 0.1091089 -1.553564e-01
[75,] 0.1091089 -1.076523e-01
[76,] 0.1091089 -8.184127e-03
[77,] 0.1091089 1.308683e-01
[78,] 0.1091089 -1.340418e-01
[79,] 0.1091089 -1.289669e-01
[80,] 0.1091089 -8.938263e-02
[81,] 0.1091089 -1.086673e-01
[82,] 0.1091089 -8.430772e-02
[83,] 0.1091089 -3.964855e-02
[84,] 0.1091089 1.065088e-01
$qraux
[1] 1.109109 1.077074
$pivot
[1] 1 2
$tol
[1] 1e-07
$rank
[1] 2
attr(,"class")
[1] "qr"
$df.residual
[1] 82
$terms
conc ~ uptake
attr(,"variables")
list(conc, uptake)
attr(,"factors")
uptake
conc 0
uptake 1
attr(,"term.labels")
[1] "uptake"
attr(,"order")
[1] 1
attr(,"intercept")
[1] 1
attr(,"response")
[1] 1
attr(,".Environment")
<environment: 0x470b3e0>
attr(,"predvars")
list(conc, uptake)
attr(,"dataClasses")
conc uptake
"numeric" "numeric"
$call
lm.gls(formula = conc ~ uptake, W = diag(nrow(CO2)))
$xlevels
named list()
attr(,"class")
[1] "lm.gls"
$coefficients
(Intercept) uptake
73.71000 13.27633
$residuals
[1] -191.131260 -302.310396 -285.726242 -217.589432 -42.364407 80.857911
[7] 399.219746 -159.268071 -261.153776 -316.261799 -278.660544 -112.728950
[13] 51.649987 338.148634 -193.786526 -328.863053 -358.746051 -282.643443
[19] -143.264506 18.459165 322.217039 -167.233868 -218.669524 -225.982763
[25] -183.070977 -5.190686 131.307960 412.496075 -102.179857 -261.153776
[31] -288.381508 -238.831558 -86.176292 103.427670 363.373658 -179.182564
[37] -177.512905 -329.538128 -175.105179 -90.159191 75.547379 376.649987
[43] -119.439085 -153.615513 -171.549814 -121.999864 16.051440 171.136947
[49] 454.980327 -138.025945 -190.789233 -229.965661 -145.897256 -3.863053
[55] 188.396174 508.085643 -128.732515 -156.270778 -166.239283 -94.119573
[61] 47.914629 228.225161 557.208059 -118.111452 -96.527299 -64.011551
[67] 25.367386 167.401589 306.555501 635.538400 -80.937731 -50.060148
[73] 12.991156 103.697726 260.335891 419.404296 735.110866 -119.439085
[79] -137.683918 -61.356285 38.643715 188.643715 350.367386 662.091057
$effects
(Intercept) uptake
-3986.8408546 1308.0368214 100.4113518 -49.5886482 -149.5886482
-225.6474908 -227.2931404 637.4708418 319.8827398 157.5888507
2.2946378 -90.2934642 -155.7638810 -196.5867059 558.0576489
229.8811212 83.4698707 -60.1770739 -151.7063333 -192.8233710
-226.2342979 495.5859372 167.4094095 -11.8259608 -155.4729054
-233.2372114 -240.4712868 -234.7050384 456.4087621 135.6441323
-53.1208211 -196.7677657 -284.0616548 -268.0011936 -242.1169364
414.0550592 121.8791789 -37.2381825 -177.7085994 -237.4725816
-238.3536016 -227.2931404 351.5833475 45.6425137 -121.9455883
-220.0623022 -363.4748477 -257.4127679 -274.9410562 340.9949217
67.8782077 -118.7690605 -270.8867457 -330.6507280 -324.1198499
-213.5281870 380.1720969 90.1139018 -54.1796637 -226.4153577
-280.8851271 -290.2368876 -265.4114730 308.1708020 0.1122832
-164.2992911 -305.8285506 -386.7693843 -378.1208211 -286.5883245
320.8769129 26.5833475 -139.9459120 -302.6520228 -352.8864220
-324.1198499 -259.0584176 369.5836712 49.8778840 -83.8272557
-253.9452645 -328.5330428 -356.9439697 -284.4706393
$rank
[1] 2
$fitted.values
[1] 286.1313 477.3104 535.7262 567.5894 542.3644 594.1421 600.7803 254.2681
[9] 436.1538 566.2618 628.6605 612.7289 623.3500 661.8514 288.7865 503.8631
[17] 608.7461 632.6434 643.2645 656.5408 677.7830 262.2339 393.6695 475.9828
[25] 533.0710 505.1907 543.6920 587.5039 197.1799 436.1538 538.3815 588.8316
[33] 586.1763 571.5723 636.6263 274.1826 352.5129 579.5381 525.1052 590.1592
[41] 599.4526 623.3500 214.4391 328.6155 421.5498 471.9999 483.9486 503.8631
[49] 545.0197 233.0259 365.7892 479.9657 495.8973 503.8631 486.6038 491.9144
[57] 223.7325 331.2708 416.2393 444.1196 452.0854 446.7748 442.7919 213.1115
[65] 271.5273 314.0116 324.6326 332.5984 368.4445 364.4616 175.9377 225.0601
[73] 237.0088 246.3023 239.6641 255.5957 264.8891 214.4391 312.6839 311.3563
[81] 311.3563 311.3563 324.6326 337.9089
$assign
NULL
$qr
$qr
(Intercept) uptake
[1,] -9.1651514 -2.494121e+02
[2,] 0.1091089 9.852398e+01
[3,] 0.1091089 8.722411e-02
[4,] 0.1091089 8.722411e-02
[5,] 0.1091089 8.722411e-02
[6,] 0.1091089 8.620913e-02
[7,] 0.1091089 1.613177e-01
[8,] 0.1091089 1.227485e-01
[9,] 0.1091089 1.298533e-01
[10,] 0.1091089 1.420331e-01
[11,] 0.1091089 1.369582e-01
[12,] 0.1091089 1.440631e-01
[13,] 0.1091089 1.531979e-01
[14,] 0.1091089 1.907522e-01
[15,] 0.1091089 4.662486e-02
[16,] 0.1091089 4.357992e-02
[17,] 0.1091089 7.098441e-02
[18,] 0.1091089 7.707430e-02
[19,] 0.1091089 8.519415e-02
[20,] 0.1091089 1.176736e-01
[21,] 0.1091089 1.623327e-01
[22,] 0.1091089 -1.325903e-02
[23,] 0.1091089 -1.630398e-02
[24,] 0.1091089 -2.036390e-02
[25,] 0.1091089 -1.427401e-02
[26,] 0.1091089 7.040592e-03
[27,] 0.1091089 7.199939e-02
[28,] 0.1091089 1.542129e-01
[29,] 0.1091089 -5.081334e-02
[30,] 0.1091089 -4.675342e-02
[31,] 0.1091089 -5.994817e-02
[32,] 0.1091089 -5.385828e-02
[33,] 0.1091089 -4.167851e-02
[34,] 0.1091089 4.560988e-02
[35,] 0.1091089 1.471080e-01
[36,] 0.1091089 -9.141259e-02
[37,] 0.1091089 -5.994817e-02
[38,] 0.1091089 -4.472345e-02
[39,] 0.1091089 -3.558862e-02
[40,] 0.1091089 2.980667e-03
[41,] 0.1091089 7.402936e-02
[42,] 0.1091089 1.613177e-01
[43,] 0.1091089 -1.512965e-01
[44,] 0.1091089 -1.330268e-01
[45,] 0.1091089 -1.259220e-01
[46,] 0.1091089 -7.618787e-02
[47,] 0.1091089 -1.178021e-01
[48,] 0.1091089 5.575969e-02
[49,] 0.1091089 1.156436e-01
[50,] 0.1091089 -1.614463e-01
[51,] 0.1091089 -1.117122e-01
[52,] 0.1091089 -1.228770e-01
[53,] 0.1091089 -1.249070e-01
[54,] 0.1091089 -8.633769e-02
[55,] 0.1091089 -8.184127e-03
[56,] 0.1091089 1.745125e-01
[57,] 0.1091089 -1.238920e-01
[58,] 0.1091089 -9.039761e-02
[59,] 0.1091089 -6.096315e-02
[60,] 0.1091089 -8.227776e-02
[61,] 0.1091089 -3.863357e-02
[62,] 0.1091089 2.429527e-02
[63,] 0.1091089 1.247784e-01
[64,] 0.1091089 -1.929107e-01
[65,] 0.1091089 -1.766710e-01
[66,] 0.1091089 -1.665212e-01
[67,] 0.1091089 -1.584014e-01
[68,] 0.1091089 -1.401317e-01
[69,] 0.1091089 -5.994817e-02
[70,] 0.1091089 1.044788e-01
[71,] 0.1091089 -1.807309e-01
[72,] 0.1091089 -1.512965e-01
[73,] 0.1091089 -1.431766e-01
[74,] 0.1091089 -1.553564e-01
[75,] 0.1091089 -1.076523e-01
[76,] 0.1091089 -8.184127e-03
[77,] 0.1091089 1.308683e-01
[78,] 0.1091089 -1.340418e-01
[79,] 0.1091089 -1.289669e-01
[80,] 0.1091089 -8.938263e-02
[81,] 0.1091089 -1.086673e-01
[82,] 0.1091089 -8.430772e-02
[83,] 0.1091089 -3.964855e-02
[84,] 0.1091089 1.065088e-01
$qraux
[1] 1.109109 1.077074
$pivot
[1] 1 2
$tol
[1] 1e-07
$rank
[1] 2
attr(,"class")
[1] "qr"
$df.residual
[1] 82
$terms
conc ~ uptake
attr(,"variables")
list(conc, uptake)
attr(,"factors")
uptake
conc 0
uptake 1
attr(,"term.labels")
[1] "uptake"
attr(,"order")
[1] 1
attr(,"intercept")
[1] 1
attr(,"response")
[1] 1
attr(,".Environment")
<environment: 0x47ab1e8>
attr(,"predvars")
list(conc, uptake)
attr(,"dataClasses")
conc uptake
"numeric" "numeric"
$call
lm.gls(formula = conc ~ uptake, W = diag(nrow(CO2)))
$xlevels
named list()
attr(,"class")
[1] "lm.gls"
GNP Unemployed Armed.Forces Population
2946.85636017 0.26352725 0.03648291 0.01116105 -1.73702984
Year Employed
-1.41879853 0.23128785
GNP Unemployed Armed.Forces Population
2946.85636017 0.26352725 0.03648291 0.01116105 -1.73702984
Year Employed
-1.41879853 0.23128785
GNP Unemployed Armed.Forces Population
2946.85636017 0.26352725 0.03648291 0.01116105 -1.73702984
Year Employed
-1.41879853 0.23128785
Call:
loglm(formula = ~Type + Origin, data = xtCars93)
Statistics:
X^2 df P(> X^2)
Likelihood Ratio 18.36179 5 0.00252554
Pearson 14.07985 5 0.01511005
Call:
loglm(formula = ~Type + Origin, data = xtCars93)
Statistics:
X^2 df P(> X^2)
Likelihood Ratio 18.36179 5 0.00252554
Pearson 14.07985 5 0.01511005
Call:
loglm(formula = ~Type + Origin, data = .)
Statistics:
X^2 df P(> X^2)
Likelihood Ratio 18.36179 5 0.00252554
Pearson 14.07985 5 0.01511005
$x
[1] 0.7500000 0.8080808 0.8661616 0.9242424 0.9823232 1.0404040 1.0984848
[8] 1.1565657 1.2146465 1.2727273 1.3308081 1.3888889 1.4469697 1.5050505
[15] 1.5631313 1.6212121 1.6792929 1.7373737 1.7954545 1.8535354 1.9116162
[22] 1.9696970 2.0277778 2.0858586 2.1439394 2.2020202 2.2601010 2.3181818
[29] 2.3762626 2.4343434 2.4924242 2.5505051 2.6085859 2.6666667 2.7247475
[36] 2.7828283 2.8409091 2.8989899 2.9570707 3.0151515 3.0732323 3.1313131
[43] 3.1893939 3.2474747 3.3055556 3.3636364 3.4217172 3.4797980 3.5378788
[50] 3.5959596 3.6540404 3.7121212 3.7702020 3.8282828 3.8863636 3.9444444
[57] 4.0025253 4.0606061 4.1186869 4.1767677 4.2348485 4.2929293 4.3510101
[64] 4.4090909 4.4671717 4.5252525 4.5833333 4.6414141 4.6994949 4.7575758
[71] 4.8156566 4.8737374 4.9318182 4.9898990 5.0479798 5.1060606 5.1641414
[78] 5.2222222 5.2803030 5.3383838 5.3964646 5.4545455 5.5126263 5.5707071
[85] 5.6287879 5.6868687 5.7449495 5.8030303 5.8611111 5.9191919 5.9772727
[92] 6.0353535 6.0934343 6.1515152 6.2095960 6.2676768 6.3257576 6.3838384
[99] 6.4419192 6.5000000
$y
[1] -686.8651 -686.3368 -685.8576 -685.4242 -685.0337 -684.6829 -684.3690
[8] -684.0884 -683.8376 -683.6132 -683.4116 -683.2294 -683.0643 -682.9154
[15] -682.7814 -682.6612 -682.5538 -682.4581 -682.3734 -682.2987 -682.2334
[22] -682.1765 -682.1274 -682.0855 -682.0504 -682.0217 -681.9990 -681.9817
[29] -681.9696 -681.9624 -681.9597 -681.9613 -681.9668 -681.9761 -681.9889
[36] -682.0050 -682.0243 -682.0464 -682.0713 -682.0989 -682.1289 -682.1612
[43] -682.1957 -682.2323 -682.2708 -682.3112 -682.3533 -682.3971 -682.4425
[50] -682.4893 -682.5375 -682.5870 -682.6378 -682.6898 -682.7429 -682.7971
[57] -682.8523 -682.9084 -682.9654 -683.0233 -683.0820 -683.1414 -683.2016
[64] -683.2624 -683.3239 -683.3860 -683.4487 -683.5119 -683.5756 -683.6398
[71] -683.7045 -683.7696 -683.8351 -683.9009 -683.9672 -684.0337 -684.1006
[78] -684.1677 -684.2352 -684.3029 -684.3708 -684.4389 -684.5073 -684.5758
[85] -684.6445 -684.7134 -684.7824 -684.8515 -684.9208 -684.9901 -685.0596
[92] -685.1291 -685.1988 -685.2685 -685.3382 -685.4080 -685.4779 -685.5477
[99] -685.6176 -685.6875
$x
[1] 0.7500000 0.8080808 0.8661616 0.9242424 0.9823232 1.0404040 1.0984848
[8] 1.1565657 1.2146465 1.2727273 1.3308081 1.3888889 1.4469697 1.5050505
[15] 1.5631313 1.6212121 1.6792929 1.7373737 1.7954545 1.8535354 1.9116162
[22] 1.9696970 2.0277778 2.0858586 2.1439394 2.2020202 2.2601010 2.3181818
[29] 2.3762626 2.4343434 2.4924242 2.5505051 2.6085859 2.6666667 2.7247475
[36] 2.7828283 2.8409091 2.8989899 2.9570707 3.0151515 3.0732323 3.1313131
[43] 3.1893939 3.2474747 3.3055556 3.3636364 3.4217172 3.4797980 3.5378788
[50] 3.5959596 3.6540404 3.7121212 3.7702020 3.8282828 3.8863636 3.9444444
[57] 4.0025253 4.0606061 4.1186869 4.1767677 4.2348485 4.2929293 4.3510101
[64] 4.4090909 4.4671717 4.5252525 4.5833333 4.6414141 4.6994949 4.7575758
[71] 4.8156566 4.8737374 4.9318182 4.9898990 5.0479798 5.1060606 5.1641414
[78] 5.2222222 5.2803030 5.3383838 5.3964646 5.4545455 5.5126263 5.5707071
[85] 5.6287879 5.6868687 5.7449495 5.8030303 5.8611111 5.9191919 5.9772727
[92] 6.0353535 6.0934343 6.1515152 6.2095960 6.2676768 6.3257576 6.3838384
[99] 6.4419192 6.5000000
$y
[1] -686.8651 -686.3368 -685.8576 -685.4242 -685.0337 -684.6829 -684.3690
[8] -684.0884 -683.8376 -683.6132 -683.4116 -683.2294 -683.0643 -682.9154
[15] -682.7814 -682.6612 -682.5538 -682.4581 -682.3734 -682.2987 -682.2334
[22] -682.1765 -682.1274 -682.0855 -682.0504 -682.0217 -681.9990 -681.9817
[29] -681.9696 -681.9624 -681.9597 -681.9613 -681.9668 -681.9761 -681.9889
[36] -682.0050 -682.0243 -682.0464 -682.0713 -682.0989 -682.1289 -682.1612
[43] -682.1957 -682.2323 -682.2708 -682.3112 -682.3533 -682.3971 -682.4425
[50] -682.4893 -682.5375 -682.5870 -682.6378 -682.6898 -682.7429 -682.7971
[57] -682.8523 -682.9084 -682.9654 -683.0233 -683.0820 -683.1414 -683.2016
[64] -683.2624 -683.3239 -683.3860 -683.4487 -683.5119 -683.5756 -683.6398
[71] -683.7045 -683.7696 -683.8351 -683.9009 -683.9672 -684.0337 -684.1006
[78] -684.1677 -684.2352 -684.3029 -684.3708 -684.4389 -684.5073 -684.5758
[85] -684.6445 -684.7134 -684.7824 -684.8515 -684.9208 -684.9901 -685.0596
[92] -685.1291 -685.1988 -685.2685 -685.3382 -685.4080 -685.4779 -685.5477
[99] -685.6176 -685.6875
Call:
polr(formula = Sat ~ Infl + Type + Cont, data = housing)
Coefficients:
InflMedium InflHigh TypeApartment TypeAtrium TypeTerrace
1.438398e-12 1.442159e-12 1.081215e-12 1.080646e-12 1.077125e-12
ContHigh
2.158385e-12
Intercepts:
Low|Medium Medium|High
-0.6931472 0.6931472
Residual Deviance: 158.2002
AIC: 174.2002
Call:
polr(formula = Sat ~ Infl + Type + Cont)
Coefficients:
InflMedium InflHigh TypeApartment TypeAtrium TypeTerrace
1.438398e-12 1.442159e-12 1.081215e-12 1.080646e-12 1.077125e-12
ContHigh
2.158385e-12
Intercepts:
Low|Medium Medium|High
-0.6931472 0.6931472
Residual Deviance: 158.2002
AIC: 174.2002
Call:
polr(formula = Sat ~ Infl + Type + Cont)
Coefficients:
InflMedium InflHigh TypeApartment TypeAtrium TypeTerrace
1.438398e-12 1.442159e-12 1.081215e-12 1.080646e-12 1.077125e-12
ContHigh
2.158385e-12
Intercepts:
Low|Medium Medium|High
-0.6931472 0.6931472
Residual Deviance: 158.2002
AIC: 174.2002
Call:
qda(cl ~ ., data = iris3df)
Prior probabilities of groups:
c s v
0.3333333 0.3333333 0.3333333
Group means:
train.Sepal.L. train.Sepal.W. train.Petal.L. train.Petal.W.
c 5.892 2.708 4.140 1.276
s 4.996 3.324 1.484 0.272
v 6.628 2.980 5.504 2.000
Call:
qda(cl ~ ., data = iris3df)
Prior probabilities of groups:
c s v
0.3333333 0.3333333 0.3333333
Group means:
train.Sepal.L. train.Sepal.W. train.Petal.L. train.Petal.W.
c 5.892 2.708 4.140 1.276
s 4.996 3.324 1.484 0.272
v 6.628 2.980 5.504 2.000
Call:
qda(cl ~ ., data = .)
Prior probabilities of groups:
c s v
0.3333333 0.3333333 0.3333333
Group means:
train.Sepal.L. train.Sepal.W. train.Petal.L. train.Petal.W.
c 5.892 2.708 4.140 1.276
s 4.996 3.324 1.484 0.272
v 6.628 2.980 5.504 2.000
Call:
rlm(formula = stack.loss ~ ., data = stackloss)
Converged in 9 iterations
Coefficients:
(Intercept) Air.Flow Water.Temp Acid.Conc.
-41.0265311 0.8293739 0.9261082 -0.1278492
Degrees of freedom: 21 total; 17 residual
Scale estimate: 2.44
Call:
rlm(formula = stack.loss ~ ., data = stackloss)
Converged in 9 iterations
Coefficients:
(Intercept) Air.Flow Water.Temp Acid.Conc.
-41.0265311 0.8293739 0.9261082 -0.1278492
Degrees of freedom: 21 total; 17 residual
Scale estimate: 2.44
Call: rlm(formula = stack.loss ~ ., data = .)
Residuals:
Min 1Q Median 3Q Max
-8.91753 -1.73127 0.06187 1.54306 6.50163
Coefficients:
Value Std. Error t value
(Intercept) -41.0265 9.8073 -4.1832
Air.Flow 0.8294 0.1112 7.4597
Water.Temp 0.9261 0.3034 3.0524
Acid.Conc. -0.1278 0.1289 -0.9922
Residual standard error: 2.441 on 17 degrees of freedom
Call: rlm(formula = stack.loss ~ ., data = ., psi = psi.hampel, init = "lts")
Residuals:
Min 1Q Median 3Q Max
-7.6599 -1.7635 -0.3594 2.3209 5.7864
Coefficients:
Value Std. Error t value
(Intercept) -40.4748 11.8929 -3.4033
Air.Flow 0.7411 0.1348 5.4967
Water.Temp 1.2251 0.3679 3.3296
Acid.Conc. -0.1455 0.1563 -0.9313
Residual standard error: 3.088 on 17 degrees of freedom
Warning message:
In is.na(y) : is.na() applied to non-(list or vector) of type 'language'
Call: rlm(formula = stack.loss ~ ., data = ., psi = psi.bisquare)
Residuals:
Min 1Q Median 3Q Max
-10.4356 -1.7065 -0.2392 0.8797 6.9326
Coefficients:
Value Std. Error t value
(Intercept) -42.2853 9.5316 -4.4363
Air.Flow 0.9275 0.1081 8.5841
Water.Temp 0.6507 0.2949 2.2068
Acid.Conc. -0.1123 0.1252 -0.8970
Residual standard error: 2.282 on 17 degrees of freedom
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.