আর-তে লেভেন টেস্ট ফাংশনটি কীভাবে ব্যবহার করবেন?

Question

আমি পরিসংখ্যানের এক নবাগত এবং আর লেভেন ফাংশনটি ব্যবহার করতে আমার সমস্যা হয় (আমি দুটি নমুনার বৈচিত্রের সাম্যতা পরীক্ষা করতে চাই)। ডকুমেন্টেশন বলে যে আমার চালানো উচিত:

levene.test (y, গ্রুপ)

তবে y এবং গ্রুপ হিসাবে আমার কী রাখা উচিত তা আমার কোনও ধারণা নেই? আমার দুটি পৃথক নমুনা রয়েছে যেগুলির মধ্যে আমি ভিন্নতার সমতা পরীক্ষা করতে চাই। আমি কি নমুনার মানগুলির এক হিসাবে y এবং দ্বিতীয়টি প্যারামিটার হিসাবে রেখে দেব?

কোন ইঙ্গিত?

Answer 1

ধরা যাক যে, আর-তে, আপনার 1 ম নমুনা নামের একটি ভেক্টর sample1এবং আপনার 2 য় নমুনা নামের একটি ভেক্টরে সঞ্চিত রয়েছে sample2।

আপনাকে প্রথমে আপনার দুটি নমুনা একক ভেক্টরে একত্রিত করতে হবে এবং দুটি গ্রুপ সংজ্ঞায়িত করে অন্য ভেক্টর তৈরি করতে হবে:

y <- c(sample1, sample2)

এবং

group <- as.factor(c(rep(1, length(sample1)), rep(2, length(sample2))))

এখন, আপনি কল করতে পারেন

library(car)
levene.test(y, group)

সম্পাদনা

আর এ চেষ্টা করার সময় আমি নিম্নলিখিত সতর্কতা পেয়েছিলাম:

'levene.test' has now been removed. Use 'leveneTest' instead...

এটি অনুসারে, আপনার leveneTestপরিবর্তে এর দিকে নজর দেওয়া উচিত ...

Comments

ধন্যবাদ! তবে আপনি কি দয়া করে দয়া করে ব্যাখ্যা করবেন যে কেন এই পথে চলতে হবে? আমি এটি বুঝতে পছন্দ করব যাতে পরের বার আমাকে প্রশ্ন জিজ্ঞাসা করতে না হয় এবং অন্যকে সাহায্য করতে পারি।

— জাকুব

@ জাকুব: আচ্ছা, এটি এই পথে চলে কারণ এটি সেই কাঠামোটি ব্যবহার করে কার্যকর করা হয়েছিল। সাহায্যে বলা হয়েছে যে প্রথম যুক্তিটি প্রতিক্রিয়া পরিবর্তনশীল হতে হবে যখন দ্বিতীয় যুক্তিকে গ্রুপ ভেরিয়েবল হতে হবে।

— ocram

বেশিরভাগ ক্ষেত্রে আর এই ধরণের ডেটা ফর্ম্যাটকে পছন্দ করে বলে মনে হয়, প্রায়শই "দীর্ঘ" হিসাবে পরিচিত। পুনরায় আকার প্যাকেজটি গলানো এবং castালাই নামক ফাংশন সরবরাহ করে যা আপনার ডেটা পুনরায় আকার দেওয়ার জন্য ব্যবহার করা যেতে পারে তবে সাধারণ দুটি ভেরিয়েবল কেসের জন্য আপনার যা প্রয়োজন তা তুলনায় এগুলি আরও জটিল।

— রাসেলপিয়ের্স

শুধু নিশ্চিত করতে, এটি নমুনা 1 এবং নমুনা 2 এর ফ্রিকোয়েন্সি বর্ণালী পরীক্ষা করবে না, সঠিক? সুতরাং, উদাহরণস্বরূপ, নমুনা 1 হ'ল: 1,0,2,1,0 এবং নমুনা 2 হ'ল: 1,1,3,0,0। নমুনা 1 এর বিতরণ তৈরি করতে এটি 1 ও 1 স্যাম্পলটির 1s বিনটি করবে না, সঠিক? আমি আশা করি যে আমার ফলো-আপ প্রশ্নটি বোঝা যায়?

— অ্যাটিকাস 29

Answer 2

ওকরামের উত্তরে সমস্ত গুরুত্বপূর্ণ অংশ রয়েছে। তবে, আপনি না চাইলে আপনাকে সমস্ত আরসিএমডিআর লোড করার দরকার নেই। সম্পর্কিত গ্রন্থাগারটি "গাড়ি" car তবে যেমন ওকরাম ইঙ্গিত করে, levene.test হ্রাস করা হয়েছে। দ্রষ্টব্য যে অবচয় হ'ল কার্যকারিতা বা কোডের পরিবর্তন নয় (এই মুহুর্তে, 09/18/2011)। এটি কেবল ফাংশনটির নাম পরিবর্তন। সুতরাং লেভেন.েস্ট এবং লেভেনটেষ্ট একই কাজ করবে। রেকর্ডের জন্য আমি ভেবেছিলাম যে আমি এই সাধারণ ক্ষেত্রে লেভিনিস্টেস্ট এবং পুনরায় ব্যবহারযোগ্য পুনরায় আকারের কোড ব্যবহার করে একটি উদাহরণ সরবরাহ করব:

#Creating example code
sample1 <- rnorm(20)
sample2 <- rnorm(20)

#General code to reshape two vectors into a long data.frame
twoVarWideToLong <- function(sample1,sample2) {
    res <- data.frame(
        GroupID=as.factor(c(rep(1, length(sample1)), rep(2, length(sample2)))),
        DV=c(sample1, sample2)
    )   
}   

#Reshaping the example data
long.data <- twoVarWideToLong(sample1,sample2)

#There are many different calls here that will work... but here is an example
leveneTest(DV~GroupID,long.data)

Answer 3

ডেটা প্রস্তুত করার সহজতম উপায় (আমার মতে) পুনরায় আকার 2 প্যাকেজটি ব্যবহার করা হচ্ছে:

#Load packages
library(reshape2)
library(car)

#Creating example data
sample1 <- rnorm(20)
sample2 <- rnorm(20)

#Combine data
sample <- as.data.frame(cbind(sample1, sample2))

#Melt data
dataset <- melt(sample)

#Compute test
leveneTest(value ~ variable, dataset)