- R-br - DINF UFPR

Distribuição de probabilidade
by Felipe Felix Costa 10 Jan '18

10 Jan '18

olá, gostaria de saber se existe alguma forma para rodar todas as possíveis distribuições sobre um determinado conjunto de dados com variáveis continuas!? pois estou procurando e não achei nada pela internet falando dessa possibilidade. meu objetivo é verificar qual melhor a distribuição para os meus dados e, a partir dessa "melhor", estabelecer um modelo geral com todos os dados e comparar com modelos individuais de áreas isoladas. HELP Att., -- *Felipe Felix Costa* Tel. 96-98122-3077

5 8

Re: [R-br] Salvando Gráficos num Loop
by sznelwar＠uol.com.br 05 Jan '18

05 Jan '18

1 0

Salvando Gráficos num Loop
by Yury Duarte 05 Jan '18

05 Jan '18

Bom dia colegas listeiros! Estou tendo um pequeno problema para automatizar a criação e salvamento de gráficos temáticos utilizando um loop. Quando rodo o script lendo um arquivo de cada vez (fora do loop), o gráfico é salvo sem prejuízo na pasta de destino, entretanto, quando peço para rodar o script dentro do loop, ele salva os gráficos em branco. Gostaria de entender o porque dessa peculiaridade, já que programa não me retorna nenhum aviso de falha ou problema. O código utilizado segue abaixo e também anexei dois arquivos distintos dos bancos de dados que estou trabalhando, caso alguém queira/precise rodar o script. Desde já, agradeço pela ajuda de todos! rm(list=ls()) library(hydroTSM) library(lattice) library(stringr) #----------Verificando Diretorio----------# setwd('C:\\Users\\Yury\\Desktop\\Fernandes') length(dir('C:\\Users\\Yury\\Desktop\\Fernandes')) a = grep('.txt', dir(getwd()), value = T) a = str_sub(a, end = 5) #----------Loop de todas as EM----------# for(i in 1:length(a)){ #----------Selecao Arquivo INMET----------# b = read.table(paste0('C:\\Users\\Yury\\Desktop\\Fernandes\\', a[i], '.txt'), header=TRUE) b.cor = data.frame(b$Data, b$Precipitacao, b$TempMaxima, b$TempMinima, b$Insolacao, b$Umidade.Relativa.Media) colnames(b.cor) = c("Data","RAIN","TMAX","TMIN","n","RHUM"); str(b.cor) b.cor$Data = as.Date(levels(b$Data))[b$Data] #----------Criar Vetor Completo de Datas----------# vetor.data = seq(as.Date('01/07/1961', "%d/%m/%Y"), as.Date('30/06/2015', "%d/%m/%Y"), 1) Data = data.frame(Data=vetor.data) #----------Corrigir Coluna de Datas com Falha----------# duff = merge(b.cor, Data, all = TRUE) str(duff) #----------Criando Graficos----------# pdf(paste0('C:\\Users\\Yury\\Desktop\\Fernandes\\Graphs\\', a[i], '.pdf'), width = 20, height = 14) matrixplot(dwi(duff, var.type="Days"),ColorRamp="Days", main = "Estação_b") dev.off() } Yury Duarte Engenheiro Agrônomo - ESALQ/USP

2 2

Dúvida em componentes principais - gráfico
by Ana Julia Righetto 04 Jan '18

04 Jan '18

Prezad@s, Realizei análise de componentes principais e fiquei com 3 componentes. Estou com dificuldade para fazer um gráfico similar ao biplot (quando tem apenas 2CP). Alguém poderia me ajudar? Obrigada desde já -- Ana Julia Righetto Estaticista - Universidade Estadual Paulista "Júlio de Mesquita Filho" - UNESP Mestre em Ciências (Estatística e Experimentação Agronômica) - ESALQ/USP - Piracicaba - SP Doutora em Ciências (Estatística e Experimentação Agronômica) - ESALQ/USP - Piracicaba - SP Telefone:(43) 3422-8700 E-mail: ajrighetto(a)gmail.com ajrighetto(a)usp.br Skype: anajrighetto Lattes: http://lattes.cnpq.br/5933607596600115

3 2

entrada de senhas por teclado
by Antonio Silva 04 Jan '18

04 Jan '18

Car@s Estou preparando um script e em uma de suas partes há o acesso a um banco de dados. Para tal é necessária a entrada de uma senha. Com a função readline consigo pedir a senha ao usuário e atribuir a uma variável que depois de usada é apagada mas o que é digitado é mostrado no console. Eu não queria que a senha aparecesse. Há alguma função ou forma de solicitar para que o usuário digite algum dado mas este não apareça na tela? Uso o R 3.4.3 64 bits no Ubuntu 16.04 Obrigado! --- Antônio Olinto Ávila da Silva Instituto de Pesca São Paulo, Brasil

2 2

Especialização em Data Science & Big Data em Curitiba - PR
by Elias T. Krainski 28 Dez '17

28 Dez '17

Caro(a)s, Alguns professores dos Departamentos de Estatística e de Ciências da Computação da Universidade Federal do Paraná estão unindo forças para abrir uma turma de Especialização em *Data Science & Big Data.* Maiores detalhes estão disponíveis em http://dsbd.leg.ufpr.br/ Att. Elias PS.: Justifico este post nesta lista por estar indiretamente associado ao R, uma vez que será uma das linguagens lecionadas :)

2 1

R Day - Encontro nacional de usuários do R
by Fernando Mayer 20 Dez '17

20 Dez '17

Car@s redistas, é com grande satisfação que anunciamos o R Day - Encontro nacional de usuários do R, um evento para congregar usuários e desenvolvedores de R em um ambiente para troca de experiências, aprendizado e contatos. O R Day acontecerá na UFPR, Curitiba (PR), no dia 22 de maio de 2018 como um evento satélite da 63ª RBras. Sendo este o primeiro evento de R no Brasil endossado oficialmente pela The R Foundation, o objetivo do R Day é reunir a ativa e grande, embora esparsa, comunidade de usuários de R no Brasil. Para atender a todos os públicos, estão programadas 3 conferências, 9 apresentações orais (short-talks) e 8 tutoriais. As apresentações orais serão voltadas para demonstração de pacotes (vinheta falada) ou estudos de caso "out of the box" relatando desafios superados com uso do R, tanto no escopo acadêmico quanto profissional. As apresentações terão duração de 15 min + 5 min para discussão. Os tutoriais terão 1 hora de duração e foram separados em duas linhas que correrão simultaneamente: uma para usuários iniciantes e outra para usuários experientes ou desenvolvedores. Alguns temas de tutorias dentro de cada linha foram sugeridos pela comissão organizadora apenas para orientação dos interessados. No entanto, serão aceitas propostas em temas diferentes. As conferências enfatizarão temas de ampla abrangência sobre R a sua comunidade de usuários. Diferentemente das demais modalidades, os conferencistas serão definidos pela comissão organizadora. Sugestões são bem vindas. Sintam-se todos convidados e estimulados a submeter propostas para apresentações orais e tutoriais. Almejamos a colaboração de todos para que o evento reuna o maior número possível de usuários do R com a intenção de promover a troca de conhecimento, recursos e contatos, além de ampliar e fortalecer nossa comunidade. O respaudo da R Foundation para realização desse evento demonstra o prestígio e reconhecimento da comunidade brasileira de usuários. Esse evento configura uma abertura de campanha para a realização de uma conferência userR! no Brasil nos próximos anos. Por isso, seu apoio e participação são fundamentais. A submissão de propostas está aberta até o dia 08/04/2018 via preenchimento de formulário online. Como R Day é um evento sem fins lucrativos e de inscrição gratuita, sua realização depende de recursos providos por patrocinadores. Aqueles que podem prestar auxílio, seja financeiro, logístico, publicitário, etc, são convidados a entrar em contato com a comissão organizadora. Todas as informações necessárias estão disponíveis no site: http://rday.leg.ufpr.br Para atualizações, acompanhe o Twitter do LEG: @LEG_UFPR e siga a #rdayCWB. Mensagens de apoio e sugestões para o evento podem ser feitos pelo Twitter ou email. Atenciosamente, Fernando Mayer P/ Comissão organizadora do R Day --- Fernando de Pol Mayer Laboratório de Estatística e Geoinformação - LEG Departamento de Estatística - DEST Universidade Federal do Paraná - UFPR URL: http://leg.ufpr.br/~fernandomayer e-mail: fernando.mayer [@] {gmail.com, ufpr.br}

1 0

leitura de arquivo de foto em formato dng (raw)
by Cleber N.Borges 19 Dez '17

19 Dez '17

Olá a todos... Bom dia... Gostaria de ler arquivos de fotos em DNG pelo R. Verifiquei que usando a função readTiff ( pacote tiff ) é possível carregar a matriz com os dados... Porém pela documentação da função, eu fiquei em dúvida se realmente não há nenhum processamento das intensidades capturados pelo sensor de luz da câmera... Alguém tem alguma de dica de como acessar essa informação contida nos arquivos DNG com o R? Pode ser alguma dica que use algum outro software para passos intermediários... mas no final gostaria de produzir matrizes com o R... Obrigado Cleber Borges ps.: imagens serão produzidas smartphone samsung s7 edge --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

1 0

Sigularidade backsolve, o que significa?
by Fernando Souza 15 Dez '17

15 Dez '17

Caros Amigos o que pode significar essa mensagem em modelos não linear mixto (nlme)? Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

3 2

FlexParamCurve
by Marcelo Laia 14 Dez '17

14 Dez '17

Alguém usa ou já usou? https://cran.r-project.org/web/packages/FlexParamCurve/index.html -- Marcelo

1 0

Dúvida gráfico lattice da função plot.nlme
by Fernando Souza 13 Dez '17

13 Dez '17

Caros colegas, Devido ao elevado número de individuos avaliados há uma sobreposição dos identificados no eixo y do gráfico, o que o torna incompreensível. Como posso organizar esses indicadores? Ou como posso removê-los. library(nlme) library(lattice) GALINHA<-structure(list(ID = structure(c(128L, 48L, 109L, 52L, 42L, 16L, 210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L, 66L, 293L, 268L, 205L, 41L, 43L, 2L, 34L, 85L, 193L, 63L, 55L, 131L, 54L, 67L, 178L, 103L, 93L, 56L, 79L, 133L, 80L, 264L, 161L, 114L, 288L, 28L, 113L, 88L, 97L, 160L, 68L, 58L, 246L, 25L, 141L, 94L, 150L, 284L, 60L, 278L, 26L, 29L, 72L, 61L, 49L, 40L, 287L, 98L, 75L, 136L, 45L, 144L, 5L, 11L, 124L, 233L, 273L, 3L, 125L, 152L, 137L, 71L, 46L, 27L, 107L, 59L, 92L, 44L, 51L, 157L, 89L, 134L, 126L, 84L, 82L, 145L, 53L, 111L, 214L, 153L, 74L, 99L, 257L, 306L, 105L, 108L, 116L, 277L, 143L, 230L, 64L, 171L, 122L, 91L, 282L, 135L, 203L, 30L, 147L, 156L, 155L, 35L, 76L, 23L, 185L, 262L, 247L, 123L, 227L, 258L, 238L, 223L, 73L, 130L, 81L, 154L, 104L, 291L, 216L, 112L, 333L, 165L, 212L, 308L, 201L, 313L, 324L, 148L, 321L, 170L, 211L, 102L, 198L, 286L, 229L, 249L, 167L, 159L, 191L, 283L, 320L, 305L, 317L, 355L, 330L, 182L, 221L, 231L, 241L, 166L, 234L, 142L, 199L, 217L, 117L, 121L, 132L, 174L, 266L, 215L, 325L, 65L, 301L, 106L, 179L, 50L, 33L, 115L, 8L, 304L, 345L, 363L, 169L, 36L, 292L, 110L, 47L, 192L, 95L, 101L, 69L, 202L, 146L, 200L, 173L, 100L, 39L, 302L, 336L, 290L, 332L, 86L, 127L, 197L, 310L, 242L, 118L, 90L, 120L, 13L, 285L, 183L, 188L, 158L, 252L, 164L, 177L, 190L, 186L, 237L, 62L, 78L, 240L, 250L, 196L, 228L, 149L, 87L, 243L, 24L, 14L, 204L, 140L, 15L, 244L, 339L, 180L, 303L, 195L, 265L, 289L, 269L, 272L, 350L, 162L, 254L, 344L, 139L, 151L, 226L, 270L, 280L, 327L, 57L, 361L, 176L, 251L, 235L, 351L, 119L, 163L, 206L, 271L, 365L, 314L, 329L, 181L, 295L, 239L, 316L, 70L, 207L, 362L, 10L, 208L, 343L, 342L, 219L, 253L, 298L, 349L, 299L, 222L, 307L, 245L, 334L, 248L, 315L, 366L, 17L, 296L, 369L, 341L, 260L, 309L, 297L, 328L, 261L, 83L, 367L, 38L, 318L, 259L, 358L, 346L, 20L, 331L, 348L, 256L, 275L, 168L, 218L, 194L, 232L, 22L, 368L, 209L, 357L, 267L, 19L, 129L, 326L, 236L, 353L, 359L, 279L, 352L, 360L, 6L, 9L, 340L, 31L, 18L, 274L, 322L, 335L, 300L, 370L, 225L, 364L, 311L, 312L, 337L, 338L, 356L, 354L, 213L, 21L, 220L, 294L, 224L, 319L, 276L, 347L, 263L, 175L, 172L, 138L, 187L, 7L, 4L, 32L, 12L, 128L, 48L, 109L, 52L, 42L, 16L, 210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L, 66L, 293L, 268L, 205L, 41L, 43L, 2L, 34L, 85L, 193L, 63L, 55L, 131L, 54L, 67L, 178L, 103L, 93L, 56L, 79L, 133L, 80L, 264L, 161L, 114L, 288L, 28L, 113L, 88L, 97L, 160L, 68L, 58L, 246L, 25L, 141L, 94L, 150L, 284L, 60L, 278L, 26L, 29L, 72L, 61L, 49L, 40L, 287L, 98L, 75L, 136L, 45L, 144L, 5L, 11L, 124L, 233L, 273L, 3L, 125L, 152L, 137L, 71L, 46L, 27L, 107L, 59L, 92L, 44L, 51L, 157L, 89L, 134L, 126L, 84L, 82L, 145L, 53L, 111L, 214L, 153L, 74L, 99L, 257L, 306L, 105L, 108L, 116L, 277L, 143L, 230L, 64L, 171L, 122L, 91L, 282L, 135L, 203L, 30L, 147L, 156L, 155L, 35L, 76L, 23L, 185L, 262L, 247L, 123L, 227L, 258L, 238L, 223L, 73L, 130L, 81L, 154L, 104L, 291L, 216L, 112L, 333L, 165L, 212L, 308L, 201L, 313L, 324L, 148L, 321L, 170L, 211L, 102L, 198L, 286L, 229L, 249L, 167L, 159L, 191L, 283L, 320L, 305L, 317L, 355L, 330L, 182L, 221L, 231L, 241L, 166L, 234L, 142L, 199L, 217L, 117L, 121L, 132L, 174L, 266L, 215L, 325L, 65L, 301L, 106L, 179L, 50L, 33L, 115L, 8L, 304L, 345L, 363L, 169L, 36L, 292L, 110L, 47L, 192L, 95L, 101L, 69L, 202L, 146L, 200L, 173L, 100L, 39L, 302L, 336L, 290L, 332L, 86L, 127L, 197L, 310L, 242L, 118L, 90L, 120L, 13L, 285L, 183L, 188L, 158L, 252L, 164L, 177L, 190L, 186L, 237L, 62L, 78L, 240L, 250L, 196L, 228L, 149L, 87L, 243L, 24L, 14L, 204L, 140L, 15L, 244L, 339L, 180L, 303L, 195L, 265L, 289L, 269L, 272L, 350L, 162L, 254L, 344L, 139L, 151L, 226L, 270L, 280L, 327L, 57L, 361L, 176L, 251L, 235L, 351L, 119L, 163L, 206L, 271L, 365L, 314L, 329L, 181L, 295L, 239L, 316L, 70L, 207L, 362L, 10L, 208L, 343L, 342L, 219L, 253L, 298L, 349L, 299L, 222L, 307L, 245L, 334L, 248L, 315L, 366L, 17L, 296L, 369L, 341L, 260L, 309L, 297L, 328L, 261L, 83L, 367L, 38L, 318L, 259L, 358L, 346L, 20L, 331L, 348L, 256L, 275L, 168L, 218L, 194L, 232L, 22L, 368L, 209L, 357L, 267L, 19L, 129L, 326L, 236L, 353L, 359L, 279L, 352L, 360L, 6L, 9L, 340L, 31L, 18L, 274L, 322L, 335L, 300L, 370L, 225L, 364L, 311L, 312L, 337L, 338L, 356L, 354L, 213L, 21L, 220L, 294L, 224L, 319L, 276L, 347L, 263L, 175L, 172L, 138L, 187L, 7L, 4L, 32L, 12L, 128L, 48L, 109L, 52L, 42L, 16L, 210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L, 66L, 293L, 268L, 205L, 41L, 43L, 2L, 34L, 85L, 193L, 63L, 55L, 131L, 54L, 67L, 178L, 103L, 93L, 56L, 79L, 133L, 80L, 264L, 161L, 114L, 288L, 28L, 113L, 88L, 97L, 160L, 68L, 58L, 246L, 25L, 141L, 94L, 150L, 284L, 60L, 278L, 26L, 29L, 72L, 61L, 49L, 40L, 287L, 98L, 75L, 136L, 45L, 144L, 5L, 11L, 124L, 233L, 273L, 3L, 125L, 152L, 137L, 71L, 46L, 27L, 107L, 59L, 92L, 44L, 51L, 157L, 89L, 134L, 126L, 84L, 82L, 145L, 53L, 111L, 214L, 153L, 74L, 99L, 257L, 306L, 105L, 108L, 116L, 277L, 143L, 230L, 64L, 171L, 122L, 91L, 282L, 135L, 203L, 30L, 147L, 156L, 155L, 35L, 76L, 23L, 185L, 262L, 247L, 123L, 227L, 258L, 238L, 223L, 73L, 130L, 81L, 154L, 104L, 291L, 216L, 112L, 333L, 165L, 212L, 308L, 201L, 313L, 324L, 148L, 321L, 170L, 211L, 102L, 198L, 286L, 229L, 249L, 167L, 159L, 191L, 283L, 320L, 305L, 317L, 355L, 330L, 182L, 221L, 231L, 241L, 166L, 234L, 142L, 199L, 217L, 117L, 121L, 132L, 174L, 266L, 215L, 325L, 65L, 301L, 106L, 179L, 50L, 33L, 115L, 8L, 304L, 345L, 363L, 169L, 36L, 292L, 110L, 47L, 192L, 95L, 101L, 69L, 202L, 146L, 200L, 173L, 100L, 39L, 302L, 336L, 290L, 332L, 86L, 127L, 197L, 310L, 242L, 118L, 90L, 120L, 13L, 285L, 183L, 188L, 158L, 252L, 164L, 177L, 190L, 186L, 237L, 62L, 78L, 240L, 250L, 196L, 228L, 149L, 87L, 243L, 24L, 14L, 204L, 140L, 15L, 244L, 339L, 180L, 303L, 195L, 265L, 289L, 269L, 272L, 350L, 162L, 254L, 344L, 139L, 151L, 226L, 270L, 280L, 327L, 57L, 361L, 176L, 251L, 235L, 351L, 119L, 163L, 206L, 271L, 365L, 314L, 329L, 181L, 295L, 239L, 316L, 70L, 207L, 362L, 10L, 208L, 343L, 342L, 219L, 253L, 298L, 349L, 299L, 222L, 307L, 245L, 334L, 248L, 315L, 366L, 17L, 296L, 369L, 341L, 260L, 309L, 297L, 328L, 261L, 83L, 367L, 38L, 318L, 259L, 358L, 346L, 20L, 331L, 348L, 256L, 275L, 168L, 218L, 194L, 232L, 22L, 368L, 209L, 357L, 267L, 19L, 129L, 326L, 236L, 353L, 359L, 279L, 352L, 360L, 6L, 9L, 340L, 31L, 18L, 274L, 322L, 335L, 300L, 370L, 225L, 364L, 311L, 312L, 337L, 338L, 356L, 354L, 213L, 21L, 220L, 294L, 224L, 319L, 276L, 347L, 263L, 175L, 172L, 138L, 187L, 7L, 4L, 32L, 12L, 128L, 48L, 109L, 52L, 42L, 16L, 210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L, 66L, 293L, 268L, 205L, 41L, 43L, 2L, 34L, 85L, 193L, 63L, 55L, 131L, 54L, 67L, 178L, 103L, 93L, 56L, 79L, 133L, 80L, 264L, 161L, 114L, 288L, 28L, 113L, 88L, 97L, 160L, 68L, 58L, 246L, 25L, 141L, 94L, 150L, 284L, 60L, 278L, 26L, 29L, 72L, 61L, 49L, 40L, 287L, 98L, 75L, 136L, 45L, 144L, 5L, 11L, 124L, 233L, 273L, 3L, 125L, 152L, 137L, 71L, 46L, 27L, 107L, 59L, 92L, 44L, 51L, 157L, 89L, 134L, 126L, 84L, 82L, 145L, 53L, 111L, 214L, 153L, 74L, 99L, 257L, 306L, 105L, 108L, 116L, 277L, 143L, 230L, 64L, 171L, 122L, 91L, 282L, 135L, 203L, 30L, 147L, 156L, 155L, 35L, 76L, 23L, 185L, 262L, 247L, 123L, 227L, 258L, 238L, 223L, 73L, 130L, 81L, 154L, 104L, 291L, 216L, 112L, 333L, 165L, 212L, 308L, 201L, 313L, 324L, 148L, 321L, 170L, 211L, 102L, 198L, 286L, 229L, 249L, 167L, 159L, 191L, 283L, 320L, 305L, 317L, 355L, 330L, 182L, 221L, 231L, 241L, 166L, 234L, 142L, 199L, 217L, 117L, 121L, 132L, 174L, 266L, 215L, 325L, 65L, 301L, 106L, 179L, 50L, 33L, 115L, 8L, 304L, 345L, 363L, 169L, 36L, 292L, 110L, 47L, 192L, 95L, 101L, 69L, 202L, 146L, 200L, 173L, 100L, 39L, 302L, 336L, 290L, 332L, 86L, 127L, 197L, 310L, 242L, 118L, 90L, 120L, 13L, 285L, 183L, 188L, 158L, 252L, 164L, 177L, 190L, 186L, 237L, 62L, 78L, 240L, 250L, 196L, 228L, 149L, 87L, 243L, 24L, 14L, 204L, 140L, 15L, 244L, 339L, 180L, 303L, 195L, 265L, 289L, 269L, 272L, 350L, 162L, 254L, 344L, 139L, 151L, 226L, 270L, 280L, 327L, 57L, 361L, 176L, 251L, 235L, 351L, 119L, 163L, 206L, 271L, 365L, 314L, 329L, 181L, 295L, 239L, 316L, 70L, 207L, 362L, 10L, 208L, 343L, 342L, 219L, 253L, 298L, 349L, 299L, 222L, 307L, 245L, 334L, 248L, 315L, 366L, 17L, 296L, 369L, 341L, 260L, 309L, 297L, 328L, 261L, 83L, 367L, 38L, 318L, 259L, 358L, 346L, 20L, 331L, 348L, 256L, 275L, 168L, 218L, 194L, 232L, 22L, 368L, 209L, 357L, 267L, 19L, 129L, 326L, 236L, 353L, 359L, 279L, 352L, 360L, 6L, 9L, 340L, 31L, 18L, 274L, 322L, 335L, 300L, 370L, 225L, 364L, 311L, 312L, 337L, 338L, 356L, 354L, 213L, 21L, 220L, 294L, 224L, 319L, 276L, 347L, 263L, 175L, 172L, 138L, 187L, 7L, 4L, 32L, 12L, 128L, 48L, 109L, 52L, 42L, 16L, 210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L, 66L, 293L, 268L, 205L, 41L, 43L, 2L, 34L, 85L, 193L, 63L, 55L, 131L, 54L, 67L, 178L, 103L, 93L, 56L, 79L, 133L, 80L, 264L, 161L, 114L, 288L, 28L, 113L, 88L, 97L, 160L, 68L, 58L, 246L, 25L, 141L, 94L, 150L, 284L, 60L, 278L, 26L, 29L, 72L, 61L, 49L, 40L, 287L, 98L, 75L, 136L, 45L, 144L, 5L, 11L, 124L, 233L, 273L, 3L, 125L, 152L, 137L, 71L, 46L, 27L, 107L, 59L, 92L, 44L, 51L, 157L, 89L, 134L, 126L, 84L, 82L, 145L, 53L, 111L, 214L, 153L, 74L, 99L, 257L, 306L, 105L, 108L, 116L, 277L, 143L, 230L, 64L, 171L, 122L, 91L, 282L, 135L, 203L, 30L, 147L, 156L, 155L, 35L, 76L, 23L, 185L, 262L, 247L, 123L, 227L, 258L, 238L, 223L, 73L, 130L, 81L, 154L, 104L, 291L, 216L, 112L, 333L, 165L, 212L, 308L, 201L, 313L, 324L, 148L, 321L, 170L, 211L, 102L, 198L, 286L, 229L, 249L, 167L, 159L, 191L, 283L, 320L, 305L, 317L, 355L, 330L, 182L, 221L, 231L, 241L, 166L, 234L, 142L, 199L, 217L, 117L, 121L, 132L, 174L, 266L, 215L, 325L, 65L, 301L, 106L, 179L, 50L, 33L, 115L, 8L, 304L, 345L, 363L, 169L, 36L, 292L, 110L, 47L, 192L, 95L, 101L, 69L, 202L, 146L, 200L, 173L, 100L, 39L, 302L, 336L, 290L, 332L, 86L, 127L, 197L, 310L, 242L, 118L, 90L, 120L, 13L, 285L, 183L, 188L, 158L, 252L, 164L, 177L, 190L, 186L, 237L, 62L, 78L, 240L, 250L, 196L, 228L, 149L, 87L, 243L, 24L, 14L, 204L, 140L, 15L, 244L, 339L, 180L, 303L, 195L, 265L, 289L, 269L, 272L, 350L, 162L, 254L, 344L, 139L, 151L, 226L, 270L, 280L, 327L, 57L, 361L, 176L, 251L, 235L, 351L, 119L, 163L, 206L, 271L, 365L, 314L, 329L, 181L, 295L, 239L, 316L, 70L, 207L, 362L, 10L, 208L, 343L, 342L, 219L, 253L, 298L, 349L, 299L, 222L, 307L, 245L, 334L, 248L, 315L, 366L, 17L, 296L, 369L, 341L, 260L, 309L, 297L, 328L, 261L, 83L, 367L, 38L, 318L, 259L, 358L, 346L, 20L, 331L, 348L, 256L, 275L, 168L, 218L, 194L, 232L, 22L, 368L, 209L, 357L, 267L, 19L, 129L, 326L, 236L, 353L, 359L, 279L, 352L, 360L, 6L, 9L, 340L, 31L, 18L, 274L, 322L, 335L, 300L, 370L, 225L, 364L, 311L, 312L, 337L, 338L, 356L, 354L, 213L, 21L, 220L, 294L, 224L, 319L, 276L, 347L, 263L, 175L, 172L, 138L, 187L, 7L, 4L, 32L, 12L, 128L, 48L, 109L, 52L, 42L, 16L, 210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L, 66L, 293L, 268L, 205L, 41L, 43L, 2L, 34L, 85L, 193L, 63L, 55L, 131L, 54L, 67L, 178L, 103L, 93L, 56L, 79L, 133L, 80L, 264L, 161L, 114L, 288L, 28L, 113L, 88L, 97L, 160L, 68L, 58L, 246L, 25L, 141L, 94L, 150L, 284L, 60L, 278L, 26L, 29L, 72L, 61L, 49L, 40L, 287L, 98L, 75L, 136L, 45L, 144L, 5L, 11L, 124L, 233L, 273L, 3L, 125L, 152L, 137L, 71L, 46L, 27L, 107L, 59L, 92L, 44L, 51L, 157L, 89L, 134L, 126L, 84L, 82L, 145L, 53L, 111L, 214L, 153L, 74L, 99L, 257L, 306L, 105L, 108L, 116L, 277L, 143L, 230L, 64L, 171L, 122L, 91L, 282L, 135L, 203L, 30L, 147L, 156L, 155L, 35L, 76L, 23L, 185L, 262L, 247L, 123L, 227L, 258L, 238L, 223L, 73L, 130L, 81L, 154L, 104L, 291L, 216L, 112L, 333L, 165L, 212L, 308L, 201L, 313L, 324L, 148L, 321L, 170L, 211L, 102L, 198L, 286L, 229L, 249L, 167L, 159L, 191L, 283L, 320L, 305L, 317L, 355L, 330L, 182L, 221L, 231L, 241L, 166L, 234L, 142L, 199L, 217L, 117L, 121L, 132L, 174L, 266L, 215L, 325L, 65L, 301L, 106L, 179L, 50L, 33L, 115L, 8L, 304L, 345L, 363L, 169L, 36L, 292L, 110L, 47L, 192L, 95L, 101L, 69L, 202L, 146L, 200L, 173L, 100L, 39L, 302L, 336L, 290L, 332L, 86L, 127L, 197L, 310L, 242L, 118L, 90L, 120L, 13L, 285L, 183L, 188L, 158L, 252L, 164L, 177L, 190L, 186L, 237L, 62L, 78L, 240L, 250L, 196L, 228L, 149L, 87L, 243L, 24L, 14L, 204L, 140L, 15L, 244L, 339L, 180L, 303L, 195L, 265L, 289L, 269L, 272L, 350L, 162L, 254L, 344L, 139L, 151L, 226L, 270L, 280L, 327L, 57L, 361L, 176L, 251L, 235L, 351L, 119L, 163L, 206L, 271L, 365L, 314L, 329L, 181L, 295L, 239L, 316L, 70L, 207L, 362L, 10L, 208L, 343L, 342L, 219L, 253L, 298L, 349L, 299L, 222L, 307L, 245L, 334L, 248L, 315L, 366L, 17L, 296L, 369L, 341L, 260L, 309L, 297L, 328L, 261L, 83L, 367L, 38L, 318L, 259L, 358L, 346L, 20L, 331L, 348L, 256L, 275L, 168L, 218L, 194L, 232L, 22L, 368L, 209L, 357L, 267L, 19L, 129L, 326L, 236L, 353L, 359L, 279L, 352L, 360L, 6L, 9L, 340L, 31L, 18L, 274L, 322L, 335L, 300L, 370L, 225L, 364L, 311L, 312L, 337L, 338L, 356L, 354L, 213L, 21L, 220L, 294L, 224L, 319L, 276L, 347L, 263L, 175L, 172L, 138L, 187L, 7L, 4L, 32L, 12L, 128L, 48L, 109L, 52L, 42L, 16L, 210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L, 66L, 293L, 268L, 205L, 41L, 43L, 2L, 34L, 85L, 193L, 63L, 55L, 131L, 54L, 67L, 178L, 103L, 93L, 56L, 79L, 133L, 80L, 264L, 161L, 114L, 288L, 28L, 113L, 88L, 97L, 160L, 68L, 58L, 246L, 25L, 141L, 94L, 150L, 284L, 60L, 278L, 26L, 29L, 72L, 61L, 49L, 40L, 287L, 98L, 75L, 136L, 45L, 144L, 5L, 11L, 124L, 233L, 273L, 3L, 125L, 152L, 137L, 71L, 46L, 27L, 107L, 59L, 92L, 44L, 51L, 157L, 89L, 134L, 126L, 84L, 82L, 145L, 53L, 111L, 214L, 153L, 74L, 99L, 257L, 306L, 105L, 108L, 116L, 277L, 143L, 230L, 64L, 171L, 122L, 91L, 282L, 135L, 203L, 30L, 147L, 156L, 155L, 35L, 76L, 23L, 185L, 262L, 247L, 123L, 227L, 258L, 238L, 223L, 73L, 130L, 81L, 154L, 104L, 291L, 216L, 112L, 333L, 165L, 212L, 308L, 201L, 313L, 324L, 148L, 321L, 170L, 211L, 102L, 198L, 286L, 229L, 249L, 167L, 159L, 191L, 283L, 320L, 305L, 317L, 355L, 330L, 182L, 221L, 231L, 241L, 166L, 234L, 142L, 199L, 217L, 117L, 121L, 132L, 174L, 266L, 215L, 325L, 65L, 301L, 106L, 179L, 50L, 33L, 115L, 8L, 304L, 345L, 363L, 169L, 36L, 292L, 110L, 47L, 192L, 95L, 101L, 69L, 202L, 146L, 200L, 173L, 100L, 39L, 302L, 336L, 290L, 332L, 86L, 127L, 197L, 310L, 242L, 118L, 90L, 120L, 13L, 285L, 183L, 188L, 158L, 252L, 164L, 177L, 190L, 186L, 237L, 62L, 78L, 240L, 250L, 196L, 228L, 149L, 87L, 243L, 24L, 14L, 204L, 140L, 15L, 244L, 339L, 180L, 303L, 195L, 265L, 289L, 269L, 272L, 350L, 162L, 254L, 344L, 139L, 151L, 226L, 270L, 280L, 327L, 57L, 361L, 176L, 251L, 235L, 351L, 119L, 163L, 206L, 271L, 365L, 314L, 329L, 181L, 295L, 239L, 316L, 70L, 207L, 362L, 10L, 208L, 343L, 342L, 219L, 253L, 298L, 349L, 299L, 222L, 307L, 245L, 334L, 248L, 315L, 366L, 17L, 296L, 369L, 341L, 260L, 309L, 297L, 328L, 261L, 83L, 367L, 38L, 318L, 259L, 358L, 346L, 20L, 331L, 348L, 256L, 275L, 168L, 218L, 194L, 232L, 22L, 368L, 209L, 357L, 267L, 19L, 129L, 326L, 236L, 353L, 359L, 279L, 352L, 360L, 6L, 9L, 340L, 31L, 18L, 274L, 322L, 335L, 300L, 370L, 225L, 364L, 311L, 312L, 337L, 338L, 356L, 354L, 213L, 21L, 220L, 294L, 224L, 319L, 276L, 347L, 263L, 175L, 172L, 138L, 187L, 7L, 4L, 32L, 12L), .Label = c("12", "23", "75", "368", "70", "336", "367", "186", "337", "281", "71", "370", "217", "238", "241", "6", "297", "340", "327", "313", "355", "322", "121", "237", "51", "58", "81", "43", "59", "115", "339", "369", "184", "24", "119", "191", "8", "308", "204", "63", "21", "5", "22", "85", "68", "80", "194", "2", "62", "183", "86", "4", "94", "30", "28", "35", "261", "49", "83", "56", "61", "228", "27", "108", "179", "17", "31", "48", "198", "278", "79", "60", "130", "98", "66", "120", "15", "229", "36", "38", "132", "92", "306", "91", "25", "209", "235", "45", "88", "215", "111", "84", "34", "53", "196", "10", "46", "65", "99", "203", "197", "149", "33", "134", "102", "181", "82", "103", "3", "193", "95", "137", "44", "41", "185", "104", "172", "214", "267", "216", "173", "110", "125", "72", "76", "90", "210", "1", "328", "131", "29", "174", "37", "89", "113", "67", "78", "365", "255", "240", "52", "169", "106", "69", "93", "200", "116", "145", "234", "54", "256", "77", "97", "133", "118", "117", "87", "221", "155", "47", "40", "252", "268", "223", "139", "167", "154", "318", "190", "147", "109", "364", "202", "175", "363", "263", "224", "32", "182", "244", "274", "163", "219", "16", "122", "226", "366", "220", "13", "225", "156", "195", "26", "320", "246", "232", "211", "150", "170", "201", "142", "199", "114", "239", "20", "269", "279", "282", "324", "7", "148", "140", "354", "96", "177", "136", "171", "319", "285", "356", "164", "290", "129", "358", "346", "257", "126", "233", "152", "107", "165", "321", "73", "168", "265", "330", "227", "128", "276", "230", "166", "213", "236", "242", "292", "50", "124", "294", "153", "231", "264", "222", "286", "253", "11", "316", "100", "127", "310", "301", "305", "123", "362", "39", "247", "176", "326", "19", "249", "258", "270", "250", "74", "341", "317", "360", "105", "57", "333", "259", "9", "112", "157", "55", "218", "151", "64", "42", "248", "207", "135", "192", "18", "357", "275", "298", "303", "287", "289", "344", "180", "205", "245", "187", "159", "101", "291", "141", "302", "212", "348", "349", "143", "272", "295", "277", "160", "309", "359", "158", "146", "342", "14", "144", "178", "329", "260", "304", "273", "162", "314", "208", "138", "293", "343", "206", "350", "351", "243", "338", "300", "284", "283", "254", "188", "312", "361", "315", "288", "251", "266", "334", "331", "353", "161", "352", "325", "311", "332", "335", "262", "280", "189", "347", "271", "296", "307", "323", "299", "345"), class = c("ordered", "factor")), TRAT = structure(c(3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("DSM", "NUT", "POL"), class = "factor"), SEMANA = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6), PESO = c(0.0395, 0.0431, 0.046, 0.043, 0.043, 0.0436, 0.042, 0.0405, 0.0417, 0.045, 0.043, 0.0416, 0.041, 0.042, 0.0445, 0.044, 0.047, 0.043, 0.049, 0.0435, 0.044, 0.0421, 0.0444, 0.041, 0.042, 0.0415, 0.0411, 0.0422, 0.042, 0.042, 0.045, 0.045, 0.0416, 0.044, 0.044, 0.041, 0.04, 0.039, 0.0408, 0.037, 0.044, 0.042, 0.041, 0.04, 0.0435, 0.039, 0.0404, 0.0466, 0.0461, 0.0435, 0.043, 0.043, 0.041, 0.0432, 0.0393, 0.0415, 0.041, 0.041, 0.042, 0.0404, 0.0385, 0.0435, 0.034, 0.0409, 0.044, 0.045, 0.045, 0.046, 0.045, 0.0396, 0.0406, 0.039, 0.038, 0.039, 0.045, 0.042, 0.042, 0.0413, 0.043, 0.042, 0.0385, 0.0415, 0.0391, 0.048, 0.043, 0.039, 0.048, 0.0388, 0.039, 0.043, 0.0436, 0.043, 0.0435, 0.0395, 0.0424, 0.0414, 0.0422, 0.0394, 0.0405, 0.041, 0.0373, 0.0466, 0.045, 0.0431, 0.0435, 0.047, 0.044, 0.042, 0.041, 0.0409, 0.042, 0.042, 0.043, 0.0375, 0.043, 0.038, 0.0402, 0.0375, 0.043, 0.044, 0.043, 0.0415, 0.0445, 0.044, 0.0445, 0.0413, 0.041, 0.043, 0.032, 0.044, 0.0392, 0.043, 0.0437, 0.034, 0.038, 0.0386, 0.037, 0.044, 0.046, 0.0445, 0.0425, 0.045, 0.0388, 0.0387, 0.039, 0.0386, 0.0421, 0.0405, 0.0412, 0.042, 0.0471, 0.042, 0.047, 0.039, 0.039, 0.043, 0.041, 0.044, 0.044, 0.042, 0.042, 0.042, 0.0422, 0.047, 0.047, 0.0415, 0.045, 0.045, 0.044, 0.0426, 0.0405, 0.044, 0.0446, 0.04, 0.0445, 0.0395, 0.0444, 0.046, 0.0465, 0.0415, 0.042, 0.042, 0.0481, 0.047, 0.042, 0.041, 0.041, 0.04, 0.0451, 0.043, 0.0415, 0.0432, 0.0414, 0.0445, 0.0427, 0.041, 0.04, 0.047, 0.043, 0.043, 0.0404, 0.0394, 0.043, 0.043, 0.04, 0.042, 0.0424, 0.042, 0.046, 0.042, 0.04, 0.0404, 0.0417, 0.0404, 0.0442, 0.041, 0.0405, 0.0426, 0.0433, 0.044, 0.0465, 0.0405, 0.0451, 0.0446, 0.043, 0.04, 0.041, 0.0454, 0.0449, 0.0433, 0.043, 0.0424, 0.0439, 0.044, 0.0452, 0.04, 0.04, 0.0416, 0.0405, 0.04, 0.044, 0.043, 0.045, 0.043, 0.043, 0.0424, 0.041, 0.0464, 0.0447, 0.0426, 0.0426, 0.042, 0.0421, 0.045, 0.0379, 0.044, 0.042, 0.042, 0.042, 0.0433, 0.034, 0.037, 0.037, 0.0445, 0.042, 0.042, 0.044, 0.045, 0.042, 0.0416, 0.0436, 0.0405, 0.0426, 0.044, 0.044, 0.0404, 0.039, 0.0435, 0.0469, 0.0435, 0.045, 0.043, 0.043, 0.038, 0.0438, 0.0433, 0.0404, 0.0435, 0.043, 0.043, 0.038, 0.041, 0.042, 0.04, 0.037, 0.038, 0.0417, 0.043, 0.047, 0.0449, 0.043, 0.04, 0.04, 0.04, 0.0421, 0.0442, 0.043, 0.041, 0.041, 0.0414, 0.0435, 0.045, 0.044, 0.044, 0.0425, 0.0465, 0.043, 0.0446, 0.045, 0.043, 0.046, 0.0423, 0.045, 0.047, 0.0464, 0.0444, 0.0446, 0.044, 0.0438, 0.044, 0.0484, 0.043, 0.0434, 0.0395, 0.0445, 0.037, 0.0395, 0.039, 0.039, 0.043, 0.0409, 0.0441, 0.0452, 0.044, 0.044, 0.0458, 0.043, 0.044, 0.04, 0.0408, 0.044, 0.044, 0.041, 0.0439, 0.045, 0.0416, 0.0447, 0.0443, 0.044, 0.042, 0.0447, 0.043, 0.044, 0.0408, 0.044, 0.044, 0.046, 0.0435, 0.047, 0.0417, 0.177, 0.1757, 0.189, 0.185, 0.1836, 0.1969, 0.1746, 0.1561, 0.1757, 0.176, 0.1957, 0.1555, 0.1684, 0.171, 0.1771, 0.1929, 0.1927, 0.185, 0.1935, 0.1883, 0.1837, 0.1763, 0.1669, 0.172, 0.172, 0.1699, 0.1722, 0.1868, 0.1836, 0.137, 0.1565, 0.1765, 0.1655, 0.18, 0.1682, 0.17, 0.1897, 0.1746, 0.1869, 0.17, 0.195, 0.207, 0.1388, 0.1815, 0.1868, 0.154, 0.1514, 0.1862, 0.178, 0.182, 0.166, 0.181, 0.146, 0.1704, 0.1541, 0.1749, 0.1976, 0.1671, 0.1706, 0.1708, 0.1597, 0.1731, 0.1465, 0.1764, 0.1647, 0.1647, 0.2056, 0.1812, 0.1746, 0.1651, 0.1705, 0.18, 0.163, 0.18, 0.183, 0.163, 0.1859, 0.1697, 0.1895, 0.1737, 0.1577, 0.1716, 0.1605, 0.1879, 0.1814, 0.1578, 0.1684, 0.157, 0.1619, 0.1793, 0.1828, 0.1705, 0.184, 0.1516, 0.1754, 0.1818, 0.1984, 0.1711, 0.1765, 0.2006, 0.1728, 0.1843, 0.178, 0.1707, 0.1761, 0.175, 0.1882, 0.1915, 0.1719, 0.1618, 0.165, 0.153, 0.167, 0.184, 0.165, 0.17, 0.1745, 0.1825, 0.1461, 0.2103, 0.1687, 0.17, 0.183, 0.18, 0.1683, 0.1898, 0.1731, 0.1827, 0.1693, 0.1859, 0.1723, 0.2012, 0.1988, 0.1447, 0.1898, 0.1522, 0.1659, 0.178, 0.18, 0.1824, 0.1778, 0.1824, 0.1854, 0.181, 0.1684, 0.1825, 0.1898, 0.186, 0.1813, 0.1852, 0.189, 0.1812, 0.1847, 0.167, 0.1523, 0.1881, 0.1753, 0.2212, 0.204, 0.205, 0.203, 0.203, 0.1786, 0.1906, 0.1907, 0.1894, 0.186, 0.18, 0.197, 0.1349, 0.1805, 0.1968, 0.1885, 0.1794, 0.1953, 0.1835, 0.2037, 0.21, 0.1988, 0.1806, 0.1847, 0.1941, 0.1857, 0.182, 0.1814, 0.178, 0.19, 0.179, 0.1911, 0.1925, 0.1763, 0.1927, 0.207, 0.183, 0.1831, 0.188, 0.1753, 0.2071, 0.1976, 0.1941, 0.1688, 0.1778, 0.1941, 0.1653, 0.19, 0.194, 0.1889, 0.1835, 0.186, 0.1812, 0.1779, 0.1744, 0.1843, 0.1808, 0.1828, 0.1701, 0.193, 0.1987, 0.1915, 0.192, 0.1966, 0.1907, 0.1963, 0.1907, 0.1855, 0.1682, 0.165, 0.1871, 0.1908, 0.183, 0.182, 0.1809, 0.1806, 0.1765, 0.1814, 0.1913, 0.1877, 0.1902, 0.1867, 0.1755, 0.196, 0.1882, 0.1984, 0.1874, 0.1835, 0.1774, 0.181, 0.188, 0.19, 0.176, 0.1825, 0.1844, 0.1873, 0.1843, 0.1694, 0.1871, 0.1706, 0.1729, 0.1794, 0.1836, 0.1459, 0.1565, 0.1565, 0.178, 0.176, 0.17, 0.1682, 0.179, 0.173, 0.1812, 0.1915, 0.172, 0.195, 0.2, 0.19, 0.1972, 0.17, 0.1772, 0.2012, 0.1857, 0.1959, 0.1834, 0.197, 0.1647, 0.1796, 0.1798, 0.1818, 0.1986, 0.2089, 0.142, 0.1882, 0.2059, 0.1882, 0.1871, 0.1706, 0.182, 0.2036, 0.196, 0.21, 0.2142, 0.1871, 0.1781, 0.1835, 0.174, 0.1831, 0.17, 0.174, 0.198, 0.186, 0.1808, 0.1983, 0.192, 0.1894, 0.1859, 0.184, 0.1847, 0.1694, 0.1734, 0.184, 0.1847, 0.1682, 0.1843, 0.1848, 0.181, 0.1915, 0.1823, 0.1796, 0.1953, 0.2031, 0.1824, 0.1882, 0.18, 0.1837, 0.1747, 0.1928, 0.1765, 0.1714, 0.1666, 0.1724, 0.1918, 0.1894, 0.1864, 0.1936, 0.1813, 0.199, 0.1896, 0.1812, 0.186, 0.1714, 0.167, 0.1824, 0.194, 0.185, 0.189, 0.2, 0.1997, 0.186, 0.1961, 0.1827, 0.1806, 0.1825, 0.18, 0.1789, 0.173, 0.1616, 0.1776, 0.191, 0.1788, 0.1976, 0.1767, 0.4056, 0.4401, 0.4615, 0.44, 0.4698, 0.4991, 0.441, 0.3404, 0.4144, 0.3682, 0.5132, 0.4332, 0.4246, 0.4535, 0.4288, 0.4961, 0.4296, 0.4883, 0.4836, 0.4675, 0.4524, 0.4419, 0.4568, 0.451, 0.473, 0.4515, 0.4702, 0.4663, 0.4743, 0.3246, 0.3918, 0.4387, 0.3816, 0.4083, 0.4273, 0.444, 0.4702, 0.441, 0.4843, 0.414, 0.4659, 0.5444, 0.31, 0.4602, 0.428, 0.4058, 0.4367, 0.4797, 0.4287, 0.4827, 0.332, 0.45, 0.334, 0.3776, 0.3914, 0.3915, 0.477, 0.348, 0.3536, 0.4312, 0.4031, 0.4253, 0.3677, 0.4695, 0.3738, 0.4018, 0.5217, 0.44, 0.4714, 0.3999, 0.4478, 0.4517, 0.404, 0.424, 0.3885, 0.4125, 0.5118, 0.3975, 0.4465, 0.4021, 0.3555, 0.3804, 0.3547, 0.4518, 0.391, 0.3592, 0.4066, 0.3966, 0.3888, 0.4716, 0.419, 0.4565, 0.4658, 0.403, 0.4524, 0.509, 0.516, 0.434, 0.433, 0.5141, 0.4718, 0.46, 0.491, 0.4407, 0.4504, 0.437, 0.4587, 0.4484, 0.4457, 0.4063, 0.461, 0.395, 0.391, 0.4757, 0.394, 0.4572, 0.4405, 0.4794, 0.3574, 0.5332, 0.3565, 0.465, 0.4525, 0.473, 0.4234, 0.5006, 0.4651, 0.4479, 0.431, 0.4939, 0.4354, 0.4916, 0.499, 0.4018, 0.4931, 0.4217, 0.4523, 0.4505, 0.44, 0.4748, 0.5106, 0.4714, 0.4307, 0.4692, 0.42, 0.5005, 0.4712, 0.4823, 0.4393, 0.4855, 0.4759, 0.4692, 0.4613, 0.4336, 0.4354, 0.4852, 0.4164, 0.5534, 0.5085, 0.5197, 0.5264, 0.5085, 0.4542, 0.4702, 0.47, 0.4701, 0.4565, 0.4972, 0.456, 0.3752, 0.449, 0.4517, 0.4309, 0.4299, 0.4839, 0.5006, 0.4957, 0.5399, 0.4951, 0.4752, 0.483, 0.496, 0.42, 0.445, 0.4378, 0.44, 0.534, 0.483, 0.5379, 0.5157, 0.4049, 0.534, 0.5041, 0.3965, 0.4214, 0.4895, 0.3962, 0.4961, 0.4786, 0.5021, 0.4374, 0.4956, 0.4731, 0.3408, 0.537, 0.4961, 0.4773, 0.5309, 0.4827, 0.416, 0.47, 0.4505, 0.4698, 0.46, 0.4632, 0.4726, 0.467, 0.5012, 0.4719, 0.445, 0.4853, 0.4754, 0.5109, 0.469, 0.4735, 0.416, 0.408, 0.4676, 0.4745, 0.4644, 0.4546, 0.5077, 0.4571, 0.4419, 0.425, 0.4886, 0.4662, 0.4973, 0.4842, 0.447, 0.522, 0.4971, 0.5032, 0.4757, 0.461, 0.467, 0.458, 0.5102, 0.4778, 0.4804, 0.5032, 0.4724, 0.4802, 0.5519, 0.4094, 0.4445, 0.449, 0.4647, 0.4809, 0.4684, 0.3682, 0.4456, 0.3842, 0.48, 0.4378, 0.4489, 0.4164, 0.461, 0.4201, 0.5157, 0.4727, 0.4375, 0.5038, 0.515, 0.499, 0.4936, 0.42, 0.4748, 0.5001, 0.4787, 0.4331, 0.4426, 0.482, 0.431, 0.4612, 0.4678, 0.4659, 0.48, 0.5247, 0.3512, 0.4568, 0.4905, 0.4669, 0.4714, 0.452, 0.4653, 0.4992, 0.492, 0.489, 0.5304, 0.501, 0.4964, 0.4658, 0.4897, 0.4783, 0.3947, 0.481, 0.501, 0.471, 0.4449, 0.4865, 0.4872, 0.4961, 0.4624, 0.4775, 0.4414, 0.44, 0.4343, 0.449, 0.4601, 0.44, 0.4997, 0.4974, 0.486, 0.509, 0.4702, 0.4749, 0.4434, 0.5167, 0.4438, 0.4847, 0.4539, 0.4285, 0.4773, 0.4846, 0.4197, 0.458, 0.453, 0.438, 0.4939, 0.4839, 0.4589, 0.473, 0.4731, 0.506, 0.4954, 0.4714, 0.4536, 0.4439, 0.4431, 0.4818, 0.496, 0.5108, 0.4134, 0.5051, 0.5247, 0.4821, 0.4583, 0.475, 0.4485, 0.4692, 0.4501, 0.4639, 0.44, 0.4004, 0.4501, 0.4843, 0.4606, 0.5253, 0.4429, 0.8378, 0.8471, 0.8886, 0.81, 0.9366, 0.9909, 0.8597, 0.6937, 0.8395, 0.832, 1.009, 0.8833, 0.7511, 0.8965, 0.8922, 0.96, 0.7807, 1.002, 0.963, 0.962, 0.8641, 0.8847, 0.8782, 0.81, 0.902, 0.7963, 0.8896, 0.8852, 0.9848, 0.737, 0.8543, 0.9315, 0.8005, 0.8761, 0.9185, 0.85, 0.9096, 0.815, 0.9353, 0.913, 0.8832, 1.0239, 0.634, 0.8885, 0.734, 0.8875, 0.8878, 0.9387, 0.8309, 0.9772, 0.718, 0.8574, 0.727, 0.843, 0.8819, 0.8765, 0.9533, 0.6587, 0.6609, 0.8953, 0.8873, 0.7675, 0.7517, 0.9561, 0.8424, 0.7728, 1.0431, 0.805, 0.9166, 0.8258, 0.9308, 0.8518, 0.8546, 0.8185, 0.7579, 0.8505, 0.9554, 0.8116, 0.8611, 0.7879, 0.775, 0.771, 0.7706, 0.8675, 0.7561, 0.7585, 0.8196, 0.7993, 0.842, 0.9372, 0.8446, 0.9328, 0.9495, 0.7644, 0.9246, 0.9592, 0.9594, 0.842, 0.9085, 0.986, 0.9085, 0.8866, 0.7695, 0.9027, 0.9104, 0.9304, 0.8761, 0.8988, 0.9277, 0.8332, 0.861, 0.776, 0.844, 0.97, 0.8026, 0.895, 0.8928, 0.918, 0.7535, 0.9975, 0.7307, 0.9997, 0.9169, 0.9444, 0.8713, 0.9523, 0.8918, 0.9314, 0.9233, 0.9, 0.8646, 0.9181, 0.9478, 0.7761, 1.0358, 0.8917, 0.8, 0.888, 0.8283, 0.9174, 0.9131, 0.927, 0.8629, 0.9495, 0.7656, 1.0016, 0.9123, 0.915, 0.8805, 0.9568, 0.9769, 0.9674, 0.9783, 0.901, 0.8097, 0.9665, 0.7767, 1.0761, 0.998, 0.988, 1.004, 0.987, 0.9206, 0.945, 0.91, 0.9553, 0.961, 0.975, 0.95, 0.8542, 0.9019, 0.9331, 0.8521, 0.8897, 0.9847, 1.0218, 0.9577, 1.0783, 0.9665, 0.9721, 0.937, 0.91, 0.8823, 0.86, 0.8604, 0.838, 1, 1, 1.0924, 0.991, 0.816, 1.0692, 0.9754, 0.847, 0.8324, 0.9739, 0.8891, 0.9609, 0.9596, 1.013, 0.8791, 0.9591, 1.0097, 0.7858, 0.994, 0.9978, 1.0076, 1.0848, 0.9174, 0.8936, 0.9353, 0.94, 0.9573, 0.8629, 0.9516, 0.903, 0.93, 0.979, 0.9293, 0.812, 0.9387, 0.9252, 0.9476, 0.9151, 0.9239, 0.8873, 0.889, 0.9308, 0.9984, 0.8824, 0.936, 0.9941, 0.962, 0.8414, 0.854, 0.9367, 0.9601, 0.9998, 0.9538, 0.904, 1.016, 0.9036, 0.977, 0.9141, 0.8887, 0.9664, 0.9467, 1.0024, 0.9383, 1.0171, 0.9536, 0.9195, 0.9381, 0.9906, 0.8393, 0.8997, 0.95, 0.9793, 1.0044, 1.0127, 0.7932, 0.9554, 0.7272, 0.9685, 0.8804, 0.9098, 0.865, 0.896, 0.8813, 0.985, 0.9001, 0.8505, 0.973, 1.0228, 1.0326, 0.9605, 0.9089, 0.969, 0.9468, 0.9253, 0.8846, 0.9015, 0.9577, 0.8712, 0.9422, 0.8845, 0.8767, 0.9465, 1.0106, 0.7609, 0.8837, 0.963, 0.8652, 0.9033, 0.958, 0.9544, 0.9772, 1.007, 0.991, 1.0612, 0.967, 0.9619, 0.8641, 0.976, 0.897, 0.8328, 0.98, 1.024, 0.95, 0.9373, 0.9992, 1.0294, 1.0174, 1, 0.9663, 0.9136, 0.875, 0.8647, 0.8879, 0.925, 0.9098, 1.0355, 1.0423, 0.974, 1.0475, 0.9684, 0.9961, 0.8761, 1.043, 0.906, 1.0113, 0.9554, 0.8868, 0.9672, 0.9689, 0.7913, 0.9557, 0.95, 0.9241, 0.977, 0.9954, 0.9414, 0.9768, 0.9361, 1.012, 0.9797, 0.935, 0.9789, 0.9177, 0.9212, 0.9467, 1.0138, 0.996, 0.844, 0.9182, 0.988, 0.9479, 0.9296, 1.0091, 0.9145, 0.9084, 0.9838, 0.931, 0.8625, 0.8116, 0.975, 0.965, 0.8768, 0.996, 0.8978, 1.3968, 1.3513, 1.3796, 1.37, 1.511, 1.5374, 1.3402, 1.2161, 1.45, 1.3467, 1.6014, 1.4685, 1.4193, 1.5277, 1.404, 1.5426, 1.3524, 1.619, 1.5212, 1.4904, 1.3892, 1.3864, 1.4216, 1.298, 1.353, 1.3753, 1.4451, 1.411, 1.5489, 1.384, 1.4021, 1.3734, 1.3706, 1.2957, 1.3287, 1.3739, 1.4868, 1.3398, 1.4879, 1.376, 1.4362, 1.6883, 1.169, 1.438, 1.1338, 1.4425, 1.461, 1.4077, 1.4075, 1.5675, 1.2991, 1.361, 1.2582, 1.3556, 1.4687, 1.3709, 1.479, 1.2, 1.171, 1.4051, 1.3693, 1.3885, 1.3117, 1.4925, 1.382, 1.307, 1.5479, 1.3287, 1.463, 1.3607, 1.4868, 1.4895, 1.471, 1.478, 1.322, 1.3986, 1.5283, 1.3912, 1.3905, 1.2882, 1.2838, 1.3064, 1.2657, 1.4752, 1.33, 1.2807, 1.3757, 1.3521, 1.3615, 1.4647, 1.3665, 1.3926, 1.5306, 1.3466, 1.4485, 1.5256, 1.4598, 1.392, 1.4577, 1.569, 1.5389, 1.4609, 1.366, 1.4558, 1.4793, 1.483, 1.4778, 1.3583, 1.4369, 1.3851, 1.3002, 1.343, 1.396, 1.6054, 1.2581, 1.4681, 1.4902, 1.516, 1.2972, 1.5183, 1.1763, 1.5713, 1.5658, 1.5176, 1.3958, 1.5345, 1.4831, 1.5042, 1.475, 1.4894, 1.3109, 1.4947, 1.4008, 1.3924, 1.6255, 1.4301, 1.3511, 1.547, 1.417, 1.4545, 1.427, 1.575, 1.4619, 1.5714, 1.395, 1.6044, 1.5129, 1.4793, 1.4363, 1.5361, 1.5828, 1.4649, 1.5128, 1.4554, 1.4241, 1.5933, 1.487, 1.6123, 1.515, 1.661, 1.614, 1.604, 1.5179, 1.5463, 1.5053, 1.5174, 1.563, 1.5646, 1.529, 1.4507, 1.521, 1.3927, 1.4156, 1.4218, 1.4786, 1.4873, 1.528, 1.5404, 1.6069, 1.529, 1.524, 1.532, 1.3764, 1.365, 1.417, 1.42, 1.5018, 1.5171, 1.6861, 1.4826, 1.3477, 1.5042, 1.497, 1.3495, 1.5278, 1.5128, 1.4574, 1.5181, 1.5579, 1.55, 1.5066, 1.5624, 1.5097, 1.2975, 1.574, 1.6447, 1.5551, 1.6489, 1.575, 1.3699, 1.4815, 1.5128, 1.52, 1.4599, 1.5155, 1.5009, 1.525, 1.5682, 1.5374, 1.4216, 1.4754, 1.4921, 1.6045, 1.5672, 1.5415, 1.4147, 1.5007, 1.3903, 1.4811, 1.4684, 1.583, 1.5141, 1.5613, 1.4138, 1.3942, 1.5426, 1.5353, 1.5059, 1.5185, 1.4718, 1.557, 1.5503, 1.5809, 1.4571, 1.3876, 1.4999, 1.567, 1.6213, 1.4722, 1.6299, 1.5003, 1.5016, 1.3979, 1.5145, 1.368, 1.511, 1.419, 1.499, 1.4829, 1.5283, 1.3489, 1.5275, 1.2543, 1.5618, 1.4681, 1.475, 1.3351, 1.382, 1.5259, 1.5143, 1.5396, 1.4435, 1.6217, 1.5953, 1.45, 1.6084, 1.516, 1.4597, 1.3616, 1.595, 1.4594, 1.5034, 1.5266, 1.5228, 1.512, 1.5051, 1.4695, 1.5701, 1.6411, 1.3285, 1.4851, 1.483, 1.5085, 1.4766, 1.6053, 1.5442, 1.5624, 1.481, 1.56, 1.6188, 1.6365, 1.4982, 1.3872, 1.5787, 1.5001, 1.4118, 1.5, 1.664, 1.538, 1.5796, 1.6087, 1.6537, 1.5894, 1.5787, 1.5931, 1.5452, 1.4897, 1.4553, 1.5442, 1.5851, 1.5819, 1.5899, 1.6965, 1.6103, 1.6574, 1.5121, 1.5663, 1.474, 1.6411, 1.545, 1.68, 1.616, 1.5187, 1.5035, 1.5558, 1.453, 1.5854, 1.607, 1.6486, 1.6085, 1.5506, 1.5697, 1.5791, 1.5083, 1.6237, 1.5324, 1.6, 1.5703, 1.5518, 1.4557, 1.5142, 1.6592, 1.5913, 1.3739, 1.5366, 1.5575, 1.5205, 1.4785, 1.6111, 1.5861, 1.4798, 1.4806, 1.4898, 1.4108, 1.3855, 1.5362, 1.5827, 1.4556, 1.612, 1.4867, 2, 1.9875, 2.0456, 1.842, 2.1236, 2.1727, 2.0817, 1.7932, 2.101, 2.0199, 2.2014, 1.9665, 2.0714, 2.1465, 2.0089, 2.1927, 1.9458, 2.202, 2.2086, 2.0236, 1.9341, 1.9841, 1.9989, 1.9223, 1.9379, 2.0421, 2.0058, 1.9416, 2.137, 1.882, 1.898, 1.963, 1.9807, 1.834, 1.898, 2.0066, 2.1234, 1.915, 2.1505, 1.977, 2.0534, 2.228, 1.6904, 2.0385, 1.8062, 2.0204, 2.0615, 1.9617, 2.021, 2.165, 1.834, 2.049, 1.909, 2.0091, 2.1625, 2.0017, 2.242, 1.705, 1.8096, 1.9774, 1.9484, 1.9542, 1.9006, 2.1364, 1.9749, 1.9247, 2.1806, 1.8766, 2.0787, 2.0609, 2.1479, 2.072, 2.193, 2.141, 2, 2.0962, 2.2043, 2.0212, 2.0448, 1.7999, 1.7583, 1.9555, 1.9188, 2.1024, 1.9143, 1.9394, 2.0909, 1.9117, 1.9921, 2.0652, 1.9214, 1.9288, 2.1423, 1.9185, 2.061, 2.1143, 2.081, 1.9681, 2.0064, 2.1208, 2.1756, 2.0744, 2.0146, 2.0347, 2.0973, 2.1789, 2.0805, 2.016, 2.0443, 2.0145, 2.0485, 1.9404, 1.9174, 2.0835, 1.8177, 2.069, 2.0467, 2.1035, 1.9164, 2.132, 1.6446, 2.0955, 2.1768, 2.1103, 2.0296, 2.1719, 2.1433, 2.0995, 2.0743, 2.0314, 1.9498, 2.0753, 2.0269, 1.9466, 2.2078, 2.1043, 2.023, 2.1835, 2.1077, 2.0485, 2.0745, 2.107, 2.1725, 2.197, 2.073, 2.2474, 2.1259, 2.147, 2.0411, 2.179, 2.2463, 2.1625, 2.2131, 2.1112, 2.0768, 2.115, 2.1757, 2.1993, 2.215, 2.238, 2.302, 2.201, 2.0725, 2.1555, 2.1632, 2.1028, 2.109, 2.1559, 2, 2.0812, 2.0636, 1.9518, 2.0405, 1.9671, 2.0485, 2.1215, 2.1124, 2.3113, 2.1994, 2.2434, 2.091, 2.116, 2.0227, 1.8781, 1.8379, 2.122, 2.2, 2.167, 2.3945, 2.0258, 1.8949, 2.1335, 2.0854, 1.9336, 2.134, 2.0199, 2.0615, 2.0121, 2.2198, 2.1395, 2.2243, 2.1077, 2.1057, 1.7483, 2.18, 2.2803, 2.2239, 2.391, 2.0472, 2.0262, 2.1252, 2.1381, 2.1761, 2.0419, 2.0137, 2.118, 2.1554, 2.138, 2.1339, 2.0934, 2.043, 2.1349, 2.1876, 2.1754, 2.14, 2.1462, 2.134, 1.9915, 2.1273, 2.0927, 2.2647, 2.2311, 2.1182, 2.015, 1.9269, 2.1666, 2.2637, 2.1654, 2.1667, 2.0613, 2.166, 2.176, 2.1652, 2.1003, 2.06, 2.0848, 2.2029, 2.2304, 2.0783, 2.2003, 2.2126, 2.0313, 2.1385, 2.2015, 2.0062, 1.9707, 2.0471, 2.0847, 2.0912, 2.1259, 2.0132, 2.2313, 1.932, 2.1384, 2.0722, 2.1778, 1.9665, 2.0136, 2.2123, 2.1932, 2.3374, 2.159, 2.2233, 2.1087, 2.179, 2.216, 2.1935, 1.9536, 2.0363, 2.1864, 2.1038, 2.1941, 2.245, 2.16, 2.1348, 2.1572, 2.1454, 2.1742, 2.2023, 2.0934, 2.1621, 2.0701, 2.16, 2.1475, 2.2273, 2.2448, 2.1836, 2.175, 2.3232, 2.2565, 2.284, 2.1939, 2.0941, 2.2709, 2.1029, 1.9158, 2.19, 2.379, 2.173, 2.2437, 2.2624, 2.2093, 2.2531, 2.25, 2.248, 2.2174, 2.1297, 2.0173, 2.1263, 2.208, 2.0977, 2.2609, 2.3667, 2.201, 2.2743, 2.1338, 2.2251, 2.102, 2.2577, 2.1854, 2.3493, 2.159, 2.168, 2.2375, 2.2589, 2.0628, 2.132, 2.1699, 2.3285, 2.192, 2.2133, 2.2179, 2.3002, 2.208, 2.2856, 2.0855, 2.25, 2.1735, 2.1595, 2.0881, 2.1715, 2.2986, 2.2804, 2.0636, 2.2535, 2.0939, 2.208, 2.1907, 2.2991, 2.217, 2.1691, 2.11, 2.0292, 2.011, 2.0168, 2.1483, 2.111, 2.0586, 2.3473, 2.1486, 2.5865, 2.4405, 2.5577, 2.456, 2.41, NA, 2.6811, 2.3779, 2.7572, 2.5393, 2.7288, NA, 2.6591, 2.8193, 2.5132, 2.654, 2.4863, 2.7741, 2.7377, 2.6767, 2.4071, 2.4123, NA, 2.356, 2.52, 2.6631, 2.4839, 2.4638, 2.59, 2.4628, 2.49, 2.643, 2.5499, 2.533, 2.466, 2.5145, 2.5911, 2.515, 2.7347, 2.626, 2.5641, 2.7665, 2.298, 2.563, 2.5209, 2.5408, 2.6245, 2.4907, 2.4728, 2.7138, 2.272, 2.6, 2.537, 2.6142, 2.7606, 2.4786, 2.7479, 2.2775, 2.3125, 2.4968, 2.4808, 2.4446, 2.4036, 2.7649, 2.543, 2.508, 2.594, 2.4322, 2.6081, NA, NA, 2.5826, 2.7014, 2.7417, NA, 2.584, 2.6155, 2.5948, 2.4967, 2.4392, 2.2845, 2.553, 2.4747, 2.5312, 2.4126, 2.4536, 2.6231, 2.5226, 2.5916, 2.585, 2.519, 2.5179, 2.6081, 2.4572, 2.5614, 2.6855, 2.6177, 2.5052, 2.5431, 2.7309, 2.7998, 2.5526, 2.554, 2.5689, 2.7456, 2.6054, 2.6992, 2.4852, 2.6347, 2.582, 2.5281, 2.76, 2.5921, 2.6748, 2.317, 2.6102, 2.6204, 2.6192, 2.3701, 2.5123, 2.2612, 2.6544, 2.7331, 2.7146, 2.5822, 2.6972, 2.7309, 2.7028, 2.6928, 2.5, 2.5893, 2.5169, 2.6186, 2.55, 2.7712, 2.6865, 2.562, 2.8347, 2.6291, 2.6815, 2.8005, 2.6727, 2.8052, 2.8221, 2.6114, 2.8143, 2.6318, 2.6814, 2.5498, 2.6703, 2.7644, 2.6983, 2.7169, 2.6292, 2.6242, 2.6612, 2.7603, 2.814, 2.7986, 2.811, 2.9153, 2.8306, 2.6517, 2.6886, 2.7007, 2.706, 2.6291, 2.7014, 2.6002, 2.6714, 2.6875, 2.5713, 2.5788, 2.5904, 2.6389, 2.7355, 2.6855, 2.8236, 2.486, 2.7896, 2.5527, 2.6436, 2.4531, 2.3533, 2.5682, NA, 2.795, 2.8655, 2.9534, 2.6312, 2.3748, 2.7729, 2.5586, 2.44, 2.6624, 2.5382, 2.5475, 2.4938, 2.6734, 2.6095, 2.6717, 2.6366, 2.5472, 2.3864, 2.7903, 2.841, 2.7699, 2.8326, 2.52, 2.5861, 2.669, 2.8013, 2.7076, 2.5727, 2.5239, 2.5757, NA, 2.7622, 2.6532, 2.6589, 2.6235, 2.7259, 2.6283, 2.6429, 2.66, 2.6561, 2.702, 2.4816, 2.5144, 2.7054, 2.7225, 2.667, 2.6976, 2.6131, 2.5205, 2.7091, NA, NA, 2.6766, 2.5982, NA, 2.7091, 2.8495, 2.6441, 2.7945, 2.6637, 2.7352, 2.7684, 2.7396, 2.7413, 2.8885, 2.6273, 2.7283, 2.8645, 2.5978, 2.6144, 2.6942, 2.7407, 2.7534, 2.8245, 2.4704, 2.9521, 2.642, 2.724, 2.7014, 2.8919, 2.575, 2.6281, 2.677, 2.7407, 2.973, 2.8089, 2.83, 2.6453, 2.7818, 2.7043, 2.8093, 2.4955, 2.6775, 2.9532, 2.1434, 2.6778, 2.8586, 2.8543, 2.6881, 2.7275, 2.7866, 2.8864, 2.7892, 2.6896, 2.8, 2.7107, 2.836, 2.7149, 2.8089, 2.9778, NA, 2.7841, 3.0114, 2.8524, 2.732, 2.8007, 2.7841, 2.8254, 2.7324, 2.5185, 2.98, NA, 2.813, 2.7309, 2.9303, 2.8687, NA, 2.8316, 2.8766, 2.7308, 2.7438, 2.6303, 2.6875, 2.6631, 2.7013, NA, 3.0026, 2.681, 2.9199, 2.7356, NA, 2.5868, 2.8237, 2.7014, 2.9015, 2.9353, 2.7509, 2.8947, 2.9439, NA, NA, 2.852, NA, NA, 2.7421, 2.8189, 2.8385, 2.7892, 3.0434, 2.6936, 2.966, 2.8017, 2.8051, 2.8421, 2.8422, 2.9182, 2.9074, 2.6816, NA, 2.6885, 2.7776, 2.6932, 2.8139, 2.7443, 2.8721, 2.7336, 2.6389, 2.6348, 2.5962, 2.6568, NA, NA, NA, NA)), .Names = c("ID", "TRAT", "SEMANA", "PESO"), row.names = c(NA, -2590L), class = c("nfnGroupedData", "nfGroupedData", "groupedData", "data.frame"), formula = PESO ~ SEMANA | ID, FUN = function (x) max(x, na.rm = TRUE), order.groups = TRUE) #-------------------------------ESTIMATIVA DOS PARAMETROS----------- modelo1.list <- nlsList(PESO ~ a * exp(-b * exp( -c * SEMANA)),start = c(a = 4.10,b = 4.53,c = 0.37 ),na.action=na.omit,data =GALINHA) #--------------------------GRÁFICO------------------------------------ plot(intervals(modelo1.list),layout = c(3,1)) -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

2 2

Panel xyplot Lattice
by Fernando Souza 05 Dez '17

05 Dez '17

Caros colegas Criei esses gráficos utilizando o xyplot {lattice}, utilizei o argumento aspect="xy" os gráficos ficaram como eu queria, porém a relação a área da imagem eles ficaram muito pequenos, como posso reduzir essa diferença? Segue o CMR, anexo inserir a imagem para vcs verem o problema. Como resolver? #install.packages("tidyr") #install.packages("lattice") #install.packages("RColorBrewer") library(tidyr) library(lattice) library(RColorBrewer) DADOS.CURTO <- read.csv(" https://www.dropbox.com/s/b4cckybgrcg86me/galinhas.csv?raw=1",head = TRUE) head(DADOS.CURTO) colnames(DADOS.CURTO) <- c("0","1","2","3","4","5","6","TRAT") DADOS<- DADOS.CURTO %>% gather('0','1','2','3','4','5','6',key ="SEMANA",value = "PESO", -TRAT) head(DADOS) str(DADOS) ID <-rep(seq(1,370), times=7) DADOS<-cbind(ID,DADOS) DADOS <- transform(DADOS, SEMANA = as.numeric(SEMANA),ID=factor(ID)) cores = grDevices::colors()[grep('gr(a|e)y', grDevices::colors(), invert = T)] parset <-list(plot.symbol = list(col = "black",fill = "black"),strip.background=list(col="lightgrey"),superpose.line = list(col = cores[1:375]),superpose.symbol = list(col = cores[1:375])) xyplot(PESO~SEMANA|TRAT,groups=ID,data=DADOS,ylab = "PESO (Kg)",pch = 16,par.settings =parset ,type=c("l","p","g"),layout = c(3,1),aspect = "xy") -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

2 3

Re: [R-br] Sobreposição de BoxPlots
by Leonard Assis 04 Dez '17

04 Dez '17

Jasmire, Acredito que faltou uma coisa básica: voce explicar o que você espera obter como resultado(visualizar) para que alguém possa propor uma solução. Leonard Em 1 de dez de 2017 6:34 PM, "Jasmine Moreira via R-br" < r-br(a)listas.c3sl.ufpr.br> escreveu: Pessoas queridas, estou encontrando um problema de sobreposição ao utilizar um boxplot com o parâmetro barwidth (usei alpha=0.6 para mostrar a sobreposição). Obrigada pela ajuda, Jasmine library(ggplot2) data(mtcars) ggplot(mtcars, aes(x=factor(cyl), y=mpg)) + geom_boxplot(aes(fill=factor(gear)), alpha=0.6, varwidth = TRUE) Resultado: Sem o parâmetro barwidth: _______________________________________________ R-br mailing list R-br(a)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.

5 6

Sobreposição de BoxPlots
by Jasmine Moreira 01 Dez '17

01 Dez '17

Pessoas queridas, estou encontrando um problema de sobreposição ao utilizar um boxplot com o parâmetro barwidth (usei alpha=0.6 para mostrar a sobreposição). Obrigada pela ajuda, Jasmine library(ggplot2) data(mtcars) ggplot(mtcars, aes(x=factor(cyl), y=mpg)) + geom_boxplot(aes(fill=factor(gear)), alpha=0.6, varwidth = TRUE) Resultado: Sem o parâmetro barwidth:

1 0

Modelo de regressão separando por sexo
by Luciane Maria Pilotto 01 Dez '17

01 Dez '17

Olá, estou tentando rodar análises de regressão separando por sexo (data=reg[reg$sexo == 1, ]) e resulta no seguinte erro: GLM.2 <- glm(dent ~ idcat2, family=binomial(logit), data=reg[reg$sexo == 1, ]) Error in model.frame.default(formula = dent ~ idcat2, data = reg[reg$sexo == : comprimentos das variáveis diferem (encontradas em 'idcat2') Alguma dica? _________________________________________ Luciane Maria Pilotto

3 2

editar eixo no ggplo2
by Ze Henrique 30 Nov '17

30 Nov '17

Boa noite pessoal, desde já agradeço a ajuda de vcs. Estou tentando fazer uma figura usando o ggplot2 mas os nomes em volta cortão, alguém poderia me ajudar, em qual função devo alterar para isso não acontecer mais, tentei varias cosias mas não consegui. Desde já obrigado. [image: Imagem inline 1] -- *José Henrique Soler Guilhen* Doutorando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

1 0

Peso espacial em X usando spreg- regressão espacial
by Edimeire Alexandra Pinto 21 Nov '17

21 Nov '17

Boa tarde, pessoal. Alguém por acaso sabe com se trabalha com modelos "parecidos" com CLIFF-ORD usando spreg do pacote sphet? Sei que o Stata faz, mas gostaria de fazer no R. Estou trabalhando com vários tipos de modelos espaciais (na realidade, quase todos do pacote sphet) e esse me gerou muita dúvida, pois preciso de que uma das variaveis explicativas com e sem peso espacial da matriz W. Por exemplo, acrescentando peso espacial a uma das variaveis explicativas do modelo. library(spdep)data(columbus)listw<-nb2listw(col.gal.nb)res1<- spreg(CRIME~HOVAL + INC+ INC, data=columbus , listw= listw, het = TRUE, verbose = FALSE,model="sarar")res2<- spreg(CRIME~HOVAL + INC+ WINC, data=columbus , listw= listw, het = TRUE, verbose = FALSE,model="sarar") ### WINC é o peso espacial multiplicado #pelo variavel INC, especie de W*INC ou mesmo res3<- spreg(CRIME~HOVAL + INC+WINC, data=columbus , listw= listw, het = TRUE, verbose = FALSE,model="error")

2 3

Reshape
by Edson Lira 16 Nov '17

16 Nov '17

Bom dia caros amigos, estou trabalhando com a base de dados abaixo(somente 6 linhas): id UF x2002 x2003 x2004 x2005 x2006 x2007 x2008 x2009 x2011 x2012 x2013 x2014 1 1 Acre 192 189 311 311 333 327 338 354 357 381 337 373 2 2 Alagoas 1251 1269 1292 1357 1413 1410 1355 1401 1344 1391 1429 1417 3 3 Amapá 204 198 239 262 252 294 317 299 302 311 332 355 4 4 Amazonas 958 1050 1433 1528 1502 1488 1531 1673 1689 1720 1771 1908 5 5 Bahia 6885 6901 7167 7336 7292 7338 7665 7792 7602 7416 7530 7879 6 6 Ceará 3705 3871 3944 4120 4143 4138 4314 4380 4121 4184 4215 4319 Estou usando o reshape para tentar transformar as linhas em colunas, ou seja, cada uma UF seria uma coluna, e assim com as demaisComo gostaria que ficasse: Acre Alagoas Amapá Amazonas Bahia Ceará ano 192 1251 3705 2002 189 1269 3871 2003 311 1293 3944 2004 ... .... Estou usando a rotina: pea<-reshape(pea1, varying=c("x2002","x2003","x2004","x2005","x2006","x2007", "x2008","x2009","x2011","x2012","x2013","x2014"), v.names="Medida", timevar="UF", times=c("Acre","Alagoas","Amapá","Amazonas","Bahia", "Ceará","Distrito Federal","Espírito Santo", "Goiás","Maranhão","Mato Grosso","Mato Grosso do Sul", "Minas Gerais","Pará","Paraíba","Paraná"," Pernambuco", "Piauí","Rio de Janeiro","Rio Grande do Norte", "Rio Grande do Sul","Rondônia","Roraima", "Santa Catarina","São Paulo","Sergipe","Tocantins"), new.row.names=1:27, direction="wide") que está me dando o erro: Erro em varying[, i] : número incorreto de dimensões Alguém tem alguma sugesstão? [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas

5 5

usar separador de milhar
by Adriele Giaretta Biase 16 Nov '17

16 Nov '17

Bom dia a todos, tenho uma tabela, em que os números não estão separados por milhar, uma tabela com muitas linhas, mas para exemplificar, tomei como exemplo essa tabela pequena abaixo: Ex: teste metodo_1 metodo_2 A 30000000 530000000 B 40000034 550003600 C 30000020 570000002 D 50000000 830000000 E 20000000 730000000 F 30000400 530000500 Eu faço graficos de barras usando cada linha de forma que os valores apareçam em cima de cada barra. A dúvida é como fazer com que esses valores (da tabela) fiquem separados por milhar. A minha intenção é fazer tanto em padrão ingles como em padrão português. Exemplo, tomando o padrão português, o teste A da tabela acima, no metodo1 queria que o R lesse 3.000.0000 e no metodo2 53.000.0000 Se tomasse o padrão ingles, o teste A no metodo1 seria 3,000,0000 e no metodo2 53,000,0000 Porém, a tabela não tem nem ponto e nem virgula para poder trocar de uma nomenclatura para outra.

2 1

Corrigir ausência de dados
by Edson Lira 15 Nov '17

15 Nov '17

Caros amigos, tenho uma série de dados de pea e pnea que têm uma ausência de dados no ano de 2010. Como faço para corrigir essa ausência? TrêsDuas dúvidas: 1. Como organizar a planilha visto que devo corrigir a ausência por estado para corrigir a ausência de pea e pnea para o ano de 2010 por estado?2. Qual a técnica que me aconselham a usar?3. E como fazê-lo com uma única rotina? Abaixo o link da planilha Download pea pnea 02 14.csv | | | Download pea pnea 02 14.csv | | | [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas

4 9

Gráfico com médias e erro agrupando var. preditora através de duas categorias
by Carlos Peixoto Neto 13 Nov '17

13 Nov '17

Boa tarde a todos. Peço um auxílio para o caso : uma espécie de peixe dividida em duas categorias de tamanho "grande" e "pequeno" e dentro dessas categorias obtenho variáveis de interação social métrica em segundos "Peixe isolado ", "peixe interação pequeno ", "peixe interação grande " estas variáveis tenho tanto para o grupo de peixes "grande" como para o "pequeno". Desejo ver o efeito na variável "Altura coluna d'água (metros)". Como representar estes dados graficamente no R, de forma que no gráfico o eixo x tenha as var. de interação social e para cada uma apareça a representação de média e erro padrão para categoria de tamanho "grande" e "pequeno" ? Segue link com médias simuladas <https://drive.google.com/open?id=1c7EkAl0s7abiHdUdbTWrmKTowsz9z0GV>. Grato.

2 1

Correlação entre variáveis categoricas
by Elias Carvalho 10 Nov '17

10 Nov '17

Boa Tarde Pessoal Tenho um banco de dados com 9 variáveis categóricas binárias e 12 variáveis categóricas ordinais. Preciso saber se há correlação entre elas, preciso de um valor que diga se há uma correlação ou associação baixa, média ou alta. O objetivo é pré-processar um banco de dados antes de criar um grafo de rede (path analysis e rede bayesiana) e eliminar os relacionamentos de variáveis com baixa corrrelação para evitar conexões causais espurias. Lendo alguns materiais foi recomendado o seguinte: - Para variáveis categoricas x categoricas usar o qui-quadrado para determinar se elas são independentes ou não.Sendo o valor P considerado a medida de correlação - Calcular também o v de Crammer (quanto menor maior é a força da correlação) Alguuém pode me confirmar se esse é o melhor método ? Ou sugere outra coisa ? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 2

Função do ExpDes.pt
by Andre Oliveira 10 Nov '17

10 Nov '17

Boa noite,não houve atualização por aqui, mas, a função faixas do pacote ExpDes.pt não funciona mais. Alguém saberia apontar uma possível causa? faixas(trat, mes, bloco, altura, quali = c(TRUE, TRUE), mcomp = "sk", fac.names = c("trat", "mes"), sigT = 0.05, sigF = 0.05) Error in faixas(trat, mes, bloco, altura, quali = c(TRUE, TRUE), mcomp = "sk", : não foi possível encontrar a função "faixas" att., André Oliveira Souza

2 1

Qual teste devo ou se dá para usar
by Felipe Felix Costa 08 Nov '17

08 Nov '17

Olá, estou com dúvidas se existe algum teste estatístico que permita eu analisar, ao mesmo tempo, o índice de área foliar em diferentes estratos de altura da vegetação e regiões diferentes. Sendo que há dependência entre os estratos nas diferentes alturas. att., abs -- *Felipe Felix Costa* Tel. 96-98122-3077

2 1

Ajuste dos graus de liberdade - método Cochran
by Euriann Yamamoto 08 Nov '17

08 Nov '17

Boa tarde, r-nautas. O R possui uma rotina para o cálculo dos ajustes dos graus de liberdade quando as variâncias são muito heterogêneas (>7)? No entanto, já fiz o cálculo à mão. Gostaria de saber se tem algum problema em continuar a análise de estabilidade e adaptabilidade com as médias originais (sem transformação), considerando essas variâncias. Att.,

2 1

Re: [R-br] Digest R-br, volume 83, assunto 2
by Wilson Cohen 08 Nov '17

08 Nov '17

Olá Pessoal Alguém poderia compartilhar um exemplo simples de uso do modelo não-linear usando o pacote nlme? Obrigado. Em 6 de novembro de 2017 11:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)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(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)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. Pacotes não carregam no R version 3.4.2 (Vanessa Felix) > 2. Re: Pacotes não carregam no R version 3.4.2 > (Fernando Antonio de souza) > 3. Re: Pacotes não carregam no R version 3.4.2 (Vanessa Felix) > 4. RES: Pacotes não carregam no R version 3.4.2 > (Jobenil Luiz Magalhães Júnior) > 5. Re: Pacotes não carregam no R version 3.4.2 > (Fernando Antonio de souza) > 6. Re: Pacotes não carregam no R version 3.4.2 (Vanessa Felix) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 6 Nov 2017 10:09:59 -0200 > From: Vanessa Felix <vanessa.felix.est(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Pacotes não carregam no R version 3.4.2 > Message-ID: > <CAJAeJpsU0JSoJZLDsFZTQReApafV-tYD+iy6ju0EASqjs+cLxQ@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > Pessoal, alguém pode me ajudar a resolver este problema?? > Isto está acontecendo com vários outros pacotes. > PS: Acabei de instalar o R na minha máquina. > > > install.packages(mapdata) > Error in install.packages(mapdata) : object 'mapdata' not found > > install.packages("mapdata") > --- Please select a CRAN mirror for use in this session --- > Warning: unable to access index for repository > https://cran-r.c3sl.ufpr.br/src/contrib: > cannot open URL 'https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES' > Warning: unable to access index for repository > https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4: > cannot open URL ' > https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4/PACKAGES' > Warning message: > package ‘mapdata’ is not available (for R version 3.4.2) > > library(mapdata) > Error: package ‘maps’ required by ‘mapdata’ could not be found > > Desde já, obrigada!! > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171106/0b524ef6/attachment-0001.html> > > ------------------------------ > > Message: 2 > Date: Mon, 06 Nov 2017 12:18:38 +0000 > From: Fernando Antonio de souza <nandodesouza(a)gmail.com> > To: Vanessa Felix <vanessa.felix.est(a)gmail.com>, a lista Brasileira > oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Pacotes não carregam no R version 3.4.2 > Message-ID: <local-ae5f2495-3148(a)nylas-mail.nylas.com> > Content-Type: text/plain; charset="utf-8" > > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171106/88a0500b/attachment-0001.html> > > ------------------------------ > > Message: 3 > Date: Mon, 6 Nov 2017 10:20:28 -0200 > From: Vanessa Felix <vanessa.felix.est(a)gmail.com> > To: Fernando Antonio de souza <nandodesouza(a)gmail.com> > Cc: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Pacotes não carregam no R version 3.4.2 > Message-ID: > <CAJAeJpssUAv_sBQg_UVYaP8ZsppdrdHiO6uOqMZgVSQT4qg > 8rw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Tentei esse Fernando, e não funcionou. > > Em 6 de novembro de 2017 10:18, Fernando Antonio de souza < > nandodesouza(a)gmail.com> escreveu: > > > pode ser algum problema com o repositório! Você ja tentou com outro > > espelho? > > > > Eu utilizo esse aqui e deu certo https://vps.fmvz.usp.br/CRAN/ > > > > options(repos = "https://vps.fmvz.usp.br/CRAN/") > > > > att > > > > ========================================= > > Fernando Souza > > Celular: (31)99796-8781 <(31)%2099796-8781> (Vivo) > > E-mail:nandodesouza@gmail.com > > <https://n1.nylas.com/link/6593b52914d1aab4e1f96df3e30975 > 198aa7fa918e747a036a54cfff041669a0/0?redirect=mailto%3Ae- > mail%253Anandodesouza%40gmail.com&recipient=vanessa.felix.est%40gmail.com> > > ========================================== > > > > On Nov 6 2017, at 10:09 am, Vanessa Felix via R-br < > > r-br(a)listas.c3sl.ufpr.br> wrote: > > > >> Pessoal, alguém pode me ajudar a resolver este problema?? > >> Isto está acontecendo com vários outros pacotes. > >> PS: Acabei de instalar o R na minha máquina. > >> > >> > install.packages(mapdata) > >> Error in install.packages(mapdata) : object 'mapdata' not found > >> > install.packages("mapdata") > >> --- Please select a CRAN mirror for use in this session --- > >> Warning: unable to access index for repository > >> https://cran-r.c3sl.ufpr.br/src/contrib > >> <https://cran-r.c3sl.ufpr.br/src/contrib?recipient=vanessa. > felix.est%40gmail.com> > >> : > >> cannot open URL 'https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES > >> <https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES? > recipient=vanessa.felix.est%40gmail.com> > >> ' > >> Warning: unable to access index for repository > >> https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4 > >> <https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4? > recipient=vanessa.felix.est%40gmail.com> > >> : > >> cannot open URL 'https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4/ > >> PACKAGES > >> <https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4/ > PACKAGES?recipient=vanessa.felix.est%40gmail.com> > >> ' > >> Warning message: > >> package ‘mapdata’ is not available (for R version 3.4.2) > >> > library(mapdata) > >> Error: package ‘maps’ required by ‘mapdata’ could not be found > >> > >> Desde já, obrigada!! > >> > >> _______________________________________________ > >> R-br mailing listR-br@listas.c3sl.ufpr.brhttps://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. > >> > >> > > > -- > > > ------------------------------------------------------------ > ------------------------------ > Vanessa Felix do Nascimento Sergio > Graduada em Estatística - ICE/UFJF > Mestre em Saúde Coletiva - IMS/UERJ > http://lattes.cnpq.br/2516291371035474 > (21) 97557-3956 > ------------------------------------------------------------ > ------------------------------- > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171106/eab9a5ac/attachment-0001.html> > > ------------------------------ > > Message: 4 > Date: Mon, 6 Nov 2017 10:53:34 -0200 > From: Jobenil Luiz Magalhães Júnior <pjobenil(a)gmail.com> > To: Vanessa Felix <vanessa.felix.est(a)gmail.com>, a lista Brasileira > oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] RES: Pacotes não carregam no R version 3.4.2 > Message-ID: <5a005b4d.6232c80a.fa547.21bf(a)mx.google.com> > Content-Type: text/plain; charset="utf-8" > > Vanessa > > Deve ser problema com o repositório. Uma solução que você pode tentar é > baixar o pacote de qualquer um dos repositórios e instalá-lo diretamente > pelo menu. > > Atenciosamente > Jobenil Júnior > > Enviado do Email para Windows 10 > > De: Vanessa Felix via R-br > Enviado:segunda-feira, 6 de novembro de 2017 10:10 > Para: r-br(a)listas.c3sl.ufpr.br > Assunto: [R-br] Pacotes não carregam no R version 3.4.2 > > Pessoal, alguém pode me ajudar a resolver este problema?? > Isto está acontecendo com vários outros pacotes. > PS: Acabei de instalar o R na minha máquina. > > > install.packages(mapdata) > Error in install.packages(mapdata) : object 'mapdata' not found > > install.packages("mapdata") > --- Please select a CRAN mirror for use in this session --- > Warning: unable to access index for repository > https://cran-r.c3sl.ufpr.br/src/contrib: > cannot open URL 'https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES' > Warning: unable to access index for repository > https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4: > cannot open URL 'https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4/ > PACKAGES' > Warning message: > package ‘mapdata’ is not available (for R version 3.4.2) > > library(mapdata) > Error: package ‘maps’ required by ‘mapdata’ could not be found > > Desde já, obrigada!! > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171106/7f48670d/attachment-0001.html> > > ------------------------------ > > Message: 5 > Date: Mon, 06 Nov 2017 13:02:42 +0000 > From: Fernando Antonio de souza <nandodesouza(a)gmail.com> > To: Vanessa Felix <vanessa.felix.est(a)gmail.com> > Cc: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Pacotes não carregam no R version 3.4.2 > Message-ID: <local-44011900-4fa9(a)nylas-mail.nylas.com> > Content-Type: text/plain; charset="utf-8" > > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171106/f3a36933/attachment-0001.html> > > ------------------------------ > > Message: 6 > Date: Mon, 6 Nov 2017 11:04:27 -0200 > From: Vanessa Felix <vanessa.felix.est(a)gmail.com> > To: Fernando Antonio de souza <nandodesouza(a)gmail.com> > Cc: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Pacotes não carregam no R version 3.4.2 > Message-ID: > <CAJAeJpsK0PSQ5F7S=dbF=HczqdFM51bYH7nQfcHL3pzjL9tPaQ@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > Eu não uso Rstudio. > Vou tentar a dica do Jobenil tbm. > > Em 6 de novembro de 2017 11:02, Fernando Antonio de souza < > nandodesouza(a)gmail.com> escreveu: > > > Você utiliza o Rstudio? Se sim, tenta fazer a instalação manual > > > > Entre nesse diretório e baixa o arquivo zip do pacote maps e mapdata > > https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4/ > > > > Na janela do Rstudo no canto inferior direito (Janela History), clique na > > aba packages >> install >> mude install from para: "Package Archive file" > > >> navegue até onde salvou os arquivos zip maps e mapdata e selecione>> > > clique em install. > > veja se assim instala > > > > ========================================= > > Fernando Souza > > Celular: (31)99796-8781 <(31)%2099796-8781> (Vivo) > > E-mail:nandodesouza@gmail.com > > <https://n1.nylas.com/link/e4be6b11068dee176cc7e8325461f6 > 612bca66bb056034be806c868a54955594/0?redirect=mailto%3Ae- > mail%253Anandodesouza%40gmail.com&recipient=vanessa.felix.est%40gmail.com> > > ========================================== > > > > On Nov 6 2017, at 10:20 am, Vanessa Felix <vanessa.felix.est(a)gmail.com> > > wrote: > > > >> Tentei esse Fernando, e não funcionou. > >> > >> Em 6 de novembro de 2017 10:18, Fernando Antonio de souza < > >> nandodesouza(a)gmail.com> escreveu: > >> > >> pode ser algum problema com o repositório! Você ja tentou com outro > >> espelho? > >> > >> Eu utilizo esse aqui e deu certo https://vps.fmvz.usp.br/CRAN/ > >> <https://vps.fmvz.usp.br/CRAN/?recipient=vanessa.felix.est%40gmail.com> > >> > >> options(repos = "https://vps.fmvz.usp.br/CRAN/ > >> <https://vps.fmvz.usp.br/CRAN/?recipient=vanessa.felix.est%40gmail.com > >") > >> > >> att > >> > >> ========================================= > >> Fernando Souza > >> Celular: (31)99796-8781 (Vivo) > >> E-mail:nandodesouza@gmail.com > >> <https://n1.nylas.com/link/6593b52914d1aab4e1f96df3e30975 > 198aa7fa918e747a036a54cfff041669a0/0?redirect=mailto%3Ae- > mail%253Anandodesouza%40gmail.com&amp%3Brecipient=vanessa. > felix.est%40gmail.com&recipient=vanessa.felix.est%40gmail.com> > >> ========================================== > >> > >> On Nov 6 2017, at 10:09 am, Vanessa Felix via R-br < > >> r-br(a)listas.c3sl.ufpr.br> wrote: > >> > >> Pessoal, alguém pode me ajudar a resolver este problema?? > >> Isto está acontecendo com vários outros pacotes. > >> PS: Acabei de instalar o R na minha máquina. > >> > >> > install.packages(mapdata) > >> Error in install.packages(mapdata) : object 'mapdata' not found > >> > install.packages("mapdata") > >> --- Please select a CRAN mirror for use in this session --- > >> Warning: unable to access index for repository > >> https://cran-r.c3sl.ufpr.br/src/contrib > >> <https://cran-r.c3sl.ufpr.br/src/contrib?recipient=vanessa. > felix.est%40gmail.com> > >> : > >> cannot open URL 'https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES > >> <https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES? > recipient=vanessa.felix.est%40gmail.com> > >> ' > >> Warning: unable to access index for repository > >> https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4 > >> <https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4? > recipient=vanessa.felix.est%40gmail.com> > >> : > >> cannot open URL 'https://cran-r.c3sl.ufpr.br/b > >> in/windows/contrib/3.4/PACKAGES > >> <https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4/ > PACKAGES?recipient=vanessa.felix.est%40gmail.com> > >> ' > >> Warning message: > >> package ‘mapdata’ is not available (for R version 3.4.2) > >> > library(mapdata) > >> Error: package ‘maps’ required by ‘mapdata’ could not be found > >> > >> Desde já, obrigada!! > >> > >> _______________________________________________ > >> R-br mailing listR-br@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/ > cgi-bin/mailman/listinfo/r-br <https://listas.inf.ufpr.br/ > cgi-bin/mailman/listinfo/r-br?recipient=vanessa.felix.est%40gmail.com> > >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia < > http://www.leg.ufpr.br/r-br-guia?recipient=vanessa.felix.est%40gmail.com>) > e forne�a c�digo m�nimo reproduz�vel. > >> > >> > >> > >> > >> -- > >> > >> > >> ------------------------------------------------------------ > >> ------------------------------ > >> Vanessa Felix do Nascimento Sergio > >> Graduada em Estatística - ICE/UFJF > >> Mestre em Saúde Coletiva - IMS/UERJ > >> http://lattes.cnpq.br/2516291371035474 > >> <http://lattes.cnpq.br/2516291371035474?recipient= > vanessa.felix.est%40gmail.com> > >> (21) 97557-3956 > >> ------------------------------------------------------------ > >> ------------------------------- > >> > >> > > > -- > > > ------------------------------------------------------------ > ------------------------------ > Vanessa Felix do Nascimento Sergio > Graduada em Estatística - ICE/UFJF > Mestre em Saúde Coletiva - IMS/UERJ > http://lattes.cnpq.br/2516291371035474 > (21) 97557-3956 > ------------------------------------------------------------ > ------------------------------- > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171106/61cfccac/attachment-0001.html> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 83, assunto 2 > **************************************** > -- *Prof. MSc. Wilson Martins da Silva* *Mestre em Ecologia Aquática e Pesca - PPGEAP/ICB/UFPA* *Especialista em Biodiversidade Amazônica - UFPA/ATMDepto de Morfologia e Ciências Fisiológicas* *Centro de Ciências Biológicas e da Saúde* *Universidade do Estado do Pará* *Campus IX - Altamira - Pará* http://lattes.cnpq.br/2440617911389759 wilsoncelula(a)gmail.com, wilsoncelula(a)hotmail.com <wilsonms(a)ufpa.br>

3 2

Comparação de modelos não lineares via teste de Vuong [nonnest2]
by Diego Vieira 07 Nov '17

07 Nov '17

Companheiros recorro a ajuda de vocês novamente. Meu objetivo é comparar dois modelos, ambos exponenciais negativo, porém para dois diferentes tratamentos. Como os modelos não são aninhados eu estou utilizando o teste de Vuong do pacote "nonnest2" para avaliar se as curvas obtidas a partir desses modelos diferem entre si. O problema é que eu consigo entender a saída do R. O "Variance test" diz que os modelos são indistinguíveis, porém o teste de razão de verossimilhança não aninhada diz que o modelo 1 e diferente do 2, isso é correto ? Segue a saída e os comandos. Agradeço a todos que me ajudarem. Abraços. SAÍDA: Model 1 Class: nls Call: nls(formula = diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == ... Model 2 Class: nls Call: nls(formula = diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == ... Variance test H0: Model 1 and Model 2 are indistinguishable H1: Model 1 and Model 2 are distinguishable w2 = 0.048, p = 0.997 Non-nested likelihood ratio test H0: Model fits are equal for the focal population H1A: Model 1 fits better than Model 2 z = 6.840, p = 3.97e-12 H1B: Model 2 fits better than Model 1 z = 6.840, p = 1 COMANDOS: dados.ajuste <- expand.grid(parA = c(0.0,0.99,2,3,4,5,6)+1, sistema = factor(c("MAFA01","MAFA04"))); perfilMAFA01; perfilMAFA04 dados.ajuste$diversidade <- c(55.00000,19.572140, 10.230871, 7.224990, 6.004212, 5.402218, 5.056630, 62.00000, 31.69523, 21.41784, 17.61259, 15.76542, 14.70719, 14.03339); dados.ajuste expMAFA01 <- nls(diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == "MAFA01",], start = list(a = 50, b = -1)) summary(expMAFA01) lines(fitted.values(expMAFA01)) expMAFA04 <- nls(diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == "MAFA04",], start = list(a = 50, b = -1)) summary(expMAFA04) lines(fitted.values(expMAFA04)) ggplot(dados.ajuste,aes(parA,diversidade))+ stat_function(fun=function(parA){coef(expMAFA01)[[1]] * parA^coef(expMAFA01)[[2]]}, aes(colour="MAFA01", linetype="MAFA01")) + stat_function(fun=function(parA){coef(expMAFA04)[[1]] * parA^coef(expMAFA04)[[2]]}, aes(colour="MAFA04", linetype="MAFA04")) AIC(expMAFA01,expMAFA04) vuongtest(expMAFA01,expMAFA04) -- Diego dos Santos Vieira Engenheiro Florestal - UFRA Mestre em Ciência Florestal - UFVJM Doutorando em Ciência Florestal - UFVJM

2 1

Pacotes não carregam no R version 3.4.2
by Vanessa Felix 06 Nov '17

06 Nov '17

Pessoal, alguém pode me ajudar a resolver este problema?? Isto está acontecendo com vários outros pacotes. PS: Acabei de instalar o R na minha máquina. > install.packages(mapdata) Error in install.packages(mapdata) : object 'mapdata' not found > install.packages("mapdata") --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository https://cran-r.c3sl.ufpr.br/src/contrib: cannot open URL 'https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES' Warning: unable to access index for repository https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4: cannot open URL ' https://cran-r.c3sl.ufpr.br/bin/windows/contrib/3.4/PACKAGES' Warning message: package ‘mapdata’ is not available (for R version 3.4.2) > library(mapdata) Error: package ‘maps’ required by ‘mapdata’ could not be found Desde já, obrigada!!

3 6

Excluir valores da matriz de correlação
by Elias Carvalho 31 Out '17

31 Out '17

Boa noite pessoal Como posso excluir de uma matriz de correlação valores muito baixo como *esse em negrito* como abaixo G001 G002 G008 C051 C071 C101 C151 R001 R02A R02C R2B1 G001 1.00000000 0.14842780 -0.114680730 0.302947649 0.32931462 0.37740045 0.30480517 -0.34884430 -0.39049631 -0.34041558 0.056306588 G002 0.14842780 1.00000000 0.070276175 0.170691108 0.09668893 *0.02550203* 0.13058268 -0.15025183 -0.15148344 -0.18614711 0.118041910 G008 -0.11468073 0.07027618 1.000000000 *0.007789258* 0.06977162 0.01623666 -0.14927575 0.01368586 *0.03586060 0.00202448* 0.486717281 C051 0.30294765 0.17069111 0.007789258 1.000000000 0.37174418 0.26435392 0.22513174 -0.47614134 -0.48203009 -0.42861271 0.372495377 C071 0.32931462 0.09668893 0.069771618 0.371744175 1.00000000 0.45945064 0.40883404 -0.24972945 -0.22757192 -0.17468398 0.312496407 -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1

Re: [R-br] Digest R-br, volume 82, assunto 19
by Magnun Rochel 29 Out '17

29 Out '17

Oi Glenio, Vc ja tentou com a função bwplot() do lattice? Ficaria algo assim: bwplot(weight ~ Litter | Treatment + sex, data= RatPupWeight) []s, Magnun 2017-10-28 16:12 GMT-02:00 <r-br-request(a)listas.c3sl.ufpr.br>: > Enviar submissões para a lista de discussão R-br para > r-br(a)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(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)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. Ajuda para elaboração de boxplot em painel > (Gilenio Borges Fernandes) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 28 Oct 2017 15:11:59 -0300 > From: Gilenio Borges Fernandes <gilenio(a)ufba.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Ajuda para elaboração de boxplot em painel > Message-ID: > <CA+7i11aPG1Oi69uA8-4uz0aL9ZjbYQpOMkWou_uuUp_+ > GR3Ung(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Prezad@s > > Estou precisando de auxilio para elaborar boxplots em painel pelo menos > aproximadamente iguais aos dos arquivos pdf anexos. > > O banco de dados é do pacote nlme do R; pode ser acessado com o CMR > > library(nlme) > > data(RatPupWeight) > > head(RatPupWeight) > > names(RatPupWeight) > > Grato > > > Gilênio Borges Fernandes > Professor Associado IV (Aposentado) > Professor Adjunto A (Substituto, 03/2015 a 03/2017) > Universidade Federal da Bahia > Instituto de Matemática e Estatística > Departamento de Estatística > Av. Adhemar de Barros, s/n – Ondina. > 40.170-110 - Salvador - BA, Brasil > Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 > Skype: gilenio.fernandes > Lattes: http://lattes.cnpq.br/6764860618464860 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171028/e953defd/attachment.html> > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: rat-pup-grafico1.pdf > Tipo: application/pdf > Tamanho: 80064 bytes > Descrição: não disponível > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171028/e953defd/attachment.pdf> > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: rat-pup-Grafico2.pdf > Tipo: application/pdf > Tamanho: 87655 bytes > Descrição: não disponível > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20171028/e953defd/attachment-0001.pdf> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 82, assunto 19 > ***************************************** > -- Magnun R. Madruga Linguista +55 11 953956550

1 0

Ajuda para elaboração de boxplot em painel
by Gilenio Borges Fernandes 28 Out '17

28 Out '17

Prezad@s Estou precisando de auxilio para elaborar boxplots em painel pelo menos aproximadamente iguais aos dos arquivos pdf anexos. O banco de dados é do pacote nlme do R; pode ser acessado com o CMR library(nlme) data(RatPupWeight) head(RatPupWeight) names(RatPupWeight) Grato Gilênio Borges Fernandes Professor Associado IV (Aposentado) Professor Adjunto A (Substituto, 03/2015 a 03/2017) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 Skype: gilenio.fernandes Lattes: http://lattes.cnpq.br/6764860618464860

1 0

Identificação de variáveis
by Elias Carvalho 25 Out '17

25 Out '17

Como posso identificar se uma variável é 'ordered factor nas seguintes condições: - o read.csv leu ela como int: - o read.csv leu ela como factor e não ordered factor: Obrigado

2 1

ajuda
by Edmar Caldas 24 Out '17

24 Out '17

Turma tem hora que funciona e hora não. alguém pode me ajudar? a variável datanasc no Rstudio está como <dttm> o que significa isso? Edmar

4 3

ANOVA in R using the function lapply ans terms
by Euriann Yamamoto 24 Out '17

24 Out '17

Hi everyone. In my dataset I´ve got 36 genotypes in 5 locations in partially balanced incomplete design 6x6, simple. We are analyzing productivity. I need to use adjusted means to do the orders analysis. TO obtain ANOVA individual to the 5 locations I run command: mi <- lapply(split(da, f=da$local), aov, formula=producao ~ repeticao/bloco + cultivar) lapply(mi, summary)#not adjusted for blocks But to transform cultivar adjusted for blocks (to obtain adjusted means) I need use command terms(formula,keep.order=TRUE) I DO NOT KNOW HOW TO USE THIS IN THE CORRECT ORDER. Something like this? m0 <- lapply(split(da, f=da$local), aov(terms(producao ~ repeticao/bloco + cultivar,keep.order = TRUE),da)) lapply(m0, summary)#adjusted to blocks Thanks Euriann

2 1

ENC: Látice parcialmente balanceado em vários locais
by Euriann Yamamoto 24 Out '17

24 Out '17

Walmes, eu consegui fazer pelo comando lapply a anova individual dos 5 locais. Só que agora não estou conseguindo colocasr os termos da fonte de variação em ordem: # ajuste do modelo de efeitos fixos para bloco dentro de repeticao #Anova Individual mi <- lapply(split(da, f=da$local), aov, formula=producao ~ repeticao/bloco + cultivar) lapply(mi, summary)#nao ajustado pra bloco #Arrumar ordem dos termos e tornar cultivar ajustado pra bloco mi <- lapply(split(da, f=da$local), aov, terms(formula=producao ~ repeticao/bloco + cultivar, keep.order = TRUE)) Ocorre um erro: Error in terms.default(formula, "Error", data = data) : no terms component nor attribute Creio que seja na fórmula. Qual seria a ordenação correta para poder utilizar esse comando terms? Att., Euriann L. M. Yamamoto Doutoranda em Agronomia Universidade Federal da Grande Dourados ________________________________ De: Euriann Yamamoto Enviado: quinta-feira, 19 de outubro de 2017 17:05 Para: Walmes Zeviani Assunto: RE: [R-br] Látice parcialmente balanceado em vários locais Walmes, segue abaixo o script que alterei para o látice parc. balanc.. Consta nele a ANOVA com tratamentos não-ajustados e bloco/repetição não-ajustados. Acontece que eu tenho 5 experimentos em látice 6x6, cada um em 5 municípios na safra agrícola 2008-09. Como automatizar a análise para que eu utilize um mesmo .txt contendo os locais, tratamentos, blocos, repetições e a variável? Script: rm(list=ls()) ls() library(lattice) library(agricolae) library(contrast) library(multcomp) library(doBy) library(latticeExtra) library(nlme) #----------------------------------------------------------------------- # lendo arquivos de dados da <- read.table("C:/ADSR/Tese/teste.txt", header = TRUE, colClasses = rep(c("factor", "numeric"), c(3, 1))) str(da) #----------------------------------------------------------------------- # ver os dados xtabs(~repeticao + bloco, da) xtabs(~repeticao + cultivar, da) xyplot(producao ~ cultivar, groups = bloco, data = da) xyplot(producao ~ cultivar | repeticao, groups = bloco, data = da) #----------------------------------------------------------------------- # ajuste do modelo de efeitos fixos para bloco dentro de repeticao m0 <- lm(producao ~ repeticao/bloco + cultivar, data = da) par(mfrow = c(2, 2)) plot(m0) layout(1) # #Teste de normalidade # shapiro.test(residuals(m0)) # #Teste de Bartlett # bartlett.test(residuals(m0)~cultivar, data=da) anova(m0) # cultivar nao ajustado para blocos #----------------------------------------------------------------------- #Só arrumar ordem dos termos m0 <- lm(terms(producao ~ repeticao/bloco + cultivar, keep.order = TRUE), data = da) anova(m0) Euriann L. M. Yamamoto Doutoranda em Agronomia Universidade Federal da Grande Dourados ________________________________ De: Walmes Zeviani <walmeszeviani(a)gmail.com> Enviado: terça-feira, 17 de outubro de 2017 19:53 Para: Euriann Yamamoto Assunto: Re: [R-br] Látice parcialmente balanceado em vários locais Não sei te responder com segurança. Testa e vê no que dá, rs. É bem provável que precise de adaptações, ainda que mínimas. Walmes.

3 2

Erro em script de apostila Embrapa Arroz e Feijão
by Euriann Yamamoto 24 Out '17

24 Out '17

Bom dia, Walmes. Estou rodando a análise 20 Análise de experimentos em vários locais (p.114) da apostila Embrapa Arroz e Feijão - Estatística experimental Aplicada 2011 e ocorreu o seguinte erro em anexo. Quando executo o comando mmqs, surge o seguinte erro: Error in qt((1 + conf.int)/2, dfdm) : Argumento não-numérico para função matemática Como resolver este problema visto que apenas copiei e colei todo o script. Att., Euriann

3 2

[OFF TOPIC] Dúvida diferente número de plantas por UA, mas, mesma área da UA
by Marcelo Laia 23 Out '17

23 Out '17

Caros colegas, Nesse final de semana me surgiu uma dúvida. Quando da análise de um experimento cuja unidade amostral possui o mesmo número de plantas, essas entram na análise normalmente. Por outro lado, tenho uma situação que, devido aos diferentes espaçamentos, a quantidade de plantas por UA é diferente. Para se fazer um teste de média eu posso utilizar as UA com números de plantas diferentes ou devo usar a média das parcelas? O experimento visa identificar qual é o melhor espaçamento de plantio para cada genótipo. Assim, a área amostral é sempre a mesma, mas, o N difere. O experimento conta com 6 blocos e os tratamentos foram sorteados dentro de cada bloco. Não sei me fiz entender, mas, lá vai minha dúvida que nada tem a ver com R. Se puder ajudar, agradeço muito. -- Marcelo

4 4

Perda de parcela - Lista R
by ana paula coelho madeira 20 Out '17

20 Out '17

Bom dia Walmes! A um tempo atrás você postou na lista de discussão do R uma sugestão acerca de experimentos em blocos com parcela perdida. Sugeriu, na ocasião, que a pessoa desse uma olhada no script http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap08dbc-iso.R a partir do chunck number 6. Estou com o mesmo problema... tenho um DBC com perda de duas parcelas. No entanto, tentei seguir o script citado acima e aparece a seguinte msg de erro: > sim <- expand.grid(blc=gl(4,1, labels="bl"), trt=gl(5,1, labels="tr")) Error in as.character.factor(1:4) : malformed fator Você poderia me ajudar? Desde já Agradeço. PS: Peço desculpas por direcionar o email a você, porém como foi uma sugestão (script) dada por você, achei mais viável.

2 1

Preencher área entre duas linhas
by Delcio Rudinei Bortolanza 19 Out '17

19 Out '17

Olá! Segue um exemplo hipotético: #=========================================== x<-1:5 set.seed(4) y0<-sample(x) set.seed(2) y1<-sample(x)+5 matplot(x,data.frame(y0,y1),pch = c(0,1),col = c(1,1),las=1) xspline(x,y0,-0.5) xspline(x,y1,-0.5) #=============================================== Como preencher com cor entre as linhas suavizadas? Obs: a função polygon() não gera o resultado pretendido, pois une os pontos com linhas retas. Att., Delcio Bortolanza.

2 3

Efeitos diretos e indiretos no lavaan
by Elias Carvalho 19 Out '17

19 Out '17

Alguem sabe me dizer se é possível extrair os efeitos diretos e indiretos de uma variável específica de um grafo de Path Analysis gerado pelo Lavaan ? Explicando melhor, tennho um estudo avaliando o relacionamento de várias variáveis com um ou mais desfechos, então gostaria de criar uma tabela com esses valores e depois um grafico que permita analisar melhor os efeitos diretos e indiretos sobre a variável de desfecho. -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

1 0

Tratamento não ajustado x tratamentos ajustados no látice e DBC
by Euriann Yamamoto 19 Out '17

19 Out '17

Boa tarde. Alguém possui algum material que fale um pouco a respeito de Tratamento não ajustado x tratamentos ajustados e quais as maneiras de utilização? Estou confusão sobre qual tipo utilizar na minha análise.

1 0

Pacotes
by Edmar Caldas 18 Out '17

18 Out '17

Boa Tarde! é possível saber o nome de todos os pacotes da versão 3.4.1?

3 3

Alocar elementos em lista...
by andrebvs＠bol.com.br 18 Out '17

18 Out '17

2 1

Merge tabelas não padronizados.
by Ronaldo Câmara de Araújo 17 Out '17

17 Out '17

Estou com um bando de dados que possui um texto ao qual, em algum local do campo tem o texto que necessito para fazer um merge com outra tabela que tem essa informação. Exemplo: SERVICO DE LOCACAO DE 02 GALPOES ONDE ABRIGARA AS INSTALACOES DO CDD ATALAIA, CONFORME DESCRICOES CONSTANTES DO DOCUMENTO REQUISITANTE. SERVICO DE LOCACAO DE IMOVEL TIPO GALPAO COM NO MINIMO 500 M², PARA FUNCIONAMENTO DO ARQUIVO GERAL DA DR/SE TRATA-SE DE PROCESSO DE RENOVACAO DO CONTRATO DE LOCACAO DE IMOVEL NAO RESIDENCIAL, LOCALIZADO NA AV. COLETORA A, S/N - CJ. MARCOS FREIRE I - N. SRª DO SOCORRO/SE, CEP 49160-000, DESTINADO AO FUNCIONAMENTO DA AC SHOPPING PREMIO. TRATA-SE DE PROCESSO DE RENOVACAO DO CONTRATO DE LOCACAO DE IMOVEL NAO RESIDENCIAL, LOCALIZADO NA AV. RUA F, Nº 70 E 70A, QD. 04, LOT. CHALES, BAIRRO AEROPORTO, ARACAJU/SE, CEP 49037-130, DESTINADO AO FUNCIONAMENTO DA CDD ATALAIA. Na outra tabela tenho o campo CDD ATALAIA 1015 DR/SE 0003 N. SRª DO SOCORRO 2105 CDD ATALAIA 1015 A minha chave vai passar a ser o nome da unidade. Alguém tem alguma sugestão de como fazer isso no R? Grato. Ronaldo Câmara de Araújo [https://mail.google.com/mail/?ui=2&ik=a38f3a957b&view=att&th=1335ff14b71079…]

2 1

Látice parcialmente balanceado em vários locais
by Euriann Yamamoto 17 Out '17

17 Out '17

Bom dia. Eu tenho 5 experimentos em látice simples 6x6 em 5 locais. Gostaria de saber se já existe alguma rotina implementada para rodar ANOVA individual e conjunta. Meu estudo é sobre adaptabilidade e e estabilidade em milho. Alguém poderia me auxiliar? Att., Euriann L. M. Yamamoto Doutoranda em Agronomia Universidade Federal da Grande Dourados

2 2

Ler arquivo exportado do SQL 2016
by Elias Carvalho 17 Out '17

17 Out '17

Pessoal Recebi um arquivo exportado do Microsoft SQL 2016 que tem extensão .txt e o mesmo com .csv e tem essa cara: Procurei, mas ainda não achei nada, alguém tem alguma expriência em importar esse arquivo para o R? Um simples read.csv por exemplo não funciona -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

1 0

Re: [R-br] Identificação das variáveis resposta em análise de correlação canônica
by sznelwar＠uol.com.br 16 Out '17

16 Out '17

2 1

Criar função CENT, complementar ao HEAD & TAIL...
by andrebvs＠bol.com.br 12 Out '17

12 Out '17

2 2

Identificação das variáveis resposta em análise de correlação canônica
by ASANTOS 11 Out '17

11 Out '17

Boa tarde Membros, Estou trabalhando com análise de correlação canônica com a função cancor() do pacote candisc, por exemplo e segundo a função, se eu faço: data(Rohwer, package="heplots") X <- as.matrix(Rohwer[,6:10]) # the PA tests Y <- as.matrix(Rohwer[,3:5]) # the aptitude/ability variables (cc <- cancor(X, Y, set.names=c("PA", "Ability")) ## 5 PA variables: n, s, ns, na, ss ## with 3 Ability variables: SAT, PPVT, Raven ## ## CanR CanRSQ Eigen percent cum scree ## 1 0.6703 0.44934 0.81599 77.30 77.30 ****************************** ## 2 0.3837 0.14719 0.17260 16.35 93.65 ****** ## 3 0.2506 0.06282 0.06704 6.35 100.00 ** Eu tenho a variância explicada com o aumento das dimensões das variáveis resposta, mas a função não identifica quem esta em 1 e contribuí com 77,30%, quais estão em 2 e contribuem com 93,65%, etc. Gostaria de saber como faço para descobrir são as 3 variáveis respostas envolvidas em cada valor dado pelo teste. Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

1 0

Como referenciar o pacote Parallel
by Elias Carvalho 10 Out '17

10 Out '17

Boa noite pessoal Estou usando o pacote parallel em um artigo, alguém sabe como devo referenciá-lo ? Procurei no google e nao encontrei nada. -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 2