Tamanho ótimo de parcela para avaliação da massa e diâmetro de cabeças de brócolis

Optimum plot size for evaluating the mass and diameter of broccoli heads

RESUMO:

O objetivo deste trabalho foi determinar o tamanho ótimo de parcela (Xo), por amostragem na área total e nas frações de áreas, a fim de avaliar a massa e o diâmetro de cabeças de brócolis. Em campo, foram demarcadas 27 filas do híbrido BRO 68^(r). Cada planta foi considerada como sendo uma unidade experimental básica (UEB). Para fins de análise, cada fila foi considerada como sendo um ensaio em branco com 100 UEB. Para os resultados das 27 filas, obtiveram-se os valores: mínimo, máximo, média, desvio padrão, CV, intervalo de confiança, o tamanho de amostra (número de filas) para estimar o Xo e para estimar a média. O tamanho ótimo de parcela para avaliar a massa e diâmetro de cabeças de brócolis é igual a cinco unidades experimentais básicas (plantas).

Palavras-chave: Brassica oleracea var. italica; variação espacial; amostragem; planejamento experimental.

ABSTRACT:

The objective of this study was to determine the optimum plot size by sampling the total area and the fractions' areas in order to evaluate the mass and diameter of broccoli heads. In the field, 27 rows of hybrid BRO 68^(r) were demarcated. Each plant was considered as an experimental basic unit (UEB). For analysis purposes, each row was regarded as being one blank test on 100 UEB. For the results of 27 rows yielded values: minimum, maximum, average, standard deviation, CV, confidence interval, the sample size (number of rows) to estimate Xo and to estimate the mean. The optimum plot size to evaluate the mass and diameter of broccoli heads equals five basic experimental units (plants).

Key words: Brassica oleracea var. italica; spatial variation; sampling; experimental design.

INTRODUÇÃO:

O brócolis (Brassica oleracea L. var. italica Plenck) é uma ótima fonte de nutrientes, além disso, o consumo regular previne alguns tipos de câncer, doenças cardíacas e purifica o sangue. Mesmo com todos os nutrientes que oferece, é um alimento pobre em calorias, sendo indicado em dietas de perda de peso (ROSA & RODRIGUES, 2001). Além disso, é uma hortaliça rica em Boro (B), que atua no metabolismo de substratos energéticos, no funcionamento do cérebro e na atividade psicomotora e cognitiva (PIZETTA et al., 2005).

Para obter melhor produtividade e qualidade de brócolis, vêm sendo realizadas pesquisas referentes às adubações, aos manejos, incidência de doenças (DINIZ, 2008; CAMPAGNOL et al., 2009; TREVISAN et al., 2003), entre outros, sem o conhecimento correto do planejamento experimental. Nesse sentido, o conhecimento da variabilidade da produção, numa área experimental, é importante para o planejamento adequado de experimentos. O planejamento define, para um determinado número de tratamentos, o delineamento experimental, o tamanho da parcela e o número de repetições (RAMALHO et al., 2000; STORCK et al., 2011).

Para conhecer a variabilidade na área experimental, podem ser conduzidos ensaios em branco, específicos para a identificação da variabilidade ou do índice de heterogeneidade da característica avaliada. Nessa metodologia, são coletados os dados da produção de pequenas parcelas denominadas de unidades experimentais básicas (UEB) e, com estas UEB, são planejadas parcelas de diferentes tamanhos, pelo agrupamento de UEB adjacentes. Para cada tamanho de parcela (X), são calculados os coeficientes de variação (CVx). Da relação CVx = A/X^B é estimado o tamanho ótimo de parcela (Xo) como sendo o ponto de máxima curvatura dessa relação (MEIER & LESSMAN, 1971), sendo A e B os parâmetros da equação. No entanto, a dificuldade e o trabalho em planejar diferentes tamanhos de parcelas pelo agrupamento de UEB adjacentes, no método de máxima curvatura, foram superados com a aplicação do método da curvatura máxima do modelo do coeficiente de variação (PARANAÍBA et al., 2009a). Nesse método, a partir dos dados das UEB de um ensaio em branco, é estimado o coeficiente de autocorrelação espacial de primeira ordem, a variância e a média, usados para o cálculo do tamanho ótimo de parcela. Este método tem sido adequado para a culturas de milho (CARGNELUTTI FILHO et al., 2011a), arroz (PARANAÍBA et al., 2009a), trigo e mandioca (PARANAÍBA et al., 2009b).

A possibilidade, para algumas espécies cultivadas em filas quase independentes (mandioca, batata, couve, brócolis, etc), de se usar os dados de algumas linhas ou grupos de linhas com determinado comprimento, de forma aleatória sobre toda a área experimental, é uma alternativa viável, segundo estudos na cultura da batata (STORCK, 2011), resultando em economia de recursos financeiros e humanos. A determinação do tamanho de parcela por amostragem de quadros menores na área experimental é favorecida pela facilidade de aplicação do método de PARANAÍBA et al. (2009a).

Usualmente, em experimentos com brócolis, foram usadas parcelas de 10 plantas na comparação da produtividade de cultivares (LALLA et al., 2010), 14 plantas na avaliação de cultivares em plantio direto sobre diferentes culturas de cobertura (MELO et al., 2010) e 14 plantas (10 úteis) na avaliação de cultivares (TREVISAN et al., 2003). No entanto, o tamanho ótimo de parcela e o número de ensaios (quadros) necessários para estimativa do tamanho ótimo de parcela em couve e em brócolis não são conhecidos. Assim, os objetivos deste trabalho foram determinar o tamanho ótimo de parcela por amostragem na área total e nas frações de áreas, a fim de avaliar a massa e o diâmetro de cabeças de brócolis.

MATERIAL E MÉTODOS:

O experimento foi conduzido na área experimental da UTFPR em Pato Branco, PR. A área está situada na latitude 26º 06' 59"S, longitude 52º 40' 59"W, e com altitude de 721,80m. O município de Pato Branco apresenta clima tipo "cfa" subtropical úmido, com chuvas bem distribuídas ao longo do ano, temperatura média no mês mais frio inferior a 18°C, temperatura média no mês mais quente acima de 22°C, com verões relativamente quentes. Foram demarcadas 27 filas de 50 metros de comprimento, espaçadas em 0,80m, resultando em 1080m² de área. A bordadura da área foi de, aproximadamente, quatro metros em todos os lados. Nessa área, foram plantadas mudas de brócolis com 0,50m entre plantas. As mudas de brócolis foram produzidas em estufa usando sementes do híbrido BRO 68^(r) e substrato Tecnomax^(r), preparadas em bandejas de poliestireno expandido de 128 células. Aos 46 dias após a semeadura, foram transplantadas para a área experimental demarcada.

Antes do plantio, a área foi preparada com subsolagem e enxada rotativa. A adubação de base foi realizada nas filas com 87,5kg ha^-1 de N (ureia), 1000kg ha^-1 de P₂O₅ (superfosfato simples) e 25kg ha^-1 de KCl (cloreto de potássio). Aos 20 e aos 44 dias após o transplante, foram realizadas adubações de cobertura com 87,5kg ha^-1 de ureia e 25kg ha^-1 de KCl. Na primeira adubação de cobertura, o adubo foi incorporado superficialmente, realizando-se também a amontoa para melhorar, principalmente, a sustentação das plantas.

Foram realizadas três aplicações foliares dos micronutrientes Boro (B) e Molibdênio (Mo), em quantidades suficientes para completo molhamento foliar e concentração de ácido bórico (1g L^-1) e molibdato de sódio (0,5g L^-1). A primeira aplicação foi realizada na fase de mudas, e, as outras, realizadas aos 15 e 30 dias após o transplantio.

No ponto de colheita, as cabeças foram cortadas e imediatamente avaliadas para evitar a perda de água por desidratação. Foram determinados diâmetro (cm) e massa (g), anotando a posição em que estavam no campo (fila e coluna). A posição de campo foi identificada pelo número de ordem da fila (1 a 27) e do número de ordem da coluna (1 a 100), formando uma matriz de 2.700 observações.

Cada planta foi considerada como sendo uma unidade experimental básica (UEB). Para fins de análise, cada fila foi considerada como sendo um ensaio em branco (ensaio de uniformidade) com 100 UEB. Para cada ensaio e para cada um dos caracteres (massa e diâmetro), foi estimada a média (m), a variância (s²⁾ e o coeficiente de variação (CV) entre as plantas. Também foi estimado o coeficiente de autocorrelação espacial de primeira ordem (), usando a sequência de numeração das plantas dentro de cada fila. Sendo X_i igual ao valor observado na cabeça i (i=1, 2, ... n=100), tem-se: D_i=X_i-m e segundo PARANAÍBA et al. (2009a). Conforme os mesmos autores, o tamanho ótimo da parcela (Xo, número de plantas adjacentes) por fila foi estimado por meio do método da curvatura máxima do modelo do coeficiente de variação (CMMCV) pela expressão , sendo que o coeficiente de variação entre as parcelas de tamanho Xo foi estimado pela expressão .

Para os resultados obtidos nas 27 filas (ensaios), obtiveram-se os valores: mínimo, máximo, média, desvio padrão, CV, intervalo de confiança (1-a=0,95), o tamanho de amostra (número de filas) para estimar o Xo (com erro de estimação equivalente a amplitude de uma UEB) e para estimar a média (com erro de estimação igual a 10% da média do caractere), com base na distribuição t (a=0,05). A normalidade dos dados observados de Xo e da média (m) foi verificada pelo aplicativo WinSTAT(r) for Microsoft(r)Excel. Para os demais procedimentos, foi usado o software R (R DEVELOPMENT CORE TEAM, 2013) e o aplicativo Microsoft Office Excel.

A área de teste (2.700 UEB) foi dividida em três quadros, com as filas variando entre 1 a 27 e as colunas entre 01 a 33 (quadro 1), entre 34 e 67 (quadro 2) e entre 68 e 100 (quadro 3). Para cada fila, dentro de cada quadro com 33 (891 UEB) ou 34 (918 UEB) colunas, foram estimadas as estatísticas: m, s2, CV, , Xo e CVXo.

Dentro de cada quadro, foi estimado o tamanho de amostra (número de filas) para estimar o Xo (com erro de estimação igual a uma planta), com base na distribuição t (α=0,05).

RESULTADOS E DISCUSSÃO:

A média geral da massa por cabeça de brócolis foi igual a 1113g. Considerando o espaçamento (0,5 x 0,8m), a produtividade foi igual a 27,8t ha-1. Este resultado é superior ao obtido em experimento de comparação de cultivares, máximo de 15,4t ha-1 para a cultivar AF649 (LALLA et al., 2010), máximo de 12,5t ha-1 para adubação orgânica de 25t ha-1 (DINIZ et al., 2008), 13,2t ha-1 (MELO et al., 2010) e 16,3t ha-1 (TREVISAN et al., 2003).

As plantas com falhas (4,7% das plantas), ou seja, com ausência de produção de cabeça de brócolis, motivado por danos mecânicos, doenças e pragas, foram substituídas pela média das duas plantas vizinhas. Normalmente, em experimentação, havendo perda de plantas não causada pelo efeito de tratamento, solução semelhante é adotada, ou seja, é obtida a média das unidades presentes (STORCK et al., 2011). Portanto, pela produtividade obtida e a baixa frequência de falhas, os dados obtidos são adequados para o estudo proposto.

A área experimental não é homogênea, fato que pode ser observado na figura 1 (massa de cabeça e diâmetro de cabeça). Nesta figura, quando LS de uma fila é menor que LI de outra fila, estas filas diferem significativamente (α<0,05) pelo teste t. Dessa forma, observam-se diferenças entre diversas filas quanto aos dois caracteres (massa e diâmetro). A justificativa para tais diferenças está nas condições iniciais de fertilidade da área e na introdução de variação, devido à adubação e manejo que foram executados por fila. Constatação semelhante foi relatada, anteriormente, com a cultura da batata (OLIVEIRA et al., 2005), possivelmente pelas mesmas razões. Nesses casos, a solução é adotar o delineamento blocos ao acaso, em que os blocos são formados pelas filas, ou seja, cada bloco é um grupo de unidades experimentais homogêneas (ZIMMERMANN, 2004; RAMALHO et al., 2005).

Figura 1
Variação espacial do limite inferior (LI) e superior (LS) do intervalo de confiança (1-p = 0,95) para a média da massa e do diâmetro de cabeças de brócolis, em 27 filas de 100 plantas.

As estimativas da média (m), s2, CV, , Xo e CVXo estão relatadas na tabela 1 (massa de cabeça) e tabela 2 (diâmetro de cabeça). O tamanho ótimo de parcela para a estimativa da massa de cabeça de brócolis é igual a 5,24 UEB (cinco plantas), variando entre os limites de 5,09 e 5,38 plantas, com 95% de confiança. O tamanho ótimo de amostra para a variável diâmetro de cabeça é menor do que o necessário para massa de cabeça, mas é igualmente homogêneo entre as 27 filas (Tabela 2). Esses resultados são inferiores ao tamanho de parcela tradicionalmente utilizado na avaliação da produção de brócolis (TREVISAN et al., 2003; LALLA et al., 2010; MELO et al., 2010).

Tabela 1
Média da massa (g cabeça-1) de brócolis, variância (s2), coeficiente de variação (CV), coeficiente de autocorrelação espacial (), tamanho ótimo de parcela (Xo) e CV para Xo (CVXo), para as diferentes filas da área de teste.

• Limite inferior (LI) e superior (LS) do intervalo de confiança (1-(=0,95) e tamanho de amostra (no, (=0,05) para erro de estimação de 10% da média (da massa) e erro de uma UEB em Xo.

Tabela 2
Média do diâmetro (cm) de cabeça de brócolis, variância (s2), coeficiente de variação (CV), coeficiente de autocorrelação espacial (), tamanho ótimo de parcela (Xo) e CV para Xo (CVXo) para as diferentes filas da área de teste.

• *Limite inferior (LI) e superior (LS) do intervalo de confiança (1-(=0,95) e tamanho de amostra (no, (=0,05) para erro de estimação de 10% da média (do diâmetro) e erro de uma UEB em Xo.

O coeficiente de autocorrelação espacial de primeira ordem (), no caso do caractere massa, é baixo e não difere significativamente (1-a=0,95) de zero, o intervalo de confiança contém o valor zero. Com isso, existe uma independência espacial entre as UEB da área experimental, propriedade importante para os pressupostos dos modelos dos delineamentos experimentais. Para o caractere diâmetro de cabeça, apesar do ser significativamente maior do que zero, a amplitude de variação (0,054< ρ <0,167; 1-α= 0,95) é pequena (Tabela 2).

Observando os valores de LI e LS de , constata-se que, para o caractere diâmetro de cabeça, a autocorrelação () é significativamente maior que zero, pois os limites entre LI (=0,054) e LS (=0,167) não contém o valor zero. Considerando que o método para estimar Xo não depende desta significância, não vemos restrições de validade dos resultados. As estimativas de m e Xo, referentes as 27 filas, se ajustam à distribuição normal, constatada pelos testes de Kolmogorov-Smirnov (α>0,60). Com isso, é possível determinar o tamanho de amostra (número de filas) usando a distribuição t. Para estimar o caractere massa (g) por cabeça de brócolis, é suficiente avaliar três fileiras de 100 plantas, considerando um erro de estimação de 10% da média. Com aproximadamente o mesmo erro de estimação, para avaliar o diâmetro da cabeça, é necessário avaliar apenas uma fila (Tabela 2). O tamanho de amostra (número de filas) para estimar o tamanho ótimo de parcela, usando o valor máximo dos dois caracteres (massa e diâmetro) é igual a três filas com erro de estimação equivalente a uma amplitude de uma UEB (menor erro possível) (Tabela 2). Dessa forma, numa área experimental, amostrar três filas de 100 plantas é o suficiente para estimar Xo (com erro de estimação equivalente a uma UEB), visto que existe um padrão de variabilidade dentro das filas.

A prática de amostragem de partes da área experimental (pequenos lotes) para a determinação do tamanho ótimo de parcela também foi relatada para ensaios de batata (STORCK, 2011), como indicativo para a redução dos recursos humanos. Dessa forma, é possível sugerir que a estimativa do tamanho ótimo de parcela, referente a uma área experimental grande, pode ser realizada a partir da amostragem de alguns lotes menores dentro da área maior. A média dos tamanhos ótimos desses lotes é válida para a área total.

Ao dividir a área experimental em três quadros (Tabela 3), as médias das estimativas doe de Xo referentes aos diferentes quadros também são homogêneos e semelhantes aos obtidos com 27 filas de 100 plantas. O tamanho de amostra (número de filas) aumentou, praticamente, na proporção inversa da redução do tamanho do ensaio em branco. Assim, ao usar 27 filas de 100 plantas, são necessárias 2,3 filas ou 230 plantas (230/2700 ou 8,5% da área total) e, ao usar 27 filas de 33 ou 34 plantas, são necessárias 9,12 filas ou 297 plantas (297/2700 = 11,0% da área total). Na condição de ensaios de 33 ou 34 plantas, observa-se maior variabilidade entre as estimativas de Xo e por isso o tamanho de amostra é maior. Este fato corrobora com o reportado por CARGNELUTTI FILHO et al. (2011b), que observaram relação entre o tamanho dos ensaios e a variação do tamanho ótimo de parcelas em estudo de simulação com dados de massa fresca de nabo forrageiro. Portanto, o tamanho do ensaio em branco não pode ser tão pequeno, porque aumenta a amplitude de variação do tamanho de parcela e de outras estatísticas.

Tabela 3
Média do coeficiente de autocorrelação espacial () e do tamanho ótimo de parcela (Xo) e tamanho de amostra (número de filas, n_o(Xo)) para estimar Xo na área de teste de produção de brócolis.

Em geral, a estratégia de amostrar "quadros" menores, dentro de uma futura área experimental, é apropriada para a determinação do plano experimental, quanto ao tamanho de parcela, número de repetições e delineamento experimental. Essa afirmação está amparada nos resultados deste estudo e nos resultados com ensaio de uniformidade de batata (STORCK, 2011), em que 39% dos quadros menores, dentro da área total, foram suficientes para determinar o tamanho ótimo de parcela.

CONCLUSÃO:

O tamanho ótimo de parcela para avaliar a massa e diâmetro de cabeças de brócolis é igual a cinco unidades experimentais básicas (plantas).

AGRADECIMENTOS

Ao Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (bolsa de produtividade em pesquisa e IC) e à Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) (Edital 6/2012).

REFERÊNCIAS:

• CAMPAGNOL, R. et al. Boro e nitrogênio na incidência de hastes ocas e no rendimento de brócolis. Ciência e Agrotecnologia, v.33, p.1477-1485, 2009. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S1413-70542009000600004 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S1413-70542009000600004
• CARGNELUTTI FILHO, A. et al. Tamanho ótimo de parcela em milho com comparação de dois métodos. Ciência Rural, v.41, p.1890-1898, 2011a. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S0103-84782011001100007 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S0103-84782011001100007
• CARGNELUTTI FILHO, A. et al. Tamanhos de parcela e de ensaio de uniformidade em nabo forrageiro. Ciência Rural, v.41, p.1517-1525, 2011b. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S0103-84782011005000119 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S0103-84782011005000119
• DINIZ, E.R. et al. Crescimento e produção de brócolis em sistema orgânico em função de doses de composto. Ciência e Agrotecnologia, v.32, p.1428-1434, 2008. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S1413-70542008000500011 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S1413-70542008000500011
• LALLA, J.G. et al. Competição de cultivares de brócolos tipo cabeça única em Campo Grande. Horticultura Brasileira, v.28, p.360-363, 2010. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S0102-05362010000300020 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S0102-05362010000300020
• MEIER, V.D.; LESSMAN, K.J. Estimation of optimum field plot shape and size testing yield in Crambe abyssinica hordnt. Crop Science, v.11, p.648-650, 1971.
• MELO, R.A.C. et al. Cultivo de brócolos de inflorescência única no verão em plantio direto. Horticultura Brasileira, v.28, p.23-28, 2010. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S0102-05362010000100005 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S0102-05362010000100005
• OLIVEIRA, S.J.R. et al. Plot size and experimental unit relationship in exploratory experiments. Scientia Agricola, v.62, p.585-589, 2005. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S0103-90162005000600012 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S0103-90162005000600012
• PARANAÍBA, P.F. et al. Tamanho ótimo de parcelas experimentais: proposição de métodos de estimação. Revista Brasileira de Biometria, v.27, p.255-268, 2009a. Disponível em: <Disponível em: http://jaguar.fcav.unesp.br/RME/biometria/inicial.php >. Acesso em: 24 jun. 2014.
  » http://jaguar.fcav.unesp.br/RME/biometria/inicial.php
• PARANAÍBA, P.F. et al. Tamanho ótimo de parcelas experimentais: comparação de métodos em experimentos de trigo e mandioca. Revista Brasileira de Biometria, v.27, p.81-90, 2009b. Disponível em: <Disponível em: http://jaguar.fcav.unesp.br/RME/biometria/inicial.php >. Acesso em: 24 jun. 2014.
  » http://jaguar.fcav.unesp.br/RME/biometria/inicial.php
• PIZETTA, L.C. et al. Resposta de brócolis, couve-flor e repolho à adubação com boro em solo arenoso. Horticultura Brasileira, v.23, p.51-56, 2005. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S0102-05362005000100011 >. Acesso em: 10 jun. 2015. doi: 10.1590/S0102-05362005000100011.
  » https://doi.org/10.1590/S0102-05362005000100011.» http://dx.doi.org/10.1590/S0102-05362005000100011
• R DEVELOPMENT CORE TEAM. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2013. Disponível em: <http://www.R-project.org>.
  » http://www.R-project.org
• RAMALHO, M.A.P. et al. Experimentação em genética e melhoramento de plantas. 2.ed. Lavras: UFLA, 2005. 322p.
• RAMALHO, M.A.P. et al. A. Experimentação em genética e melhoramento de plantas Lavras: EDUFLA, 2000. 322p.
• ROSA E.A.S; RODRIGUES, A.S. Total and individual glucosinalate content in 11 broccoli cultivars grown in early and late seasons. HortScience, v.36, p.56-59, 2001. Disponível em: <Disponível em: http://hortsci.ashspublications.org/content/36/1/56.full.pdf+html >. Acesso em: 10 jun. 2015.
  » http://hortsci.ashspublications.org/content/36/1/56.full.pdf+html
• STORCK, L. Partial collection of data on potato yield for experimental planning. Field Crops Research, v.121, p.286-290, 2011. Disponível em: <Disponível em: http://dx.doi.org/10.1016/jfcr2010.12.018 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1016/jfcr2010.12.018
• STORCK, L. et al. Experimentação vegetal. 3.ed. Santa Maria: UFSM, 2011. 198p.
• TREVISAN, J.N. et al. Rendimento de cultivares de brócolis semeadas em outubro na região centro do Rio Grande do Sul. Ciência Rural, v.33, p.233-239, 2003. Disponível em: <Disponível em: http://dx.doi.org/10.1590/S0103-84782003000200009 >. Acesso em: 10 jun. 2015.
  » http://dx.doi.org/10.1590/S0103-84782003000200009
• ZIMMERMANN, F.J.P. (Ed.). Estatística aplicada à pesquisa agrícola. Santo Antonio de Goiás: Embrapa Arroz e Feijão, 2004. 402p.

• 1
  *CR-2015-0236.R2

Datas de Publicação

• Publicação nesta coleção
  30 Nov 2015
• Data do Fascículo
  Mar 2016

Histórico

• Recebido
  22 Fev 2015
• Aceito
  23 Jun 2015
• Revisado
  23 Out 2015