Nothing
## ----setup, include = FALSE---------------------------------------------------
library("smurf")
## -----------------------------------------------------------------------------
data("rent", package = "catdata")
# Urban district in Munich
rent$area <- as.factor(rent$area)
# Decade of construction
rent$year <- as.factor(floor(rent$year / 10) * 10)
# Number of rooms
rent$rooms <- as.factor(rent$rooms)
# Quality of the house with levels "fair", "good" and "excellent"
rent$quality <- as.factor(rent$good + 2 * rent$best)
levels(rent$quality) <- c("fair", "good", "excellent")
# Floor space divided in categories (0, 30), [30, 40), ..., [130, 140)
sizeClasses <- c(0, seq(30, 140, 10))
rent$size <- as.factor(sizeClasses[findInterval(rent$size, sizeClasses)])
# Is warm water present?
rent$warm <- factor(rent$warm, labels = c("yes", "no"))
# Is central heating present?
rent$central <- factor(rent$central, labels = c("yes", "no"))
# Does the bathroom have tiles?
rent$tiles <- factor(rent$tiles, labels = c("yes", "no"))
# Is there special furniture in the bathroom?
rent$bathextra <- factor(rent$bathextra, labels = c("no", "yes"))
# Is the kitchen well-equipped?
rent$kitchen <- factor(rent$kitchen, labels = c("no", "yes"))
## ---- eval = FALSE------------------------------------------------------------
# rentm ~ p(area, pen = "gflasso") +
# p(year, pen = "flasso", refcat = 2000) + p(rooms, pen = "flasso") +
# p(quality, pen = "flasso") + p(size, pen = "flasso") +
# p(warm, pen = "grouplasso", group = 1) + p(central, pen = "grouplasso", group = 1) +
# p(tiles, pen = "none") + bathextra +
# p(kitchen, pen = "lasso")
## ---- eval = FALSE------------------------------------------------------------
# p(pred1, pred2, pen = "2dflasso")
## -----------------------------------------------------------------------------
formu <- rentm ~ p(area, pen = "gflasso") +
p(year, pen = "flasso") + p(rooms, pen = "flasso") +
p(quality, pen = "flasso") + p(size, pen = "flasso") +
p(warm, pen = "lasso") + p(central, pen = "lasso") +
p(tiles, pen = "lasso") + p(bathextra, pen = "lasso") +
p(kitchen, pen = "lasso")
## ---- warning = FALSE---------------------------------------------------------
munich.fit <- glmsmurf(formula = formu, family = gaussian(), data = rent,
pen.weights = "glm.stand", lambda = 0.01404071)
## ---- warning = FALSE---------------------------------------------------------
plot(munich.fit)
plot_reest(munich.fit)
## -----------------------------------------------------------------------------
summary(munich.fit)
## ---- eval = FALSE------------------------------------------------------------
# munich.fit.cv <- glmsmurf(formula = formu, family = gaussian(), data = rent,
# pen.weights = "glm.stand", lambda = "cv1se.dev",
# control = list(lambda.length = 25L))
## ---- eval = FALSE------------------------------------------------------------
# munich.fit.cv$lambda
## ---- eval = FALSE------------------------------------------------------------
# plot_lambda(munich.fit.cv)
## ----fig1, echo = FALSE, out.width = "545px"----------------------------------
knitr::include_graphics("Plot_lambda1.png")
## ---- eval = FALSE------------------------------------------------------------
# # Zoomed plot
# plot_lambda(munich.fit.cv, xlim = c(-7, -3.5), ylim = c(1575, 1750))
## ----fig2, echo = FALSE, out.width = "545px"----------------------------------
knitr::include_graphics("Plot_lambda2.png")
## ---- eval = FALSE------------------------------------------------------------
# munich_adj <- matrix(0, 25, 25)
# colnames(munich_adj) <- rownames(munich_adj) <- 1:25
# munich_adj[1, c(2, 3, 5, 12, 13)] <- 1
# munich_adj[2, c(1, 3, 5, 6, 8, 18)] <- 1
# munich_adj[3, c(1, 2, 4, 8, 9, 12)] <- 1
# munich_adj[4, c(3, 9, 11, 12)] <- 1
# munich_adj[5, c(1, 2, 13, 14, 16, 17, 18)] <- 1
# munich_adj[6, c(2, 7, 8, 18, 19)] <- 1
# munich_adj[7, c(6, 8, 19, 20, 25)] <- 1
# munich_adj[8, c(2, 3, 6, 7, 9, 25)] <- 1
# munich_adj[9, c(3, 4, 8, 10, 11, 21, 25)] <- 1
# munich_adj[10, c(9, 11, 21, 23, 24)] <- 1
# munich_adj[11, c(4, 9, 10, 12, 24)] <- 1
# munich_adj[12, c(1, 3, 4, 11, 13)] <- 1
# munich_adj[13, c(1, 5, 12, 14, 15)] <- 1
# munich_adj[14, c(5, 13, 15, 16)] <- 1
# munich_adj[15, c(13, 14, 16)] <- 1
# munich_adj[16, c(5, 14, 15, 17)] <- 1
# munich_adj[17, c(5, 16, 18)] <- 1
# munich_adj[18, c(2, 5, 6, 17, 19)] <- 1
# munich_adj[19, c(6, 7, 18, 20)] <- 1
# munich_adj[20, c(7, 19, 21, 25)] <- 1
# munich_adj[21, c(9, 10, 20, 22, 23, 25)] <- 1
# munich_adj[22, c(21, 23)] <- 1
# munich_adj[23, c(10, 21, 22, 24)] <- 1
# munich_adj[24, c(10, 11, 23)] <- 1
# munich_adj[25, c(7, 8, 9, 20, 21)] <- 1
## ---- eval = FALSE------------------------------------------------------------
# formu2 <- rentm ~ p(area, pen = "ggflasso") +
# p(year, pen = "flasso") + p(size, pen = "flasso") +
# p(rooms, pen = "flasso") + p(quality, pen = "flasso") +
# p(warm, pen = "lasso") + p(central, pen = "lasso") +
# p(tiles, pen = "lasso") + p(bathextra, pen = "lasso") +
# p(kitchen, pen = "lasso")
#
# munich.fit2 <- glmsmurf(formula = formu2, family = gaussian(), data = rent,
# pen.weights = "glm.stand", lambda = 0.048423,
# adj.matrix = list(area = munich_adj))
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.