[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