我使用rpy2
在python中嵌入一些R,並在python中調用lm
。我的目標是提取參數t-stat或p值。我不知道該怎麼做。一般在R我使用summary(model)$coefficients[1,4]
。我如何在python環境中調用它?提取係數p值rpy2
1
A
回答
-1
0
有一些(但不是很多)在the rpy2 docs和the Pandas docs有用的信息。
看起來,檢索R
命令的詳細結果的最佳方式如下。
我們開始與通常的進口:
import pandas as pd
from rpy2.robjects import r as R
from rpy2.robjects.packages import importr
from rpy2.robjects import pandas2ri
stats = importr('stats')
base = importr('base')
現在運行在R上的線性模型和檢索係數:
# Equivalent of lm(Sepal.Length ~ Sepal.Width, data='iris')
lm = stats.lm("Sepal.Length ~ Sepal.Width", data=R['iris'])
# Equivalent of summary(lm)
summary = base.summary(lm)
# Extract the coefficients
coeffs = summary.rx2('coefficients')
然後,我們可以創建一個從coeffs
對象熊貓數據幀:
# Build a DataFrame from the coefficients tables
df = pd.DataFrame(pandas2ri.ri2py(coeffs),
index=coeffs.names[0], columns=coeffs.names[1])
這現在使我們能夠使用係數,因爲我們將n Python中正式使用:
In [11]: df['Pr(>|t|)'] # p-values!
Out[11]:
(Intercept) 6.469702e-28
Sepal.Width 1.518983e-01
Name: Pr(>|t|), dtype: float64
In [12]: df.loc['Sepal.Width', 'Pr(>|t|)']
Out[12]: 0.15189826071144744
相關問題
- 1. 從fit.mult.impute提取係數的p值
- 2. 提取p值
- 3. R Package quantreg:提取p值
- 4. 提取回歸係數值
- 5. 相關係數的p值
- 6. 僅從邏輯模型中提取其p值有意義的係數
- 7. 提取係數
- 8. R中提取係數G值
- 9. 解釋Lasso迴歸p值與係數
- 10. lmPerm P值取決於係數的順序
- 11. Rpy2 pandas2ri.ri2py()被轉換NA值以整數
- 12. 提取係數和從LME
- 13. p-400數據庫系統?
- 14. 提取參數值
- 15. lmer從anova獲取p值
- 16. 從p獲取值:SelectOneMenu
- 17. 聯繫表單不提取值
- 18. mgcv:如何提取自適應平滑的P樣條的結,基,係數和預測值?
- 19. Bash:採取參數或讀取-p提示它
- 20. 多義串聯關係的P值
- 21. 在rpy2中按名稱提取列表元素?
- 22. 如何在Rpy2中提取因子向量中的索引
- 23. 如何從sumurca對象中提取p值?
- 24. 從nlme/lme4輸出中提取固定效果的p值
- 25. 使用COEF和summary.lm與R中robcov(提取p值)
- 26. R:獲取p值在多次線性迴歸的所有係數
- 27. 從命令提示符運行rpy2
- 28. 提取文本<p></p>與BeautifulSoup
- 29. 用Python/Scrapy在h1內提取p-1
- 30. 使用熊貓系列提取熊貓數據框的值
顯然不是。 (至少不是我能找到的。) – LondonRob 2017-09-27 10:46:03