[R-br] teste wilcoxon
alanarocha em sapo.pt
alanarocha em sapo.pt
Terça Junho 18 13:29:28 BRT 2013
Olá a todos,
executei a base de dados original:
DoseT-9 DoseT-6 DoseT-3 DoseT+3 DoseT+6 DoseT+9
1 54020104 930,0 900,0 900,0 390,0 1.170,0
2 54020178 1.860,0 930,0 900,0 420,0 840,0
3 54020303 1.860,0 2.790,0 2.700,0 1.620,0 1.440,0
4 54020492 930,0 930,0 900,0 420,0 840,0
5 54020637 1.080,0 480,0 870,0 2.970,0 2.700,0
6 54020698 2.790,0 1.590,0 930,0 390,0 720,0
NA
1 2.700,0
2 630,0
3 1.530,0
4 720,0
5 2.520,0
6 2.340,0
depois tentei executar só duas colunas da base de dados
onde quero aplicar o teste de wilcoxon para elas o resultado não entendo:
> setwd("f:/Ana_trabalho/Cinacalset")
> dados<-read.table("Cinacalset.csv",header=TRUE,sep=";",dec=",")
> colnames(dados)<-c("DoseT-9", "DoseT-6", "DoseT-3", "DoseT+3",
> "DoseT+6", "DoseT+9")
> dados[,-c(1,2,5,6,7)]
DoseT-3 DoseT+3
1 900,0 900,0
2 930,0 900,0
3 2.790,0 2.700,0
4 930,0 900,0
5 480,0 870,0
6 1.590,0 930,0
7 930,0 900,0
8 750,0 780,0
9 3.720,0 3.600,0
10 420,0 60,0
11 2.790,0 2.790,0
12 930,0 900,0
13 3.720,0 3.720,0
14 930,0 900,0
15 930,0 1.800,0
16 900,0 930,0
17 1.860,0 1.860,0
18 1.800,0 1.800,0
19 390,0 270,0
20 750,0 780,0
21 930,0 900,0
22 930,0 930,0
23 540,0 450,0
24 390,0 390,0
25 1.860,0 1.860,0
26 2.340,0 2.340,0
27 210,0 390,0
28 390,0 930,0
29 1.440,0 1.350,0
30 930,0 900,0
31 1.410,0 1.560,0
32 840,0 930,0
33 900,0 930,0
34 1.860,0 1.710,0
35 900,0 600,0
36 900,0 1.860,0
37 900,0 930,0
38 900,0 930,0
39 900,0 930,0
40 900,0 930,0
41 390,0 390,0
42 900,0 930,0
43 930,0 900,0
44 2.700,0 2.790,0
45 900,0 930,0
46 900,0 930,0
47 900,0 930,0
48 900,0 930,0
49 900,0 1.860,0
50 1.800,0 930,0
51 1.800,0 930,0
52 3.600,0 3.720,0
53 1.800,0 1.860,0
54 1.860,0 1.860,0
55 1.860,0 1.860,0
56 930,0 900,0
57 360,0 930,0
58 390,0 390,0
59 1.800,0 1.860,0
60 900,0 1.860,0
61 900,0 930,0
62 930,0 930,0
63 930,0 930,0
64 930,0 900,0
65 930,0 900,0
66 900,0 930,0
67 930,0 900,0
68 2.790,0 2.790,0
69 1.860,0 1.800,0
70 390,0 900,0
71 930,0 900,0
72 360,0 390,0
73 900,0 930,0
74 930,0 930,0
75 2.520,0 2.040,0
76 900,0 840,0
77 1.680,0 2.820,0
78 1.680,0 1.860,0
79 840,0 930,0
80 930,0 930,0
81 1.680,0 1.860,0
82 1.860,0 1.800,0
83 900,0 840,0
84 930,0 1.860,0
85 930,0 900,0
86 930,0 390,0
87 360,0 360,0
88 930,0 930,0
89 930,0 900,0
90 930,0 900,0
91 840,0 1.860,0
92 360,0 360,0
93 1.800,0 1.680,0
wilcox.test(dados[,-c(1,2,5,6,7)]$DoseT-3,
+ dados[,-c(1,2,5,6,7)]$DoseT+3, paired=TRUE)
Error in wilcox.test.default(dados[, -c(1, 2, 5, 6, 7)]$DoseT - 3, dados[, :
not enough (finite) 'x' observations
Cumprimentos Ana
Quoting r-br-request em listas.c3sl.ufpr.br:
> Enviar submissões para a lista de discussão R-br para
> r-br em listas.c3sl.ufpr.br
>
> Para se cadastrar ou descadastrar via WWW, visite o endereço
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> ou, via email, envie uma mensagem com a palavra 'help' no assunto ou
> corpo da mensagem para
> r-br-request em listas.c3sl.ufpr.br
>
> Você poderá entrar em contato com a pessoa que gerencia a lista pelo
> endereço
> r-br-owner em listas.c3sl.ufpr.br
>
> Quando responder, por favor edite sua linha Assunto assim ela será
> mais específica que "Re: Contents of R-br digest..."
>
>
> Tópicos de Hoje:
>
> 1. trabalho em R test wilcox ou regressão (alanarocha em sapo.pt)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 18 Jun 2013 11:06:30 +0100
> From: alanarocha em sapo.pt
> To: r-br <r-br em listas.c3sl.ufpr.br>, fjsequeira em fc.ul.pt,
> sandro.matos em vodafone.com, rcoster em gmail.com,
> cezar-fonseca em hotmail.com
> Subject: [R-br] trabalho em R test wilcox ou regressão
> Message-ID:
> <20130618110630.Horde.PzJaJaeWCxtRwDEmmikGXKA em mail.sapo.pt>
> Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes
>
>
> Bom dia tenho a seguinte questão a resolver, alguém me sabe dizer se
> devo utilizar o teste de wilcoxon emparelhado ou um modelo de
> regressão linear
> 1. Avaliar as médias para cada caso, considerando 2 períodos
> a. Período antes é igual aos Doset-9 ao Doset-3, que passará a ser
> designado como H
> b. Período depois é igual aos Doset+3 ao Doset+9, que passará a ser
> designado como C
> DoseT-9 DoseT-6 DoseT-3 DoseT+3 DoseT+6 DoseT+9
> 930,0 900,0 900,0 390,0 1.170,0 2.700,0
> 1.860,0 930,0 900,0 420,0 840,0 630,0
> 1.860,0 2.790,0 2.700,0 1.620,0 1.440,0 1.530,0
> 930,0 930,0 900,0 420,0 840,0 720,0
> 1.080,0 480,0 870,0 2.970,0 2.700,0 2.520,0
> 2.790,0 1.590,0 930,0 390,0 720,0 2.340,0
> 930,0 930,0 900,0 390,0 390,0 390,0
> 1.860,0 750,0 780,0 720,0 780,0 720,0
> 3.720,0 3.720,0 3.600,0 840,0 840,0 840,0
> 780,0 420,0 60,0 390,0 390,0 420,0
> 2.700,0 2.790,0 2.790,0 720,0 780,0 360,0
> 930,0 930,0 900,0 420,0 420,0 360,0
> 2.700,0 3.720,0 3.720,0 1.170,0 1.170,0 540,0
> 1.860,0 930,0 900,0 420,0 420,0 720,0
> .....
> no total são 93 dados para cada coluna.
>
> a minha ideia é fazer o calculo das médias entre doset-9 e doset+9
> entre -6 e +6 e entre -3 e +3.
> Ana
>
>
> ------------------------------
>
> Subject: Legenda do Digest
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
> ------------------------------
>
> Fim da Digest R-br, volume 30, assunto 24
> *****************************************
Mais detalhes sobre a lista de discussão R-br