R語(yǔ)言RcppEigen計(jì)算點(diǎn)乘與矩陣乘法連乘算法錯(cuò)誤解決
計(jì)算點(diǎn)乘與矩陣乘法連乘計(jì)算錯(cuò)誤
當(dāng)我們想將 R 中的連乘(如下公式所示)修改成 Rcpp 代碼時(shí),
t(X)^2 %*% X
理論上我們只用在 .cpp 代碼中輸入下述語(yǔ)句即可(默認(rèn)使用了 RcppEigen 庫(kù)):
X.adjoint().array().square() * X.array().square();
但實(shí)際上這樣會(huì)會(huì)出現(xiàn)問(wèn)題,原因是 X.adjoint().array().square() 與 X.array().square() 沒(méi)有成功轉(zhuǎn)化成 Eigen::MatrixXd 形式。雖然單獨(dú)顯示兩部分都沒(méi)問(wèn)題(可以直接輸出看到),但是連乘的時(shí)候就會(huì)出錯(cuò)。這時(shí)我們分開(kāi)進(jìn)行一個(gè)默認(rèn)轉(zhuǎn)換,然后再進(jìn)行矩陣乘法即可,如下所示:
Eigen::MatrixXd A, B; A = X.adjoint().array().square(); B = X.array().square(); A * B;
這樣就可以解決之前的計(jì)算錯(cuò)誤問(wèn)題了。
以上就是RcppEigen計(jì)算點(diǎn)乘與矩陣乘法連乘計(jì)算錯(cuò)誤解決的詳細(xì)內(nèi)容,更多關(guān)于RcppEigen計(jì)算點(diǎn)乘與矩陣乘法的資料請(qǐng)關(guān)注腳本之家其它相關(guān)文章!
相關(guān)文章
R語(yǔ)言繪圖數(shù)據(jù)可視化Ridgeline plot山脊圖畫(huà)法
這篇文章主要為大家介紹了R語(yǔ)言繪圖數(shù)據(jù)可視化Ridgeline plot山脊圖畫(huà)法的示例詳解,有需要的朋友可以借鑒參考下,希望能夠有所幫助2022-02-02
R語(yǔ)言-如何循環(huán)讀取excel并保存為RData
這篇文章主要介紹了R語(yǔ)言循環(huán)讀取excel并保存為RData的操作,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2021-04-04
R語(yǔ)言 實(shí)現(xiàn)將factor轉(zhuǎn)換成numeric方法
這篇文章主要介紹了R語(yǔ)言 實(shí)現(xiàn)將factor轉(zhuǎn)換成numeric方法,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2021-03-03
R語(yǔ)言 data.frame()命令的使用說(shuō)明
這篇文章主要介紹了R語(yǔ)言 data.frame()命令的使用說(shuō)明,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2021-04-04
R語(yǔ)言ComplexHeatmap繪制復(fù)雜熱圖heatmap
ComplexHeatmap可以繪制很復(fù)雜的熱圖,能滿足日常以及文章所需,本文介紹一下文章常見(jiàn)的復(fù)雜熱圖的繪制方式,含代碼,文章比較長(zhǎng),建議耐心看一下,也許就有你需要的熱圖2022-06-06

