- R-br - DINF UFPR

Função com vários if e else if que não funciona
by ASANTOS 08 Fev '16

08 Fev '16

Caros Listeiros, Estou tentando criar uma função abaixo para converter um SpatialGridDataFrame originado de uma imagem de satélite para reflectância, mas a função que eu criei (radconvL()) pelo excesso de else if não esta funcionando. Meu problema é que não dá erro, o que normalmente ajuda a detectar o problema, mas no meu caso a função nem funciona, segue meu CRM abaixo, alguém poderia me ajudar? require(raster) require(sp) # RasterLayer inventado r <- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) band2<- as(r, 'SpatialGridDataFrame') ### transforma em SpatialGridDataFrame ##Função para conversão DN para reflectância radconvL<-function(x, band = 2) { BAND <- c("2", "3", "4", "5") band <- pmatch(band, BAND) Lmin = 0 Qmax = 127 if (is.na(band)) stop("invalid band") if (band == 2) { Lmax = 120.64 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 3) { Lmax = 151.31 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if (band == 4) { Lmax = 157.57 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } else if(band == 5) { Lmax = 69.03 results <- x x <- as.vector(as.matrix(x)) x <- Lmin + ((Lmax-Lmin)*x)/Qmax if (class(results) == "SpatialGridDataFrame") results@data[, 1] <- x else if (is.data.frame(x)) results <- data.frame(matrix(x, nrow = nrow(results), ncol = ncol(results))) else results <- x results } } ###### teste<-radconvL(band2, band=2) ## Aqui não funciona e não sei porque teste Muito 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

7 19

Função de alisamento de Henderson
by Daniel Marcelino 07 Fev '16

07 Fev '16

Olá, eu estou querendo fazer um teste com a função de alisamento para médias móveis de Henderson, mas pesquisando no R e Google eu não a encontrei. Na verdade, eu enontrei no pacote `pracma`, mas não funciona na verdade. Alguém sabe de algum pacote ou implementação avulsa dessa função? Eu econtrei apenas notas e demonstrações: http://www.ausstats.abs.gov.au/ausstats/free.nsf/0/D829449095630207CA256D78… Daniel

1 0

Grafico Open High Loss Close (OHLC) de barras
by regis barros 05 Fev '16

05 Fev '16

Bom dia PessoalGostaria de saber se o R consegue plotar o gráfico de barras OHLC para cotações de valores de ações.E Como poderia fazer este gráfico no R. Regis

4 4

(sem assunto)
by Rudiney S Pereira 04 Fev '16

04 Fev '16

Bom dia! Preciso fazer a carga de imagens formato Geotiff e, não consigo pois um dos pacotes rgdal, não está disponível na lista para a instalação. Uso Linux Distribuição DEBIAN, SID. A versão do R é a última. Agradeço desde já! Rudiney -- Prof. Dr. Rudiney Soares Pereira UFSM - CCR - DER - Núcleo de Geotecnologias Prédio 44 J, Sala 214 RS 509 KM 9 - Campus Universitário 97.105-900 Santa Maria RS Brasil site: http://www.ufsm.br/geosere skype: rudiney.s.pereira fone: 55 3220-9468 55 9118-4142

2 2

Gráfico Open High Low Close
by regis barros 04 Fev '16

04 Fev '16

Bom dia PessoalGostaria de saber como posso plotar um gráfico de abertura, maximo, minimo e fechamento, referente a cotação de ações na bolsa de valores.Grato Regis Godoy BarrosLicenciatura em Física

2 1

Vetor com vários formatos de datas
by ASANTOS 03 Fev '16

03 Fev '16

Caros Listeiros, Tenho um vetor com diferentes formatos de datas abaixo, que se tento tirar a média não funciona e se tento explicitar um format=%d/%m/%Y", acabo deixando de fora o %d-%m-%Y e %d.%m.%Y. Aonde pergunto é possível trabalhar com um vetor de datas em diferentes formatos? Segue CRM: datas<-c("12-10-2011","10.03.2013","05/06/2016" ) time_a<-range(as.Date(datas)) time_a time_b<-range(as.Date(tableFE[,11], format = "%d/%m/%Y")) time_b 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 3

I Congresso de Aplicações de R em Administração
by Rodrigo Coster 03 Fev '16

03 Fev '16

Para quem interessar (via http://cara.coppead.ufrj.br/home) O I Congresso de Aplicações de R em Administração – CARA – visa disseminar o universo de estudos na área de Administração que utilizam a linguagem R como forma de alavancar a eficiência e qualidade dos trabalhos elaborados. Profissionais dos setores público e privado, de distintas áreas como Finanças e Educação, ministrarão palestras sobre a utilização do R no seu dia a dia. Haverá espaço para exposição de trabalhos das diversas áreas de Administração, que tenham sido desenvolvidos utilizando a linguagem R. O foco não deverá ser a discussão do código desenvolvido, mas sim a análise proposta, e do poder estatístico oferecido pela linguagem. Interessados devem submeter até três artigos pelo sistema disponível neste site, observando o calendário do evento. O CARA é destinado tanto a profissionais que nunca escreveram um código em R, bem como a profissionais experientes que desejarem conhecer trabalhos que têm sido desenvolvidos utilizando-se desta poderosa linguagem. Serão também oferecidas duas atividades: um workshop Análise de Business Intelligence, para profissionais que possuem um conhecimento básico de R, e um minicurso Apresentação do R, voltado àqueles que desejarem ter o primeiro contato com a linguagem. O evento ocorrerá nos dias 13 a 15 de abril de 2016, no Instituto COPPEAD de Administração. []'s

1 0

Adicionar um vetor a uma tabela mysql
by José Jesus Filho 03 Fev '16

03 Fev '16

Olá Pessoal, Esto começando a trabalhar com SQL a partir do R, então perdoem-me se cometer algum erro estúpido. Eu exportei do R uma tabela chamada "legislacao" para uma base de dados no MySQL usando o pacote RMySQL. Olá pessoal, Acontece que eu tenho um vetor "seguranca" com o mesmo número de linhas que a tabela "legislacao" e gostaria de adicioná-lo a essa tabela no MySQL como uma nova caluna sem ter de exportá-la novamente. No R eu simplesmente faria isso: legislacao$seguranca<-seguranca Já com o RMySQL, eu não sei, pois só encontrei "statements" para adicionar linhas ,não colunas. Claro é possível criar uma coluna com o "statement" ALTER TABLE legislacao ADD seguranca, mas isso cria uma coluna vazia com este nome. Eu quero preenchê-la com os valores correspondente. Obrigado, José.

4 5

Problema com datas que se tornam numéricas
by ASANTOS 02 Fev '16

02 Fev '16

Caros Listeiros, Desenvolvi uma função que extrai dentre várias informações de uma tabela, a menor e a maior data, meu problema é que mesmo especificando a data como as.Date(), ela fica numérica após a execução da função, alguém poderia me dar uma luz? Segue CRM: ## Função para leitura da tabela HTML readFE <- function(url0){ require(RCurl); require(XML) url1 <- gsub("(^.*)(/s/.*)(\\?.*$)", "https://dl.dropboxusercontent.com\\2", url0); url1 page <- getURL(url1) # page <- getURL(url1, ssl.verifypeer = FALSE) ### SSL certificate problem classes <- c(rep("factor", 6), rep("numeric", 4), "factor") return(readHTMLTable(page, head=T, colClasses=classes)[[1]]) } tableFE<- readFE("https://www.dropbox.com/s/foeh9oqbiqlawvg/TALHOES_AGENTE2.htm?dl=1") ## Sumariza as informações do banco de dados summaryFE<-function(x, db=tableFE){ ## Função Identificadora da empresa enterprise<-function(x, db = tableFE) { x<-NULL results <- x x<-db[1,1] results <- x results } #--# #Exemplo ENTERP<-enterprise(db = tableFE) # ## Função do período do monitoramento timeFE<-function(x, db = tableFE) { x<-NULL minT<-min(as.Date(db[,11])) maxT<-max(as.Date(db[,11])) results <- x x<-cbind(as.Date(minT),as.Date(maxT)) colnames(x)<-c("Ínicio","Fim") results <- x results } #--# ##Enumera o número de grupos de pragas monitorados npragsFE<-function(x){ npragsFE<-unique(db[,6]) } # result<-list(enterprise(x),timeFE(x),npragsFE(x)) return(result) } # ## Sumariza summaryFE(tableFE) -- ====================================================================== 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

2 2

problema biblioteca netcdf
by Mateus Dias Nunes 02 Fev '16

02 Fev '16

Olá Quando tento chamar a biblioteca "library(ncdf)" resulta no erro: > library(ncdf) Error in library(ncdf) : there is no package called ‘ncdf’ ai quando vou instalar o pacote, apresenta-se esta informação: package ‘ncdf’ is not available (for R version 3.2.3) como eu poderia resolver este problema. att ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 81125154 ____________________________________________________________________________

3 2

Criar função para tirar a média entre grupos sob condições dadas
by ASANTOS 31 Jan '16

31 Jan '16

Caros listeiros, Estou tentando criar uma função aggPestFE que para uma primeira variável se for "Formiga", vai tirar a media da coluna 5, mas se for "Besouros" ou "Barata" vai tira a média da coluna 4. Sendo que quero o resultado dado em forma de list, ficando o meu resultado: [[1]] [1] "Barata" [[1]][[3]] Projeto Talhao Injuria(%) 1 Abre Campo 15 25 2 Volta Redonda 73 0 3 Vitoria 74 50 [[2]] [1] "Besouros" [[1]][[3]] Projeto Talhao Injuria(%) 1 Abre Campo 15 0 2 Volta Redonda 73 25 3 Vitoria 74 50 [[2]] [1] "Formiga" [[1]][[3]] Projeto Talhao Area (m2) 1 Abre Campo 15 12.5 2 Volta Redonda 73 1.5 3 Vitoria 74 23.8 Para isso estou tentando: ##Dados artificiais Bichos<-c("Besouros","Barata","Formiga","Barata","Formiga","Besouros","Formiga", "Besouros","Barata","Formiga","Barata","Formiga","Besouros","Formiga","Formiga") Talhao<-c(73,15,74,75,15,15,15,73,15,15,73,15,73,74,74) Projeto<-c("Abre Campo", "Volta Redonda","Vitoria","Abre Campo", "Volta Redonda", "Vitoria","Abre Campo", "Volta Redonda","Vitoria","Abre Campo", "Volta Redonda", "Vitoria","Abre Campo", "Volta Redonda","Vitoria") Injuria<-c(25,100,0,25,0,100,0,50,25,0,25,0,25,0,0) Area<-c(0,0,12.5,0,7.5,0,1.5,0,0,23.8,0,5.3,0,2.0,11.3) tableFE<-data.frame(Bichos, Talhao, Projeto, Injuria, Area) ## Banco de dados criado ##Função aggPestFE<-function(x, db=tableFE){ mylist <- list() if (tableFE[,1] != "Formigas") { for (i in length(tableFE[,1])){ GP_FE2<-tableFE[tableFE[,1]==tableFE[i],] aggdata <-aggregate(GP_FE2[,4], list(GP_FE2[,3],GP_FE2[,2]), mean) colnames(aggdata)<-c("Regional","Projeto","Talhão","Injúria média (%)") tmp <- list(aggdata) mylist[NPRAGS[i]] <- tmp result<-mylist return(result) } else if (tableFE[,1] == "Formigas") { GP_FE2<-[tableFE[,1] aggdata <-aggregate(GP_FE2[,5], list(GP_FE2[,3],GP_FE2[,2]), mean) colnames(aggdata)<-c("Regional","Projeto","Talhão","Injúria média (%)") tmp <- list(aggdata) mylist[NPRAGS[i]] <- tmp result<-mylist return(result) } } RES<-list(result,aggdata3) result<-RES return(result) } #Teste aggPestFE(tableFE) # Sem sucesso, alguém poderia me ajudar? Obrigado e abraços, -- ====================================================================== 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

5 23

Separar strings no R
by Wagner Tassinari 30 Jan '16

30 Jan '16

Saudações grupo, Tenho uma variável que representa o nome completo de uma pessoa, ex: Manuel dos Santos de Jesus. Gostaria de separar em duas variáveis o nome e o úlltimo sobrenome respectivamente, por exemplo: VarNome = Manuel VarSobrenome = Jesus. Agradeço muitíssimo a todos, -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

6 6

Erro função lme
by Victor Fioreze 29 Jan '16

29 Jan '16

Olá retorno aqui com a minha duvida anterior: Estou rodando um script com a função lme do pacote nlme e esta dando um erro.Após eu rodar o modelo usando lme aparece: Erro em MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 Revisei os dados e o script e esta tudo ok.O que pode estar acontecendo? Att, pois não sei responder no fórum.... O arquivoa ser carregado esta no anexo Segue o meu comando para que alguém possa tentar> cruza <- read.table("C:/Users/Victor/Dropbox/projeto bezerras/Dados/cernelha.txt", header=T,)> cruza Bezerro Tratamento ACd Data DN DJN Sexo ACD Idade1 423 1 81.33 41937 41930 41948 FALSE 81.0 672 423 1 81.59 41944 41930 41948 FALSE 81.0 743 423 1 82.07 41951 41930 41948 FALSE 81.0 814 423 1 83.62 41958 41930 41948 FALSE 81.0 885 423 1 83.72 41965 41930 41948 FALSE 81.0 906 425 1 81.38 41937 41930 41962 FALSE 81.0 677 425 1 81.95 41944 41930 41962 FALSE 81.0 748 425 1 82.40 41951 41930 41962 FALSE 81.0 819 425 1 83.00 41958 41930 41962 FALSE 81.0 8810 425 1 83.07 41965 41930 41962 FALSE 81.0 9011 426 1 81.09 41937 41930 41968 FALSE 81.0 6712 426 1 81.69 41944 41930 41968 FALSE 81.0 7413 426 1 82.29 41951 41930 41968 FALSE 81.0 8114 426 1 83.20 41958 41930 41968 FALSE 81.0 8815 426 1 83.28 41965 41930 41968 FALSE 81.0 9016 428 1 81.69 41937 41930 41980 FALSE 81.0 6717 428 1 82.32 41944 41930 41980 FALSE 81.0 7418 428 1 83.34 41951 41930 41980 FALSE 81.0 8119 428 1 83.65 41958 41930 41980 FALSE 81.0 8820 428 1 83.74 41965 41930 41980 FALSE 81.0 9021 421 2 82.72 41937 41930 41931 FALSE 79.5 6722 421 2 84.07 41944 41930 41931 FALSE 79.5 7423 421 2 84.98 41951 41930 41931 FALSE 79.5 8124 421 2 85.77 41958 41930 41931 FALSE 79.5 8825 421 2 85.99 41965 41930 41931 FALSE 79.5 9026 429 2 81.29 41937 41930 41989 FALSE 81.0 6727 429 2 83.21 41944 41930 41989 FALSE 81.0 7428 429 2 84.47 41951 41930 41989 FALSE 81.0 8129 429 2 85.75 41958 41930 41989 FALSE 81.0 8830 429 2 85.91 41965 41930 41989 FALSE 81.0 9031 430 2 85.61 41937 41930 42007 FALSE 83.5 6732 430 2 87.20 41944 41930 42007 FALSE 83.5 7433 430 2 87.53 41951 41930 42007 FALSE 83.5 8134 430 2 87.87 41958 41930 42007 FALSE 83.5 8835 430 2 88.02 41965 41930 42007 FALSE 83.5 9036 438 2 90.66 41937 41930 42063 FALSE 89.0 6737 438 2 91.99 41944 41930 42063 FALSE 89.0 7438 438 2 92.80 41951 41930 42063 FALSE 89.0 8139 438 2 94.07 41958 41930 42063 FALSE 89.0 8840 438 2 94.25 41965 41930 42063 FALSE 89.0 9041 422 3 83.43 41937 41930 41933 FALSE 83.0 6742 422 3 84.17 41944 41930 41933 FALSE 83.0 7443 422 3 84.65 41951 41930 41933 FALSE 83.0 8144 422 3 84.91 41958 41930 41933 FALSE 83.0 8845 422 3 84.97 41965 41930 41933 FALSE 83.0 9046 427 3 79.06 41937 41930 41978 FALSE 78.4 6747 427 3 79.86 41944 41930 41978 FALSE 78.4 7448 427 3 80.65 41951 41930 41978 FALSE 78.4 8149 427 3 81.66 41958 41930 41978 FALSE 78.4 8850 427 3 81.78 41965 41930 41978 FALSE 78.4 9051 432 3 82.79 41937 41930 42018 FALSE 80.2 6752 432 3 83.80 41944 41930 42018 FALSE 80.2 7453 432 3 84.64 41951 41930 42018 FALSE 80.2 8154 432 3 85.60 41958 41930 42018 FALSE 80.2 8855 432 3 85.79 41965 41930 42018 FALSE 80.2 9056 435 3 89.95 41937 41930 42025 FALSE 88.5 6757 435 3 90.69 41944 41930 42025 FALSE 88.5 7458 435 3 91.34 41951 41930 42025 FALSE 88.5 8159 435 3 92.30 41958 41930 42025 FALSE 88.5 8860 435 3 92.43 41965 41930 42025 FALSE 88.5 9061 424 4 82.57 41937 41930 41954 FALSE 82.5 6762 424 4 83.02 41944 41930 41954 FALSE 82.5 7463 424 4 83.17 41951 41930 41954 FALSE 82.5 8164 424 4 83.33 41958 41930 41954 FALSE 82.5 8865 424 4 83.35 41965 41930 41954 FALSE 82.5 9066 431 4 82.99 41937 41930 42013 FALSE 81.0 6767 431 4 85.85 41944 41930 42013 FALSE 81.0 7468 431 4 87.54 41951 41930 42013 FALSE 81.0 8169 431 4 89.33 41958 41930 42013 FALSE 81.0 8870 431 4 89.62 41965 41930 42013 FALSE 81.0 9071 433 4 84.94 41937 41930 42020 FALSE 81.4 6772 433 4 86.87 41944 41930 42020 FALSE 81.4 7473 433 4 88.33 41951 41930 42020 FALSE 81.4 8174 433 4 88.99 41958 41930 42020 FALSE 81.4 8875 433 4 89.25 41965 41930 42020 FALSE 81.4 9076 436 4 87.64 41937 41930 42036 FALSE 87.3 6777 436 4 88.81 41944 41930 42036 FALSE 87.3 7478 436 4 89.17 41951 41930 42036 FALSE 87.3 8179 436 4 90.48 41958 41930 42036 FALSE 87.3 8880 436 4 90.58 41965 41930 42036 FALSE 87.3 90 str(cruza)'data.frame': 80 obs. of 9 variables: $ Bezerro : int 423 423 423 423 423 425 425 425 425 425 ... $ Tratamento: int 1 1 1 1 1 1 1 1 1 1 ... $ ACd : num 81.3 81.6 82.1 83.6 83.7 ... $ Data : int 41937 41944 41951 41958 41965 41937 41944 41951 41958 41965 ... $ DN : int 41930 41930 41930 41930 41930 41930 41930 41930 41930 41930 ... $ DJN : int 41948 41948 41948 41948 41948 41962 41962 41962 41962 41962 ... $ Sexo : logi FALSE FALSE FALSE FALSE FALSE FALSE ... $ ACD : num 81 81 81 81 81 81 81 81 81 81 ... $ Idade : int 67 74 81 88 90 67 74 81 88 90 ... library(nlme)> attach(cruza)> fTrat<-as.factor(Tratamento)> class(fTrat)[1] "factor"> class(ACd)[1] "numeric"> class(ACD)[1] "numeric"> class(Data)[1] "integer"> Data<-as.factor(Data)> simcomp.ran <- lme(ACd ~ DJN + ACD + fTrat + Data + fTrat:Data + ACD*fTrat + fTrat*DJN,+ random= ~1 | Bezerro/fTrat)Erro em MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 Muito obrigado pela ajuda e desculpe não saber responder no fórum. Abraço a todos! Victor Ionatan FiorezeMédico Veterinário e ZootecnistaMestrando em Nutrição de Ruminantes, PPGZ-UFPelBolsista Capes Integrante GERUMEN (53) 81032887lattes: http://lattes.cnpq.br/5026587958713372

2 1

Aumentar a memória do R
by Wagner Tassinari 28 Jan '16

28 Jan '16

Olá pessoal, Como faço para liberar mais memoria para o R rodar ? Estou utilizando os comandos: > memory.size() [1] 3306.77 > memory.size(TRUE) [1] 3995 > memory.limit() [1] 32676 Obrigado -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

3 2

Re: [R-br] Dúvida em looping e data.frame
by Paulo Henrique Pimenta 28 Jan '16

28 Jan '16

Usei a função simulacao <- function(proporcao, moedas){ for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) a <- print(moedas(n,p)) return(a) }b <- simulacao(proporcao, moedas)No final ele possui somente um vetor com os dados gerados com n = 30 e p = 0.5 e não todos os gerados pelas combinações. Como eu poderia ter uma vetor ou alguma outra forma de armazenar estes dados? Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435. > Message: 11 > Date: Wed, 27 Jan 2016 09:33:10 +0300 > From: Paulo Henrique Pimenta <paulopimenta6(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Dúvida em looping e data.frame > Message-ID: <BAY182-W241C2B8A79DFC8287003AFF4D90(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > Olá a todos > Estou escrevendo um código que calcula o lance de uma moeda jogada e a probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, num bloco do meu código fiquei com uma dúvida relacionada a vetorização. Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > ####Laboratorio 1 de R > ###cara eh 1###coroa eh 0vetor_moedas <- NULLvetor <- NULLa <- NULL > moeda <- function(p){ > if(runif(1) < p){ > y <- 1 > }else{ > y <- 0 } return(y) } ####separando uma funcao da outra#### > moedas <- function(numeros,probabilidade){ > for (i in 1:numeros){ > vetor[i] <- moeda(probabilidade) > } > return(vetor) > } > > ####separando uma funcao da outra#### > proporcao <- function(moedas){ > numeros_1 <- sum(vetor_moedas) > numeros_totais <- length(vetor_moedas) > proporcao_final <- numeros_1/numeros_totais > return(proporcao_final) > } > ####separando uma funcao da outra#### > simulacao <- function(proporcao, moedas){ > ###variando o "n" e o "p" > for (i in 1:500){ > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como data.frame todos os dados feitos pela função moedas, como fazer com que a2 também aloque um outro data.frame que aloque informações da função moedas e assim por diante... > n <- c(10, 20, 30, 40, 50) p <- c(0.2, 0.3, 0.4, 0.5) a[i] <- data.frame(moedas(n[i],p[i])) > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > } > > return(a) > } > ###############################################################Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda###lancamento de 10 moedas - funcao moedas > ###comeco das chamas > #########################################prob <- 0.4 > moeda(prob) > #########################################n_1 <- 100 > moedas(n_1, prob) > ######################################### > vetor_moedas <- moedas(n_1, prob) > ######################################### > proporcao(vetor_moedas) > ######################################### > simulacao(proporcao,moedas) > #################################### > Um grande abraço a todos!! > > Paulo Henrique de A. S. Pimenta. > Graduando em Meteorologia (Bacharelado) - IAG/USP. > Fone: +5511981318435. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160127/683e3d06/atta…> > > ------------------------------ > > Message: 12 > Date: Wed, 27 Jan 2016 10:57:57 -0300 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Dúvida em looping e data.frame > Message-ID: > <CABmC8gmxHtjGenrn6us=VnsgRE2E1B_zCY_nBZ+Vnyp1CZX5nw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Paulo, bom dia! > > A instrução for (i in 1:500) não parece coerente, pois me parece que você > deveria variar em função da iteração n x p. Ao variar n você irá criar > vetores de tamanhos diferentes, impossibilitando o armazenamento direto em > um data.frame(). Poderia utilizar uma lista. > > Você pode tentar algo como: > simulacao <- function(proporcao, moedas){for (n in c(10, 20, 30)) for (p in > c(0.3, 0.4, 0.5)) print(moedas(n,p))} > simulacao(proporcao, moedas) > > Outro ponto é que você poderia usar rbinom() para os sorteios, enxugando > seu código. Deixo uma ideia abaixo: > > ### <code r> > args(rbinom) ### function (n, size, prob) > vetor_moedas <- rbinom(100, 1, .4); vetor_moedas > sum(vetor_moedas)/length(vetor_moedas) ### proporção final > > for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) print(rbinom(n, 1,p)) > # [1] 0 1 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 0 0 1 > # [1] 0 1 1 0 0 0 0 1 0 1 > # [1] 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 > # [1] 0 1 0 0 0 1 0 1 0 0 0 1 1 0 0 1 0 0 0 0 > # [1] 0 0 0 1 0 1 1 1 0 1 1 0 1 0 0 1 1 1 1 1 > # [1] 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 1 1 1 0 1 0 1 0 0 1 1 1 0 0 0 > # [1] 1 0 0 1 1 1 1 1 1 1 0 0 1 1 0 1 1 1 1 0 0 1 1 1 0 1 1 1 1 0 > ### </code> > > > > ================================================ > Éder Comunello > PhD Student in Agricultural Systems Engineering (USP/Esalq) > Brazilian Agricultural Research Corporation (Embrapa) > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > > > > > Em 27 de janeiro de 2016 03:33, Paulo Henrique Pimenta < > paulopimenta6(a)hotmail.com> escreveu: > > > Olá a todos > > > > Estou escrevendo um código que calcula o lance de uma moeda jogada e a > > probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, > > num bloco do meu código fiquei com uma dúvida relacionada a vetorização. > > Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > > > ####Laboratorio 1 de R > > > > ###cara eh 1 > > ###coroa eh 0 > > vetor_moedas <- NULL > > vetor <- NULL > > a <- NULL > > > > moeda <- function(p){ > > > > if(runif(1) < p){ > > > > y <- 1 > > > > }else{ > > > > y <- 0 > > > > } > > return(y) > > > > } > > > > ####separando uma funcao da outra#### > > > > moedas <- function(numeros,probabilidade){ > > > > for (i in 1:numeros){ > > > > vetor[i] <- > > moeda(probabilidade) > > > > } > > > > return(vetor) > > > > } > > > > > > ####separando uma funcao da outra#### > > > > proporcao <- function(moedas){ > > > > numeros_1 <- sum(vetor_moedas) > > > > numeros_totais <- length(vetor_moedas) > > > > proporcao_final <- numeros_1/numeros_totais > > > > return(proporcao_final) > > > > } > > > > ####separando uma funcao da outra#### > > > > simulacao <- function(proporcao, moedas){ > > > > ###variando o "n" e o "p" > > > > for (i in 1:500){ > > > > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como > > data.frame todos os dados feitos pela função moedas, como fazer com que a2 > > também aloque um outro data.frame que aloque informações da função moedas e > > assim por diante... > > > > n <- c(10, 20, 30, 40, 50) > > p <- c(0.2, 0.3, 0.4, 0.5) > > > > a[i] <- data.frame(moedas(n[i],p[i])) > > > > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e > > recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > > > > } > > > > > > return(a) > > > > } > > > > ############################## > > ############################## > > ###Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda > > ###lancamento de 10 moedas - funcao moedas > > > > ###comeco das chamas > > > > ######################################### > > prob <- 0.4 > > > > moeda(prob) > > > > ######################################### > > n_1 <- 100 > > > > moedas(n_1, prob) > > > > ######################################### > > > > vetor_moedas <- moedas(n_1, prob) > > > > ######################################### > > > > proporcao(vetor_moedas) > > > > ######################################### > > > > simulacao(proporcao,moedas) > > > > ### > > ### > > ############################## > > > > > > Um grande abraço a todos!! > > > > Paulo Henrique de A. S. Pimenta. > > > > Graduando em Meteorologia (Bacharelado) - IAG/USP. > > Fone: +5511981318435.

2 2

Re: [R-br] Digest R-br, volume 61, assunto 17
by Paulo Henrique Pimenta 27 Jan '16

27 Jan '16

Usei a função simulacao <- function(proporcao, moedas){ for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) a <- print(moedas(n,p)) return(a) } b <- simulacao(proporcao, moedas) No final ele possui somente um vetor com os dados gerados com n = 30 e p = 0.5 e não todos os gerados pelas combinações. Como eu poderia ter uma vetor ou alguma outra forma de armazenar estes dados? Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435. > > Message: 11 > Date: Wed, 27 Jan 2016 09:33:10 +0300 > From: Paulo Henrique Pimenta <paulopimenta6(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Dúvida em looping e data.frame > Message-ID: <BAY182-W241C2B8A79DFC8287003AFF4D90(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > Olá a todos > Estou escrevendo um código que calcula o lance de uma moeda jogada e a probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, num bloco do meu código fiquei com uma dúvida relacionada a vetorização. Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > ####Laboratorio 1 de R > ###cara eh 1###coroa eh 0vetor_moedas <- NULLvetor <- NULLa <- NULL > moeda <- function(p){ > if(runif(1) < p){ > y <- 1 > }else{ > y <- 0 } return(y) } ####separando uma funcao da outra#### > moedas <- function(numeros,probabilidade){ > for (i in 1:numeros){ > vetor[i] <- moeda(probabilidade) > } > return(vetor) > } > > ####separando uma funcao da outra#### > proporcao <- function(moedas){ > numeros_1 <- sum(vetor_moedas) > numeros_totais <- length(vetor_moedas) > proporcao_final <- numeros_1/numeros_totais > return(proporcao_final) > } > ####separando uma funcao da outra#### > simulacao <- function(proporcao, moedas){ > ###variando o "n" e o "p" > for (i in 1:500){ > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como data.frame todos os dados feitos pela função moedas, como fazer com que a2 também aloque um outro data.frame que aloque informações da função moedas e assim por diante... > n <- c(10, 20, 30, 40, 50) p <- c(0.2, 0.3, 0.4, 0.5) a[i] <- data.frame(moedas(n[i],p[i])) > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > } > > return(a) > } > ###############################################################Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda###lancamento de 10 moedas - funcao moedas > ###comeco das chamas > #########################################prob <- 0.4 > moeda(prob) > #########################################n_1 <- 100 > moedas(n_1, prob) > ######################################### > vetor_moedas <- moedas(n_1, prob) > ######################################### > proporcao(vetor_moedas) > ######################################### > simulacao(proporcao,moedas) > #################################### > Um grande abraço a todos!! > > Paulo Henrique de A. S. Pimenta. > Graduando em Meteorologia (Bacharelado) - IAG/USP. > Fone: +5511981318435. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160127/683e3d06/atta…> > > ------------------------------ > > Message: 12 > Date: Wed, 27 Jan 2016 10:57:57 -0300 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Dúvida em looping e data.frame > Message-ID: > <CABmC8gmxHtjGenrn6us=VnsgRE2E1B_zCY_nBZ+Vnyp1CZX5nw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Paulo, bom dia! > > A instrução for (i in 1:500) não parece coerente, pois me parece que você > deveria variar em função da iteração n x p. Ao variar n você irá criar > vetores de tamanhos diferentes, impossibilitando o armazenamento direto em > um data.frame(). Poderia utilizar uma lista. > > Você pode tentar algo como: > simulacao <- function(proporcao, moedas){for (n in c(10, 20, 30)) for (p in > c(0.3, 0.4, 0.5)) print(moedas(n,p))} > simulacao(proporcao, moedas) > > Outro ponto é que você poderia usar rbinom() para os sorteios, enxugando > seu código. Deixo uma ideia abaixo: > > ### <code r> > args(rbinom) ### function (n, size, prob) > vetor_moedas <- rbinom(100, 1, .4); vetor_moedas > sum(vetor_moedas)/length(vetor_moedas) ### proporção final > > for (n in c(10, 20, 30)) for (p in c(0.3, 0.4, 0.5)) print(rbinom(n, 1,p)) > # [1] 0 1 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 0 0 1 > # [1] 0 1 1 0 0 0 0 1 0 1 > # [1] 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 > # [1] 0 1 0 0 0 1 0 1 0 0 0 1 1 0 0 1 0 0 0 0 > # [1] 0 0 0 1 0 1 1 1 0 1 1 0 1 0 0 1 1 1 1 1 > # [1] 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 0 > # [1] 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 1 1 1 0 1 0 1 0 0 1 1 1 0 0 0 > # [1] 1 0 0 1 1 1 1 1 1 1 0 0 1 1 0 1 1 1 1 0 0 1 1 1 0 1 1 1 1 0 > ### </code> > > > > ================================================ > Éder Comunello > PhD Student in Agricultural Systems Engineering (USP/Esalq) > Brazilian Agricultural Research Corporation (Embrapa) > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > > > > > Em 27 de janeiro de 2016 03:33, Paulo Henrique Pimenta < > paulopimenta6(a)hotmail.com> escreveu: > > > Olá a todos > > > > Estou escrevendo um código que calcula o lance de uma moeda jogada e a > > probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, > > num bloco do meu código fiquei com uma dúvida relacionada a vetorização. > > Abaixo passarei todo o código e irei destacar o que me traz a dúvida: > > > > ####Laboratorio 1 de R > > > > ###cara eh 1 > > ###coroa eh 0 > > vetor_moedas <- NULL > > vetor <- NULL > > a <- NULL > > > > moeda <- function(p){ > > > > if(runif(1) < p){ > > > > y <- 1 > > > > }else{ > > > > y <- 0 > > > > } > > return(y) > > > > } > > > > ####separando uma funcao da outra#### > > > > moedas <- function(numeros,probabilidade){ > > > > for (i in 1:numeros){ > > > > vetor[i] <- > > moeda(probabilidade) > > > > } > > > > return(vetor) > > > > } > > > > > > ####separando uma funcao da outra#### > > > > proporcao <- function(moedas){ > > > > numeros_1 <- sum(vetor_moedas) > > > > numeros_totais <- length(vetor_moedas) > > > > proporcao_final <- numeros_1/numeros_totais > > > > return(proporcao_final) > > > > } > > > > ####separando uma funcao da outra#### > > > > simulacao <- function(proporcao, moedas){ > > > > ###variando o "n" e o "p" > > > > for (i in 1:500){ > > > > \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como > > data.frame todos os dados feitos pela função moedas, como fazer com que a2 > > também aloque um outro data.frame que aloque informações da função moedas e > > assim por diante... > > > > n <- c(10, 20, 30, 40, 50) > > p <- c(0.2, 0.3, 0.4, 0.5) > > > > a[i] <- data.frame(moedas(n[i],p[i])) > > > > \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e > > recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! > > > > } > > > > > > return(a) > > > > } > > > > ############################## > > ############################## > > ###Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda > > ###lancamento de 10 moedas - funcao moedas > > > > ###comeco das chamas > > > > ######################################### > > prob <- 0.4 > > > > moeda(prob) > > > > ######################################### > > n_1 <- 100 > > > > moedas(n_1, prob) > > > > ######################################### > > > > vetor_moedas <- moedas(n_1, prob) > > > > ######################################### > > > > proporcao(vetor_moedas) > > > > ######################################### > > > > simulacao(proporcao,moedas) > > > > ### > > ### > > ############################## > > > > > > Um grande abraço a todos!! > > > > Paulo Henrique de A. S. Pimenta. > > > > Graduando em Meteorologia (Bacharelado) - IAG/USP. > > Fone: +5511981318435. > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160127/515a392f/atta…> > > ------------------------------ > > 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 61, assunto 17 > *****************************************

1 1

Validando o CPF
by Wagner Tassinari 27 Jan '16

27 Jan '16

Olá pessoal, Alguém conhece algum código em R para verificar a consistência do número do CPF ? Obrigado -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

2 1

Dúvida em looping e data.frame
by Paulo Henrique Pimenta 27 Jan '16

27 Jan '16

Olá a todos Estou escrevendo um código que calcula o lance de uma moeda jogada e a probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, num bloco do meu código fiquei com uma dúvida relacionada a vetorização. Abaixo passarei todo o código e irei destacar o que me traz a dúvida: ####Laboratorio 1 de R ###cara eh 1###coroa eh 0vetor_moedas <- NULLvetor <- NULLa <- NULL moeda <- function(p){ if(runif(1) < p){ y <- 1 }else{ y <- 0 } return(y) } ####separando uma funcao da outra#### moedas <- function(numeros,probabilidade){ for (i in 1:numeros){ vetor[i] <- moeda(probabilidade) } return(vetor) } ####separando uma funcao da outra#### proporcao <- function(moedas){ numeros_1 <- sum(vetor_moedas) numeros_totais <- length(vetor_moedas) proporcao_final <- numeros_1/numeros_totais return(proporcao_final) } ####separando uma funcao da outra#### simulacao <- function(proporcao, moedas){ ###variando o "n" e o "p" for (i in 1:500){ \\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como data.frame todos os dados feitos pela função moedas, como fazer com que a2 também aloque um outro data.frame que aloque informações da função moedas e assim por diante... n <- c(10, 20, 30, 40, 50) p <- c(0.2, 0.3, 0.4, 0.5) a[i] <- data.frame(moedas(n[i],p[i])) \\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar! } return(a) } ###############################################################Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda###lancamento de 10 moedas - funcao moedas ###comeco das chamas #########################################prob <- 0.4 moeda(prob) #########################################n_1 <- 100 moedas(n_1, prob) ######################################### vetor_moedas <- moedas(n_1, prob) ######################################### proporcao(vetor_moedas) ######################################### simulacao(proporcao,moedas) #################################### Um grande abraço a todos!! Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435.

2 1

Criar um convex hull para um objeto shapefile
by ASANTOS 26 Jan '16

26 Jan '16

Prezados Listeiros, Eu gostaria de criar um convex hull que envolvesse 4 polígonos, representados por um objeto oriundo da leitura de um readShapeLines(), sem ter que transformar o objeto em SpatialPolygonsDataFrame, isso é possível? Meu CRM: #Pacotes require(rgdal) require(maptools) #------------------------------------------------------------------------------- #Crio 4 polígonos sr <- SpatialPolygons(list( Polygons(list(Polygon(cbind(c(180114, 180553, 181127, 181477, 181294, 181007, 180409, 180162, 180114), c(332349, 332057, 332342, 333250, 333558, 333676, 332618, 332413, 332349)))),'1'), Polygons(list(Polygon(cbind(c(180042, 180545, 180553, 180314, 179955, 179142, 179437, 179524, 179979, 180042), c(332373, 332026, 331426, 330889, 330683, 331133, 331623, 332152, 332357, 332373)))),'2'), Polygons(list(Polygon(cbind(c(179110, 179907, 180433, 180712, 180752, 180329, 179875, 179668, 179572, 179269, 178879, 178600, 178544, 179046, 179110), c(331086, 330620, 330494, 330265, 330075, 330233, 330336, 330004, 329783, 329665, 329720, 329933, 330478, 331062, 331086)))),'3'), Polygons(list(Polygon(cbind(c(180304, 180403,179632,179420,180304), c(332791, 333204, 333635, 333058, 332791)))),'4'))) plot(sr) #Converto em polígono espacial srdf=SpatialPolygonsDataFrame(sr, data.frame(row.names=c('1','2','3','4'), PIDS=1:4)) srdf@data #Criação do shapefile para ilustrar o problema writeOGR(srdf, getwd(), 'POLY', 'ESRI Shapefile') #Leitura do shapefile contorno_line_X <- readShapeLines ("POLY.shp") #Plot plot(contorno_line_X) #Tentativa de criar o convex hull df.data = as.data.frame(contorno_line_X) ch <- chull(df.data) lines(ch, col="red") ##Não funciona 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 5

Pacote rJava
by Gislaine Batistela 26 Jan '16

26 Jan '16

Olá! Estou com dificuldade para carregar o pacote rJava. Por favor, alguém poderia me ajudar? Vejam as mensagens: > utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- trying URL 'http://www.vps.fmvz.usp.br/CRAN/bin/windows/contrib/3.0/rJava_0.9-6.zip' Content type 'application/zip' length 758820 bytes (741 Kb) opened URL downloaded 741 Kb package ‘rJava’ successfully unpacked and MD5 sums checked The downloaded binary packages are in C:\Users\Unesp\AppData\Local\Temp\Rtmpsx6PLo\downloaded_packages > local({pkg <- select.list(sort(.packages(all.available = TRUE)),graphics=TRUE) + if(nchar(pkg)) library(pkg, character.only=TRUE)}) Error : .onLoad failed in loadNamespace() for 'rJava', details: call: fun(libname, pkgname) error: JAVA_HOME cannot be determined from the Registry Error: package or namespace load failed for ‘rJava’ -- Gislaine Cristina Batistela

3 3

Dúvidas em condicional composto
by Edson Lira 26 Jan '16

26 Jan '16

Caros listeiros, tenho essa configuração em um banco eem dpm brpm 1 Alto Alto Médio 2 Alto Alto Alto 3 Alto Alto Médio 4 Alto Alto Alto 5 Alto Alto Médio 6 Alto Alto Baixo 7 Alto Alto Baixo Se duas dessas variáveis apresentarem o nível alto, digo que há um estresse no indivíduo. Por exemplo se eem ou dpm for alto, digo que há um estresse, basta uma cominação qualquer de duas ou três variáveis apresentando alto, digo que há estresse. eem=alto e dpm= alto ->estresseeem=alto e brpm= alto ->estressedpm=alto e brpm= alto ->estressedpm=medio e brpm= alto ->sem estresse Não sei como implementar no R. Alguém poderia me sugerir alguma rotina? [ ]'s Edson Lira Estatístico Manaus-Amazonas

2 2

PNS - pesquisa nacional de saude
by Elias Teixeira Krainski 26 Jan '16

26 Jan '16

Caros(as) Achei interessante esse post sobre a pesquisa nacional de saude: http://www.r-bloggers.com/analyze-the-pesquisa-nacional-de-saude-pns-with-r/ Qualquer um pode iniciar pelos scripts disponibilizados e analizar dados da PNS ao nivel de individuos. Abs, Elias

2 1

Frequência de missings
by Wagner Tassinari 22 Jan '16

22 Jan '16

Olá pessoal, Tenho uma variável character ou factor com NAs, daí eu gostaria de tabula-lo também para ver a frequência, eu uso o comando: Table(banco$var[banco$var=='NA']) e não roda. Obrigado

3 2

arquivo .rmd - Grafico da rede neural não aparece
by Leonard de Assis 22 Jan '16

22 Jan '16

Título da msg explica > sessionInfo() R version 3.2.3 (2015-12-10) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows >= 8 x64 (build 9200) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] nnet_7.3-11 caret_6.0-64 ggplot2_2.0.0.9001 lattice_0.20-33 [5] mlbench_2.1-1 RevoUtilsMath_3.2.3 loaded via a namespace (and not attached): [1] Rcpp_0.12.3 compiler_3.2.3 nloptr_1.0.4 plyr_1.8.3 [5] iterators_1.0.8 tools_3.2.3 digest_0.6.9 lme4_1.1-10 [9] nlme_3.1-124 gtable_0.1.2 mgcv_1.8-10 Matrix_1.2-3 [13] foreach_1.4.3 yaml_2.1.13 parallel_3.2.3 SparseM_1.7 [17] stringr_1.0.0.9000 knitr_1.12 MatrixModels_0.4-2 stats4_3.2.3 [21] grid_3.2.3 rmarkdown_0.9.4 minqa_1.2.4 reshape2_1.4.1.9000 [25] car_2.1-1 magrittr_1.5 scales_0.3.0 codetools_0.2-14 [29] htmltools_0.3 MASS_7.3-45 splines_3.2.3 rsconnect_0.4.1.4 [33] pbkrtest_0.4-5 colorspace_1.2-6 quantreg_5.19 stringi_1.0-1 [37] munsell_0.4.2 > ____ código --- title: "Neural Network example" author: "Leonard Mendonça de Assis" date: "21 de janeiro de 2016" output: pdf_document: highlight: pygments number_sections: yes html_document: highlight: pygments number_sections: yes word_document: highlight: pygments --- # Session start ```{r} require(MASS) require(neuralnet) require(plyr) require(boot) require(knitr) data <- Boston kable(head(data)) kable( apply(data,2,function(x) sum(is.na(x))), caption='NAs per variables' ) ``` # Adjusting a linear model ```{r} set.seed(500) index <- sample(1:nrow(data),round(0.75*nrow(data))) train <- data[index,] test <- data[-index,] lm.fit <- glm(medv~., data=train) summary(lm.fit) pr.lm <- predict(lm.fit,test) MSE.lm <- sum((pr.lm - test$medv)^2)/nrow(test) ``` # Adjusting a Neural Net As a first step, we are going to address data preprocessing. It is good practice to normalize your data before training a neural network. I cannot emphasize enough how important this step is: depending on your dataset, avoiding normalization may lead to useless results or to a very difficult training process (most of the times the algorithm will not converge before the number of maximum iterations allowed). You can choose different methods to scale the data (z-normalization, min-max scale, etc). I chose to use the min-max method and scale the data in the interval [0,1]. Usually scaling in the intervals [0,1] or [-1,1] tends to give better results. We therefore scale and split the data before moving on: ```{r} maxs <- apply(data, 2, max) mins <- apply(data, 2, min) scaled <- as.data.frame(scale(data, center = mins, scale = maxs - mins)) train_ <- scaled[index,] test_ <- scaled[-index,] ``` Note that scale returns a matrix that needs to be coerced into a data.frame. ```{r} n <- names(train_) f <- as.formula(paste("medv ~", paste(n[!n %in% "medv"], collapse = " + "))) nn <- neuralnet(f,data=train_,hidden=c(5,3),linear.output=T) ``` And the Neural Nework looks like this: ```{r} plot(nn) ``` ```{r} sessionInfo() ```

1 0

Re: [R-br] RES: Ler um HTM
by Éder Comunello 21 Jan '16

21 Jan '16

Leonardo, boa tarde! Voltando na discussão... Acredito que não tem como fazer a conexão direta com o arquivo .zip alocado na web, mas uma vez que você baixa o arquivo com download.file(), pode usar unz() para não precisar descompactar. Tem um post antigo em < https://stat.ethz.ch/pipermail/r-help/2008-October/176266.html>. Abaixo tem um exemplo bem simples, mas que pode ser útil. É feita a leitura de uma tabela .html a partir de um arquivo .zip alocado no Google Drive. ### <code r> setwd("D:\\TEMP") url.zip <- "https://drive.google.com/open?id=0B8y1jJExW6kabHU3SXV3dXB2OEk" url.zip <- sub("^.*\\?(id.*$)", " https://drive.google.com/uc?export=download&\\1", url.zip) url.zip ### https://drive.google.com/uc?export=download&id=0B8y1jJExW6kabHU3SXV3dXB2OEk" if (!file.exists("html.zip")) download.file(url.zip, "html.zip", mode="wb") unzip("html.zip", list=T) ### conteúdo closeAllConnections() con <- unz("html.zip", "htmlTable.html") XML::readHTMLTable(readLines(con)) # $`NULL` # First Name Last Name Points # 1 Jill Smith 50 # 2 Eve Jackson 94 close(con) ### </code> ================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] Em 7 de janeiro de 2016 10:54, Leonard de Assis <assis.leonard(a)gmail.com> escreveu: > Paulo, o que achei via google de ontem pra hoje vai nessa direção mesmo. > > > > Estou rodando vários exemplos que achei e, em paralelo, adaptando meu > código. > > > > Eu já consigo ler os trecos dentro do HTML, o que agarra é esse processo > de descompactar o htm dentro. > > > > Fiz um esqueminha em C usando curl e xml que funciona. Se ficar demorando > muito, porto essa joça pra ficar no formato que o R entende e resolvo por > esse caminho. > > > > Ou então, deixo como está, que é baixar o zip, descompactar e ler. > > > > *De:* R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] *Em nome de *Paulo > Nogueira Starzynski > *Enviada em:* quinta-feira, 7 de janeiro de 2016 09:39 > *Para:* R-BR <r-br(a)listas.c3sl.ufpr.br> > *Assunto:* Re: [R-br] Ler um HTM > > > > Leonardo, > > talvez você já tenha alguma vez trilhado o caminho das pedras do > webscrapping, mas, se não, recomendo os pacotes XML, RCurl e rvest (eu uso > mais os dois primeiros). > > > Abraços, > > Paulo > > > > Em 6 de janeiro de 2016 23:45, Leonard de Assis <assis.leonard(a)gmail.com> > escreveu: > > Boa noite a todos > > > > Estou com um problema idiota em mãos. Preciso ler um arquivo .HTM dentro > de um arquivo .zip. > > Até agora, o jeito que achei mais produtivo foi: > > 1. baixar o zip (posso fazer isso com o R via download.file) > > 2. extrair o .HTM do arquivo (posso fazer isso via unzip) > > 3. ler a tabela contida no HTM utilizando readHTMLTable (lib XML) > > > > A questão que tenho é a seguinte: > > > > Existe um jeito mais rápido que este? Eu já consegui ler .csv direto de > arquivo .zip utilizando o comando unz e read.table. Será que funciona assim > também em HTML? Se sim, qual seria a func? > > > > Leonard > > > > OS: ainda tenho que pensar como fazer webscrap desse htm, kkk > > > > > > > _______________________________________________ > 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. > > > > _______________________________________________ > 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. >

1 0

extraindo dados de arquivo img
by Hélio Gallo Rocha 21 Jan '16

21 Jan '16

Boa tarde Listeiros, Alguém poderia me dar uma mão nesta tarefa: O site : http://spirits.jrc.ec.europa.eu/ dispõem arquivos climatológicos com as extenções : img e hdr Quem duipõem um CRM ou indicaria um pacote que extraia os dados e suas coordenadas? Grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Gerador automático de listagens
by Nei 20 Jan '16

20 Jan '16

Oi Pessoal, Uma ferramenta que pode ser útil: http://www.fakenamegenerator.com/order.php Um Abraço, Nei

2 3

Erro função lme
by Victor Fioreze 20 Jan '16

20 Jan '16

Olá! Estou rodando um script com a função lme do pacote nlme e esta dando um erro.Após eu rodar o modelo usando lme aparece: Erro em MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 Revisei os dados e o script e esta tudo ok.O que pode estar acontecendo? Att, Victor Ionatan FiorezeMédico Veterinário e ZootecnistaMestrando em Nutrição de Ruminantes, PPGZ-UFPelBolsista Capes Integrante GERUMEN (53) 81032887lattes: http://lattes.cnpq.br/5026587958713372 Este e-mail foi enviado por um computador sem vírus e protegido pelo Avast. www.avast.com

2 1

Conversão de matriz em raster georeferenciado
by ASANTOS 19 Jan '16

19 Jan '16

Prezados Listeiros, Gostaria de realizar a conversão de uma matriz em raster georeferenciado com duas informações, a representação de plantas em uma matriz (B) e juntar as informações geográficas conhecendo-se o espaçamento entre plantas (d1), no entanto, gostaria de saber se conforme o meu CRM abaixo essa conversão esta correta, pois não estou bem certo de que realização de um cbind() entre os objetos B e d1 correspondam a real posição desses elementos. Alguém poderia me ajudar? ##Dados artificiais injuria<-rnorm(1550, mean=12) ##Crio a matriz com os dados artificiais B = matrix( data=injuria, nrow=25, ncol=62) # Inserindo coordenadas espaciais p0<-c(680234.482625,7803361.55128)# Ponto inicial c0 <- seq(p0[1],p0[1]+182,by=2) d0 <- seq(p0[2],p0[2]+50,by=3) d1 <- expand.grid(x=c0, y=d0) ## Coordenadas do grid espaçadas 3 x 2 d1<-d1[-(1:14),] ## Removo algumas corrdenadas para ficar com o mesmo número de dados plot(d1[,1],d1[,2]) length(d1[,1]) require(raster) FORM<-cbind(d1,B) colnames(FORM)<-c('x','y','area') r <- rasterFromXYZ(as.data.frame(FORM)[, c("x", "y", "area")]) extent(r) <- c(680234.5,680416.5,7803362,7803410) projection(r) <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") plot(r) # 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 4

Múltiplas funções em uma única função
by ASANTOS 15 Jan '16

15 Jan '16

Prezados Listeiros, Gostaria de criar uma função com múltiplas funções em seu interior, no meu caso gostaria de criar uma função denominada MYsummary, com três funções no interior dela, sendo que à partir de banco de dados artificial: ##Banco de dados artificial who<-rep(c("JOSE","CARLOS"),40) dates<-rep(as.Date(c("2007-06-22", "2004-02-13")),40) var1<-rnorm(80) tableFE<-data.frame(who, dates,var1) head(tableFE) Tenho uma primeira função que identifica o objeto who, ou seja, quem são as pessoas: ## Função Identificadora da pessoa people<-function(x, db = tableFE) { x<-NULL results <- x x<-unique(db[,1]) results <- x results } PER<-people(db = tableFE) #--# print(PER) A segunda função me fornece o período da coleta de dados: ## Função do período timeFE<-function(x, db = tableFE) { x<-NULL minT<-min(db[,2]) maxT<-max(db[,2]) results <- x x<-cbind(minT,maxT) colnames(x)<-c("start","end") results <- x results } TIM<-timeFE(db = tableFE) #--# print(TIM) A terceira função calcula a média: MD<-mean(tableFE[,3]) print(MD) Porém quando tento juntar as três funções na minha função MYsummary(), para que as três funções sejam executadas concomitantemente, exibindo os três prints ao mesmo tempo não funciona, sendo: MYsummary<-function(x, db = tableFE) { c( ## Função Identificadora da pessoa people<-function(x, db = tableFE) { x<-NULL results <- x x<-unique(db[,1]) results <- x results print(results) } , ## Função do período timeFE<-function(x, db = tableFE) { x<-NULL minT<-min(db[,2]) maxT<-max(db[,2]) results <- x x<-cbind(minT,maxT) colnames(x)<-c("start","end") results <- x results print(results) } , MD=mean(tableFE[,3]) ) } MYsummary(tableFE)## Não funciona Alguém poderia me dar uma luz? Muito 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

1 1

Classificação de pontos em cada polígono em um shapefile
by ASANTOS 14 Jan '16

14 Jan '16

Caros Listeiros Eu criei um shapefile POLY e gostaria de lançar alguns pontos aleatórios (objeto pts) e depois classifica-los com a função over() para saber em qual polígono esta contido cada ponto, no entanto, todas as modificações que eu faço geram em incompatibilidade na avaliação de argumentos entre ppp e cordinates e gostaria que alguém me desse uma dica de como otimizar essa operação sem ficar fazendo inúmeras transformações de objetos, segue CRM: #Pacotes require(shapefiles) require(rgdal) require(maptools) require(spatstat) #------------------------------------------------------------------------------- #Criação de 4 polígonos sr1=Polygons(list(Polygon(cbind(c(180114, 180553, 181127, 181477, 181294, 181007, 180409, 180162, 180114), c(332349, 332057, 332342, 333250, 333558, 333676, 332618, 332413, 332349)))),'1') sr2=Polygons(list(Polygon(cbind(c(180042, 180545, 180553, 180314, 179955, 179142, 179437, 179524, 179979, 180042), c(332373, 332026, 331426, 330889, 330683, 331133, 331623, 332152, 332357, 332373)))),'2') sr3=Polygons(list(Polygon(cbind(c(179110, 179907, 180433, 180712, 180752, 180329, 179875, 179668, 179572, 179269, 178879, 178600, 178544, 179046, 179110), c(331086, 330620, 330494, 330265, 330075, 330233, 330336, 330004, 329783, 329665, 329720, 329933, 330478, 331062, 331086)))),'3') sr4=Polygons(list(Polygon(cbind(c(180304, 180403,179632,179420,180304), c(332791, 333204, 333635, 333058, 332791)))),'4') #Converte em polígonos espaciais sr=SpatialPolygons(list(sr1,sr2,sr3,sr4)) srdf=SpatialPolygonsDataFrame(sr, data.frame(row.names=c('1','2','3','4'), PIDS=1:4)) #Cria o shapefile writeOGR(srdf, getwd(), 'POLY', 'ESRI Shapefile') #Leitura do shapefile contorno_line_X <- readShapeSpatial("POLY.shp") plot(contorno_line_X) #Seleção de pontos contorno_line_spat <- as(contorno_line_X,"owin") pts <- rpoispp(lambda=0.000012, win=contorno_line_spat) points(pts, col="blue") #Classificação dos pontos em cada área pts$PIDS <- over(pts, contorno_line_X)$PIDS pontos<-as.data.frame(pts) pontos #------------------------------------------------------------------------------- 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 6

Large Merge
by Wagner Tassinari 14 Jan '16

14 Jan '16

Olá pessoal, tudo bem ? Estou com um problema, estou querendo fazer o merge entre dois bancos de dados: > dim(bancolink) [1] 1630328 8 > dim(fichaAlink) [1] 2954930 6 Minha variável - chave primária está com o formato character. Daí qunado eu faço: > CadUnico_fichaA_merge = merge(bancolink, fichaAlink, by.x=c('CPF'), by.y=c('NUMERO_CPF')) Error in vecseq(f__, len__, if (allow.cartesian || notjoin || !anyDuplicated(f__, : Join results in 150596694 rows; more than 4585258 = nrow(x)+nrow(i). Check for duplicate key values in i each of which join to the same group in x over and over again. If that's ok, try by=.EACHI to run j for each group to avoid the large allocation. If you are sure you wish to proceed, rerun with allow.cartesian=TRUE. Otherwise, please search for this error message in the FAQ, Wiki, Stack Overflow and datatable-help for advice. Tentei fazer usando o formato data.table, daí ele dá um erro também: > CadfAmerge = bancolink[fichaAlink] Error in `[.data.table`(bancolink, fichaAlink) : When i is a data.table (or character vector), x must be keyed (i.e. sorted, and, marked as sorted) so data.table knows which columns to join to and take advantage of x being sorted. Call setkey(x,...) first, see ?setkey. Desde já agradeço a todos, -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

4 3

Formatando datas
by Wagner Tassinari 14 Jan '16

14 Jan '16

Olá pessoal, Estou com um probleminha. Em meu banco de dados, minha variável data se encontra desta forma: > data [1] 4042009 4042009 27072009 27072009 27072009 6082010 6082010 6082010 [9] 6082010 6082010> library(chron)> data2 = chron(dates=data, format=c('dd/mm/yy'))> data2 [1] 24/Aug/13036 24/Aug/13036 23/Aug/76090 23/Aug/76090 23/Aug/76090 24/Dec/18621 [7] 24/Dec/18621 24/Dec/18621 24/Dec/18621 24/Dec/18621 O certo seria: [1] 04/04/2009 04/04/2009 27/07/2009 27/07/2009 06/08/2010 ... Alguém poderia me dar uma dica ? Obrigado -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis) <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…> Este e-mail foi enviado por um computador sem vírus e protegido pelo Avast. www.avast.com <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…> <#DDB4FAA8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

6 5

Seleção de dígitos dos valores de um vetor
by ASANTOS 13 Jan '16

13 Jan '16

Prezados Listeiros, Tenho um vetor de dados db<-c(201002,486003,486016,486014,486004,486020,486002,486023,486027) e gostaria de extrair somente os três últimos dígitos de cada valor. Isso de uma maneira que eu tenha um novo vetor de dados sendo newdb: 002 003 016 014 004 020 002 023 027 Para isto pergunto, se alguém conhece alguma função que faça isso? Muito 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

3 3

Critério de Seleção de Modelos
by Andressa Siroky 10 Jan '16

10 Jan '16

Boa tarde! Quais critérios de seleção são adequados quando os modelos não são encaixados? Vi que AIC e BIC não são. Obg Andressa Siroky

5 4

Ler um HTM
by Leonard de Assis 07 Jan '16

07 Jan '16

Boa noite a todos Estou com um problema idiota em mãos. Preciso ler um arquivo .HTM dentro de um arquivo .zip. Até agora, o jeito que achei mais produtivo foi: 1. baixar o zip (posso fazer isso com o R via download.file) 2. extrair o .HTM do arquivo (posso fazer isso via unzip) 3. ler a tabela contida no HTM utilizando readHTMLTable (lib XML) A questão que tenho é a seguinte: Existe um jeito mais rápido que este? Eu já consegui ler .csv direto de arquivo .zip utilizando o comando unz e read.table. Será que funciona assim também em HTML? Se sim, qual seria a func? Leonard OS: ainda tenho que pensar como fazer webscrap desse htm, kkk

2 2

Leitura de tabela em formato HTML no dropbox
by ASANTOS 29 Dez '15

29 Dez '15

Prezados Listeiros, Estou tentando ler uma tabela em HTML que hospedei no dropbox, porém a função readHTMLTable() parece não funcionar se o arquivo esta no dropbox, alguém poderia me ajudar? Segue CRM: ##Leitura de uma tabela em HTML usando httr e XML ##Pacotes require(httr) require(XML) ## Abrindo o arquivo HTML com a tabela hospedado no dropbox FILE <- GET(url="https://www.dropbox.com/s/mb316ghr4irxipr/TALHOES_AGENTES.htm?dl=0") ## Lendo a tabela tables <- getNodeSet(htmlParse(FILE), "//table") FE_tab <- readHTMLTable(tables[2], header = c("empresa","desc_projeto","desc_regiao", "cadastrador_por","cod_talhao","descricao", "formiga_area","qtd_destruido","latitude", "longitude","data_cadastro"), colClasses = c("character","character","character", "character","character","character", "character","character","character", "character","character"), trim = TRUE, stringsAsFactors = FALSE ) head(FE_tab) # 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

1 0

Como evitar recursão infinita?
by Cleber N.Borges 21 Dez '15

21 Dez '15

olá a todos, Gostaria de saber se há alguma maneira de se evitar recursão quando houver um aninhamento de funções? Exemplo: f1 <- function() { "faz alguma coisa" e; chama f2() } f2 <- function() { "faz alguma coisa" e; chama f1() } Gostaria de, alguma forma, quando f1 (ou 2) fosse chamada dentro de f2 (ou 1), a função não entrasse mais na recursão! Agradeço pela atenção. cleber --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

2 1

Loop para leitura de picos
by José Lucas Safanelli 20 Dez '15

20 Dez '15

Boa tarde, estou com o seguinte caso: DADO: amostra = c(49.936157, 48.133488, 45.775238, 43.229397, 40.021607, 35.927558, 31.365912, 26.849848, 23.233177, 20.770278, 18.572067, 15.590341, 11.936978, 8.478520, 5.572115, 3.272254, 1.642041, 0.999997, 1.473886, 2.616919, 4.169691, 6.464019, 9.790810, 13.820101, 17.978423, 21.674281, 24.065947, 24.875069, 24.590811, 23.620608, 22.240339, 20.676521, 19.088229, 17.689105, 16.720912, 16.358739, 16.792281, 18.321545, 21.053896, 24.777756, 29.121638, 33.551339, 37.628649, 41.272700, 44.405581, 46.899820, 48.602420, 49.253714, 48.816692, 47.380737, 44.920718) wavelength = seq(400, 450, by=1) plot(wavelenght, amostra, xlab = "Comp. Onda (nm)", ylab = 'Valor', type = "l") axis(side=1, at=seq(400, 450, by=5)) PROBLEMA: Estou tentando ler os valores mínimos e máximos de cada pico para ter o valor da amplitude. Consigo fazer isso delimitando o intervalo de comp. de ondas da posição do pico, mas é muito laborioso EX: max(amostra[,1:20]) e min(amostraamostra[,1:20]). Amplitude seria a diferença deles. Para várias amostras tenho utilizado a função apply para dentro da linha, PORÉM decorrente da heterogeneidade das amostras, as posições variam causando erro de leitura e modelagem. Outro fato é que existem inúmeros picos dentro de uma faixa de interesse, e delimitando visualmente torna-se inviável. Para isso, estou tentando criar uma rotina que leia todos os picos dentro da faixa de interesse, para que possa utiliza-los diferentemente nas modelagens. Tentei o código abaixo, mas me retorna o erro que o argumento é de comprimento zero! for(i in 1:length(amostra)){ if((amostra[i+1] > amostra[i])&(amostra[i-1]>amostra[i])) print(amostra[i]) } (Error in if ((amostra[i + 1] > amostra[i]) & (amostra[i - 1] > amostra[i])) print(amostra[i]) : argument is of length zero) Um resultado que observei foi de do seguinte código, somente com uma condição colocada nele: for(i in 1:length(amostra)){ if(amostra[i+1] > amostra[i]) print(amostra[i]) } Nesse código, consegue-se ler os valores crescentes dos picos, desconsiderando os valores do objeto quando os valores são decrescentes. A partir dele que tentei colocar duas condições depois do IF para que ele lê-se somente os valores mínimos, mas sem sucesso. Portanto, peço ajuda para confecção de um código que faça a leitura dos min e max, ou amplitude, de todos os picos da amostra, retornando em um objeto/dataframe para modelagem. Grato pela atenção! José Lucas.

2 1

Re: [R-br] Digest R-br, volume 60, assunto 14
by José Lucas Safanelli 20 Dez '15

20 Dez '15

Pessoal, consegui resolver aqui com uma pesquisa mais aprofundada. Segue o código de leitura dos picos e vales para a amostra utilizada: maximo = data.frame() for(i in 1:length(amostra)){ if( max(amostra[i+1], amostra[i-1]) < amostra[i]) maximo = rbind(maximo, data.frame(print(amostra[i]))) } #LEITURA DOS MINIMOS minimo = data.frame() for(i in 1:length(amostra)){ if( min(amostra[i+1], amostra[i-1]) > amostra[i]) minimo = rbind(minimo, data.frame(print(amostra[i]))) } Atenciosamente, José Lucas. Message: 1 > Date: Sat, 19 Dec 2015 19:04:52 -0200 > From: José Lucas Safanelli <zecojls(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Loop para leitura de picos > Message-ID: > < > CADni8G5B-N9M4LqpZTVkAOgoSkTO56aK0icdg-PjVf0Us8PcFQ(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Boa tarde, estou com o seguinte caso: > > DADO: > amostra = c(49.936157, 48.133488, 45.775238, 43.229397, 40.021607, > 35.927558, 31.365912, 26.849848, 23.233177, 20.770278, 18.572067, > 15.590341, 11.936978, 8.478520, 5.572115, 3.272254, 1.642041, 0.999997, > 1.473886, 2.616919, 4.169691, 6.464019, 9.790810, 13.820101, 17.978423, > 21.674281, 24.065947, 24.875069, 24.590811, 23.620608, 22.240339, > 20.676521, 19.088229, 17.689105, 16.720912, 16.358739, 16.792281, > 18.321545, 21.053896, 24.777756, 29.121638, 33.551339, 37.628649, > 41.272700, 44.405581, 46.899820, 48.602420, 49.253714, 48.816692, > 47.380737, 44.920718) > wavelength = seq(400, 450, by=1) > plot(wavelenght, amostra, xlab = "Comp. Onda (nm)", ylab = 'Valor', type = > "l") > axis(side=1, at=seq(400, 450, by=5)) > > PROBLEMA: > Estou tentando ler os valores mínimos e máximos de cada pico para ter o > valor da amplitude. > Consigo fazer isso delimitando o intervalo de comp. de ondas da posição do > pico, mas é muito laborioso > EX: max(amostra[,1:20]) e min(amostraamostra[,1:20]). Amplitude seria a > diferença deles. > > Para várias amostras tenho utilizado a função apply para dentro da linha, > PORÉM decorrente da heterogeneidade das amostras, as posições variam > causando erro de leitura e modelagem. > > Outro fato é que existem inúmeros picos dentro de uma faixa de interesse, e > delimitando visualmente torna-se inviável. Para isso, estou tentando criar > uma rotina que leia todos os picos dentro da faixa de interesse, para que > possa utiliza-los diferentemente nas modelagens. > > Tentei o código abaixo, mas me retorna o erro que o argumento é de > comprimento zero! > > for(i in 1:length(amostra)){ > if((amostra[i+1] > amostra[i])&(amostra[i-1]>amostra[i])) > print(amostra[i]) > } > > (Error in if ((amostra[i + 1] > amostra[i]) & (amostra[i - 1] > > amostra[i])) print(amostra[i]) : argument is of length zero) > > Um resultado que observei foi de do seguinte código, somente com uma > condição colocada nele: > > for(i in 1:length(amostra)){ > if(amostra[i+1] > amostra[i]) > print(amostra[i]) > } > > Nesse código, consegue-se ler os valores crescentes dos picos, > desconsiderando os valores do objeto quando os valores são decrescentes. A > partir dele que tentei colocar duas condições depois do IF para que ele > lê-se somente os valores mínimos, mas sem sucesso. > > Portanto, peço ajuda para confecção de um código que faça a leitura dos min > e max, ou amplitude, de todos os picos da amostra, retornando em um > objeto/dataframe para modelagem. > > Grato pela atenção! José Lucas. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151219/add2cc3b/atta… > > >

1 0

Erro installar r DEBIAN
by Fernando Antonio de souza 18 Dez '15

18 Dez '15

Caros colegas Uso um pc Debian jessie 64 bits. Meu sistema estava funcionando plenamente. após fazer a atualização do r para a recente versão 3.2.3 passei a receber a seguinte mensagem em meu sistema ao rodar o camando dpkg --configure -a no shell : Erros foram encontrados durante o processamento de: r-base-core r-cran-boot r-cran-cluster r-cran-foreign r-cran-mass r-cran-kernsmooth r-cran-lattice r-cran-nlme r-cran-matrix r-cran-mgcv r-cran-survival r-cran-rpart r-cran-class r-cran-nnet r-cran-spatial r-cran-codetools r-recommended r-base r-base-dev r-base-html E: Sub-process /usr/bin/dpkg returned an error code (1) Isto é apenas o final, pois a saída é muito grande. Mas basicamente parece ser uma questão de dependências. Pelo que entendi está ocorrendo um problema ao instalar o pacote r-base-core que impede a instalação do r-base e assim dos demais: f:~$ sudo dpkg --configure r-base-core Configurando r-base-core (3.2.3-1~jessiecran.0) ... dpkg: erro ao processar o pacote r-base-core (--configure): sub-processo script post-installation instalado retornou estado de saída de erro 10 A processar 'triggers' para menu (2.1.47) ... Erros foram encontrados durante o processamento de: r-base-core :~$ sudo dpkg --configure r-base dpkg: problemas com dependências impedem a configuração de r-base: r-base depende de r-base-core (>= 3.2.3-1~jessiecran.0); porém: Pacote r-base-core não está configurado ainda. r-base depende de r-recommended (= 3.2.3-1~jessiecran.0); porém: Pacote r-recommended não está configurado ainda. dpkg: erro ao processar o pacote r-base (--configure): problemas de dependência - deixando desconfigurado Erros foram encontrados durante o processamento de: r-base Eu ja fiz a remoção completa do R utilizando o synaptic e tentei reinstalar e sempre volta a mesma coisa. Após reinstalar e digitar R no shell o sofware abriu com as seguintes messagens Warning message: package "methods" in options("defaultPackages") was not found Error : .onLoad falhou em loadNamespace() para 'utils', detalhes: chamada: options(op.utils[toset]) erro: valor inválido para 'editor' Error : .onLoad falhou em loadNamespace() para 'utils', detalhes: chamada: options(op.utils[toset]) erro: valor inválido para 'editor' Além disso: Warning message: package ‘utils’ in options("defaultPackages") was not found Error : .onLoad falhou em loadNamespace() para 'utils', detalhes: chamada: options(op.utils[toset]) erro: valor inválido para 'editor' Além disso: Warning message: package ‘stats’ in options("defaultPackages") was not found Durante a inicialização - Warning message: package ‘methods’ in options("defaultPackages") was not found Alguém sabe como corrigir este problema. Eu ja procurei em vários foruns debian/ubuntu e ja tentei comandos como apt-get install -f apt-get clean apt-get autoremove etc O que eu faço? -- ======================================================================= Fernando Souza Zootecnista, DSc. Produção Animal e-mail:nandodesouza@gmail.com https://producaoanimalcomr.wordpress.com/ ========================================================================

2 5

Criar um spatial points e shapefile com coordenadas e seus atributos
by ASANTOS 18 Dez '15

18 Dez '15

Caros Listeiros, Estou tentando criar um spatial points e posteriormente um shapefile com quatro variáveis, longitude (long), latitude (lat), espécie(esp) e cobertura (cob). Quando faço um cbind() apenas das informações geográficas, tudo OK, mas quando tento inserir os atributos esp e cob, a função SpatialPoints() não funciona. Alguém com experiência poderia me dar uma luz? Segue CRM: #Pacotes require(shapefiles); require(sp) #Pontos e variáveis esp e cob long<-c(-54.8956314,-54.8986246,-54.8985713, -54.8987823,-54.898895,-54.8988127) lat<-c(-30.4011373,-30.4121832,-30.4121513, -30.4120925,-30.4122421,-30.4123977) esp<-c("a","a","b","a","b","b") cob<-c(74.3,56.2,96.3,66.3,21.6,55.2) final.pts<-cbind(long,lat,esp,cob) #Transforma em spatial points pts.sampling = SpatialPoints(final.pts, proj4string=CRS("+proj=longlat +datum=WGS84")) ANTcoor.df <- SpatialPointsDataFrame(pts.sampling, data.frame(id=1:length(pts.sampling))) #Salva em shapefile writeOGR(ANTcoor.df, getwd(), 'ANTcoorF', 'ESRI Shapefile') 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 1

[fora] Exemplo bem-humorado de regressão à média
by Leonardo Ferreira Fontenelle 17 Dez '15

17 Dez '15

Assim como nos anos anteriores, esta edição natalina do periódico médico britânico *The BMJ* trouxe alguns artigos originais inusitados e divertidos. Acredito que boa parte dos companheiros da lista de discussão estarão interessados neste artigo como uma ilustração do conceito de regressão à média: http://dx.doi.org/10.1136/bmj.h6304 Att, Leonardo Ferreira Fontenelle[1] Links: 1. http://lattes.cnpq.br/9234772336296638

1 0

Meetup R São Paulo Users Group - 2016
by Luis Benites Sánchez 17 Dez '15

17 Dez '15

Caros, Para o 2016 temos previstos os seguintes Meetup futuros de R para o 2016. ######################################################### *Sexta-feira, 4 de Março de 2016 14:00* - R e Julia uma comparação Parte I - Criação de pacotes - Parte I *Sexta-feira, 18 de Março de 2016 14:00* - Birnbaum-Saunders SN com aplicações em R - R e Julia uma comparação Parte II *Sexta-feira, 8 de Abril de 2016 14:00* - Técnicas de amostragem com R - Uma introdução ao dados de contagem com R *Sexta-feira, 6 de Maio de 2016 14:00* - ggplot - Ajustando modelos semi-paramétricos log-simétricos com R *Sexta-feira, 3 de Junho de 2016 14:00* - Visualizações em R *Sexta-feira, 22 de Julho de 2016 14:00* - Modelos de regressão com dados censurados SMSN - Introdução a amostragem com R *Sexta-feira, 29 de Julho de 2016 14:00* - Uma introdução a programação em paralelo em R - R Markdown - Parte 2 ######################################################### As datas e horários podem variar. As atualizações estarão disponíveis no site do Meetup. As palestras serão transmitidas ao vivo pelo YouTube. Para poder participar, pode confirmar sua participação gratuitamente no site do Meetup :) Para mais informações: Site: www.meetup.com/pt/useR-SP Canal de YouTube: www.youtube.com/channel/UCb9PdJYYaElE6Uqb1KxhE8w e-mail: meetuprsp(a)gmail.com Um abraço, Luis Benites Sánchez Rocío Paola Maehara Organizadores do São Paulo R Users Group

1 0

Re: [R-br] Matrix de contrastes em fatorial DBC com a função glht!
by Ivan Bezerra Allaman 17 Dez '15

17 Dez '15

Prezado Walmes! Como sempre, muito obrigado! A função apcMatrix é muitíssimo útil. Uma sugestão! Embora eu saiba que tu não gostes de letrinhas, talvez seria interessante criar alguma função semelhante a função "cld", já que a mesma só trabalha caso seja indicado no argumento "linfct" a função "mcp" com o nome do tratamento de interesse. Em casos mais gerais como foi neste exemplo, a função "cld" não trabalha. Abraços!

3 3

Re: [R-br] manipulação arquivos
by Rafael Tieppo 16 Dez '15

16 Dez '15

Paulo, MUITO obrigado! Tenho que estudar e aprender a usar as funções de manipulações! Salvou o dia! Abrax On Wednesday, December 16, 2015 4:03 PM, Paulo Abreu <abreups(a)gmail.com> wrote: Rafael, Porque você não tenta primeiro fazer um “pré-processamento” desse arquivo para só depois ler com o R?Por exemplo, se você usar o 'grep' com uma expressão regular para isolar somente as linhas com dados, você salva isso num outro arquivo e depois lê esse outro arquivo com o R.Salvei esse exemplo de dados que você passou num arquivo chamado teste.txt e executei o seguinte comando: grep "^ 198" < teste.txt O resultado são todas as linhas que começam com um espaço em branco e depois têm os caracteres 198 Assim: 1983 001 1 0.0 13.8 1983 002 2 0.0 13.8 1983 003 3 0.0 13.7 1984 001 1 2.2 13.8 1984 002 2 54.8 13.8 1984 003 3 4.8 13.7 Se você fizer grep "^ 198" < teste.txt > resultado.txt O arquivo resultado.txt vai ter só as linhas de dados.Coloca um cabeçalho à mão e aí fica mais fácil de ler com o R.Ajuste a expressão regular para ficar mais adequada aos dados que você tem (sei que você passou só um pedaço dos dados). Espero que ajude. Paulo Abreu Em 16 de dez de 2015, à(s) 16:14, Rafael Tieppo <rafaeltieppo(a)yahoo.com.br> escreveu: Prezados, Por favor, Alguém poderia me indicar algum pacote ou algum método para ler várias Tabelas em um único arquivo ".txt", quando há texto no meio do arquivo. Creio que o símbolo "@" permita iniciar a captura dos dados, mas não tenho ideia de como criar uma lógica. Obrigado EXEMPLO do ARQUIVO: *RUN 1 : POT_S1 @YEAR DOY DAS PRED DAYLD 1983 001 1 0.0 13.8 1983 002 2 0.0 13.8 1983 003 3 0.0 13.7 *RUN 2 : POT_S1 @YEAR DOY DAS PRED DAYLD 1984 001 1 2.2 13.8 1984 002 2 54.8 13.8 1984 003 3 4.8 13.7 Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://www2.unemat.br/rafaeltieppo blog: https://fuidebicicleta.wordpress.com"Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente". _______________________________________________ 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.

1 0

manipulação arquivos
by Rafael Tieppo 16 Dez '15

16 Dez '15

Prezados, Por favor, Alguém poderia me indicar algum pacote ou algum método para ler várias Tabelas em um único arquivo ".txt", quando há texto no meio do arquivo. Creio que o símbolo "@" permita iniciar a captura dos dados, mas não tenho ideia de como criar uma lógica. Obrigado EXEMPLO do ARQUIVO: *RUN 1 : POT_S1 @YEAR DOY DAS PRED DAYLD 1983 001 1 0.0 13.8 1983 002 2 0.0 13.8 1983 003 3 0.0 13.7 *RUN 2 : POT_S1 @YEAR DOY DAS PRED DAYLD 1984 001 1 2.2 13.8 1984 002 2 54.8 13.8 1984 003 3 4.8 13.7 Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://www2.unemat.br/rafaeltieppo blog: https://fuidebicicleta.wordpress.com"Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".

1 0

Matrix de contrastes em fatorial DBC com a função glht!
by Ivan Bezerra Allaman 15 Dez '15

15 Dez '15

Prezados, bom dia! Tenho os seguintes dados: gnexpo <- structure(list(celula = structure(c(7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 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, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 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, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("cd14il10", "cd14il12", "cd14TGFb", "cd14TNFa", "cd4il10", "cd4il17", "cd4il4", "cd4il6", "cd4infg", "cd8il10", "cd8il17", "cd8il4", "cd8il6", "cd8infg"), class = "factor"), tratamento = structure(c(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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 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, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("sangue", "sangueesp", "sangueesplei", "sanguelei"), class = "factor"), sujeito = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L), .Label = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16"), class = "factor"), ybox = c(NA, -2.69292076130798, -1.7308938232588, -2.71897413664256, -2.52195729242333, -2.07143460646478, -2.28616618805722, -2.23934403380362, -1.5816386015359, -1.45531244494296, NA, -1.78763628106266, -1.78763628106266, NA, -1.67799460050211, -2.05420193782142, NA, -2.69292076130798, -2.00495057061834, NA, -2.52195729242333, -2.39220210880234, -2.02098207859771, -2.21718134741238, -2.20638335201612, NA, -1.90186608973923, -1.96640375569663, NA, NA, -2.17503854559032, -2.13546702289896, NA, -2.4374515925038, -2.37786948402795, NA, -4.34782608695652, -2.66817235203095, NA, -1.93700472646838, -2.39220210880234, NA, NA, NA, NA, NA, -2.19576462068308, -2.13546702289896, NA, NA, -2.17503854559032, -2.42197056348201, -2.22816571976642, -2.37786948402795, -1.95893989862709, -2.71897413664256, -1.67799460050211, -1.7308938232588, NA, NA, NA, -1.84893147202209, -1.84893147202209, -2.20638335201612, NA, 2.29386873816499, 1.01505310461186, 1.59806601020187, 1.70359042610938, 0.629371021625038, 1.89488135689716, 1.34173827270132, 1.97933615752285, 1.32950237522194, 1.80954049266222, 0.544292233915625, 2.58628244712855, 2.60920818919467, 1.96791029953215, 1.37792682396749, NA, 2.42730191595763, 0.918037037053, 1.43659188396698, 1.52650689488233, 0.780499053243246, 1.95064066159234, 1.35791699022796, 2.0629046336914, 1.93612722331712, 1.57334871090689, 0.402847055658559, 2.14819950257032, 2.60090083364325, 1.65670196502143, 1.43659188396698, NA, 2.10868691208822, 1.21518877022432, 1.1617541378147, 1.3376694301362, 0.610168122899425, 1.96504299434836, 1.3376694301362, 2.02991623442824, 1.15267768767035, 0.91276233665236, 0.177480055692309, 2.21231887578234, 1.79072828872115, 1.62934834431764, 1.40555319415085, NA, 2.0629046336914, 1.17527566321877, 1.36995100720572, 1.6702237928893, 0.262246858924378, 1.59455599161021, 1.44043274603715, -1.4953647932855, 1.05902858248096, 0.875334134463397, 0.114818592199316, 2.3745242548625, 2.31314212920843, 1.88592320743286, 1.42114242322365, -1.4953647932855, -2.71897413664256, -1.99707061963371, -4.34782608695652, -4.34782608695652, -2.06276388770908, -1.83626720563971, -1.80550931322309, -2.40689528336595, -2.31097837032429, -2.1853184563932, -1.7308938232588, -2.06276388770908, NA, -2.66817235203095, -2.13546702289896, NA, -2.69292076130798, -2.64458778082735, -2.40689528336595, -4.34782608695652, -2.08021732374643, -2.66817235203095, -2.71897413664256, -1.67799460050211, NA, -1.78763628106266, -1.99707061963371, NA, NA, -2.19576462068308, -2.28616618805722, NA, -2.71897413664256, -2.64458778082735, -2.69292076130798, -4.34782608695652, -4.34782608695652, -1.41703514811409, -2.06276388770908, -1.38036185998819, NA, -1.34514600593702, NA, -2.08021732374643, -1.24705599641255, -2.1853184563932, -2.46973767079137, NA, -2.45336438654387, -4.34782608695652, -4.34782608695652, -2.80686212698278, -2.40689528336595, -2.6220484374455, -2.71897413664256, -1.7308938232588, -2.07143460646478, -1.7308938232588, NA, NA, -2.08911542185436, -1.84893147202209, -2.15495397463673, NA, -2.22816571976642, -1.7308938232588, -2.66817235203095, -2.52195729242333, -1.53739362224875, -4.34782608695652, -2.01292026207036, -1.78763628106266, -2.52195729242333, -1.67799460050211, -1.7308938232588, -1.9739467612432, NA, -2.37786948402795, -1.34514600593702, NA, -1.95155330658925, -2.64458778082735, -1.78763628106266, -2.10727211374119, -1.94424216534924, -4.34782608695652, -2.69292076130798, -1.7308938232588, NA, -1.95155330658925, NA, NA, NA, NA, -1.7308938232588, NA, -2.08911542185436, -2.35020867465772, -2.06276388770908, -2.80686212698278, -2.21718134741238, -1.87492464519061, -2.08021732374643, -1.8618143278143, NA, NA, -2.22816571976642, NA, NA, -2.03739213457179, -1.83001441700635, NA, -0.675361229952281, -2.15495397463673, -1.86834047734623, -2.54052969919257, -2.21718134741238, -2.35020867465772, -2.08911542185436, -2.08021732374643, -2.10727211374119, -1.27859646831934, -1.38036185998819, -1.00235979040798, -1.00235979040798, -1.15837435544212, -2.15495397463673, NA, -4.34782608695652, -2.64458778082735, -4.34782608695652, -4.34782608695652, -4.34782608695652, -4.34782608695652, -2.4374515925038, -2.19576462068308, -2.50399662742379, NA, -2.2623152494604, -1.93700472646838, NA, -4.34782608695652, -2.25072434606231, NA, -4.34782608695652, -4.34782608695652, -2.69292076130798, -2.52195729242333, -2.39220210880234, -4.34782608695652, NA, -2.1853184563932, NA, NA, -2.37786948402795, NA, NA, -2.39220210880234, -2.71897413664256, NA, -4.34782608695652, -4.34782608695652, -2.69292076130798, -2.80686212698278, -2.66817235203095, -4.34782608695652, -4.34782608695652, -1.95893989862709, NA, NA, NA, -1.78763628106266, -1.83626720563971, -2.66817235203095, -2.27412592309772, NA, -2.40689528336595, -2.64458778082735, -2.71897413664256, -2.54052969919257, -2.19576462068308, -4.34782608695652, -4.34782608695652, -1.95893989862709, -1.95155330658925, NA, NA, NA, -4.34782608695652, -2.71897413664256, NA, NA, -0.846625250115866, -0.424146292887046, -1.05165266989135, -0.278371048130955, -1.24705599641255, -0.978598328138476, -0.978598328138476, -0.382587090747639, -0.623794941549055, NA, -0.410129619681842, -0.955384996473032, -1.5816386015359, -1.45531244494296, -0.711134555233659, NA, -1.4953647932855, -1.10354472792737, -0.558555816099048, -0.481958704693426, -0.846625250115866, -0.657903436683406, -0.542819523945636, -0.623794941549055, NA, -0.542819523945636, NA, -0.265954948927854, -0.910492044990005, -1.31126094846664, -0.846625250115866, NA, -2.32377374500676, -1.0267000794542, -0.888762381154637, -0.511985037199133, -0.290915281172533, -1.24705599641255, -1.38036185998819, -0.329348187317446, -0.182393356563792, -1.05165266989135, -0.193990047506904, -0.0932951673799442, -1.05165266989135, -1.18701626669498, -0.826177946574437, NA, -1.18701626669498, -0.867480330972702, -0.955384996473032, -0.558555816099048, -0.693101362150131, -0.806120545038559, -0.846625250115866, -0.481958704693426, -1.13056852939884, -1.07725416865755, -0.438332632355958, NA, -1.31126094846664, -1.05165266989135, -0.826177946574437, NA, -0.382587090747639, 0.0776462231901975, -0.316400590353713, NA, -0.729472182374671, -0.0825871589763347, NA, -0.125972397047986, -0.253664014288352, 0.0776462231901975, -0.424146292887046, NA, -0.316400590353713, NA, -0.0719684003695356, NA, -0.278371048130955, -0.159518997859581, 0.0198477924053798, 0.417635857120927, -0.607125627348296, 0.105622563120547, NA, 0.278588642396027, -0.675361229952281, 0.482353675529317, -0.496873898568359, NA, -0.438332632355958, -1.27859646831934, -0.410129619681842, NA, 0.150982355897241, -0.205696150863317, 0.212021636625242, 0.660835179222585, -0.826177946574437, -0.574513982772485, NA, -0.205696150863317, 0.410260691983402, -0.278371048130955, -0.303590729063761, NA, -0.369052612827179, -1.13056852939884, 0.114818592199316, NA, -0.452693434554865, -0.265954948927854, -0.217514069018959, 0.395394483699403, -0.193990047506904, 0.0198477924053798, NA, 0.0393981466797574, 0.15987295146349, 0.979995806626161, 0.721848977768072, 0.949319687828156, -0.496873898568359, -0.846625250115866, 0.310704279174891, -0.205696150863317, -1.67799460050211, -0.511985037199133, -1.41703514811409, -1.00235979040798, -0.452693434554865, -1.45531244494296, -0.806120545038559, -0.342436972562116, -0.382587090747639, NA, -1.24705599641255, -1.07725416865755, -1.21655441108309, -1.13056852939884, -1.00235979040798, -0.329348187317446, -1.21655441108309, -1.10354472792737, -1.7308938232588, -0.0825871589763347, -0.329348187317446, -0.955384996473032, -0.978598328138476, -0.748126323001209, -1.00235979040798, -0.303590729063761, -0.290915281172533, 0.245708807412923, -1.38036185998819, -1.24705599641255, -0.846625250115866, -0.786436366916452, -0.826177946574437, -0.786436366916452, -0.675361229952281, -0.978598328138476, -0.452693434554865, -1.0267000794542, -0.846625250115866, -0.496873898568359, -0.438332632355958, -0.527298140744417, -0.396278019492523, -0.0201558427538927, -1.13056852939884, -1.34514600593702, -1.27859646831934, -0.396278019492523, -0.657903436683406, -0.867480330972702, -1.13056852939884, -0.786436366916452, -0.693101362150131, -0.910492044990005, -1.5816386015359, -0.806120545038559, -0.978598328138476, NA, -0.452693434554865, -0.729472182374671, -1.07725416865755, -1.4953647932855, -1.07725416865755, -0.0614371998817261, -1.34514600593702, -1.07725416865755, -1.27859646831934, -1.10354472792737, -0.527298140744417, -1.27859646831934, -0.438332632355958, -0.867480330972702, -0.355670536085818, 0.114818592199316, -0.424146292887046, -0.452693434554865, -1.5816386015359, -1.67799460050211, -1.41703514811409, NA, -1.13056852939884, -1.13056852939884, -1.07725416865755, -1.45531244494296, -0.748126323001209, -1.07725416865755, -1.45531244494296, -0.711134555233659, -0.767109823860092, -0.657903436683406, -1.07725416865755, -0.826177946574437, -0.978598328138476, -1.24705599641255, -0.888762381154637, -0.0825871589763347, -1.15837435544212, -1.18701626669498, -0.607125627348296, -0.711134555233659, -0.290915281172533, -0.558555816099048, -0.867480330972702, -0.452693434554865, -0.452693434554865, -0.806120545038559, -1.21655441108309, -0.316400590353713, -0.888762381154637, -1.88156803086095, NA, NA, -0.826177946574437, -1.13056852939884, -0.496873898568359, -0.542819523945636, -0.623794941549055, -0.826177946574437, -1.53739362224875, -0.888762381154637, -0.438332632355958, -0.623794941549055, -0.558555816099048, -0.159518997859581, -0.978598328138476, -0.888762381154637, -0.303590729063761, -0.675361229952281, -0.693101362150131, -1.0267000794542, -1.07725416865755, -1.15837435544212, -1.45531244494296, -0.767109823860092, -0.303590729063761, -0.607125627348296, -0.693101362150131, NA, -0.438332632355958, -0.382587090747639, -1.31126094846664, -1.7308938232588, -1.45531244494296, -0.438332632355958, -1.34514600593702, -1.78763628106266, -0.846625250115866, -0.786436366916452, -1.05165266989135, -1.05165266989135, -1.24705599641255, -0.481958704693426, -0.511985037199133, 0.0490649462655013, -0.675361229952281, NA, -0.955384996473032, -1.4953647932855, -0.342436972562116, -0.382587090747639, -0.767109823860092, -1.34514600593702, -0.253664014288352, -0.806120545038559, -1.24705599641255, -1.41703514811409, -1.07725416865755, -0.511985037199133, NA, -0.590701349454554, -0.0201558427538927, -0.104094168973296, -1.05165266989135, -1.45531244494296, -1.62838825505761, -0.0406309528360936, -0.329348187317446, -1.45531244494296, -1.07725416865755, -1.00235979040798, -1.10354472792737, -0.748126323001209, -1.27859646831934, -0.511985037199133, 0.150982355897241, -0.382587090747639, -0.396278019492523, -1.31126094846664, -1.34514600593702, -1.15837435544212, -0.729472182374671, NA, 1.7167663207627, 2.43634147030751, 0.245708807412923, 1.15267768767035, 0.365185236456101, 0.517089544114438, 0.73376127136365, 1.59806601020187, 0.679400605803935, 0.791963538514193, 1.72004536647892, 1.79702046385978, 0.424973006905875, 0.721848977768072, 1.92150068931436, -0.640717861878449, -0.104094168973296, NA, NA, NA, -0.170903755458391, -1.24705599641255, NA, NA, 0.30274478985439, 0.00996172558875369, -0.137055389065604, 0.349834251678161, -0.253664014288352, 0.544292233915625, 0.286687590134947, NA, -0.767109823860092, 0.133022612878126, -0.424146292887046, -0.170903755458391, -0.0406309528360936, 1.34173827270132, 0.342095353591695, 0.523937993488919, 1.34173827270132, 0.133022612878126, -0.0509919154844358, -0.278371048130955, NA, 1.97648616730155, 1.88592320743286, 1.25844070072036, 0.902159500232335, 0.0963625187904701, 1.35388670327868, 0.564366845252007, -1.24705599641255, 0.757309220404252, 1.24558240333043, 0.853536835529698, -1.18701626669498, 0.286687590134947, 0.150982355897241, 0.446760797795615, -0.253664014288352, -0.826177946574437, 1.70689340655833, 0.387902502823721, 2.58628244712855, 2.49194504353092, 1.57690093891415, 2.901014602786, 0.0776462231901975, 0.310704279174891, 2.72805294137269, -0.481958704693426, -0.104094168973296, -0.0201558427538927, 0.32648803558755, 0.0296595084509125, -0.867480330972702, 1.32540405417028, -0.217514069018959, NA, -0.496873898568359, 2.00479472433388, 1.84659785469478, 1.94194605742328, -0.205696150863317, -0.748126323001209, 0.310704279174891, 0.660835179222585, 0.475305203762996, -0.329348187317446, 0.654595194602954, -0.170903755458391, -0.410129619681842, 1.81576868042314, 0.836985867234145, 0.0296595084509125, -0.511985037199133, -0.396278019492523, NA, 0.168705227406102, -0.342436972562116, -0.438332632355958, -0.452693434554865, NA, -1.10354472792737, -0.910492044990005, 0.446760797795615, 0.0198477924053798, NA, NA, NA, 0.622997622202259, -0.137055389065604, 0.294739430916844, 0.262246858924378, -0.114985961040159, 0.105622563120547, 0.203468274579351, 0.0963625187904701, -0.590701349454554, -1.7308938232588, -0.481958704693426, NA, -0.290915281172533, 0.212021636625242, NA, NA, 2.23489664779158, 2.75336054283323, 2.93328701678478, 1.82507164076158, NA, 1.75251413217418, 2.12986067381282, 2.13511807776368, 2.58837711038031, 1.70359042610938, 1.2369549109739, 0.446760797795615, 0.791963538514193, 0.954473759568006, 1.33359081568175, 1.21956483123272, 2.49851667302163, 2.48535220761645, 2.48974979559189, 3.16770417289536, 2.89377723669538, 1.83740303366685, 0.831429364687842, 2.07918795894332, -0.0719684003695356, 2.21484038563627, 1.35388670327868, 0.745580804308699, 2.26705821548107, 1.43274246477897, 0.803343047103057, 1.20640191321435, 2.33697431241692, 1.34984685630819, 2.68051900094631, 2.00479472433388, 1.94774692527104, 1.19756827033542, 0.597226509378741, 2.47873800731025, 2.29386873816499, 0.809001384770533, 1.28806529107744, 0.73376127136365, 1.49704766870529, NA, 1.07819479726929, 1.4250177788916, 1.84353799187812, NA, 1.68364466583674, 1.32540405417028, 1.93321102883982, -0.396278019492523, 0.557705833690471, 1.83432789943533, 0.212021636625242, 0.709842232833092, 1.43274246477897, NA, 1.19756827033542, 1.39770525114091, 1.65670196502143, 1.76534118794863, NA, 0.453950279688335, 2.17667216064136, 1.78757408686753, NA, 1.63965502207382, 1.56264900151517, 1.17975830926739, 1.60156910558301, 1.22828256682037, 0.32648803558755, 0.703802899501406, 1.55906797859696, 1.47085674653963, 1.30891085553319, 0.721848977768072, 2.52243371930239, 2.23240078706643, 2.04097581379485, 2.56944838544469, 1.86182177231003, 0.123951614253518, -0.125972397047986, 1.31304925442411, 1.02984564759017, 1.1572221751824, 1.46708280045183, 0.721848977768072, 1.01009176362354, 1.96504299434836, 1.26695835965958, 0.523937993488919, NA, 0.853536835529698, 2.22238551779573, 2.06834779093892, 2.78791091530525, 0.349834251678161, -0.329348187317446, 1.84659785469478, 1.27120101787992, 0.503295005047169, 1.36193776900057, 0.482353675529317, 0.584170503622459, NA, 1.44426509458991, NA, NA, 1.70028139904404, 1.67694674039951, 1.5301547136733, 2.24980568722757, 0.380370630683735, 0.372798375628128, 2.19202902290734, 0.439535041984402, -0.0509919154844358, 0.622997622202259, NA, -0.657903436683406, 1.82507164076158, -0.0932951673799442, 1.07819479726929), celtrat = structure(c(7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 34L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 54L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 41L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 55L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 38L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 39L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 53L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 42L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 56L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 44L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 45L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L, 46L), .Label = c("cd14il10sangue", "cd14il12sangue", "cd14TGFbsangue", "cd14TNFasangue", "cd4il10sangue", "cd4il17sangue", "cd4il4sangue", "cd4il6sangue", "cd4infgsangue", "cd8il10sangue", "cd8il17sangue", "cd8il4sangue", "cd8il6sangue", "cd8infgsangue", "cd14il10sangueesp", "cd14il12sangueesp", "cd14TGFbsangueesp", "cd14TNFasangueesp", "cd4il10sangueesp", "cd4il17sangueesp", "cd4il4sangueesp", "cd4il6sangueesp", "cd4infgsangueesp", "cd8il10sangueesp", "cd8il17sangueesp", "cd8il4sangueesp", "cd8il6sangueesp", "cd8infgsangueesp", "cd14il10sangueesplei", "cd14il12sangueesplei", "cd14TGFbsangueesplei", "cd14TNFasangueesplei", "cd4il10sangueesplei", "cd4il17sangueesplei", "cd4il4sangueesplei", "cd4il6sangueesplei", "cd4infgsangueesplei", "cd8il10sangueesplei", "cd8il17sangueesplei", "cd8il4sangueesplei", "cd8il6sangueesplei", "cd8infgsangueesplei", "cd14il10sanguelei", "cd14il12sanguelei", "cd14TGFbsanguelei", "cd14TNFasanguelei", "cd4il10sanguelei", "cd4il17sanguelei", "cd4il4sanguelei", "cd4il6sanguelei", "cd4infgsanguelei", "cd8il10sanguelei", "cd8il17sanguelei", "cd8il4sanguelei", "cd8il6sanguelei", "cd8infgsanguelei"), class = "factor")), .Names = c("celula", "tratamento", "sujeito", "ybox", "celtrat"), class = "data.frame", row.names = 225:1120) modnexp1 <- lm(ybox ~ sujeito + celula*tratamento, gnexpo) tabanovn <- Anova(modnexp1,type='III') library(phia) inter <- testInteractions(modnexp1, fixed='celula',across='tratamento') library(multcomp) Tukey <- contrMat(with(gnexpo,tapply(ybox,tratamento,function(x)length(na.omit(x)))), "Tukey")#Acredito que a glht esteja equivocada ao montar os contrastes com dados desbalanceados de acordo com que tem na literatura. Mais ok!!! K1 <- cbind(Tukey, matrix(0, nrow = nrow(Tukey), ncol = ncol(Tukey))) rownames(K1) <- paste(levels(gnexpo$celula)[1], rownames(K1), sep = ":") K2 <- cbind(matrix(0, nrow = nrow(Tukey), ncol = ncol(Tukey)), Tukey) rownames(K2) <- paste(levels(gnexpo$celula)[2], rownames(K2), sep = ":") k1 <- rbind(K1,K2) K <- cbind(matrix(0,nrow(k1),ncol=16), k1, matrix(0,nrow(k1),ncol=47)) colnames(K) <- c(rep('nada',16), rep(colnames(Tukey),2), rep('nada',47)) # Minha dúvida é se de fato a matrix K está montada adequadamente! comp <- lm(ybox ~ sujeito + celtrat-1, gnexpo) rescompe <- glht(comp, linfct=K) rescomp1e <- summary(rescompe) Desde já grato!

2 1

Análise de séries temporais
by Gabriela Alves 15 Dez '15

15 Dez '15

Boa Noite, gostaria te esclarecer uma dúvida que muito impede o avanço do meu projeto, para uns parecerá trivial, mas confesso que para mim está bem complicado.Estou fazendo um estudo do comportamento climático do Estado de São Paulo para os anos considerados normais em comparação aos anos considerados anômalos. Isso me permite gerar uma séria temporal com dados diários de precipitação, umidade e temperatura para os verões de 1961 à 2015. Ao iniciar o processo pelo software R me deparei com a dificuldade falhas serem incorporadas (já que alguns dias, ou meses e até verões completos podem não ter sido aferidos), tendo em vista essa dificuldade fui orientada a utilizar a metodologia de Fourier - ARIMA pois assim seria possível ter uma série completa com dados e falhas para todos os verões previsto. Todavia estou com muita dificuldade em avançar no desenvolvimento dessa série temporal modelagem ARIMA, é algo novo para mim e sendo bem sincera não sei por onde começar. Fazendo buscas descobri que usando o filtro de Kalmam da função auto.arima, é possível a imputação dos dados. Peço gentilmente uma orientação para esclarecer essa minha dúvida. Desde já agradeço. Att. Gabriela Alves

2 2

Dúvida no output da função lme
by Karen Castillioni 14 Dez '15

14 Dez '15

Pessoal, rodei a função lme para testar o efeito de 4 tratamentos (Controle, Corte, H e H+Ct) sobre a cobertura de uma vegetação. Porém ao olhar o output, fiquei em dúvida sobre qual seria o valor de p que devo usar para indicar se houve ou não diferença significativa entre os tratamentos. Alguém sabe qual seria o valor que devo usar? O do "Intercept"? Segue o output: Fixed effects: Cobertura ~ Tratamento Value Std.Error DF t-value p-value (Intercept) 44.79167 3.306008 84 13.548567 0.0000 TratamentoCorte 6.04167 4.675401 8 1.292224 0.2324 TratamentoH -41.91667 4.675401 8 -8.965363 0.0000 TratamentoH+Ct -35.75000 4.675401 8 -7.646403 0.0001 Correlation: (Intr) TrtmnC TrtmnH TratamentoCorte -0.707 TratamentoH -0.707 0.500 TratamentoH+Ct -0.707 0.500 0.500 Obrigada! Karen

2 1