নিম্নলিখিত পরিস্থিতি তদন্তকারী (আই), পর্যালোচক / সম্পাদক (আর, সিআরএএন সম্পর্কিত নয়) এবং প্লট স্রষ্টা হিসাবে আমার (এম) এর ত্রয়ীর মধ্যে সর্বাধিক জিজ্ঞাসিত প্রশ্নে পরিণত হয়েছে। আমরা ধরে নিতে পারি যে (আর) টিপিকাল মেডিকেল বিগ বস রিভিউর, যিনি কেবল জানেন যে প্রতিটি প্লটের অবশ্যই ত্রুটি বার থাকা উচিত, অন্যথায় এটি ভুল। কোনও পরিসংখ্যান পর্যালোচক জড়িত থাকলে, সমস্যাগুলি খুব কম সমালোচিত হয়।
দৃশ্যপট
একটি সাধারণ ফার্মাকোলজিকাল ক্রস-ওভার গবেষণায়, দুটি ওষুধ এ এবং বি গ্লুকোজ স্তরের প্রভাবের জন্য পরীক্ষা করা হয়। প্রতিটি রোগীর দু'বার এলোমেলো ক্রমে এবং কোনও বহন-ধারণা গ্রহণের অধীনে পরীক্ষা করা হয়। প্রাথমিক এন্ডপয়েন্টটি হ'ল গ্লুকোজ (বিএ) এর মধ্যে পার্থক্য এবং আমরা ধরে নিই যে একটি জোড়াযুক্ত টি-পরীক্ষা পর্যাপ্ত।
(আমি) এমন একটি প্লট চাই যা উভয় ক্ষেত্রেই পরম গ্লুকোজ স্তর প্রদর্শন করে। তিনি ত্রুটি বারগুলির জন্য (আর) এর আকাঙ্ক্ষাকে ভয় করেন এবং বার গ্রাফগুলিতে স্ট্যান্ডার্ড ত্রুটিগুলির জন্য বলেন। আসুন এখানে বার গ্রাফ যুদ্ধ শুরু করা যাক ._)
(আমি): এটি সত্য হতে পারে না। বারগুলি ওভারল্যাপ করে, এবং আমাদের পি = 0.03 রয়েছে? আমি হাই স্কুলে যা শিখেছি তা নয়।
(এম): আমাদের এখানে একটি জোড়াযুক্ত নকশা রয়েছে। অনুরোধ করা ত্রুটি বারগুলি সম্পূর্ণ অপ্রাসঙ্গিক, জুটিযুক্ত পার্থক্যের এসই / সিআই কী গণনা করা হয়, যা প্লটে প্রদর্শিত হয় না। আমার যদি পছন্দ হয় এবং খুব বেশি ডেটা না থাকে তবে আমি নীচের প্লটটি পছন্দ করব
যোগ করা হয়েছে 1: এটি বেশ কয়েকটি প্রতিক্রিয়াতে উল্লিখিত সমান্তরাল স্থানাঙ্ক প্লট
(এম): লাইনগুলি জুটিবদ্ধ হওয়া দেখায় এবং বেশিরভাগ লাইন উঠে যায় এবং এটি সঠিক ধারণা, কারণ opeাল যা গণনা করে (ঠিক আছে, এটি শ্রেণিবদ্ধ, তবে তবুও)।
(আমি): এই চিত্রটি বিভ্রান্তিকর। এটি কেউ বুঝতে পারে না এবং এর কোনও ত্রুটি বার নেই (আর লুকিয়ে আছে)।
(এম): আমরা আরও একটি প্লট যুক্ত করতে পারি যা পার্থক্যটির প্রাসঙ্গিক আত্মবিশ্বাসের ব্যবধানটি দেখায়। শূন্য-রেখা থেকে দূরত্ব প্রভাবের আকারের একটি ছাপ দেয়।
(আমি): কেউ এটা করে না
(আর): এবং এটি মূল্যবান গাছগুলি অপচয় করে
(এম): (একটি ভাল জার্মান হিসাবে): হ্যাঁ, গাছগুলিতে পয়েন্ট নেওয়া হয়েছে। তবে তবুও আমি এটি ব্যবহার করি (এবং এটি কখনই প্রকাশিত হবে না) যখন আমাদের একাধিক চিকিত্সা এবং একাধিক বৈপরীত্য হয়।
কোন পরামর্শ ? আর-কোডটি নীচে রয়েছে, যদি আপনি কোনও প্লট তৈরি করতে চান।
# Graphics for Crossover experiments
library(ggplot2)
library(plyr)
theme_set(theme_bw()+theme(panel.margin=grid::unit(0,"lines")))
n = 20
effect = 5
set.seed(4711)
glu0 = rnorm(n,120,30)
glu1 = glu0 + rnorm(n,effect,7)
dt = data.frame(patient = rep(paste0("P",10:(9+n))),
treatment = rep(c("A","B"), each=n),glucose = c(glu0,glu1))
dt1 = ddply(dt,.(treatment), function(x){
data.frame(glucose = mean(x$glucose), se = sqrt(var(x$glucose)/nrow(x)) )})
tt = t.test(glucose~treatment,paired=TRUE,data=dt,conf.int=TRUE)
dt2 = data.frame(diff = -tt$estimate,low=-tt$conf.int[2], up=-tt$conf.int[1])
p = paste("p =",signif(tt$p.value,2))
png(height=300,width=300)
ggplot(dt1, aes(x=treatment, y=glucose, fill=treatment))+
geom_bar(stat="identity")+
geom_errorbar(aes(ymin=glucose-se, ymax=glucose+se),size=1., width=0.3)+
geom_text(aes(1.5,150),label=p,size=6)
ggplot(dt,aes(x=treatment,y=glucose, group=patient))+ylim(0,190)+
geom_line()+geom_point(size=4.5)+
geom_text(aes(1.5,60),label=p,size=6)
ggplot(dt2,aes(x="",y=diff))+
geom_errorbar(aes(ymin=low,ymax=up),size=1.5,width=0.2)+
geom_text(aes(1,-0.8),label=p,size=6)+
ylab("95% CI of difference glucose B-A")+ ylim(-10,10)+
theme(panel.border=element_blank(), panel.grid.major.x=element_blank(),
panel.grid.major.y=element_line(size=1,colour="grey88"))
dev.off()