使用mgcv将样条线添加到多项式logit模型

时间:2018-12-05 16:32:18

标签: r logistic-regression gam multinomial mgcv

我正在尝试训练多项式logit模型,而在此过程中,我不妨将其设置为GAM并在组合中添加样条曲线。

我尝试使用<svg width="100%" height="100%" viewBox="0 0 1000 1000" preserveAspectRatio="xMidYMin slice" x="0" y="0" tabindex="1"> <g transform="translate(100 100)"> <g transform="translate(100 100)"> <circle id="thecircle" r="50" cx="25" cy="25" fill="yellow" /> </g> </g> <svg>,但到目前为止,我仅设法产生了错误。以下是使用mgcv数据集的一些示例。

您是否知道如何使其收敛?
谢谢

iris

模型2(花瓣宽度上有1个样条线)崩溃,并显示以下警告:

library(tidyverse)
library(mgcv)
my_iris <- iris %>% 
  mutate(y =    as.integer(as.factor(Species)) -1) 
set.seed(8)

# this works (no spline)
mod1 <- gam(
  data = my_iris,
  formula = list(
    y ~ Sepal.Length + Petal.Length,
      ~ Sepal.Length + Petal.Length),
  family=multinom(K=2))

Warning message:
 In newton(lsp = lsp, X = G$X, y = G$y, Eb = G$Eb, UrS = G$UrS, L = G$L,  :
Fitting terminated with step failure - check results carefully

模型3(花瓣长度上有1个样条线)崩溃,并显示以下警告:

mod2 <- gam(
  data = my_iris,
  formula = list(
    y ~ Sepal.Length + s(Petal.Width),
      ~ Sepal.Length + Petal.Length),
  family=multinom(K=2))

Error in eigen(hess1, symmetric = TRUE) : 
infinite or missing values in 'x'

模型4(在Sepal.Length上有2个样条线)崩溃,并显示以下警告:

mod3 <- gam(
  data = my_iris,
  formula = list(
    y ~ Sepal.Length + s(Petal.Length),
      ~ Sepal.Length + Petal.Length),
  family=multinom(K=2))

Error in gam.fit5(x, y, sp, Sl = Sl, weights = weights, offset = offset,  : 
non finite values in Hessian

0 个答案:

没有答案