library(tidyverse)
rwd <- read_csv("https://denvirlab.marshall.edu/BMR617-2023/data/RenalWesternData.csv")
mTor <- filter(rwd, Protein=="mTOR")
mTor <- pivot_wider(mTor, names_from=TissueType, values_from=Expression)
mTor <- mutate(mTor, Diff=Tumor-Normal)

t.test(pull(mTor, Tumor), pull(mTor, Normal))
t.test(pull(mTor, Tumor), pull(mTor, Normal), paired=TRUE)

t.test(pull(mTor, Diff), mu=0)

mTor <- filter(rwd, Protein=="mTOR")
ggplot(mTor, aes(x=TissueType, y=Expression)) + geom_point()
ggplot(mTor, aes(x=TissueType, y=Expression)) + geom_point() + 
  geom_line(aes(group=Patient))