【R】 dataframe を群で分けてt検定

dataframe を群で分けてt検定

今回は、deep learningで作成した頭蓋骨CTデータのPSNR, SSIMを近位側、尾側に分けて比較。

データはこんな感じにずらっと並べてある。boxplotなどへそのままいける状態

group分けは"proximal"と”distal”

R t検定

R t検定

手順
subsetでgroup分け
↓
t検定

〜

dat <- read.csv("/Users/Dai_Oura/Desktop/proximal_SSIM.csv",header=TRUE)
dat <- na.omit(dat)
attach(dat)

library(dplyr)
library(psych)


#library for plot

library(ggplot2)
library(ggsignif)
library(gridExtra)
library(ggpubr)
library(rstatix)

#show plot in 1 figure
library(tidyverse)
library(patchwork)

# save figure as editable
library("mschart")

describeBy(dat, group)


proximal <-subset(dat,group=="proximal")
distal     <-subset(dat,group=="distal")

t_test_PSNR <- t.test(x=proximal$PSNR, y=distal$PSNR,paired=F)

t_test_SSIM <- t.test(x=proximal$SSIM, y=distal$SSIM,paired=F) 

t_test_NMSE <- t.test(x=proximal$NMSE, y=distal$NMSE,paired=F) 

〜

subsetでデータをわけるのがポイント！！