Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
## ----setup--------------------------------------------------------------------
library('genset')
## -----------------------------------------------------------------------------
library('genset')
## -----------------------------------------------------------------------------
y <- mtcars$mpg
x1 <- mtcars$hp
x2 <- mtcars$wt
## -----------------------------------------------------------------------------
set1 <- data.frame(y, x1, x2)
## -----------------------------------------------------------------------------
multi.fun <- function(x) {
c(mean = mean(x), media=median(x), sd=sd(x))
}
round(multi.fun(set1$y), 1)
round(multi.fun(set1$x1), 0)
round(multi.fun(set1$x2), 3)
## -----------------------------------------------------------------------------
summary(lm(y ~ x1, x2, data=set1))
## -----------------------------------------------------------------------------
set.seed(101)
set2 <- genset(y, x1, x2, method=1, option="x1")
## -----------------------------------------------------------------------------
round(multi.fun(set2$y), 1)
round(multi.fun(set2$x1), 0)
round(multi.fun(set2$x2), 3)
## -----------------------------------------------------------------------------
summary(lm(y ~ x1 + x2, data=set2))
## -----------------------------------------------------------------------------
y <- mtcars$mpg
x1 <- mtcars$wt
x2 <- mtcars$vs
## -----------------------------------------------------------------------------
set3 <- data.frame(y, x1, x2)
## -----------------------------------------------------------------------------
v.shape <- subset(set3, x2==0)
straight <- subset(set3, x2==1)
## -----------------------------------------------------------------------------
multi.fun <- function(x) {
c(mean = mean(x), media=median(x), sd=sd(x))
}
round(multi.fun(v.shape$y), 1)
round(multi.fun(v.shape$x1), 3)
round(multi.fun(straight$y), 1)
round(multi.fun(straight$x1), 3)
## -----------------------------------------------------------------------------
summary(lm(y ~ x1 + factor(x2), data=set3))
## -----------------------------------------------------------------------------
set.seed(123)
set4 <- genset(y, x1, factor(x2), method=2, option="x2")
## -----------------------------------------------------------------------------
v.shape <- subset(set4, x2==0)
straight <- subset(set4, x2==1)
## -----------------------------------------------------------------------------
multi.fun <- function(x) {
c(mean = mean(x), media=median(x), sd=sd(x))
}
round(multi.fun(v.shape$y), 1)
round(multi.fun(v.shape$x1), 3)
round(multi.fun(straight$y), 1)
round(multi.fun(straight$x1), 3)
## -----------------------------------------------------------------------------
summary(lm(y ~ x1 + factor(x2), data=set4))
## ---- echo=FALSE, include=FALSE, results='asis'-------------------------------
knitr::kable(head(set1, 10))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.