Erro no teste de Scott-Knott
Caro Colegas, Boa tarde. Estou tentando fazer um teste de Scott-Knott componente 1 (Can1) de uma discriminante canonica, entretanto aparece o erro "There were 14 warnings (use warnings() to see them)", usando "warnings()" aparece as mensagens abaixo. Alguém saberia como solucionar esse problema? Obrigado. -- Alisson Lucrecio da Costa ###########------------------------------------------------------------------------########### str(can1_scores) 'data.frame': 40 obs. of 2 variables: $ Tempo: Factor w/ 7 levels "0","2","4","6",..: 1 1 1 1 1 1 1 1 1 1 ... $ Can1 : num -3.98 -2.08 -3.15 -2.86 -2.43 ... can1_scores Tempo Can1 1 0 -3.97903121 2 0 -2.08094673 3 0 -3.14865243 4 0 -2.86189110 5 0 -2.43198428 6 0 -2.61786107 7 0 -2.43352998 8 0 -2.40806466 9 0 -2.46111031 10 0 -3.71884011 11 2 -2.05034264 12 2 -2.00522155 13 2 -3.04292682 14 2 -1.07030352 15 2 -1.59189582 16 4 -0.48869443 17 4 -0.38234727 18 4 -0.28869622 19 4 -1.63394988 20 4 -0.38818468 21 6 1.04914435 22 6 1.35412327 23 6 0.06221095 24 6 1.72914975 25 6 1.49862908 26 8 1.88308179 27 8 3.76305713 28 8 2.61399465 29 8 3.24530432 30 8 3.28546337 31 10 2.19512933 32 10 2.71337547 33 10 2.47481491 34 10 1.49630701 35 10 1.79375735 36 12 2.24198697 37 12 2.26729170 38 12 1.57615671 39 12 2.00861523 40 12 1.83288137 av1 <- aov(Can1~Tempo, data=can1_scores) summary(av1) Df Sum Sq Mean SqLevels Means SK(5%) 8 2.9581803 a 10 2.1346768 b 12 1.9853864 b 6 1.1386515 c 4 -0.6363745 d 2 -1.9521381 e 0 -2.8141912 f There were 14 warnings (use warnings() to see them) F value Pr(>F) Tempo 6 193.01 32.17 88.92 <2e-16 *** Residuals 33 11.94 0.36 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 sk1 <- summary(SK(av1, which ='Tempo')) Levels Means SK(5%) 8 2.9581803 a 10 2.1346768 b 12 1.9853864 b 6 1.1386515 c 4 -0.6363745 d 2 -1.9521381 e 0 -2.8141912 f There were 14 warnings (use warnings() to see them) warnings() Warning messages: 1: In if ((lam < valchisq) | (ord1 == k)) { ... : the condition has length > 1 and only the first element will be used 2: In if ((lam < valchisq) | (ord1 == k)) { ... : the condition has length > 1 and only the first element will be used 3: In if ((lam < valchisq) | (ord1 == k)) { ... : the condition has length > 1 and only the first element will be used 4: In if (lam > valchisq) { ... : the condition has length > 1 and only the first element will be used 5: In if (lam < valchisq) { ... : the condition has length > 1 and only the first element will be used 6: In if ((lam < valchisq) | (ord1 == k)) { ... : the condition has length > 1 and only the first element will be used 7: In if (lam > valchisq) { ... : the condition has length > 1 and only the first element will be used 8: In if (lam < valchisq) { ... : the condition has length > 1 and only the first element will be used 9: In if ((lam < valchisq) | (ord1 == k)) { ... : the condition has length > 1 and only the first element will be used 10: In if (lam > valchisq) { ... : the condition has length > 1 and only the first element will be used 11: In if (lam < valchisq) { ... : the condition has length > 1 and only the first element will be used 12: In if ((lam < valchisq) | (ord1 == k)) { ... : the condition has length > 1 and only the first element will be used 13: In if (lam > valchisq) { ... : the condition has length > 1 and only the first element will be used 14: In if (lam < valchisq) { ... : the condition has length > 1 and only the first element will be used ###########------------------------------------------------------------------------###########
Apesar de ter input e output do R, não é reproduzível porque os dados não estão em um formato de imediata leitura pelo R. Procure fornecer o dput() dos dados ao invés ou junto do print(). O dput() permite leitura imediata via copia/cola, isso queima etapas e convida mais, além do mais você não menciona os pacotes carregados na função, presumo que seja o pacote do colega JC Faria, mas não exija adivinhação dos leitores. À disposição. Walmes.
Caro Walmes Boa noite. O pacote é o ScottKnott
dput(can1_scores)structure(list(Tempo = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L), .Label = c("0", "2", "4", "6", "8", "10", "12"), class = "factor"), Can1 = c(-3.97903120666332, -2.08094673287043, -3.14865243445412, -2.86189109535914, -2.43198427698181, -2.6178610707061, -2.43352997837457, -2.40806466202522, -2.46111031174995, -3.71884011193399, -2.05034263703473, -2.0052215528943, -3.0429268229614, -1.07030351926172, -1.59189581930415, -0.488694426925713, -0.382347271019961, -0.28869621735626, -1.63394987567049, -0.388184676668998, 1.0491443457085, 1.35412326627344, 0.06221094779304, 1.72914975430105, 1.49862908063138, 1.88308178933803, 3.76305713116495, 2.6139946529437, 3.24530431502092, 3.28546336805719, 2.19512933319006, 2.71337547039165, 2.4748149110824, 1.49630700969963, 1.79375734523888, 2.24198697232919, 2.26729170076295, 1.57615670948704, 2.00861522820103, 1.83288136860133 )), .Names = c("Tempo", "Can1"), row.names = c(NA, 40L), class = "data.frame")
str(can1_scores)'data.frame': 40 obs. of 2 variables: $ Tempo: Factor w/ 7 levels "0","2","4","6",..: 1 1 1 1 1 1 1 1 1 1 ... $ Can1 : num -3.98 -2.08 -3.15 -2.86 -2.43 ...
av1 <- aov(Can1~Tempo, data=can1_scores)> summary(av1) Df Sum Sq Mean Sq F value Pr(>F) Tempo 6 193.01 32.17 88.92 <2e-16 *** Residuals 33 11.94 0.36
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1> sk1 <- summary(SK(av1, which ='Tempo')) Levels Means SK(5%) 8 2.9581803 a 10 2.1346768 b 12 1.9853864 b 6 1.1386515 c 4 -0.6363745 d 2 -1.9521381 e 0 -2.8141912 fThere were 14 warnings (use warnings() to see them) 2014-05-13 17:39 GMT-03:00 walmes . <walmeszeviani@gmail.com>:
Apesar de ter input e output do R, não é reproduzível porque os dados não estão em um formato de imediata leitura pelo R. Procure fornecer o dput() dos dados ao invés ou junto do print(). O dput() permite leitura imediata via copia/cola, isso queima etapas e convida mais, além do mais você não menciona os pacotes carregados na função, presumo que seja o pacote do colega JC Faria, mas não exija adivinhação dos leitores.
À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Alisson Lucrecio da Costa
participantes (2)
-
Alisson Lucrecio -
walmes .