誰會在 %T>% 上給我建議,
代碼是
c(1031, 8, 0, 968, 66, 12) %>%
matrix(ncol = 3) %>%
'colnames<-'(c("Improved", "Hospitalized", "Death")) %>%
'rownames<-'(c("Treated", "Placebo")) %>%
as.table %T>% fisher.test %>% chisq.test
然而,它總是繞過fisher.test,只進行chisq.test。我的計劃是將資料分別輸出到 Fisher 和 chisq 測驗中。
非常感謝您的建議!
查爾斯
uj5u.com熱心網友回復:
你不需要三通管。
library(magrittr)
c(1031, 8, 0, 968, 66, 12) %>%
matrix(ncol = 3) %>%
'colnames<-'(c("Improved", "Hospitalized", "Death")) %>%
'rownames<-'(c("Treated", "Placebo")) %>%
as.table %>%
{fisher.test(.) %>% print()
chisq.test(.) %>% print()}
#>
#> Fisher's Exact Test for Count Data
#>
#> data: .
#> p-value < 2.2e-16
#> alternative hypothesis: two.sided
#>
#>
#> Pearson's Chi-squared test
#>
#> data: .
#> X-squared = 2012.4, df = 2, p-value < 2.2e-16
由reprex 包于 2022-01-22 創建(v2.0.1)
uj5u.com熱心網友回復:
問題是雖然資料在您的測驗后回傳,但您的測驗結果并沒有列印出來。這可以通過定義一個Print強制列印每個步驟的功能來克服。然后,您可以%T>%根據需要進行多次。Print這樣做的一個好處是,您還可以通過剝離包裝“關閉”那些您不想列印到控制臺的步驟。
library(magrittr)
Print <- function(f) function(x) print(f(x))
c(1031, 8, 0, 968, 66, 12) %>%
matrix(ncol = 3) %>%
'colnames<-'(c("Improved", "Hospitalized", "Death")) %>%
'rownames<-'(c("Treated", "Placebo")) %>%
as.table %T>% Print(fisher.test)() %T>% Print(chisq.test)() %T>% Print(t.test)()
請注意,()后面Print(...)是強制性的。
輸出
Fisher's Exact Test for Count Data
data: x
p-value < 2.2e-16
alternative hypothesis: two.sided
Pearson's Chi-squared test
data: x
X-squared = 2012.4, df = 2, p-value < 2.2e-16
One Sample t-test
data: x
t = 1.6823, df = 5, p-value = 0.1533
alternative hypothesis: true mean is not equal to 0
95 percent confidence interval:
-183.4795 878.4795
sample estimates:
mean of x
347.5
Improved Hospitalized Death
Treated 1031 0 66
Placebo 8 968 12
轉載請註明出處,本文鏈接:https://www.uj5u.com/gongcheng/419369.html
標籤:
