GLHT對中的R

多重回歸係數我已經運行一個GLM二項式模型GLHT對中的R

fit <- glm(highlow ~ V1 + V2 + V3 + V4 + V5 + V6 + V7 + V8 + V9 + V10, 
     family="binomial")

爲了測試V1 = V2的零假設我已經使用下面的代碼。

glht.mod <- glht(fit, linfct = c("V1 - V2 = 0")) 
summary(glht.mod)

我的問題是我可以測試是否V1 = V2 = V3（所有三個係數相等的零假設 - 注意，這是不一樣的檢測無論是在單獨的迭代V1 = V2和V2 = V3）？

如果是任何幫助，我可以使用下面的代碼在SAS實現這一

proc logistic; 
    model highlow = V1 V2 V3 V4 V5 V6 V7 V8 V9 V10; 
    test1: V1 = V2; 
    test2: V1 = V2 = V3; 
run;

來源

2017-05-05 Jorvik77

縱觀SAS手冊的' test'語句中的邏輯'（https://support.sas.com/documentation/cdl/en/statug/63347/HTML/default/viewer.htm#statug_logistic_sect025.htm），你可以找到斷言「test3 ：test a1 = a2 = a3; test4：test a1 = a2，a2 = a3;「是等同的。 –

這是正確的，但在SAS手冊的test3和test4中，測試都在同一步驟中運行。這是我不知道如何在R中完成的。我能夠在Test1 V1 = V2的SAS中執行相同的操作。 test2 V2 = V3;我無法做的是Test1 V1 = V2 = V3;或者如果您更喜歡Test1 V1 = V2，V2 = V3; – Jorvik77

謝謝。看起來這是我以後的解決方案。 – Jorvik77

你的問題的不可能性的解決方案：

set.seed(1) 
n <- 1000 
highlow <- factor(runif(n)>0.5) 
X <- matrix(rnorm(n*10),nrow=n) 
df <- data.frame(highlow, X) 
names(df) <- c("highlow", paste("V",1:10,sep="")) 

fit <- glm(highlow ~ V1 + V2 + V3 + V4 + V5 + V6 + V7 + V8 + V9 + V10, 
     family="binomial", data=df) 

library(car) 
linearHypothesis(fit, c("V1-V2", "V2-V3"), c(0,0)) 


################ 
Linear hypothesis test 

Hypothesis: 
V1 - V2 = 0 
V2 - V3 = 0 

Model 1: restricted model 
Model 2: highlow ~ V1 + V2 + V3 + V4 + V5 + V6 + V7 + V8 + V9 + V10 

    Res.Df Df Chisq Pr(>Chisq) 
1 991      
2 989 2 0.2761  0.871

來源

2017-05-05 16:54:12

回答

相關問題