2014-04-18 42 views
0

我有一個關於Smacof包中的smacofSym函數的問題。我通過RStudio版本0.98.501使用R版本3.1.0。smacofSym的困難 - 多維縮放

我使用下面的命令:以下數據

MDSdata <- smacofSym(DJaccardMatrix, ndim=2, metric=FALSE, verbose=TRUE). 

我已經包含細節我使用(DJaccardMatrix)。每次運行smacofSym時,我都會得到一個最終配置正確的配置。下面是結果的樣本:

MDSdata$conf 
     D1   D2 
1 0.06259624 -0.01494732 
2 0.06276541 -0.01480409 
3 0.06266933 -0.01492375 
4 0.06262438 -0.01496111 
5 0.06243336 -0.01496193 
6 0.06258047 -0.01502270 
7 0.06247747 -0.01500037 ....... 

檢查結果我跑XLSTAT相同的矩陣,並得到我所期待的,一個更分佈的一系列點。在看過其他一些幫助請求之後,我嘗試將smacofSym同時作爲矩陣和dist,但都沒有影響結果。

這是我在DJaccardMatrix信息爲矩陣:

num [1:121, 1:121] 0 0.969 0.679 0.704 0.939 ... 
attr(*, "dimnames")=List of 2 
..$ : chr [1:121] "1" "2" "3" "4" ... 
..$ : chr [1:121] "1" "2" "3" "4" ... 

這是我在DJaccardMatrix信息作爲DIST對象:

Class 'dist' atomic [1:7260] 0.969 0.679 0.704 0.939 0.8 ... 
    ..- attr(*, "Size")= int 121 
    ..- attr(*, "call")= language as.dist.default(m = dissmat) 
    ..- attr(*, "Diag")= logi FALSE 
    ..- attr(*, "Upper")= logi FALSE 

我感謝的人有任何建議。我認爲這是非常基本的東西,但我絕對沒有找到它。 (在附註中 - 可以自由地忽略它,因爲它涉及的是解釋 - smacof報告和克魯斯卡爾壓力之間的非對稱應力之間的關係是什麼?)

回答

0

這個答案涉及到括號內的括號問題結束語:「smacof報告與Kruskal應力之間的非應力應力之間的關係是什麼」

克魯斯卡應力(或應力-1)是由smacof報告的非應力應力的平方根(stress.nm)。

所以,如果你有運行smacofSym得到一個名爲mod模型:

Stress-1 = mod$stress.nm^.5