擬合R中的重複測量

Question

擬合R中的重複測量，收斂問題。我有以下適合這是許多數據集之一，它不會收斂。我做了其他的套。這個數據集和模型在SAS工作...我可以得到一些方向在做什麼來做這項工作在R？看的東西（矩陣，選項設置，關於r/splus的這個主題的參考...）。

fit.gls <- gls(resp~grpnum+grpnum/day,data=long, corr=cormat,na.action=na.omit) 
Error in glsEstimate(object, control = control) : 
    computed "gls" fit is singular, rank 62 

我已閱讀以下內容，仍然試圖通它的工作... Converting Repeated Measures mixed model formula from SAS to R

Answer 1

的問題是數據。 gls需要反轉矩陣才能正常工作（有關估計協變量的公式，請參閱Wikipedia）。對於您特定的數據集，該矩陣是不可逆的。

可以允許奇異值與control參數被允許：

fit.gls <- gls(resp~grpnum+grpnum/day,data=long, corr=cormat,na.action=na.omit, control = list(singular.ok = TRUE))

與此小心，因爲你可能會得到不好的結果！隨後檢查型號。

關於選項的更多詳細信息，請參閱gls和glsConrol的幫助。