问题描述
我是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 (将#修改为@)