使用gam模型分析浮游植物丰度和环境参数的错误mgcv软件包

问题描述

我是R的新手,也是这个论坛的新手。我一直在尝试使用gam模型来根据环境预测因子对浮游植物物种计数数据进行建模,但是我陷入了错误
我的代码如下:

file <- read.csv("sg1.csv",header = TRUE,sep = ";",dec = ".",check.names = FALSE,na.strings=c("","NA")) #my dataset contains empty cells that I substitute with NA
data.selected <- file[,c(5,6,14:19)] #I select only the columns on which I am interested
data.no_na <- na.omit(data.selected)
colnames(data.no_na) <- c("T","S","P","Si","DIN","DIN_P","Si_DIN","Diato")
set.seed(123)
training.samples <- data.no_na$Diato %>% createDataPartition(p =0.8,list = FALSE) #to use Diatom as outcome variable
train.data <- data.no_na[training.samples,]
test.data <- data.no.na[-training.samples,]
model <- gam(Diato ~ s(T) + s(S) + s(P) + s(Si),data = train.data)

运行代码时,出现以下错误 smooth.construct.tp.smooth.spec(object,dk $ data,dk $ knots)中的错误: NA / NaN / Inf in chiamata a funzione esterna(arg 1) Inoltre:警告消息: 1:在mean.default(xx)中:参数不是数字或逻辑:返回NA 2:在Ops.factor(xx,shift [i])中:“-”对因子没有意义

我看到只有在将T参数放在命令行中并且使用data.no_na $ TI分析值时,才会发生这种情况。no_na$ TI获取值列表和'3011级别:1.25321 10 10.001 10.0043 10.0094 10.025 10.0304 .. S'结尾。

有人可以帮助我了解发生了什么以及我做错了什么吗?先感谢您! 如果您需要更多信息,请告诉我。

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)

相关问答

Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其...
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。...
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbc...