পটভূমি
আমি একটি মেটা-বিশ্লেষণ পরিচালনা করছি যা আগে প্রকাশিত ডেটা অন্তর্ভুক্ত। প্রায়শই, চিকিত্সার মধ্যে পার্থক্যগুলি পি-মানগুলি, কমপক্ষে উল্লেখযোগ্য পার্থক্য (এলএসডি) এবং অন্যান্য পরিসংখ্যানগুলির সাথে রিপোর্ট করা হয় তবে ভিন্নতার কোনও সরাসরি অনুমান সরবরাহ করে না।
আমি যে মডেলটি ব্যবহার করছি তার প্রসঙ্গে, বৈকল্পিকতার একটি অতিমাত্রায় ধারণা ঠিক আছে।
সমস্যা
এখানে রূপান্তরগুলির একটি তালিকা রয়েছে যেখানে (সাভিল 2003) যা আমি বিবেচনা করছি, প্রতিক্রিয়া প্রশংসিত; নীচে, আমি ধরে নিলাম যে সুতরাং এবং ভেরিয়েবলগুলি সাধারণত বিতরণ করা হয় যদি না অন্যথায় বলা হয়:এস ই = √ 1 - α / 2 = 0.975
প্রশ্নাবলী:
প্রদত্ত , , এবং চিকিত্সার অর্থ এবংএন ˉ এক্স 1 ˉ এক্স 2 এস ই = ˉ এক্স 1 - ˉ এক্স 2
প্রদত্ত এলএসডি (রোজেনবার্গদ 2004) , , , যেখানে ব্লক সংখ্যা, এবং RCBD জন্য ডিফল্ট n বি বি এন = বি এস ই = এল এস ডি
প্রদত্ত এমএসডি (ন্যূনতম উল্লেখযোগ্য পার্থক্য) (ওয়াং 2000) , , , ডিএফ =α 2 এন - 2 এস ই = এম এস ডি
একটি 95% আত্মবিশ্বাস অন্তর্বর্তী (সাভিল 2003) দেওয়া হয়েছে (গড় থেকে উচ্চ বা নিম্ন আত্মবিশ্বাসের সীমা পর্যন্ত পরিমাপ করা হয়েছে), এবংn এস ই = সি আই
টুকির এইচএসডি, , যেখানে হল 'স্টাটিজাইজড রেঞ্জ স্ট্যাটিস্টিকস',q এস E = এইচ এস ডি
এই সমীকরণগুলি encapsulate করতে একটি আর ফাংশন:
উদাহরণ ডেটা:
data <- data.frame(Y=rep(1,5), stat=rep(1,5), n=rep(4,5), statname=c('SD', 'MSE', 'LSD', 'HSD', 'MSD')
উদাহরণ ব্যবহার:
transformstats(data)
transformstats
ফাংশন:transformstats <- function(data) { ## Transformation of stats to SE ## transform SD to SE if ("SD" %in% data$statname) { sdi <- which(data$statname == "SD") data$stat[sdi] <- data$stat[sdi] / sqrt(data$n[sdi]) data$statname[sdi] <- "SE" } ## transform MSE to SE if ("MSE" %in% data$statname) { msei <- which(data$statname == "MSE") data$stat[msei] <- sqrt (data$stat[msei]/data$n[msei]) data$statname[msei] <- "SE" } ## 95%CI measured from mean to upper or lower CI ## SE = CI/t if ("95%CI" %in% data$statname) { cii <- which(data$statname == '95%CI') data$stat[cii] <- data$stat[cii]/qt(0.975,data$n[cii]) data$statname[cii] <- "SE" } ## Fisher's Least Significant Difference (LSD) ## conservatively assume no within block replication if ("LSD" %in% data$statname) { lsdi <- which(data$statname == "LSD") data$stat[lsdi] <- data$stat[lsdi] / (qt(0.975,data$n[lsdi]) * sqrt( (2 * data$n[lsdi]))) data$statname[lsdi] <- "SE" } ## Tukey's Honestly Significant Difference (HSD), ## conservatively assuming 3 groups being tested so df =2 if ("HSD" %in% data$statname) { hsdi <- which(data$statname == "HSD" & data$n > 1) data$stat[hsdi] <- data$stat[hsdi] / (qtukey(0.975, data$n[lsdi], df = 2)) data$statname[hsdi] <- "SE" } ## MSD Minimum Squared Difference ## MSD = t_{\alpha/2, 2n-2}*SD*sqrt(2/n) ## SE = MSD*n/(t*sqrt(2)) if ("MSD" %in% data$statname) { msdi <- which(data$statname == "MSD") data$stat[msdi] <- data$stat[msdi] * data$n[msdi] / (qt(0.975,2*data$n[lsdi]-2)*sqrt(2)) data$statname[msdi] <- "SE" } if (FALSE %in% c('SE','none') %in% data$statname) { print(paste(trait, ': ERROR!!! data contains untransformed statistics')) } return(data) }
তথ্যসূত্র
স্যাভিল 2003 ক্যান জে এক্সপটল সাইক। (PDF)
ওয়াং এট আল। 2000 এনভি। Tox ডাউনলোড। এবং কেম 19 (1): 113-117 (লিঙ্ক)