Maxichamps 1/43
940015202 AUDI 100 1990 SILVER METALLIC
940030110 MERCEDES-BENZ C-CLASS T-MODEL (S203) 2001 GREY METALLIC
940030111 MERCEDES-BENZ C-CLASS T-MODEL (S203) 2001 RED
940036200 MERCEDES-BENZ S-CLASS (W220) 1998 DARK RED METALLIC
940036201 MERCEDES-BENZ S-CLASS (W220) 1998 SILVER METALLIC
940103320 LAMBORGHINI URRACO 1974 GREEN
940103321 LAMBORGHINI URRACO 1974 RED METALLIC
940161500 HONDA CIVIC 1990 WHITE
940161501 HONDA CIVIC 1990 BLACK
940171710 VOLVO 740 BREAK 1986 WHITE
940171711 VOLVO 740 BREAK 1986 GOLD
940015202 AUDI 100 1990 SILVER METALLIC
940030110 MERCEDES-BENZ C-CLASS T-MODEL (S203) 2001 GREY METALLIC
940030111 MERCEDES-BENZ C-CLASS T-MODEL (S203) 2001 RED
940036200 MERCEDES-BENZ S-CLASS (W220) 1998 DARK RED METALLIC
940036201 MERCEDES-BENZ S-CLASS (W220) 1998 SILVER METALLIC
940103320 LAMBORGHINI URRACO 1974 GREEN
940103321 LAMBORGHINI URRACO 1974 RED METALLIC
940161500 HONDA CIVIC 1990 WHITE
940161501 HONDA CIVIC 1990 BLACK
940171710 VOLVO 740 BREAK 1986 WHITE
940171711 VOLVO 740 BREAK 1986 GOLD
#logistic regression mode##r##用图像表示确定y估计x
mod.fit<-glm(formula = kill/total ~ picloram, weight = total, family =
binomial(link = logit), data = data4)
summary(mod.fit)
plot(x = data4$picloram, y = data4$kill/data4$total, xlab = "Picloram", ylab =
"Estimated kill rate", col = rep(x = c("black", "blue", "orange"), each = 4),
pch = rep(x = 1:3, each = 4))
curve(expr = plogis(q = mod.fit$coefficients[1] + mod.fit$coefficients[2]*x), add
= TRUE, col = "red", xlim = c(min(data4$picloram), max(data4$picloram)))
legend(x = 3, y = 0.4, legend = c("Rep 1", "Rep 2", "Rep 3"), col = c("black", "blue", "orange"),
pch = c(1,2,3), bty = "n")
pi0<-0.9
LD.x<-(log(pi0/(1-pi0)) - mod.fit$coefficients[1])/mod.fit$coefficients[2]
as.numeric(LD.x)
segments(x0 = -1, y0 = pi0, x1 = LD.x, y1 = pi0, lty = "dashed")
segments(x0 = LD.x, y0 = -1, x1 = LD.x, y1 = pi0, lty = "dashed")
mod.fit<-glm(formula = kill/total ~ picloram, weight = total, family =
binomial(link = logit), data = data4)
summary(mod.fit)
plot(x = data4$picloram, y = data4$kill/data4$total, xlab = "Picloram", ylab =
"Estimated kill rate", col = rep(x = c("black", "blue", "orange"), each = 4),
pch = rep(x = 1:3, each = 4))
curve(expr = plogis(q = mod.fit$coefficients[1] + mod.fit$coefficients[2]*x), add
= TRUE, col = "red", xlim = c(min(data4$picloram), max(data4$picloram)))
legend(x = 3, y = 0.4, legend = c("Rep 1", "Rep 2", "Rep 3"), col = c("black", "blue", "orange"),
pch = c(1,2,3), bty = "n")
pi0<-0.9
LD.x<-(log(pi0/(1-pi0)) - mod.fit$coefficients[1])/mod.fit$coefficients[2]
as.numeric(LD.x)
segments(x0 = -1, y0 = pi0, x1 = LD.x, y1 = pi0, lty = "dashed")
segments(x0 = LD.x, y0 = -1, x1 = LD.x, y1 = pi0, lty = "dashed")
#顶级东方神颜##red velvet唯一官方c位裴珠泫# 如果爱有回应那一定是节气&裴珠泫浪漫天才小兔带着节气送的花在舞台上谱写出属于你我的爱意之章是偷偷去打卡粉丝应援的害羞株宝是在凌晨还记录着粉丝给的意见的真挚小泫上帝说世界上要有善神才能有物质守恒所以创造了最好的姐姐Always Stay With Irene
✋热门推荐