- R-br - DINF UFPR

Adicionar os desvios padrão num gráfico
by Wenceslau Teixeira 12 Dez '12

12 Dez '12

Caro Walnes, Muito obrigado pela atenção. Funcionou tudo aqui ...apenas não aparece no grafico no eixo x as datas. Se voce ou outro membro da lista puder me ajudar em mais isto Segue o script que você corrigiu x <- c("30/04/2010","13/05/2010","28/05/2010", "15/06/2010","30/06/2010","15/07/2010", "30/07/2010","15/08/2010","30/08/2010", "15/09/2010","30/09/2010","15/10/2010", "30/10/2010","15/11/2010","30/11/2010", "15/12/2010","30/12/2010","15/01/2011", "30/01/2011","15/02/2011","28/02/2011", "15/03/2011","30/03/2011","15/04/2011", "30/04/2011","15/05/2011","30/05/2011", "15/06/2011","30/06/2011","15/07/2011") y <- c(0.278, 0.068, 0.064, 0.028, 0.030, 0.024, NA, 0.025, 0.039, 0.077, 0.081, NA, 0.043, NA, 0.343, NA, 0.195, 0.090, 0.027, 0.157, 0.122, 0.069, 0.197, 0.107, 0.103, 0.154,0.164, 0.051,0.243,0.145) SDFE10 <- c(0.294,0.033,0.058,0.013,0.006,0.006,NA,0.018,0.031,0,0, NA,0.021,NA,0.435,NA,0.057,0.049,0.005,0.220,0.159,0.029, 0.194,0.114,0.109,0.116,0.187,0,0.266,0.121) # colocar num data.frame, eliminar NA, converter para data da <- data.frame(x=as.POSIXlt(x, format="%d/%m/%Y"), y=y, sd=SDFE10) da <- na.omit(da) da$x plot(y~x, data=da, ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})),xlab="Dia", pch=4, type="b", xaxt="n",ylim=c(-0.1,0.8)) with(da, arrows(x, y-sd, x, y+sd, code=3, angle=90,length=0.03)) #axis.POSIXct(side=1, at=da$x, format="%d/%m") axis.POSIXct(side=1, at=seq(min(da$x), max(da$x), by="month"), format="%d/%m") Cordialmente, Wenceslau Teixeira Embrapa Solos -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

1 0

Forma "correta" e/ou mais "elegante"...
by Mauricio Morgado Oliveira 12 Dez '12

12 Dez '12

Gostaria de alguma dica (sugestão) para escrever de forma "correta" e/ou mais "elegante". Segue abaixo minha antiga programação no sas e no R (tentei anteriomente dentro de um for (loop), mas fica muito demorado). Att, Mauricio M. Oliveira Programação no SAS: if TP_PRODUTO='RM' then pai=COMPRESS('RM'||Cod_Rebanho||racapai||ano||nomepai,' '); else if livropai^=' ' & registropa^=' ' & registropa^='0' & registropa^='-' then pai=COMPRESS(livropai||'.'||registropa,' '); else if nomepai=' ' & tatuagempa=' ' then do; pai=COMPRESS(Cod_Rebanho||raca||ano||"DESC",' '); TP_PRODUTO='RM'; end; else pai=COMPRESS(Cod_Rebanho||racapai||tatuagempa||nomepai,' '); Script no R: Base[,"PAI"][Base[,"TP_PRODUTO"] == "RM"] <- paste("RM",Base[,"Cod_Rebanho"][Base[,"TP_PRODUTO"] == "RM"], Base[,"RACAPAI"] [Base[,"TP_PRODUTO"] == "RM"], Base[,"anoN"][Base[,"TP_PRODUTO"] == "RM"], Base[,"NOMEPAI"][Base[,"TP_PRODUTO"] == "RM"], sep="") Base[,"PAI"][Base[,"TP_PRODUTO"] != "RM" & Base[,"LIVROPAI"]!="" & Base[,"REGISTROPA"]!="" & Base[,"REGISTROPA"]!="0" & Base[,"REGISTRO"]!= "-"] <- paste(Base[,"LIVROPAI"][Base[,"TP_PRODUTO"] != "RM" & Base[,"LIVROPAI"]!="" & Base[,"REGISTROPA"]!="" & Base[,"REGISTROPA"]!="0" & Base[,"REGISTRO"]!="-"], ".", Base[,"REGISTROPA"][Base[,"TP_PRODUTO"] != "RM" & Base[,"LIVROPAI"]!="" & Base[,"REGISTROPA"]!="" & Base[,"REGISTROPA"]!="0" & Base[,"REGISTRO"]!="-"], sep="") Base[,"PAI"][Base[,"TP_PRODUTO"] != "RM" & Base[,"NOMEPAI"]=="" & Base[,"TATUAGEMPA"]==""] <- paste(Base[,"Cod_Rebanho"] [Base[,"TP_PRODUTO"] != "RM" & Base[,"NOMEPAI"]=="" & Base[,"TATUAGEMPA"]==""], Base[,"RACA"][Base[,"TP_PRODUTO"] != "RM" & Base[,"NOMEPAI"]=="" & Base[,"TATUAGEMPA"]==""],"DESC", sep="")

2 2

Adicionar os desvios padrão num gráfico
by Wenceslau 12 Dez '12

12 Dez '12

Caros membros da lista, Estou fazendo um gráfico simples (datas vs teores de ferro). Eu queria adicionar os desvios padrão das avaliações dos teores em cada data. Que foram calculados separadamente. Mas não estou conseguindo. Talvez alguém tenha experiência com estes gráficos e possa me ajudar. Segue o script ########################### # SCRIPT PARA GRAFICOS DE ELEMENTOS NO SOLO #DATAS DE COLETAS x<-1:30 x<-factor x<- c("30/04/2010","13/05/2010","28/05/2010","15/06/2010","30/06/2010", "15/07/2010", "30/07/2010" ,"15/08/2010", "30/08/2010","15/09/2010","30/09/2010","15/10/2010","30/10/2010","15/11/2010 ","30/11/2010","15/12/2010", "30/12/2010", "15/01/2011", "30/01/2011", "15/02/2011", "28/02/2011" ,"15/03/2011","30/03/2011","15/04/2011","30/04/2011","15/05/2011","30/05/201 1","15/06/2011","30/06/2011", "15/07/2011") #DADOS DE FE AVALIADOS NO LABORATORIO y<-1:30 y <- c(0.278, 0.068, 0.064, 0.028, 0.030, 0.024, NA, 0.025, 0.039, 0.077, 0.081, NA, 0.043, NA, 0.343, NA, 0.195, 0.090, 0.027, 0.157, 0.122, 0.069, 0.197, 0.107, 0.103, 0.154,0.164, 0.051,0.243,0.145) # DESVIOS PADRAO DOS TEORES DE FE SDFE10 <- c(0.294, 0.033, 0.058,0.013,0.006,0.006,NA,0.018,0.031,0,0,NA,0.021,NA,0.435,NA,0.057,0.049, 0.005,0.220,0.159,0.029,0.194,0.114,0.109,0.116,0.187,0,0.266,0.121) #PLOTAR O GRAFICO plot(y, ylim=c(-0.1,0.7), axes=F, ylab="Concentração de Ferro na solução do solo (mg L-1)", xlab="Dia", pch=4, type="b") #Adicionar os desvios padrão (SDFE10) arrows(y,y-SDFE10,y, y+SDFE10, code = 3, angle = 90,length = 0.03) Wenceslau Geraldes Teixeira Embrapa Solos (www.cnps.embrapa.br) Rua Jardim Botânico, 1024 CEP 22460-000 Rio de Janeiro - RJ - Brasil Tel: 55 21 21794583 Skype: wgt007 E-mail: wenceslau(a)cnps.embrapa.br / wenceslau.teixeira(a)pq.cnpq.br CV: http://lattes.cnpq.br/4136381865367906 ________________________________ II Encontro Brasileiro de Física do Solo 2sd Brazilian Soil Physics Meeting Rio de Janeiro - 06 a 10 de Maio de 2013 http://www.bspm2013.org/

2 1

Imputação múltipla: como fazer operações com os dados imputados.
by Lucas Petri Damiani 12 Dez '12

12 Dez '12

Boa noite, sou leitor assíduo da lista, embora nunca tenha participado das discussões. Porém, hoje vou me arriscar e passar um problema pra vcs. Espero que possam me ajudar. Estou realizando uma análise num banco de dados de cerca de 10 mil pacientes de Trauma. O objetivo da pesquisadora é testar a validade de um escore de gravidade (preditor de sobrevida) em pacientes que entram no pronto socorro e eventualmente promover alterações nesse escore. Porém, algumas medidas que ela deseja testar possuem cerca de 20% de missing. E a exclusão desses pacientes faz com que o novo escore predito fique um bocado diferente do original, o que me fez ir atrás de alternativas para tratar esses dados faltantes e assim me arrisquei nesse pacote *mice. *Em testes mais simples que fiz, fiquei bastante empolgado com o método. Enfim, tentei reproduzir abaixo o problema que encontrei. ##### Banco de dados de exemplo: dados <- matrix(c( 4, 6, 5, 34, 120, 5, 1, 0, NA, 1, 1, 3, 80, 6, 0, 1, 2, 4, 2, 24, 90, 5, 1, 1, 3, 5, 5, 5, 98, 5, 1, 1, NA, 6, 5, 40, 94, 4, 0, 0 , NA, 6, 5, 12, 87, 3, 0, 1, 1, 5, NA, 32, 70, 2, 0, 1, 1, 5, 1, 7, 60, 3, 1, 0, 2, NA, 2, 25, 80, 5, 1, 0, 4, 6, 5, 24, 101, 7, 1, 0, 1, 5, 1, 1, 120, 8, 1, 0, 3, 5, 3, 12, 200, 6, 0, 1, 3, 3, 1, 3, 40, 5, 1, 0, NA, 5, 1, 35, 60, 4, 0, 1, 1, 4, 1, 32, 87, 3, 1, 0, 4, 5, 1, 23, 30, 2, 1, 1, 2, 1, 1, 5, 46, 1, 0,1), byrow=T, ncol=8) dados <- as.data.frame(dados) colnames(dados) <- c("ECGL_AO", "ECGL_MRM", "ECGL_MRV", "FR", "PAS", "GRAU", "IOT", "SOBREVIDA") dados$IOT <- as.factor(dados$IOT) dados$ECGL <- apply(dados[, c("ECGL_AO", "ECGL_MRV", "ECGL_MRM")],1,sum) dados$ECGL_cut <- cut(dados$ECGL, breaks=c(3,4,6,9,13,16), right=FALSE) # Notem que a medida RTS é um escore derivado de outras medidas, que pretendo imputar. dados$RTS <- (0.9368*(dados$ECGL_cut=="[4,6)") + 2*0.9368*(dados$ECGL_cut=="[6,9)") + 3*0.9368*(dados$ECGL_cut=="[9,13)") + 4*0.9368*(dados$ECGL_cut=="[13,16)") + 0.7326*(dados$PAS>=1 & dados$PAS<=49) + 2*0.7326*(dados$PAS>=50 & dados$PAS<=75) + 3*0.7326*(dados$PAS>=76 & dados$PAS<=89) + 4*0.7326*(dados$PAS>=90) + 0.2908*(dados$FR>=1 & dados$FR<=5) + 2*0.2908*(dados$FR>=6 & dados$FR<=9) + 4*0.2908*(dados$FR>=10 & dados$FR<=29) + 3*0.2908*(dados$FR>=30)) #Imputando dados nas variáveis "ECLG": require(mice) PRED <- (1 - diag(1, ncol(dados[,1:8]))) PRED[8,1:8] <- PRED[1:8,8] <- 0 teste <- mice(dados[,1:8], method=c("pmm","pmm", "pmm", "", "", "", "", ""), seed=15, predictorMatrix = PRED) #O meu problema começa agora. Eu não sei como calcular medidas derivadas no objeto teste. E a fórmula abaixo retorna erro, já conferi várias vezes muito dos parênteses e etc, mas creio que o problema seja que a fórmula é muito grande entro da função I(x), pois quando eu reduzo o modelo, eu consigo resultados satisfatórios. Me corrijam se estiver errado, por favor. fit_error <- with(teste, glm(SOBREVIDA ~ I(0.9368*((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 4) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 5) + 2*0.9368*(((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 6) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 8)) + 3*0.9368*I(((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 9) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 12)) + 4*0.9368*I(((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 13) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 15)) + 0.7326*I(PAS>=1 & PAS<=49) + 2*0.7326*I(PAS>=50 & PAS<=75) + 3*0.7326*I(PAS>=76 & PAS<=89) + 4*0.7326*I(PAS>=90) + 0.2908*I(FR>=1 & FR<=5) + 2*0.2908*I(FR>=6 & FR<=9) + 4*0.2908*I(FR>=10 & FR<=29) + 3*0.2908*I(FR>=30)) + IOT, family="binomial")) #Enfim, gostaria de rodar algo como isso escrito abaixo, mas não sei como definir essa variável RTS dentro do arquivo teste. fit <- with(teste, glm(SOBREVIDA ~ RTS + IOT, family="binomial")) fit2 <- glm(SOBREVIDA ~ RTS + IOT, family="binomial", data=dados) pool(fit) # Uma saída que imagino é gerar os cinco bancos imputados através da função "complete(teste, m)" e ajustar a regressão. Mas dessa forma não posso utilizar a função *pool. *Teria que reprogramá-la. De qualquer forma, estou perdido. Qualquer ideia é bem vinda. Abraços, Lucas Petri Damiani

2 2

Re: [R-br] Trabalhar com model.matrix
by Walmes Zeviani 11 Dez '12

11 Dez '12

A model.matrix() te dá a matriz sob a restrição paramétrica considerada no ajuste, que por padrão é a contraste tratamento (anula 1 nível). Para ter a matriz sem restrições você vai ter que fazer a matriz para cada termo do seu modelo, removendo o intercepto e colando lado à lado, por exemplo mA <- model.matrix(~-1+A) mB <- model.matrix(~-1+B) mAB <- model.matrix(~-1+A:B) m <- cbind(A, B, AB) Submeta suas dúvidas para lista de discussão de R, mais pessoas podem contribuir, podem haver diferentes soluções, respostas em menos tempo, as mensagens ficam arquivadas para referência/consulta futura, as mensagens podem ser úteis para outros usuários. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Teste no R
by Márcia Varela da Silva 11 Dez '12

11 Dez '12

Não consigo + criar um novo tópico na R-br. Qual comando para aplicação do Teste de Bonferroni no R (para MANOVA - Comparações múltiplas multivariadas)? Desde já agradeço Att. Márcia Varela da SilvaMestranda em Produção Vegetal (Fitossanidade) - Universidade Federal do Espírito Santo(28)9964-3488(27)9812-1832

2 1

Grupo de estudo de Geoestatistica
by Hélio Gallo Rocha 11 Dez '12

11 Dez '12

Prezados, Não sei é prática desta lista, mas haveriam interessados em formar um grupo de estudo especifico em geoestatística? Hélio

15 17

Alterar source code
by Davi 10 Dez '12

10 Dez '12

Boa tarde! Estou precisando alterar o source code da função readline(), do base package. A questão é que preciso de mais de 256 caracteres e nos detalhes da ajuda da função (abaixo) parece que isto é possível de ser feito. Alguém sabe como? Obrigado! Details The prompt string will be truncated to a maximum allowed length, normally 256 chars (but can be changed in the source code). ___________________________________________________ Davi Butturi-Gomes Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP - Piracicaba Mestre em Biometria - IB/UNESP - Botucatu Ecólogo - IB/UNESP - Rio Claro Skype: davibg87 (19)8206.6565

2 1

Re: [R-br] Shapefile
by Thiago V. dos Santos 10 Dez '12

10 Dez '12

Taynãna, Estou respondendo com cópia para a lista, pois esse assunto tem despertado interesse de muitos usuários que podem se beneficiar com a solução da sua dúvida. Não é possível entender completamente o seu problema ser ter amostra dos dados. Se você tem um objeto da classe SpatialPolygon que precisa ser convertido para SpatialPolygonDataFrame para ser escrito em shapefile, você pode usar a função 'SpatialPolygonsDataFrame' do pacote sp. Você precisa conhecer a natureza dos seus dados para saber o que deve fazer. Vou postar abaixo um exemplo de conversão entre classes SpatialPolygon* que encontrei no manual do pacote maptools. Esse exemplo deve te ajudar: library(maptools) xx <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1], IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66")) SPxx <- as(xx, "SpatialPolygons") class(SPxx) IDs <- sapply(slot(SPxx, "polygons"), function(x) slot(x, "ID")) df <- data.frame(rep(0, length(IDs)), row.names=IDs) SPDFxx <- SpatialPolygonsDataFrame(SPxx, df) class(SPDFxx) summary(SPDFxx) tf <- tempfile() writePolyShape(SPDFxx, tf) getinfo.shape(tf) Espero que ajude, Thiago. ________________________________ From: Taynãna César Simões <tay_estatistica(a)yahoo.com.br> To: Thiago V. dos Santos <thi_veloso(a)yahoo.com.br> Sent: Monday, December 10, 2012 8:59 AM Subject: Re: [R-br] Shapefile Thiago e Paulo, Ainda não consegui transformar o objeto SpatialPolygons em um shapefile. Podem me dar um exemplo de como colocar os argumentos? O help não ajudou muito. Obrigada, Taynãna ________________________________ De: Thiago V. dos Santos <thi_veloso(a)yahoo.com.br> Para: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>; Taynãna César Simões <tay_estatistica(a)yahoo.com.br> Enviadas: Sexta-feira, 7 de Dezembro de 2012 15:21 Assunto: Re: [R-br] Shapefile Taynãna, Há funções apropriadas para escrever shapefiles nos pacotes maptools(writePolyShape) e rgdal (writeOGR), mas você provavelmente terá que converter seus dados em SpatialPolygonsDataFrame antes. Thiago. ________________________________ From: Taynãna César Simões <tay_estatistica(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Sent: Friday, December 7, 2012 2:45 PM Subject: [R-br] Shapefile Pessoal, Alguém sabe como exportar um objeto da classe SpatialPolygons como um arquivo shapefile? Att, Taynãna Fiocruz/RJ _______________________________________________ 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

Implementação Função
by FHRB Toledo 10 Dez '12

10 Dez '12

Senhores, Dia ou outro me vejo com a necessidade de aplicar o teste de Bartlett (função bartllet.test(), do pacote stats). Em um desses usos regulares notei certo "comportamento" anômalo na mesma. O uso que dou a ela era para testar a homogeneidade de variâncias entre modelos ajustados, ou seja, aplicava-a sobre uma lista de saídas da função lm(), da forma: bartlett.test(list(x = lm(...), z = lm(...), w = lm(...), y = lm(...))), De acordo com a sintaxe apresentada no help() da mesma ( http://stat.ethz.ch/R-manual/R-patched/library/stats/html/bartlett.test.html ). Para investigar o porque desse comportamento anômalo, além do help() indicado no link, busquei a citação do help, de autoria próprio Bartlett ( http://www.jstor.org/stable/96803) E a fonte da função disponível no tar.gz do R 2.15.2, que baixei do CRAN ( http://cran.r-project.org/src/base/R-2/R-2.15.2.tar.gz, acesso hoje 10 de dezembro de 2012). Transcrevo no final da mensagem um trecho da fonte, que se encontra em: ~/R-2.15.2/src/library/stats/R/bartlett.test.R. Pois bem, lendo a código fonte, o artigo do Bartlett e consultando o Livro Steel, Torrie and Dickey "Principles and Procedures of Statistics: A biometrial Approach", terceira edição, capítulo 19, página 480. Pude perceber que todo o fundamento do teste aplica-se a variâncias. Entretanto ao observar a implementação da função (linha destacada) observei que ao "capturar" sum(obj$residuals^2), a função estaria atendo-se a soma de quadrados do resíduo, e não a variância. Obviamente, uma soma de quadrados, dividido pelo seu respectivo grau de liberdade é uma quadrado médio e por sua vez uma variância, entretanto, isso não é realizado! Pelos testes que fiz, ao utilizar a função com diferentes modelos ajustados com o mesmo grau de liberdade dos resíduos, pela proporcionalidade das somas de quadrados o resultado do teste é satisfatório, entretanto, para diferentes modelos, com diferentes graus de liberdades o resultado do teste não confere com o obtido calculado a mão! Diante do exposto, gostaria de ouvir a opinião dos senhores quanto a esse fato, estou certo? Em caso positivo, qual seria um meio elegante de entrar em contato com o R CORE TEAM no sentido de apresentar esse fato? Agradeço qualquer explicação, sugestão, etc. att, FH bartlett.test.default <- function(x, g, ...) { LM <- FALSE if (is.list(x)) { if (length(x) < 2L) stop("'x' must be a list with at least 2 elements") DNAME <- deparse(substitute(x)) if (all(sapply(x, function(obj) inherits(obj, "lm")))) LM <- TRUE else x <- lapply(x, function(x) x <- x[is.finite(x)]) k <- length(x) } else { if (length(x) != length(g)) stop("'x' and 'g' must have the same length") DNAME <- paste(deparse(substitute(x)), "and", deparse(substitute(g))) OK <- complete.cases(x, g) x <- x[OK] g <- factor(g[OK]) k <- nlevels(g) if (k < 2) stop("all observations are in the same group") x <- split(x, g) } if (LM) { n <- sapply(x, function(obj) obj$df.resid) v <- sapply(x, function(obj) sum(obj$residuals^2)) ***** DESTAQUE ***** DESTAQUE ***** DESTAQUE ***** } else { n <- sapply(x, "length") - 1 if (any(n <= 0)) stop("there must be at least 2 observations in each group") v <- sapply(x, "var") } n.total <- sum(n) v.total <- sum(n * v) / n.total STATISTIC <- ((n.total * log(v.total) - sum(n * log(v))) / (1 + (sum(1 / n) - 1 / n.total) / (3 * (k - 1)))) PARAMETER <- k - 1 PVAL <- pchisq(STATISTIC, PARAMETER, lower.tail = FALSE) names(STATISTIC) <- "Bartlett's K-squared" names(PARAMETER) <- "df" RVAL <- list(statistic = STATISTIC, parameter = PARAMETER, p.value = PVAL, data.name = DNAME, method = "Bartlett test of homogeneity of variances") class(RVAL) <- "htest" return(RVAL) }

1 0

Re: [R-br] Falha em ANAVA
by Ivan Bezerra Allaman 10 Dez '12

10 Dez '12

Jovem, aconselho que leia o help da função glht novamente! Não tem sentido usar covariate_average=TRUE, pois não há covariável no modelo. Não tem sentido usar interaction_average=TRUE, pois não foi significativo. Não tem sentido utilizar procedimentos de comparação múltipla quando se tem apenas dois grupos. Ao fazer, summary(glht(analise, linfct=mcp(Manej="Tukey",interaction_average=FALSE))) Tem-se o mesmo resultado da Anova como era de se esperar, pois a estatística t ao quadrado é equivalente ao estatística F. Antes de partir para vias de fato (inferência), faça uma boa análise exploratória de dados. Existem três outliers na sua base. Avalie se eles são influentes ou não. Sabemos que valores extremos influenciam a média. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

Grupo de geoestatistica - script s100
by Hélio Gallo Rocha 10 Dez '12

10 Dez '12

Caros da lista. Anexei um script e um txt (recorte) Foi criado com objetivo de analisar uma situação de forma rápida, depois de "acostumado" com ele, em 10 min. tem os cinco mapas krigados É bastante útil quando se vai analisar uma tabela com muitos dados, como uma análise de solo. Tentei dar algumas explicações no inicio do script. A ideia é dar uma enxugada nele, tirar linhas desnecessárias. Está aberto para considerações, análises, correções, enriquecimento .... Acho interessante que qualquer mudança nele, poste e indique as mudanças e principalmente o porque para que possamos aprender quaisquer dúvidas, posso tentar esclarecer. Hélio

3 3

Falha em ANAVA
by Fernando Antonio de souza 09 Dez '12

09 Dez '12

Caros colegas, Estou realizando uma análise de variânica do peso vazio inicial de ovelhas em função da idade gestacional e manejo. Após rodar a análise a análise de variância tipo 3 não indicou haver diferenças para o nível do manejo (0,05) entretanto ao utilizar a função glht (multcomp) a análise mostrou haver diferenças entre os pesos. Em qual análise confiar? segue CMR library(multcomp) library(agricolae) library(doBy) library(phia) library(car) analise<-lm(PVINICIAL~Manej,data=GANHO1) Anova(analise,type="II") anova(analise) summary(glht(analise,linfct=mcp(Manej="Tukey",interaction_average=TRUE,covariate_average=TRUE))) dput(GANHO1)structure(list(Animal = c(25L, 38L, 42L, 50L, 53L, 18L, 22L, 34L, 41L, 44L, 48L, 2L, 12L, 14L, 39L, 46L, 4L, 5L, 21L, 26L, 49L, 3L, 32L, 37L, 43L, 47L, 8L, 15L, 30L, 40L, 1L, 17L, 23L, 28L, 7L, 9L, 19L, 35L), Gest = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("0", "100", "130", "140"), class = "factor"), Manej = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = "factor"), PVINICIAL = c(36.5, 37, 25.6, 30.3, 30.1, 27.3, 31, 25.1, 32.5, 29.6, 35, 40, 34.3, 33.8, 33.7, 31.4, 40, 32.3, 34.5, 29.5, 39.7, 38, 33.2, 34.5, 34.8, 31.3, 29.2, 32, 30.5, 27.5, 37, 32, 36.5, 40.6, 35.1, 29.8, 36.2, 28.6)), .Names = c("Animal", "Gest", "Manej", "PVINICIAL" ), class = "data.frame", row.names = c(NA, -38L))

4 6

Mensagem de Boas Festas postada no ano passado com código R
by Mauro Sznelwar 09 Dez '12

09 Dez '12

No ano passado alguém postou mensagem de Boas Festas para rodar no R. Queria saber onde poderia pegá-la novamente.

3 2

Grupo de estudo de Geoestatistica
by Jose Vilton Costa 08 Dez '12

08 Dez '12

Colegas,Achei a idéia ótima e eu também tenho interesse em participar. Se a idéia seguir a diante, a minha sugestão é que o nome do grupo mudasse para "Grupo de estudos em Estatística Espacial". Penso que este nome é mais amplo e cobriria de uma forma mais adequada as diversas áreas de estudo/aplicação dos dados de natureza espacial. Att,José Vilton > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 22, assunto 10 > To: r-br(a)listas.c3sl.ufpr.br > Date: Sat, 8 Dec 2012 12:00:02 -0200 > > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Grupo de estudo de Geoestatistica (Eder Comunello) > 2. Re: Grupo de estudo de Geoestatistica (Hélio Gallo Rocha) > 3. Re: Fwd: 2013 IBS Board Election Results > (Guilherme Moraes Ferraudo) > 4. Re: Grupo de estudo de Geoestatistica (Guilherme Moraes Ferraudo) > 5. Re: Grupo de estudo de Geoestatistica (Wagner Tassinari) > 6. Re: Grupo de estudo de Geoestatistica (Luis Iván Ortiz Valencia) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 8 Dec 2012 06:43:44 -0300 > From: Eder Comunello <ecomunel(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Grupo de estudo de Geoestatistica > Message-ID: > <CAOwgnzK0JjJcmCB61CRs_q8+Aw6AhA2fuLNGV73sbp8ttDzyMA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Hélio e colegas, bom dia! > > Tenho interesse no tema, mas também tenho dúvidas em relação à forma de > implementar este grupo. Pra mim é bastante cômodo que seja aqui na lista > mesmo, pois aproveito pra aprender sobre o uso do R em outras vertentes que > também me interessam. > > ================================================ > Éder Comunello > > Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) > Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] > > Researcher at Embrapa Western Region Agriculture > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > ================================================ > UTC-03:00 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121208/810d188c/atta…> > > ------------------------------ > > Message: 2 > Date: Sat, 8 Dec 2012 09:10:05 -0200 > From: Hélio Gallo Rocha <heliogallorocha(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Grupo de estudo de Geoestatistica > Message-ID: > <CA+ppPW5e7sGW-yfLVOoO1EwF7NKqgah5NJeHU=H8oPMaCiBJJw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Mauro, e colegas. > Não sabia desta lista, mas eu particularmente tenho dificuldade co m > a língua inglesa, > > Como sou professor de topografia e fazendo mestrado penso da seguinte forma > > > Objetivo: > Ajudar na resolução de problemas de possíveis interessados como > professores, agrônomos, ambientalistas. > > Metas, pensando em AP > Preparar mapas precisos > Calcular dosagens através dos mapas > Trabalhar co shape files > etc > > Para ajudar na organização, começaríamos em cima de um script (dados=S100), > tipo passo a passo, assim até quem nunca trabalhou com geoest conseguiria > segui > > Assim construiríamos em equipe um script adaptável a qualquer situação, com > um moderador e se possível que seja o que tem mais conhecimento do > assunto, > Quem gostaria de ser o moderador? > > > Tenho um script quase pronto (que ainda está longe de acabar...) > > > Disponibilizo este script inicial, em cima dos dados S100, fácil de adaptar > , para auxiliar a > Estudo da determinação da hipótese da Envelopado > Estudo de distancia máxima: auxilia a escolha da máxima > Cinco krigagens (asentimento e 4 validações) na mesma página > Relatório com AIC (FALTA VALIDAR) > Exportado para XLS > > > Hélio > > > > Em 8 de dezembro de 2012 07:43, Eder Comunello <ecomunel(a)gmail.com>escreveu: > > > Hélio e colegas, bom dia! > > > > Tenho interesse no tema, mas também tenho dúvidas em relação à forma de > > implementar este grupo. Pra mim é bastante cômodo que seja aqui na lista > > mesmo, pois aproveito pra aprender sobre o uso do R em outras vertentes que > > também me interessam. > > > > ================================================ > > Éder Comunello > > > > Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) > > Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] > > > > Researcher at Embrapa Western Region Agriculture > > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > > ================================================ > > UTC-03:00 > > > > > > _______________________________________________ > > 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/20121208/d8f4206e/atta…> > > ------------------------------ > > Message: 3 > Date: Sat, 8 Dec 2012 10:08:47 -0200 > From: Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Fwd: 2013 IBS Board Election Results > Message-ID: > <CA++u5OhL=A7WiF_2PZiKCC2QETy7Xyfu4RY7-NKacwWVFM1R3w(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Parabens Prof. Paulo Justiniano por essa merecida conquista. > > Abs > > 2012/12/8 Daniel Marcelino <dmsilva.br(a)gmail.com> > > > Parabéns Paulo Justiniano. > > > > Daniel > > > > On Fri, Dec 7, 2012 at 6:01 PM, Leandro Marino < > > leandromarino(a)leandromarino.com.br> wrote: > > > >> Gostaria de parabenizar um de nossos membros o Paulo Justiniano por esta > >> conquista que vem junto com muito trabalho. > >> > >> Paulo, parabéns. > >> > >> Abraços, > >> Leandro Marino > >> > >> Enviado de um dispositivo móvel. > >> Sent from a mobile device. > >> ---------- Forwarded message ---------- > >> From: <ibs(a)tibs.org> > >> Date: 7 Dec 2012 18:17 > >> Subject: 2013 IBS Board Election Results > >> To: <leandromarino(a)leandromarino.com.br> > >> Cc: > >> > >> Click here<http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856870&m=2417987&u=Bostr…>for web version. > >> > >> > >> 2013 IBS EXECUTIVE BOARD ELECTION RESULTS > >> Dear IBS Members, > >> > >> Thank you for your outstanding participation in the recent Society-wide > >> election for directors of the new IBS Executive Board. More than 25% of the > >> eligible voters participated in the election and 999 ballots were cast! > >> This is outstanding! > >> > >> We are pleased to announce that the following members have been elected > >> for either a term of two or four years beginning 1 January 2013. > >> > >> Wesley Johnson (2)WNAR James Carpenter (4)BIR José Pinheiro (4)ENAR Andreas > >> Ziegler (4)DR Sharon-Lise Normand (4)ENAR Paulo J. Ribeiro (4)RBras Alan > >> Welsh (4)AR Kerri Mengersen (2)AR Michael Schemper (2)RoeS Peter Njuho > >> (2)GSAf Vicente Núñez-Antón (2)Resp Lisa Lavange (2)ENAR > >> > >> The remaining members of the Executive Board include President Clarice > >> Demétrio (RBras), President-Elect John Hinde (BIR) and Secretary-Treasurer > >> Linda Young (ENAR). > >> > >> Please join me in congratulating all of our newly elected Executive Board > >> members. If you have ideas or wish to contribute to IBS as a volunteer in > >> some manner, please contact a member of the Executive Board and let them > >> know. The first in-person meeting of the Board will be 16-17 March in > >> Freiburg, Germany. > >> > >> All of us as the International Biometric Office wish everyone a very > >> happy Holiday Season and a most prosperous New Year! > >> > >> Dee Ann Walker > >> Executive Director > >> Melissa Mays > >> Membership Services Coordinator > >> > >> Al Baggett > >> Director, Education & Certification > >> Lisette Werbowetzki > >> Director, Membership & Marketing > >> > >> > >> > >> > >> > >> > >> *International Biometric Society* > >> 1444 I Street, NW Suite 700 > >> Washington, DC 20005 > >> Ph: 202-712-9049 > >> www.biometricsociety.org<http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856871&m=2417987&u=Bostr…> > >> Unsubscribe here.<http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856872&m=2417987&u=Bostr…> > >> <http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856873&m=2417987&u=Bostr…> > >> > >> 1444 I Street, NW, Suite 700, Washington, DC, 20005 > >> > >> > >> _______________________________________________ > >> 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. > >> > > > > > > > > -- > > "Small steps toward a much better world" > > > > \begin{signature} > > Daniel Marcelino > > Land Phone 1+514 343 6111 #3799 > > 3200 Jean Brillant, Office C5071 > > Montreal, QC; H3T 1N8 > > Canada > > \end{signature} > > > > _______________________________________________ > > 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. > > > > > > -- > Guilherme Moraes Ferraudo > Jaboticabal/Campinas - SP > http://lattes.cnpq.br/2096118558794430 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121208/55b8481c/atta…> > > ------------------------------ > > Message: 4 > Date: Sat, 8 Dec 2012 10:09:58 -0200 > From: Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Grupo de estudo de Geoestatistica > Message-ID: > <CA++u5Oh6yS9U8c2DhhSNs_=uvXLUTqnE-8fAJdnVOZHVGSBvyA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Eu tenho muito interesse. > > Em 8 de dezembro de 2012 09:10, Hélio Gallo Rocha <heliogallorocha(a)gmail.com > > escreveu: > > > Mauro, e colegas. > > Não sabia desta lista, mas eu particularmente tenho dificuldade co m > > a língua inglesa, > > > > Como sou professor de topografia e fazendo mestrado penso da seguinte forma > > > > > > Objetivo: > > Ajudar na resolução de problemas de possíveis interessados como > > professores, agrônomos, ambientalistas. > > > > Metas, pensando em AP > > Preparar mapas precisos > > Calcular dosagens através dos mapas > > Trabalhar co shape files > > etc > > > > Para ajudar na organização, começaríamos em cima de um script > > (dados=S100), tipo passo a passo, assim até quem nunca trabalhou com geoest > > conseguiria segui > > > > Assim construiríamos em equipe um script adaptável a qualquer situação, > > com um moderador e se possível que seja o que tem mais conhecimento do > > assunto, > > Quem gostaria de ser o moderador? > > > > > > Tenho um script quase pronto (que ainda está longe de acabar...) > > > > > > Disponibilizo este script inicial, em cima dos dados S100, fácil de > > adaptar , para auxiliar a > > Estudo da determinação da hipótese da Envelopado > > Estudo de distancia máxima: auxilia a escolha da máxima > > Cinco krigagens (asentimento e 4 validações) na mesma página > > Relatório com AIC (FALTA VALIDAR) > > Exportado para XLS > > > > > > Hélio > > > > > > > > Em 8 de dezembro de 2012 07:43, Eder Comunello <ecomunel(a)gmail.com>escreveu: > > > >> Hélio e colegas, bom dia! > >> > >> Tenho interesse no tema, mas também tenho dúvidas em relação à forma de > >> implementar este grupo. Pra mim é bastante cômodo que seja aqui na lista > >> mesmo, pois aproveito pra aprender sobre o uso do R em outras vertentes que > >> também me interessam. > >> > >> ================================================ > >> Éder Comunello > >> > >> Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) > >> Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] > >> > >> Researcher at Embrapa Western Region Agriculture > >> Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > >> ================================================ > >> UTC-03:00 > >> > >> > >> _______________________________________________ > >> 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. > > > > > > -- > Guilherme Moraes Ferraudo > Jaboticabal/Campinas - SP > http://lattes.cnpq.br/2096118558794430 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121208/41734424/atta…> > > ------------------------------ > > Message: 5 > Date: Sat, 8 Dec 2012 04:38:05 -0800 (PST) > From: Wagner Tassinari <wtassinari(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, Gilbert > Queiroz <gilbert_queiroz(a)yahoo.com.br> > Subject: Re: [R-br] Grupo de estudo de Geoestatistica > Message-ID: <1354970285.6880.YahooMailNeo(a)web39301.mail.mud.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Estou dentro !!!! > > > Wagner Tassinari - UFRRJ > > > > > ________________________________ > De: Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > Para: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Sexta-feira, 7 de Dezembro de 2012 20:34 > Assunto: Re: [R-br] Grupo de estudo de Geoestatistica > > > Estou interessado... > Gilbert-IME > > > > ________________________________ > De: Hélio Gallo Rocha <heliogallorocha(a)gmail.com> > Para: r-br <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Sexta-feira, 7 de Dezembro de 2012 17:53 > Assunto: [R-br] Grupo de estudo de Geoestatistica > > > Prezados, > > Não sei é prática desta lista, mas haveriam interessados em formar um grupo de estudo especifico em geoestatística? > > > Hélio > _______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121208/e6f6d548/atta…> > > ------------------------------ > > Message: 6 > Date: Sat, 8 Dec 2012 11:12:49 -0200 > From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Grupo de estudo de Geoestatistica > Message-ID: > <CAJ+S-MGo54ENwjZdBddR9scNGemuw1-P7R4Z9aGJ-rGjSOgqYQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Considero muito boa iniciativa. gostaria de participar. > > Iván > > Em 8 de dezembro de 2012 10:38, Wagner Tassinari > <wtassinari(a)yahoo.com.br>escreveu: > > > Estou dentro !!!! > > > > > > Wagner Tassinari - UFRRJ > > > > > > ------------------------------ > > *De:* Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > > *Para:* "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > > *Enviadas:* Sexta-feira, 7 de Dezembro de 2012 20:34 > > *Assunto:* Re: [R-br] Grupo de estudo de Geoestatistica > > > > Estou interessado... > > Gilbert-IME > > > > ------------------------------ > > *De:* Hélio Gallo Rocha <heliogallorocha(a)gmail.com> > > *Para:* r-br <r-br(a)listas.c3sl.ufpr.br> > > *Enviadas:* Sexta-feira, 7 de Dezembro de 2012 17:53 > > *Assunto:* [R-br] Grupo de estudo de Geoestatistica > > > > Prezados, > > > > Não sei é prática desta lista, mas haveriam interessados em formar um > > grupo de estudo especifico em geoestatística? > > > > > > Hélio > > > > _______________________________________________ > > 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. > > > > > > _______________________________________________ > > 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. > > > > > > -- > Luis Iván Ortiz Valencia > Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ > Estatístico Msc. > Spatial Analyst Msc. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121208/c24c05ac/atta…> > > ------------------------------ > > _______________________________________________ > 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 22, assunto 10 > *****************************************

1 0

Re: [R-br] Grupo de estudo de Geoestatistica
by Eder Comunello 08 Dez '12

08 Dez '12

Hélio e colegas, bom dia! Tenho interesse no tema, mas também tenho dúvidas em relação à forma de implementar este grupo. Pra mim é bastante cômodo que seja aqui na lista mesmo, pois aproveito pra aprender sobre o uso do R em outras vertentes que também me interessam. ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

3 2

Fwd: 2013 IBS Board Election Results
by Leandro Marino 08 Dez '12

08 Dez '12

Gostaria de parabenizar um de nossos membros o Paulo Justiniano por esta conquista que vem junto com muito trabalho. Paulo, parabéns. Abraços, Leandro Marino Enviado de um dispositivo móvel. Sent from a mobile device. ---------- Forwarded message ---------- From: <ibs(a)tibs.org> Date: 7 Dec 2012 18:17 Subject: 2013 IBS Board Election Results To: <leandromarino(a)leandromarino.com.br> Cc: Click here<http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856870&m=2417987&u=Bostr…>for web version. 2013 IBS EXECUTIVE BOARD ELECTION RESULTS Dear IBS Members, Thank you for your outstanding participation in the recent Society-wide election for directors of the new IBS Executive Board. More than 25% of the eligible voters participated in the election and 999 ballots were cast! This is outstanding! We are pleased to announce that the following members have been elected for either a term of two or four years beginning 1 January 2013. Wesley Johnson (2)WNARJames Carpenter (4)BIRJosé Pinheiro (4)ENARAndreas Ziegler (4)DRSharon-Lise Normand (4)ENARPaulo J. Ribeiro (4)RBrasAlan Welsh (4)ARKerri Mengersen (2)ARMichael Schemper (2)RoeSPeter Njuho (2)GSAfVicente Núñez-Antón (2)RespLisa Lavange (2)ENAR The remaining members of the Executive Board include President Clarice Demétrio (RBras), President-Elect John Hinde (BIR) and Secretary-Treasurer Linda Young (ENAR). Please join me in congratulating all of our newly elected Executive Board members. If you have ideas or wish to contribute to IBS as a volunteer in some manner, please contact a member of the Executive Board and let them know. The first in-person meeting of the Board will be 16-17 March in Freiburg, Germany. All of us as the International Biometric Office wish everyone a very happy Holiday Season and a most prosperous New Year! Dee Ann Walker Executive Director Melissa Mays Membership Services Coordinator Al Baggett Director, Education & Certification Lisette Werbowetzki Director, Membership & Marketing *International Biometric Society* 1444 I Street, NW Suite 700 Washington, DC 20005 Ph: 202-712-9049 www.biometricsociety.org<http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856871&m=2417987&u=Bostr…> Unsubscribe here.<http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856872&m=2417987&u=Bostr…> <http://www.mmsend1.com/link.cfm?r=1080888716&sid=21856873&m=2417987&u=Bostr…> 1444 I Street, NW, Suite 700, Washington, DC, 20005

3 2

Re: [R-br] Intersecção entre dois shapefiles no R
by Eder Comunello 08 Dez '12

08 Dez '12

Senhores, bom dia! Além das dicas dos colegas, pode ser interessante dar uma olhada no pacote {aRT} em <http://www.leg.ufpr.br/doku.php/software:art>. O pacote {aRT} faz a integração do R com a biblioteca TerraLib, permitindo assim integrar também o TerraView (INPE). Ainda estou para estudar adequadamente este pacote para saber do seu potencial na operação com 'shapes', mas já na vi na página um exemplo para união de polígonos (em 'Datasets and scripts'). ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

1 0

Correlação Espacial
by Ricardo Dobrovolski 07 Dez '12

07 Dez '12

Caros colegas, Alguém conhece pacotes do R que permitam realizar correlações com correção dos graus de liberdade devido à autocorrelação espacial? Os métodos mais comuns para essa correção estão descritos em: Dutilleul, P. 1993. Modifying the t-test for assessing the correlation between two spatial processes. Biometrics 49: 305-314. Clifford, P., Richardson, S. & Hémon, D. (1989) Assessing the significance of the correlation between two spatial processes. Biometrics, 45, 123–134. Atenciosamente, Ricardo Dobrovolski

2 1

Semivariogramas experimentais com envelopes simulados
by Hélio Gallo Rocha 07 Dez '12

07 Dez '12

Prezados, boa tarde. Lendo uma dissertação sobre geoestatística aplicada ao levantamento de bicho-mineiro em café. Ela usa um recurso para verificar se o atributo possui dependência espacial, chamado de envelope simulado, interessante. Minha dúvida é sobre o parâmetro nsim, assim descrito: nsim - número de simulações utilizadas para calcular os envelopes. O padrão é 99. Seria o nível de significância? Se não for, teria como verificar e mudar esse valor? Exemplos s100.ml <- likfit (s100, ini c = (0,5, 0,5), fix.nugget = TRUE) s100.vario <- variog (s100, max.dist = 1) s100.env <- variog.model.env (s100, obj.v = s100.vario, model.pars = s100.ml) parcela (s100.vario, env = s100.env) Grato Hélio

2 2

ajuda com gráfico
by D'Apolito 07 Dez '12

07 Dez '12

Desculpem não saber o nome desse tipo de gráfico para colocar no título do e-mail. Alguém sabe uma função/pacote que faça esse tipo de gráfico (os da parte superior)?? http://ars.els-cdn.com/content/image/1-s2.0-S0367253011000703-gr2.jpg grato Carlos

8 8

Erro: Intersecção entre dois shapefiles no R
by Ricardo Dobrovolski 07 Dez '12

07 Dez '12

Oi Pessoal, Tentei realizar a intersecção entre dois shapes por meio da função gIntersect do pacote rgeos. No entanto, ocorreu o seguinte erro, conforme mostrado abaixo: > wdpa<- readShapePoly(file.choose()) > grid_05<- readShapePoly(file.choose()) > inter_wdpa_grid<-gIntersection(wdpa,grid_05) Erro em createPolygonsComment(p) : rgeos_PolyCreateComment: orphaned hole, cannot find containing polygon for hole at index 2 Em 6 de dezembro de 2012 20:22, Ricardo Dobrovolski <rdobrovolski(a)gmail.com>escreveu: > Oi Pessoal, > > Eu gostaria de saber se vocês conhecem pacotes no R nos quais eu possa > fazer operações entre shapefiles como intersecções, clips, etc, de maneira > semelhante às operações que se podem realizar no ArcGis. > > Atenciosamente, > > Ricardo Dobrovolski >

2 1

Tinn-R 2.4.1.5
by Jose Claudio Faria 07 Dez '12

07 Dez '12

Caros usuários do Tinn-R, Uma nova versão estável do editor Tinn-R está disponível no servidor SourceForge: http://sourceforge.net/projects/tinn-r/ As novidades: http://sourceforge.net/p/tinn-r/news/?source=navbar Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Derivar equações com R
by Fernando Antonio de souza 07 Dez '12

07 Dez '12

Caros colegas Gerei uma equação da forma = log10(PB)~a+b*log10(PCVZ) o qual relaciona o conteúdo de proteína corporal animal em função de seu peso corporal vazio. Segundo a metodologia que estou seguindo esta equação deve ser derivada para se estimar o conteúdo de PB por Kg de PCVZ. Há uma forma de fazer esta derivada? Alguém poderia me ensinar como fazer? grato pela atenção. modelo5<-lm(log10(PB)~log10(PCVZ),data=PB) summary(modelo5) dput(PB)structure(list(Animal = c(25L, 38L, 42L, 50L, 53L, 18L, 22L, 34L, 41L, 44L, 48L, 2L, 12L, 14L, 39L, 46L, 4L, 5L, 21L, 26L, 49L, 3L, 32L, 37L, 43L, 47L, 8L, 15L, 30L, 40L, 1L, 17L, 23L, 28L, 7L, 9L, 19L, 35L), Gest = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L), Manej = c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), rep = c(4L, 5L, 3L, 1L, 2L, 1L, 2L, 6L, 4L, 3L, 5L, 3L, 2L, 5L, 1L, 4L, 2L, 4L, 1L, 3L, 5L, 1L, 5L, 2L, 3L, 4L, 1L, 4L, 3L, 2L, 2L, 3L, 4L, 1L, 3L, 2L, 4L, 1L), PV = c(40.5, 40.9, 35.5, 34.5, 36.7, 26.9, 35.5, 31.2, 31.5, 30, 35.5, 40.8, 36.6, 36.6, 36.5, 37.9, 41.9, 35, 35.5, 32.5, 42, 51.8, 50, 46.5, 38.7, 42.7, 35.4, 38.4, 36.5, 31.2, 49, 46.4, 48.4, 42.8, 40.5, 42.3, 50.2, 33.5), PCVZ = c(34.7811, 33.6882, 22.8156, 26.4014, 26.7112, 19.1745, 26.6094, 22.4639, 25.6466, 22.1737, 25.773, 35.4571, 30.8157, 28.1626, 31.8387, 30.3712, 32.9897, 26.5665, 27.9632, 23.4118, 38.1716, 42.3432, 38.2408, 36.4839, 31.3696, 36.7408, 25.9655, 28.1035, 28.9001, 26.0481, 39.8412, 37.9058, 40.8529, 43.4908, 30.3399, 28.4752, 34.3194, 24.6319), PB = c(9.18910573, 9.11864608, 5.653976965, 6.870682232, 6.838822605, 5.119182736, 6.932858931, 6.047701201, 5.593529478, 5.379896679, 6.611844309, 9.230423109, 8.184234256, 7.20898144, 8.703687287, 8.645435859, 7.017519902, 5.976945686, 7.366941986, 4.644961891, 9.085415538, 11.06095042, 9.461960295, 8.85340357, 8.008399724, 7.058114299, 5.333367984, 5.670736687, 7.363494857, 6.014652251, 7.077788459, 7.890982758, 8.935434043, 9.484189031, 6.804250745, 7.034024007, 7.628728687, 5.124003503)), .Names = c("Animal", "Gest", "Manej", "rep", "PV", "PCVZ", "PB"), class = "data.frame", row.names = c(NA, -38L))

3 6

Intersecção entre dois shapefiles no R
by Ricardo Dobrovolski 06 Dez '12

06 Dez '12

Oi Pessoal, Eu gostaria de saber se vocês conhecem pacotes no R nos quais eu possa fazer operações entre shapefiles como intersecções, clips, etc, de maneira semelhante às operações que se podem realizar no ArcGis. Atenciosamente, Ricardo Dobrovolski

3 2

Ler um arquivo grande
by Daniel Marcelino 06 Dez '12

06 Dez '12

Caros, Estou tentando ler um arquivo de texto com 930 mb, mas sem sucesso desde ontem. Meu computador é multicore Macbook 8 gb memória, mas não consegue finalizar de ler o arquivo e "congela". Alguém tem alguma ideia de como superar isso sem precisar abrir o arquivo e dividir manualmente o banco? Daniel

4 5

Problemas com Encoding
by Marcos Silva 06 Dez '12

06 Dez '12

Caros amigos, Estou tentando usar o pacote XML para ler alguns dados em tabelas existentes em páginas html. Especificamente, estou tentando ler páginas da previdência social contendo informações sobre o CRP (Certificado de Regularidade Previdenciária) O problema é que estou obtendo uns caracteres estranhos, coisas do tipo "SituaÃ§Ã£o". Estou fazendo da seguinte forma: library(XML) html <- readHTMLTable("F:\\Regime Próprio\\html\\ExtratoPrevidenciario.htm") html <- html[[3]][, c('V1', 'V2')] Já tentei o iconv() sem resultado. Já tentei ler após mudar a opção "encoding" em options(). Ah, estou usando o windows velho de guerra... sorry... :) A página pode ser vista no seguinte link: http://www1.previdencia.gov.br/sps/app/crp/ExtratoRegularidadeRegimes.asp?C… Se alguém tiver alguma idéia de corrigir... Desde já obrigado. Abs. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006

2 3

Ajuda Mapa
by Sérgio Henrique almeida da silva ju 06 Dez '12

06 Dez '12

Amigos Estou plotando um mapa, no qual, eu quero destacar que em um município qualquer tem uma doença x. no meu banco a categoria da variável muniprio é "sim" e "não", "sim" significa que ali tem a doença x. Fiz dessa forma, porém não ficou legal, como posso fazer isso? brks <- c("não", "sim") cols <- c("blue", "white") plot(mapa, col=cols[findInterval(mapa.dt$muniprio, brks)], border="white") Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 2

Poisson multivariada
by Leonard de Assis 06 Dez '12

06 Dez '12

Existe algum pacote que lida com Poisson multivariada? Eu achei alguma coisa no rseek: http://www.rseek.org/?cx=010923144343702598753%3Aboaz1reyxd4&q=bivariate+po… mas o pacote mencionado por eles (bivpois) não está mais no CRAN e está dando erro de instalação quando baixa a ultima versão disponível. Ahhh, meu interesse é uma Poisson bivariada inflada por zeros. -- []s Leonard de Assis http://about.me/ldeassis

2 1

Re: [R-br] Como criar data.frame com linhas contendo NA's
by Eder Comunello 05 Dez '12

05 Dez '12

Taynãna, boa tarde! Poderia enviar a saída do console no email? Para mim a saída do último comando aparece assim: > dados[order(dados$local),] ### ordenando dia local dado1 dado2 1 1 1 0.10 15 3 2 1 0.15 16 5 3 1 NA NA 7 4 1 0.20 20 2 1 2 0.27 23 4 2 2 NA NA 6 3 2 NA NA 8 4 2 0.15 25 > E seu pedido era: "Preciso criar um data.frame full, completando os dados com NA´s, tipo: dia local dado1 dado2 1 1 0.10 15 2 1 0.15 16 3 1 NA NA 4 1 0.20 20 1 2 0.27 23 2 2 NA NA 3 2 NA NA 4 2 0.15 25" Como aparece no teu sistema? ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

8 11

Distribuição Gamma Incompleta
by Leonardo Monteiro 05 Dez '12

05 Dez '12

Boa tarde, to trabalhando com séries de dados diários de precipitação da ANA (Hidroweb) e gostaria de saber se alguem sabe um pct ou script para distr. gamma incompleta, com o objetivo de preenchimento dessas falhas? obrigado, -- Leonardo Amaral Monteiro - Eng° Agrônomo Mestre em "Física do Ambiente Agrícola" - Agrometeorologia Doutorando do PPG em "Engenharia de Sistemas Agrícolas" Escola Superior de Agricultura "Luiz de Queiroz" - Universidade de São Paulo Piracicaba,SP. Fone: (19)34294123 - Ramal:236

3 2

Script ANCOVA
by Murilo de Lima Arantes 05 Dez '12

05 Dez '12

Olá amigos, Alguém podo me passar um script para ANCOVA ou alguma tutorial para esta análise.Ela é feita em lm também? Abraço! Murilo de LimaBiólogo - PPG Biologia de Água Doce e Pesca Interior - INPAManaus-AM (92 - 82434321/36114190 / Skype: murilodlima)

2 1

R e action
by alanarocha＠sapo.pt 05 Dez '12

05 Dez '12

Queria perguntar se alguem trabalha com o action. Ao utilizar a função somario base por factor para alem das estatísticas descritivas aparece tambem limite inferior da média e limite superior da média, o que significam estes limites? Eu penso que tenham a ver com intervalos de confiaça mas não tenho a certeza obrigada. Ana

2 1

uma pergunta
by alanarocha＠sapo.pt 04 Dez '12

04 Dez '12

Sem querer enviei o mail sem assunto. Queria perguntar se alguem trabalha com o action. Ao utilizar a função somario base por factor para alem das estatísticas descritivas aparece tambem limite inferior e limite superior, o que significam estes limites? obrigada. Ana

2 1

(sem assunto)
by alanarocha＠sapo.pt 04 Dez '12

04 Dez '12

caros usuários, não sei se posso fazer esta pergunta: oalguem trabalha com o action? ao fazer um comando sumario de base por factor o resultado para alem das estat=C3=ADsticas descritivas mostra tambem um limite inferior e um limite superior. o que significam estes limites? obrigada. Ana

1 0

multcomp/doBy comparação multipla com interação
by Fernando A. Souza 03 Dez '12

03 Dez '12

Caros colegas Estou fazendo uma análise de um experimento DIC com arranjo fatorial (Gest e Manej). A resposta medida é o peso vivo. A análise de variância mostrou que há efeito de interação entre os dois fatores avaliados. Eu gostaria de saber onde está esta diferenças. A minha dúvida é o seguinte como fazer a matrix de contrastes para a interação? estou utilizando o pacote doBy. A outra dúvida é a seguinte: o R utiliza como padrão o contraste tratamento o qual compara todos os tratamento como o grupo controle( em meu caso Gest=0 e Manej=1) como faço para alterar esta forma de contraste, para um que compare todos os tratamentos entre si ? >dput(abate[,1:5]) structure(list(Animal = c(1L, 2L, 3L, 4L, 5L, 7L, 8L, 9L, 12L, 14L, 15L, 17L, 18L, 19L, 21L, 22L, 23L, 25L, 26L, 28L, 30L, 32L, 34L, 35L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 46L, 47L, 48L, 49L, 50L, 53L, NA, NA), Gest = structure(c(4L, 2L, 3L, 2L, 2L, 4L, 3L, 4L, 2L, 2L, 3L, 4L, 1L, 4L, 2L, 1L, 4L, 1L, 2L, 4L, 3L, 3L, 1L, 4L, 3L, 1L, 2L, 3L, 1L, 1L, 3L, 1L, 2L, 3L, 1L, 2L, 1L, 1L, NA, NA), .Label = c("0", "100", "130", "140"), class = "factor"), Manej = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, NA, NA ), .Label = c("1", "2"), class = "factor"), rep = c(2L, 3L, 1L, 2L, 4L, 3L, 1L, 2L, 2L, 5L, 4L, 3L, 1L, 4L, 1L, 2L, 4L, 4L, 3L, 1L, 3L, 5L, 6L, 1L, 2L, 5L, 1L, 2L, 4L, 3L, 3L, 3L, 4L, 4L, 5L, 5L, 1L, 2L, NA, NA), PV = c(49000L, 40800L, 51800L, 41900L, 35000L, 40500L, 35400L, 42300L, 36600L, 36600L, 38400L, 46400L, 26900L, 50200L, 35500L, 35500L, 48400L, 40500L, 32500L, 42800L, 36500L, 50000L, 31200L, 33500L, 46500L, 40900L, 36500L, 31200L, 31500L, 35500L, 38700L, 30000L, 37900L, 42700L, 35500L, 42000L, 34500L, 36700L, NA, NA)), .Names = c("Animal", "Gest", "Manej", "rep", "PV"), class = "data.frame", row.names = c(NA, -40L)) abate$Gest<-factor(abate$Gest) abate$Manej<-factor(abate$Manej) C<-lm(PV~Gest*Manej,abate) Anova(C,type="III") -- Fernando Antônio de Souza Zootecnista, Msc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

3 4

Áreas correlatas - Hospitais são acusados de substituir médicos por enfermeiros em cirurgias
by Richard Santos 03 Dez '12

03 Dez '12

Se profissionais de áreas correlatas podem exercer as profissões correlatas, esse caso não precisaria ter punição, certo? Porque ciências humanas inclusive é mais intuitivo do que ciências exatas. Quanto à quantidade de conhecimento estatística é tão robusta quanto medicina. E ambas causam mortes quando feito errado. É essa opinião que eu não entendo, para mim, profissão é coisa muito séria. Abraços, Richard http://www.correiobraziliense.com.br/app/noticia/cidades/2012/11/27/interna…

5 8

Ponderar Casos
by Sérgio Henrique almeida da silva ju 03 Dez '12

03 Dez '12

Olá Como posso ponderar meu banco por uma variável x qualquer, para que todos os resultados, summary, table, etc, sejam ponderados por esse variável (0,5, por exemplo), algo equivalente ao WEIGHT BY do SPSS? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

1 0

agradecimentos
by alanarocha＠sapo.pt 03 Dez '12

03 Dez '12

presados, jagradeço a auda de todos e fiquei a perceber melhor em relação as duvidas que coloquei. Ana Rocha

1 0

agradecimentos
by alanarocha＠sapo.pt 03 Dez '12

03 Dez '12

presados, jagradeço a auda de todos e fiquei a perceber melhor em relação as duvidas que coloquei. Ana Rocha

1 0

agradecimentos
by alanarocha＠sapo.pt 03 Dez '12

03 Dez '12

presados, jagradeço a auda de todos e fiquei a perceber melhor em relação as duvidas que coloquei. Ana Rocha

1 0

Re: [R-br] Como criar data.frame com linhas contendo NA's
by Eder Comunello 03 Dez '12

03 Dez '12

Olá, pessoal! Testei as dicas passadas e esta variação também parece atender a pedida. dia <- c(1,2,4,1,4) local <- c(1,1,1,2,2) dado1 <- c(0.10,0.15,0.20,0.27,0.15) dado2 <- c(15,16,20,23,25) dados <- data.frame(dia,local,dado1,dado2); dados tmp <- data.frame(dia=c(1:4,1:4),local=c(rep(1,4),rep(2,4))); tmp dados <- merge(tmp,dados,all=T); dados ### desordenado dados[order(dados$local),] ### ordenando str(dados) ### verifica estrutura Saudações, ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

3 2

Formatação de números na legenda
by Edson Lira 03 Dez '12

03 Dez '12

Bom dia a todos. Como faço para mofidificar estes números que são de uma legenda em um formato tipo: 3.76e+03 em 3760 4.95e+03 em 4950 Meus dados para usar na legenda é: mao.dados$ano2011 <- cut(mao.dados$ano2011, 4) levels(mao.dados$ano2011) que me dá essa saída: "(3.76e+03,4.95e+03]" "(4.95e+03,6.14e+03]" "(6.14e+03,7.33e+03]" "(7.33e+03,8.52e+03]" [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 1

Concursos Ufal - Maceió
by Alejandro C. Frery 03 Dez '12

03 Dez '12

Bom dia pessoal, A Ufal está com vários concursos abertos e, em particular, o IC - Instituto de Computação, abrirá concurso para seis vagas. As inscrições para doutores estão abertas até o 17/12/2012. Seguem os detalhes de titulação exigida, disciplinas vinculadas e pontos de cada uma dessas seis vagas. Mais informações e formulários em http://www.copeve.ufal.br/index.php?opcao=concurso&idConcurso=153 Automação e Controle Graduação em Engenharia Elétrica ou Engenharia de Computação ou Áreas Afins (conforme Engenharias 4 da CAPES/CNPq), e Doutorado em Engenharia Elétrica ou Engenharia de Computação ou Áreas Afins (conforme Engenharias 4 da CAPES/CNPq). DISCIPLINAS: 1. Teoria de Controle; 2. Inteligência Artificial; 3. Sistemas de Eventos Discretos; 4. Projeto de automação e controle; 5. Sinais e sistemas; 6. Eletrônica. PONTOS: 1. Sinais e sistemas contínuos e discretos no tempo: representação e análise no domínio do tempo; 2. Transformada de Laplace e transformada Z: análise, aplicações e caracterização de sistemas; 3. Modelagem de sistemas dinâmicos utilizando função de transferência e variáveis de estado; 4. Análise de resposta transitória de sistemas dinâmicos; 5. Estabilidade de sistemas dinâmicos e métodos de análise de estabilidade: Routh-Hurwitz, Nyquist e análise tipo Lyapunov; 6. Análise e projeto de controladores por método do lugar das raízes; 7. Análise e projeto de controladores por método de resposta frequência; 8. Ações de controle básicas e projeto de controladores PID; 9. Análise e projeto de sistemas com realimentação via espaço de estados, projeto de observadores e servossistemas; 10. Introdução aos sistemas de controle digital. ---------------------------------------------------------------------------------- Banco de Dados e Estrutura de Dados Graduação em Sistema de Informação ou Ciência da Computação ou Engenharia de Computação, e Doutorado em Computação ou Áreas Afins. DISCIPLINAS: 1. Estrutura de Dados; 2. Banco de Dados; 3. Introdução à Engenharia de Computação; 4. Programação 1; 5. Introdução à Engenharia de Computação; PONTOS: 1. Mineração de Dados e Descoberta de Conhecimento; 2. Estratégias e Estruturas de Dados para Indexação em SGBDs; 3. Algoritmos e Estruturas de Dados Não Lineares; 4. Aplicações das Estruturas de Dados Não Lineares; 5. Algoritmos de Busca e Ordenação; 6. Indexação Semântica de Dados; 7. Integração e Interoperabilidade de Dados; 8. Banco de Dados Multidimensionais e Temporais; 9. Recuperação de Dados e Controle de Concorrência; 10. Processamento e Otimização de Consulta. ---------------------------------------------------------------------------------- Circuitos Elétricos e Eletrônicos Graduação em Engenharia Elétrica ou Engenharia de Computação ou Áreas Afins, e Doutorado em Engenharia Elétrica, ou Engenharia de Computação, ou Áreas Afins (conforme Engenharias 4 da CAPES/CNPq). DISCIPLINAS: 1. Introdução à Engenharia de Computação; 2. Circuitos Elétricos; 3. Eletrônica; 4. Circuitos digitais; 5. Sistemas Digitais; 6. Organização e Arquitetura de Computadores; 7. Disciplinas Eletivas na Área de Circuitos Elétricos e Eletrônicos. PONTOS: 1. Teoria de circuitos elétricos; 2. Circuitos analógicos; 3. Dispositivos eletrônicos; 4. Amplificadores operacionais e suas aplicações; 5. Instrumentação eletrônica; 6. Sensores, Atuadores e condicionadores de sinais; 7. Circuitos digitais; 8. Projeto de sistemas digitais; 9. Organização básica de processadores; 10. Hardware configurável. ---------------------------------------------------------------------------------- Comunicação Analógica e Digital Graduação em Engenharia Elétrica ou Engenharia de Computação ou Matemática Aplicada e Computacional ou Física ou Áreas Afins (conforme Engenharias 4 da CAPES/CNPq), e Doutorado em Engenharia Elétrica ou Engenharia de Computação ou Matemática Aplicada e Computacional ou Física ou Áreas Afins (conforme Engenharias 4 da CAPES/CNPq). DISCIPLINAS: 1. Princípios de Comunicação; 2. Sinais e Sistemas; 3. Probabilidade e Estatística; 4. Processamento Digital de Sinais; 5. Métodos Numéricos. PONTOS: 1. Sinais e sistemas contínuos e discretos no tempo: representação e análise no domínio do tempo; 2. Representação de Fourier para sinais contínuos e discretos no tempo; 3. Modulação de onda contínua; 4. Modulação de pulso; 5. Transmissão de pulsos em banda base; 6. Análise de espaço de sinais e transmissão digital em banda passante; 7. Amostragem e recuperação de sinais (conversão A/D e D/A); 8. Transformada Z: análise, aplicações e caracterização de sistemas; 9. Análise no domínio da frequência de sistemas lineares e invariantes; 10. Realização de sistemas discretos, predição linear e filtragem linear ótima. ---------------------------------------------------------------------------------- Sistemas Embarcados Graduação em Engenharia Elétrica ou Engenharia de Computação ou Ciência da Computação, e Doutorado em Engenharia Elétrica ou Engenharia de Computação ou Ciência da Computação ou Engenharias 4 (CAPES/CNPq). DISCIPLINAS: 1. Introdução à Engenharia de Computação; 2. Sistemas Embarcados; 3. Sistemas Operacionais; 4. Sistemas de Eventos Discretos. PONTOS: 1. Desenvolvimento de software em sistemas embarcados; 2. Computação Móvel e Ubíqua; 3. Sistemas Operacionais em Sistemas Embarcados; 4. Arquitetura de Microprocessadores em Sistemas Embarcados; 5. Processador Digital de Sinais e Aplicações; 6. Projeto de Sistemas Embarcados; 7. Interfaces de Entrada e Saída em Sistemas Embarcados; 8. Redes sem Fio e seus Padrões; 9. Sistemas de Tempo Real; 10. Testes em Sistemas Embarcados. ---------------------------------------------------------------------------------- Software Básico e Arquitetura de Computadores Graduação em Sistema de Informação ou Ciência da Computação ou Engenharia de Computação, e Doutorado em Computação ou Áreas Afins. DISCIPLINAS: 1. Organização e Arquitetura de Computadores; 2. Sistemas Operacionais; 3. Compiladores; 4. Introdução à Engenharia de Computação; 5. Circuitos Digitais. PONTOS: 1. Desenvolvimento de Sistemas Operacionais para Arquitetura Intel x86; 2. Desenvolvimento de Sistemas Operacionais para Arquitetura ARM; 3. Máquinas Virtuais - Conceitos, Projeto e Implementação; 4. Desenvolvimento de Sistemas Operacionais; 5. Sistemas Operacionais de Tempo Real; 6. Sistemas Operacionais Fidedignos; 7. “Internet das Coisas” no contexto de Sistemas Operacionais e Máquinas Virtuais; 8. “Computação nas Nuvens” no contexto de Sistemas Operacionais e Máquinas Virtuais; 9. Compiladores: infraestrutura LLVM; 10. Geração e Otimização de Código de Máquina. -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008

1 0

Re: [R-br] Classe Date
by Eder Comunello 02 Dez '12

02 Dez '12

Daniel, boa noite! Segue uma sugestão, alojando as conversões em novas colunas. O formato de saída é configurável e o uso de as.factor() funcionaria, mas iria cria levels para cada 'data' (n=103). a<- data.frame(Date.end = as.Date(seq(from=as.POSIXct('2012-06-26', tz="UTC"), to=as.POSIXct('2012-10-06', tz="UTC"), by="day"))) a$Date.end2 <- as.integer(format(a$Date.end,"%Y%m%d")) a$Date.end3 <- as.character(format(a$Date.end,"%Y%m%d")) str(a) Espero que ajude, ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

1 0

Classe Date
by Daniel Marcelino 02 Dez '12

02 Dez '12

Caros, tenho um data frame que tem colunas de classe as.Date. Eu tô fazendo uma função para uma tabela eu gostaria que o R considerasse a classe date como factor ou numeric quando executasse a função. a<- data.frame(Date.end = as.Date(seq(from=as.POSIXct('2012-06-26', tz="UTC"), to=as.POSIXct('2012-10-06', tz="UTC"), by="day"))) Aqui eu reproduzo parte da função que converte factor para numeric. Não consegui achar um forma de fazer o mesmo com Date. if (is.factor(x[[i]]) || is.logical(x[[i]])) { x[[i]] <- as.numeric(x[[i]]) rownames(tabela)[i] <- paste(rownames(tabela)[i], "!#", sep = " ") } Alguém tem alguma dica? -- "Small steps toward a much better world" \begin{signature} Daniel Marcelino Land Phone 1+514 343 6111 #3799 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature}

1 0

geoestatistica
by Natalia Martins 02 Dez '12

02 Dez '12

Prezados bom dia, estou cursando uma disciplina de amostragem e tenho como tema de seminario a aplicação de amostragem em geoestatistica. Para este estudo tenho as "linhas" que um trator passa em um talhao e o objetivo é verificar a cada quantas linhas o trator pode passar sem que haja prejuízo nas analises de solo. Acredito ser uma amostragem sistematica na qual devo ir variando os K's e verificando até qual K os resultados sejam iguais quando considero todas as linhas. No entanto no momento de fazer a análise geoestatistica na linha não estou conseguindo. Gostaria de ajuda. Estrutura dos dados LONG LAT NDRE NDVI linha -48.0149775 -21.3646713 0.171400000 0.435000000 1 . . . -48.0149775 -21.3646713 0.171400000 0.435000000 63 Abraços Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

3 5

Como criar data.frame com linhas contendo NA's
by Taynãna César Simões 01 Dez '12

01 Dez '12

Pessoal, Tentei fazer 'for' dentro de 'for', mas está ficando complicado. Suponham que tenho um data.frame com dados espaco-tempo irregulares: foram coletados. dia local dado1 dado2 1 1 0.10 15 2 1 0.15 16 4 1 0.20 20 1 2 0.27 23 4 2 0.15 25 Preciso criar um data.frame full, completando os dados com NA´s, tipo: dia local dado1 dado2 1 1 0.10 15 2 1 0.15 16 3 1 NA NA 4 1 0.20 20 1 2 0.27 23 2 2 NA NA 3 2 NA NA 4 2 0.15 25 Tenho dezenhas de semanas e centenas de locais observados, além de várias covariáveis. Alguma sugestão? Obrigada. Att, Taynãna Fiocruz/RJ

3 6

Cotação Histórica de Ações Bovespa
by T Branquinho 01 Dez '12

01 Dez '12

Caros Colegas da Lista, Alguém já conseguiu importar as cotações históricas das Ações da Bovespa pelo R? Se sim, como o fez? http://www.bmfbovespa.com.br/shared/iframe.aspx?idioma=pt-br&url=http://www… No mais agradeço, Tales

6 6

Re: [R-br] leitura de dados
by Eder Comunello 01 Dez '12

01 Dez '12

Olá, bom dia! Sugiro verificar se o seu arquivo "dados.r" é realmente um arquivo de dados com <file.show('dados.r')>. Com essa extensão era esperado ser o script de comandos. Se for realmente o arquivo de dados com uma extensão inapropriada, verifique os argumentos sep e header de read.table(). Sendo o script, pode ser que os dados sejam passados via linha de comando. Sendo assim, bastaria abrir o script e executar o trecho de código correspondente para criar o objeto de dados. E uma outra possibilidade, já aventada pelos colegas, é tratar-se de uma imagem do workspace (.RData) salva com o nome do script (.R). Nesse caso após ler o arquivo com <load('dados.r')>, você poderá usar o comando <ls()> para visualizar os objetos carregados pela leitura da imagem. Identifique o objeto de dados e, caso queira, é possível separá-lo ('exportá-lo') para não causar mais confusão (pode usar dput(), write.table()). De qualquer modo, procure acertar as extensões. Não é obrigatório, mas é útil pra quem tá começando. Pra ajudar mais, seria necessário seu código (CMR), mas espero que ajude, ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00 >

1 0