1

考虑一个最小的工作示例（例如二项式模型）：

test.a.tset <- rnorm(10)
test.b.tset <- rnorm(10)
c <- runif(10)
c[c < 0.5] <- 0
c[c >= 0.5] <- 1
df <- data.frame(test.a.tset,test.b.tset,c)

使用正则表达式，我想c使用结构回归所有变量test."anything".tset：

summary(glm(paste("c ~ ",paste(colnames((df[, grep("test\\.\\w+\\.tset", colnames(df))])),
        collapse = "+"), sep = ""), data = df, family=binomial))

到目前为止，没有任何问题。现在我们cbind进入发挥作用的部分。假设我想使用不同的统计模型（例如rbprobitGibbs来自bayesm包），它需要一个设计矩阵作为输入。因此，我需要将数据框转换为适当的格式。

X <- cbind(df$test.a.tset,df$test.b.tset)

或者，如果我想再次使用正则表达式（我什至添加第二个grep以确保仅选择引号内的部分）：

X2 <- cbind(grep("[^\"]+",paste(paste("df$", colnames((df[, grep("test\\.\\w+\\.tset", colnames(df))])), 
            sep = ""), collapse = ","), value = TRUE))

但是有一个区别：

> X
            [,1]         [,2]
 [1,] -0.4525601 -1.240484170
 [2,]  0.3135625  1.240519383
 [3,] -0.2883953 -0.554670224
 [4,] -1.3696994 -1.373690426
 [5,]  0.8514529 -0.063945537
 [6,] -1.1804205 -0.314132743
 [7,] -1.0161170 -0.001605679
 [8,]  1.0072168  0.938921869
 [9,] -0.8797069 -1.158626865
[10,] -0.9113297  1.641201924
> X2
     [,1]                        
[1,] "df$test.a.tset,df$test.b.tset"

从我的角度来看，问题似乎在于grep将选定的值作为引号内的字符串返回，而在某种程度上glm忽略了引号中的引号"df$test.a.tset,test.b.tset"，cbind但没有。即粘贴后对 X2 的调用实际上读作：

X2 <- cbind("df$test.a.tset,df$test.b.tset")

问题：有没有办法获得与使用正则表达式相同的X2结果X？

4

1 回答 1

1

该代码grep("test\\.\\w+\\.tset", colnames(df))将返回与您的模式匹配的列的索引。如果您想仅使用这些列构建矩阵，则可以使用：

X3 <- as.matrix(df[,grep("test\\.\\w+\\.tset", colnames(df))])
于 2014-03-25T13:58:06.637 回答