R中logistic回归的误差

R中logistic回归的误差,r,statistics,logistic-regression,R,Statistics,Logistic Regression,有人能告诉我哪里出了问题吗。我试图在R中执行逻辑回归。为了测试代码,我有一个小的数据框,一行和所有列。我运行了以下代码: Cs <- c(A[,4],A[,5],A[,6],A[,7],A[,8],A[,9],A[,10],A[,11],A[,12],A[,13],A[,14],A[,15],A[,16],A[,17],A[,18]) Ts <- c(A[,19],A[,20],A[,21],A[,22],A[,23],A[,24],A[,25],A[,26],A[,27],A[,2

有人能告诉我哪里出了问题吗。我试图在R中执行逻辑回归。为了测试代码,我有一个小的数据框,一行和所有列。我运行了以下代码:

Cs <- c(A[,4],A[,5],A[,6],A[,7],A[,8],A[,9],A[,10],A[,11],A[,12],A[,13],A[,14],A[,15],A[,16],A[,17],A[,18])
Ts <- c(A[,19],A[,20],A[,21],A[,22],A[,23],A[,24],A[,25],A[,26],A[,27],A[,28],A[,29],A[,30],A[,31],A[,32],A[,33])
regions.X <- cbind(Cs, Ts)
regions.IV <- c('s','s','s','s','s','s','s','c','c','c','c','c','c','c','c')
result <- glm(regions.X ~ regions.IV, family = 'binomial')
summary(result)$coefficients[2,4]

Cs在运行模型之前,您是否检查了变量的结构?或者只是期望它能工作?我不知道这个结构会影响相同的代码在apply中用作函数?在任何情况下,我的测试数据帧的结构如下:>str(a)'data.frame':1 obs。共33个变量:$chr:系数w/21级“chr1”、“chr10”,..:1$start:int 100027723$end:int 100028202$numCs1:num 32$numCs2:num 38$numCs3:num 30$numCs4:num 29$numCs5:num 43$numCs6:num 41$numCs7:num 34$numCs8:num 27$numCs9:num 29$numCs11:num 31$numCs12:num 30$numCs13:num 53$numCs14:num 52$numCs15:num 54$numTs1:num 10$numTs2:num 16$numTs3:num 17$numTs4:num 9。。。。等等,因为您正在执行此
apply(mtcars,1,函数(x)x[1:5])
vs此
mtcars[,1:5]
在第一段代码中。因此,请再次访问我的第一条评论,这次不要把它吹掉。我很感谢你在这方面的意见,并不是贬低你的评论。你能澄清一下你认为代码有什么问题吗?如果我弄错了,请纠正我。据我所知,当apply与MARGIN=1一起使用时(例如apply(mtcars,1,函数(x)x[1:5]),MARGIN 1告诉函数我正在迭代我的行,所以我只需要指定列索引。相反,如果我没有使用apply函数,我需要同时指定行和列索引,因此mtcars[,1:5]。
region.log <- apply(A, MARGIN = 1, function(m){Cs <- c(m[4],m[5],m[6],m[7],m[8],m[9],m[10],m[11],m[12],m[13],m[14],m[15],m[16],m[17],m[18])
Ts <- c(m[19],m[20],m[21],m[22],m[23],m[24],m[25],m[26],m[27],m[28],m[29],m[30],m[31], m[32],m[33])
regions.X <- cbind(Cs, Ts)
regions.IV <- c('s','s','s','s','s','s','s','c','c','c','c','c','c','c','c')
result <- glm(regions.X ~ regions.IV, family = 'binomial')
summary(result)$coefficients[2,4]})
Error in `[[<-.data.frame`(`*tmp*`, i, value = c(13L, 15L, 11L, 10L, 17L,  : 
replacement has 30 rows, data has 15