2015-03-25 102 views
1

我使用MatchIt根據傾向得分進行匹配,然後檢查標準化差異。我試圖手動計算一些匹配變量的標準差,並且我沒有得到與模型輸出相同的結果。MatchIT中標準化差異的公式?

#run the model 
ps<-matchit(y ~ ., data=mydata, method = "nearest", distance = "logit") 
#store the results as a data frame 
comp <- as.data.frame(summary(ps, standardize = TRUE)["sum.matched"]) 
#print the results 
comp[, "sum.matched.Std..Mean.Diff." , drop=FALSE] 

我的問題這裏:有誰知道公式手動計算是在matchIt包使用的標準化的區別是什麼?我對平均值和比例的標準化差異公式感興趣(它們是不同的嗎?也許這就是爲什麼)

+0

問一個公式真的不是一個具體的編程問題。如果您想了解更多關於統計模型或方法的信息,最好在[stats.se]處詢問。如果你想知道這個函數的功能,請查看[如何找到R函數的源代碼](http://stackoverflow.com/questions/19226816/how-can-i-view-the-source-代碼功能) – MrFlick 2015-03-25 02:25:45

+0

我很難讀取該函數的源代碼。這是我需要幫助的。將編輯以突出顯示。 – user2592729 2015-03-25 03:05:20

回答

0

我有同樣的問題,我認爲MatchIt使用的公式與最常用的公式不同。

根據Ho等的文章(http://imai.princeton.edu/research/files/matchit.pdf) 「的標準化= TRUE選項將打印出標準版本的平衡措施,其中平均差異是標準化(分)由 標準偏差在 原始治療組「。

所以包裝公式中的分母應該是s(處理),而不是s(處理後)^ 2 + s(對照)^ 2的平方根。