2011-11-28 21 views
7

Tôi chắc chắn đây là lệnh dễ trong R, nhưng vì một lý do nào đó, tôi gặp khó khăn khi tìm kiếm dung dịch.Kiểm tra ý nghĩa trong R, xác định xem tỷ lệ trong một cột có khác biệt đáng kể so với cột khác trong biến duy nhất

Tôi đang cố gắng chạy một nhóm crosstabs (sử dụng lệnh table()) trong R và mỗi tab có hai cột (điều trị và không điều trị). Tôi muốn biết liệu sự khác biệt giữa các cột có khác nhau đáng kể đối với tất cả các hàng hay không (các hàng là một số lựa chọn câu trả lời từ một khảo sát). Tôi không quan tâm đến ý nghĩa tổng thể, chỉ trong điều trị so sánh chéo so với không điều trị.

Loại phân tích này rất dễ dàng trong SPSS (liên kết dưới đây để minh họa những gì tôi đang nói), nhưng dường như tôi không thể làm việc này trong R. Bạn có biết tôi có thể làm điều này không?

http://help.vovici.net/robohelp/robohelp/server/general/projects_fhpro/survey_workbench_MX/Significance_testing.htm

EDITED: Dưới đây là một ví dụ về R về những gì tôi muốn nói:

treatmentVar <-c(0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1) # treatment is 1 or 0 
question1 <-c(1,2,2,3,1,1,2,2,3,1,1,2,2,3,1,3) #choices available are 1, 2, or 3 
Questiontab <- table(question1, treatmentVar) 
Questiontab 

Tôi có bảng như thế này^(percentaged theo cột trên treatmentVar), và tôi muốn để xem liệu có sự khác biệt đáng kể giữa mỗi lựa chọn câu hỏi (hàng) từ điều trị 0 đến điều trị hay không 1. Vì vậy, trong ví dụ trên, tôi muốn biết liệu có sự khác biệt đáng kể giữa 4 và 2 (hàng 1), 3 và 3 (hàng 2) và 1 và 3 (hàng 3). Vì vậy, trong ví dụ này, các lựa chọn cho câu hỏi 1 có thể khác biệt đáng kể cho các lựa chọn 1 và 3 (vì sự khác biệt là 2) nhưng sự khác biệt cho lựa chọn 2 không phải vì sự khác biệt là 0. Cuối cùng, tôi đang cố gắng xác định loại ý nghĩa này. Tôi hy vọng rằng sẽ giúp.

Cảm ơn!

+0

Sẽ dễ dàng hơn nếu mọi người trả lời câu hỏi của bạn nếu bạn cung cấp ví dụ có thể tái sản xuất, bao gồm cả tập dữ liệu mẫu nhỏ. – joran

+0

xin lỗi về điều đó. Tôi chỉ cần chỉnh sửa câu hỏi của tôi với một ví dụ R –

Trả lời

5

Sử dụng ví dụ của bạn, hoặc là chisq.test hoặc prop.test (tương đương trong trường hợp này):

> chisq.test(Questiontab) 

     Pearson's Chi-squared test 

data: Questiontab 
X-squared = 1.6667, df = 2, p-value = 0.4346 

Warning message: 
In chisq.test(Questiontab) : Chi-squared approximation may be incorrect 
> prop.test(Questiontab) 

     3-sample test for equality of proportions without continuity 
     correction 

data: Questiontab 
X-squared = 1.6667, df = 2, p-value = 0.4346 
alternative hypothesis: two.sided 
sample estimates: 
    prop 1 prop 2 prop 3 
0.6666667 0.5000000 0.2500000 

Warning message: 
In prop.test(Questiontab) : Chi-squared approximation may be incorrect 

Lưu ý các cảnh báo; các xét nghiệm này không nhất thiết phải phù hợp với những con số nhỏ như vậy.

9

Tôi nghĩ rằng chức năng bạn đang tìm kiếm là pairwise.prop.test(). Xem ?pairwise.prop.test để biết ví dụ.

+1

+10, nếu tôi có thể. – Andrie

+0

Cảm ơn John, nhưng đó không phải là những gì tôi yêu cầu. Lỗi của tôi, tôi nên rõ ràng hơn trong câu hỏi của tôi. Tôi đã cập nhật câu hỏi của mình với một ví dụ R. Tôi hy vọng rằng sẽ giúp làm rõ. –

Các vấn đề liên quan