2015-05-15 34 views
0

相對較新的R和第一次發佈,所以如果我的問題中缺少某些東西,我很抱歉。 我使用純素的ordistep函數進行變量選擇,使用默認的「雙向」方向方法。 我使用Hellinger轉換的物種丰度數據框作爲響應變量和一個12列獨立變量數據框。 我似乎得到不同的包含變量,如果我多次運行相同的代碼(如下)。我假設有一個僞代碼生成器涉及set.seed不能約束。有沒有辦法重現一致的結果?Ordistep的複製結果

set.seed(2000) 
step.both <- ordistep(pre_met.rda, pstep=1000) 

謝謝。

回答

0

ordistep使用標準R RNG和榮譽set.seed()。我無法重現您的問題,但在調用之前將種子設置爲相同的值時,我會得到完全相同的結果。

當前版本的ordistep沒有pstep參數。雖然我們現在有了一個新的實施方案,但我們一直很榮幸R set.seed()

+0

謝謝Dr. Oksanen。我相信我的錯誤是,我認爲在劇本開始處設置種子在會話期間保持不變。上面的問題只發生在我只重複運行ordistep行時。一旦我開始設置種子,每次運行ordistep時,結果都會返回相同的結果。再次感謝! –