R 使用gam模型分析浮游植物丰度和环境参数时出错(mgcv软件包)

R 使用gam模型分析浮游植物丰度和环境参数时出错(mgcv软件包),r,gam,mgcv,R,Gam,Mgcv,我对R和这个论坛都很陌生。我一直在尝试使用gam模型对浮游植物物种计数数据与环境预测值进行建模,但我遇到了错误。 我的代码如下: file <- read.csv("sg1.csv", header = TRUE, sep = ";", dec = ".", check.names = FALSE,na.strings=c("","NA")) #my dataset contains emp

我对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)

文件尝试添加行
数据。不,如果数据不包含数字级别,则$T不起作用,这里就是这种情况。如果假定
T
都是数值,则应找出向量中没有数值的原因。