- R-br - DINF UFPR

Derivada parcial de função "numérica"
by julio cesar araujo 06 Jun '11

06 Jun '11

Boa noite, Não sei se defini bem o "título" desse e-mail, ms minha situação é a seguinte: Tenho valores de uma variável "G" (é uma cópula, que estimei não-pramétricamente) que encontrei via simulação da combinação de duas outras variáveis "x" e "y", avaliadas num "grid" de probabilidades (quantilico) ui1 e ui2 (que estão definidos entre 0 e 1). Agora preciso de uma derivada parcial de G com relação a U1. Alguem sabe como posso fazer para obter tal valor??? Desde já agradeço. Julio Cesar Araujo Mestrando em Economia

1 0

OF TOPIC
by Edson Lira 06 Jun '11

06 Jun '11

Caros amigos da lista, comprei um Mac, alguém tem algum material que oriente o uso do R no MAC? Edson Lira Estatístico Ma-Am

3 2

Oportunidade para estatistico no Rio de Janeiro
by Pedro Emmanuel Alvarenga Americano do Brasil 06 Jun '11

06 Jun '11

Amigos de R, Desculpem o tópico fora escopo da lista, mas... Recebi um email sobre edital de inovação tecnológica da FIOCRUZ, que permite agregar recursos humanos através de bolsas de pesquisa a projetos cadastrados. http://www.castelo.fiocruz.br/vpplr/inovatec.php Assim, estava pensando em submeter um projeto para desenvolvimento e aperfeiçoamento de funções R. Como todas as bolsas e etc vindas das agencias de fomento no Brasil, ha algumas restrições por parte dos candidatos. O único quesito de minha parte é dominar R e conseguir desenvolver funções (que é quase redundante). Os que estiverem interessados ou que quiserem maiores informações favor entrar em contato com o meu email: emmanuel.brasil(a)gmail.com . Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.openoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

Teste de Tukey! Experimento Desbalanceado! Erro.
by Thiago De paula protásio 04 Jun '11

04 Jun '11

Walmes; Desculpe incomodar novamente, mas não estou conseguindo obter as médias ajustadas. Veja o que está ocorrendo. Talvez, vc possa me ajudar! lapply(levels(da$Tratamento),+ function(i){+ contrast(m0, type="average", list(Tratamento=i, Bloco=levels(da$Bloco)))+ }+ )list() > summary(glht(m0, linfct=mcp(Tratamento="Tukey")))Erro em mcp2matrix(model, linfct = linfct) : Variable(s) ‘Tratamento’ of class ‘integer’ is/are not contained as a factor in ‘model’. Muito obrigado pela ajuda! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

3 2

Matrix Modelo
by Eder David Borges da Silva 04 Jun '11

04 Jun '11

Pessoal, gostaria da ajuda de vocês para criar uma matrix do modelo, como exemplo abaixo: LOCAL <- rep(c('local1','local2'),each=4) Rep <- 1:8 CLONE <- rep(1:2,4) BLOCO <- rep(1:4,each=2) dados <- data.frame(LOCAL,BLOCO,CLONE,Rep) dados$Int <- paste(dados$LOCAL,dados$CLONE,sep='') dados$Obs <- 1 dados$TCH <- sort(rnorm(nrow(dados),10,1)) dados$IND <- as.factor(dados$IND) dados$CLONE <- as.factor(dados$CLONE) dados$BLOCO <- as.factor(dados$BLOCO) summary(dados) op <- options(contrasts=c("contr.treatment", "contr.poly")) MM <- model.matrix(~BLOCO+CLONE+LOCAL:CLONE,dados) MM (Intercept) BLOCO2 BLOCO3 BLOCO4 CLONE2 CLONE1:LOCALlocal2 CLONE2:LOCALlocal2 1 1 0 0 0 0 0 0 2 1 0 0 0 1 0 0 3 1 1 0 0 0 0 0 4 1 1 0 0 1 0 0 5 1 0 1 0 0 1 0 6 1 0 1 0 1 0 1 7 1 0 0 1 0 1 0 8 1 0 0 1 1 0 1 esta foi a matrix criado pelo model.matrix, porem gostaria de uma matrix diferente, pois nesta criada pela model.matrix o bloco1 não aparece pois é estimado por direfença, assim como o clone1, eu quero que todos os niveis aparecem tendo como nome das colunas: (Intercept) BLOCO1 BLOCO2 BLOCO3 BLOCO4 CLONE1 CLONE2 CLONE1:LOCAL1 CLONE1:LOCAL2 CLONE2:LOCAL1 CLONE2:LOCAL2 Alguem tem uma dica? Att

3 5

Teste de Tukey! Experimento Desbalanceado!
by Thiago De paula protásio 04 Jun '11

04 Jun '11

Muito obrigado Walmes..... Esse exemplo do Rídiculas será de grande valia.Abraço! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Teste de Tukey! Experimento Desbalanceado!
by Thiago De paula protásio 04 Jun '11

04 Jun '11

Pessoal; Venho novamente pedir a ajuda de vcs!Tenho que analisar um experimento em blocos, mas os tratamentos possuem diferentes números de repetições. Estou utilizando a soma de quadrados do tipo III. library(car) anava1=aov(QEA~Tratamento+Bloco)Anova(anava1,type="III") Contudo, não estou conseguindo efetuar o teste de Tukey. Alguém pode me ajudar?Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

3 2

grafico de barras
by Gilson Sanchez 04 Jun '11

04 Jun '11

Pessoal, mais uma vez na procura de vc's, é a primeira vez que quero fazer um grafico de barras e eu não sei se eu estou fazendo certo ou errado. Então eu quero fazer um grafico com as seguintes atribuições: 1. quero colocar três barras (chla, chlb, caro) em cada espécie. 2. quero colocar os nomes das espécies no eixo x. 3. no titulo do eixo y = Concentração~de~pigmentos~cloroplastídicos~/~(mu~mol~.~g^{-1}))), quero fazer quebra em "/" 4. colocar a legenda. esp chla chlb caro espécie1 0.25 0.08 0.09 espécie2 0.18 0.06 0.06 espécie3 0.21 0.06 0.07 barplot(dados$chla, density=c(20), col="black", angle=c(45), xlab="Espécies", font.lab=2, axes=TRUE, ylab=expression(bold(Concentração~de~pigmentos~cloroplastídicos~(mu~mol~.~g^{-1}))), ylim=c(0,0.5)) axis(1, at=1:3, labels=dados$esp) par(new=TRUE) barplot(dados$chlb, density=c(21), col="black", angle=c(100)) axes=TRUE, frame=T, ann=F) par(new=TRUE) barplot(dados$caro, density=c(22), col="black", angle=c(200)) Acho que esta facil mas eu não estou conseguindo. Obrigado mais uma vez pela ajuda. -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

3 6

Re: [R-br] Instalação do pacote rparallel
by Junior Beleti 03 Jun '11

03 Jun '11

Boa noite Ivan.. Eu ja vi um com o mesmo título, mas do autor Justin R. Smith. O livro que voce citou não conheço, mas procurei e consegui baixar. Ele cita modelos de computação e análise de algoritmos, e depois descreve algoritmos sendo implementados com computacao paralela. Bem interessante. Obrigado pela dica. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

1 0

[Dúvida] Gerar filmes com imagem no R
by Pedro Rafael 03 Jun '11

03 Jun '11

Pessoal tenho vários mapas de um mesmo geoespaço de uma mesma variável e gostaria de gerar um filme com esses mapas no R -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

3 2

ANOVA por medidas repetidas
by Evandro Malanski 03 Jun '11

03 Jun '11

Pessoal, tenho algumas dúvidas quanto a ANOVA. Se puderem me ajudar, ficarei muito grato. Montei o seguinte experimento com medição de organismos: 3 analisadores diferentes 3 magnificações no microscópio possíveis 2 técnicas de medição diferentes 90 organismos Dos 90 organismos, 30 foram medidos pela magnificação 1, outros 30 na magnificação 2, e os últimos 30 na magnificação 3, por cada um dos 3 analisadores, e usando as 2 técnicas. Ou seja, para cada organismo, tenho 6 medidas. Eu montei a planilha da seguinte forma, mostrando somente as 10 primeiras medições: > compara Larva Técnica Analisador Magnificação Medida 1 1 1 1 20x 3.8303 2 2 1 1 20x 4.2159 3 3 1 1 20x 4.4470 4 4 1 1 20x 4.9549 5 5 1 1 20x 4.0265 6 6 1 1 20x 3.9235 7 7 1 1 20x 3.7697 8 8 1 1 20x 3.9166 9 9 1 1 20x 3.5777 10 10 1 1 20x 3.7515 ... Como eu tenho medidas repetidas, usei a seguinte análise do R: compara.aov<- aov(Medida~Analisador+Magnificação+Técnica+ Analisador:Magnificação+Analisador:Técnica+Magnificação:Técnica+ Analisador:Magnificação:Técnica) Ou seja, até mesmo a interação entre os fatores analiso. Contudo, agora que enviei para correção o meu trabalho, questionaram que a análise está incorreta. Segue o sumário da análise que fiz, e destaco o grau de liberdade do resíduo, que foi a partir daí que me questionaram: > summary(compara.aov) Df Sum Sq Mean Sq F value Pr(>F) Analisador 2 0.43 0.21 0.2253 0.7984 Magnificação 2 1278.94 639.47 672.4928 <2e-16 *** Técnica 1 0.03 0.03 0.0288 0.8654 Analisador:Magnificação 4 0.13 0.03 0.0330 0.9979 Analisador:Técnica 2 0.44 0.22 0.2336 0.7918 Magnificação:Técnica 2 0.86 0.43 0.4521 0.6366 Analisador:Magnificação:Técnica 4 0.10 0.03 0.0269 0.9986 Residuals 522 496.37 0.95 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Esta análise eu fiz para observar se haveriam diferenças significativas entre técnicas. Falaram que minha comparação não está baseada nas 90 larvas medidas, mas sim em todas as 540 medidas que eu tenho, como verificado nos graus de liberdade dos resíduos. Alguém saberia me explicar melhor isso? E, como eu poderia solucionar esse erro? Modificando a organização da planilha dos dados? Modificando a análise? Agradeço a ajuda. ____________________ Evandro Malanski, Oc. Curso de Mestrado em Oceanografia Biológica Laboratório de Ecologia do Ictioplâncton - LEI Instituto de Oceanografia - IO Universidade Federal de Rio Grande - FURG +55(53)3233-6529 / +55(53)8439-6755 Master Course in Biological Oceanography Ichthyoplankton Ecology Laboratory - LEI Oceanography Institute - IO Federal University of Rio Grande - FURG +55(53)3233-6529 / +55(53)8439-6755

7 9

Re: [R-br] ANOVA por medidas repetida
by Ivan Bezerra Allaman 03 Jun '11

03 Jun '11

Evandro, Desculpe de novo!! Desconsidere o modelo que lhe passei, pois vc não tem as mesmas larvas dentro de cada analisador. Você teria aí um caso de efeitos aninhados. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] ANOVA por medidas repetida
by Ivan Bezerra Allaman 03 Jun '11

03 Jun '11

Desculpe Evandro, agora que fui entender!! Bão a minha opinião é o seguinte. Na verdade, a sua repetição são os analisadores, que segundo o seu objetivo, entrariam como fonte de blocagem, afinal de contas, você quer apontar a melhor técnica independentemente de que os analisa não é mesmo? Logo, se o interesse também é avaliar se existe alguma interação entre técnicas e magnificação, o modelo ficaria da seguinte forma (foi como eu enxerguei as coisas). modelo <- lm(MEDIDA ~ ANALISADOR(Bloco) + Magnificação + Técnica + Magnificação:Técnica) anova(modelo) A grosso modo seria isso. Vamos esperar os demais colegas opinarem. Abraço. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] ANOVA por medidas repetidas
by Ivan Bezerra Allaman 03 Jun '11

03 Jun '11

Evandro, Realmente tem alguma coisa estranha. Os seus fatores deveriam estar diluídos no total de observações (90 larvas) se é que me entende. Terás que enviar sua base para que possamos ajudá-lo. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

[Dúvida] Listar todos métodos de um pacote
by Pedro Rafael 03 Jun '11

03 Jun '11

Pessoal como faço para ver todos os métodos de um determinado pacote? Pedro

3 2

Re: [R-br] adicionar texto ou legenda em barchart
by Alexandre Heinemann 03 Jun '11

03 Jun '11

Olá Walmes, Ambas as soluções resolveram o problema. Muito obrigado. Alexandre

1 0

Re: [R-br] Instalação do pacote rparallel
by Ivan Bezerra Allaman 02 Jun '11

02 Jun '11

Boa noite Júnior! Estava procurando algo que não tem nada a ver com o termo "parallel", no entanto, encontrei o seguinte livro: "The Design and Analysis of Parallel Algorithms" do autor "Selim G. Aki". Conhece? Te interessa? Tem o livro? Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

teste de Homogeneidade e Normalidade
by Eduardo Vargas 02 Jun '11

02 Jun '11

Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de médias de germinação de sementes, e as varianveis são tratamentos(trat) fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em `[.data.frame`(x, complete.cases(x)) : undefined columns selected Será que alguem pode me ajudar? -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com

10 24

Re: [R-br] teste de Homogeneidade e Normalidade
by Ivan Bezerra Allaman 02 Jun '11

02 Jun '11

Adriano, Para dizer que alguém está doido você precisa ter argumentos mais forte do que o dele. E pelo o que eu conheço dos dois, você precisará estudar muitoooooooooooooooooo para ter um argumento melhor do que a do Prof. Walmes. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Fwd: [Dúvida] Operador @
by Pedro Rafael 02 Jun '11

02 Jun '11

Pessoal qual a funcionalidade do operador @ e qual a diferença entre ele e o operador $? Pedro

5 4

adicionar texto ou legenda em barchart
by Alexandre Heinemann 02 Jun '11

02 Jun '11

Caros Colegas, Gostaria de adicionar o seguinte texto "a)" em um gráfico barchart. Esse texto deve ficar acima do grafico do lado esquerdo. Entretanto, qdo insiro o comando no codigo do barchart #legend=list(top=list(fun=textGrob("a)"), x=0, just="left", gp=gpar(cex=0.8))) a legenda atual "1-Nov, 1-Dec, 31-Dec" que está acima some. Segue abaixo o codigo que estou utilizando: #----------início---------------------------------------------- library(lattice) library(grid) T1 <- expand.grid(modelo=rep(c("A", "B", "C", "D", "E"),12)) T2 <- expand.grid (data = rep(c(rep(c("1-Nov"),5),rep(c("1-Dec"),5), rep(c("31-Dec"),5)),4)) T3 <- c(rep(c("X1"),15), rep(c("X2"),15),rep(c("X3"),15),rep(c("X4"),15)) T4 <- rnorm(60, mean=10, sd=5) teste <- cbind(T1,T2,T3,T4) teste$data <- factor(teste$data, levels = c("1-Nov","1-Dec","31-Dec" )) teste$modelo = factor(teste$modelo, levels = c("A", "B", "C", "D", "E")) teste$T3 <- factor(teste$T3, levels = c("X1","X2","X3","X4" )) barchart (modelo ~ teste$T4 |teste$T3 , data = teste, xlab = "XXXXX",ylab = " XXXXX", groups = teste$data, auto.key = list( space="top", columns = 3, points = FALSE, rectangles = TRUE, cex=0.8, size=1.4, adj=1,between=0.2, between.colums=0.1),col=gray.colors(3)[1:3], box.ratio = 1.5,aspect=0.6, layout=c(1,4),par.settings = list(superpose.polygon = list(col =gray.colors(3)),fontsize=list(text=9)), par.strip.text=list(cex=0.9),scales = list(x = list(tick.number = 6)), strip=strip.custom(bg="gray95",factor.levels=c("X1","X2","X3","X4")), panel=function(x, y, ...) { panel.barchart(x, y, ...) panel.abline(v=(9.8), col="gray60", lty=2) }) #----------------fim---------------------------------------------------------------------------- Atenciosamente, Alexandre

2 2

[Dúvida] Como pintar um poligono específico de um Mapa
by Pedro Rafael 02 Jun '11

02 Jun '11

Pessoal como faço para pintar um polígono específico de um mapa que importei no R com o comando readShapePoly? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 2

Re: [R-br] [Dúvida ] Importar nome dos municípios no MAPA
by Pedro Rafael 02 Jun '11

02 Jun '11

Como posso pintar um poligono especifico do mapa que eu importei com uma cor diferente? Pedro Em 2 de junho de 2011 14:28, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Como posso pintar um poligono especifico do mapa que eu importei com uma > cor diferente? > > Pedro > > Em 2 de junho de 2011 13:17, Pedro Andrade <pedro.andrade(a)inpe.br>escreveu: > > Procure na pagina do prof. Paulo Justiniano: www.leg.ufpr.br/~paulojus >> >> >> Em 02/06/2011 12:30, Pedro Rafael escreveu: >> >>> Valeu Pedro. Tu sabe de alguma apostila sobre geoestatística no R? >>> >>> Pedro >>> >>> Em 2 de junho de 2011 11:59, Pedro Andrade <pedro.andrade(a)inpe.br >>> <mailto:pedro.andrade@inpe.br>> escreveu: >>> >>> O R não é bom para estas coisas. uma solucao é usar o comando pdf >>> para salvar o mapa vetorial, e depois abrir ele no acrobat para >>> poder fazer o zoom. >>> >>> >>> Em 02/06/2011 11:55, Pedro Rafael escreveu: >>> >>> Existe alguma forma de dar zoom no mapa sem perder a qualidade? >>> >>> Pedro >>> >>> Em 2 de junho de 2011 11:36, Pedro Andrade >>> <pedro.andrade(a)inpe.br <mailto:pedro.andrade@inpe.br> >>> <mailto:pedro.andrade@inpe.br <mailto:pedro.andrade@inpe.br>>> >>> >>> escreveu: >>> >>> >>> Existe um exemplo em >>> http://r-spatial.sourceforge.net/gallery/ >>> >>> library(sp) >>> library(maptools) >>> nc2 <- readShapePoly(system.file("shapes/sids.shp", >>> package="maptools")[1], proj4string=CRS("+proj=longlat >>> +datum=NAD27")) >>> plot(nc2) >>> invisible(text(getSpPPolygonsLabptSlots(nc2), >>> labels=as.character(nc2$NAME), cex=0.4)) >>> >>> um abraço, >>> >>> Pedro >>> >>> >>> Em 02/06/2011 11:14, Pedro Rafael escreveu: >>> >>> Pessoal como faço para colocar dentro de cada poligono o >>> nome do >>> município? Importei o shapefile do estado da paraíba e >>> queria que os >>> municípios (pologonos) ficassem com o nome. Existe >>> alguma forma? >>> Pedro >>> >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr.br> >>> <mailto:R-br@listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr.br >>> >> >>> >>> >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>> >>> >>> -- >>> Pedro Ribeiro de Andrade, Dr. >>> Earth System Science Center (CCST) >>> National Institute for Space Research (INPE) >>> Sao Jose dos Campos, Brazil >>> >>> voice +55 12 3208 7130 >>> fax +55 12 3208 6468 >>> >>> >>> >>> >>> -- >>> Saudações, >>> Pedro Rafael Diniz Marinho - Estatístico SES-PB >>> Currículo >>> Lattes: >>> http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 >>> >>> >>> -- >>> Pedro Ribeiro de Andrade, Dr. >>> Earth System Science Center (CCST) >>> National Institute for Space Research (INPE) >>> Sao Jose dos Campos, Brazil >>> >>> voice +55 12 3208 7130 >>> fax +55 12 3208 6468 >>> >>> >>> >>> >>> -- >>> Saudações, >>> Pedro Rafael Diniz Marinho - Estatístico SES-PB >>> Currículo >>> Lattes: >>> http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 >>> >>> >> -- >> Pedro Ribeiro de Andrade, Dr. >> Earth System Science Center (CCST) >> National Institute for Space Research (INPE) >> Sao Jose dos Campos, Brazil >> >> voice +55 12 3208 7130 >> fax +55 12 3208 6468 >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 1

Fwd: [Dúvida ] Importar nome dos municípios no MAPA
by Pedro Rafael 02 Jun '11

02 Jun '11

Valeu Pedro. Tu sabe de alguma apostila sobre geoestatística no R? Pedro ---------- Mensagem encaminhada ---------- De: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> Data: 2 de junho de 2011 12:30 Assunto: Re: [R-br] [Dúvida ] Importar nome dos municípios no MAPA Para: pedro.andrade(a)inpe.br Valeu Pedro. Tu sabe de alguma apostila sobre geoestatística no R? Pedro Em 2 de junho de 2011 11:59, Pedro Andrade <pedro.andrade(a)inpe.br> escreveu: O R não é bom para estas coisas. uma solucao é usar o comando pdf para > salvar o mapa vetorial, e depois abrir ele no acrobat para poder fazer o > zoom. > > > Em 02/06/2011 11:55, Pedro Rafael escreveu: > >> Existe alguma forma de dar zoom no mapa sem perder a qualidade? >> >> Pedro >> >> Em 2 de junho de 2011 11:36, Pedro Andrade <pedro.andrade(a)inpe.br >> <mailto:pedro.andrade@inpe.br>> escreveu: >> >> >> Existe um exemplo em http://r-spatial.sourceforge.net/gallery/ >> >> library(sp) >> library(maptools) >> nc2 <- readShapePoly(system.file("shapes/sids.shp", >> package="maptools")[1], proj4string=CRS("+proj=longlat +datum=NAD27")) >> plot(nc2) >> invisible(text(getSpPPolygonsLabptSlots(nc2), >> labels=as.character(nc2$NAME), cex=0.4)) >> >> um abraço, >> >> Pedro >> >> >> Em 02/06/2011 11:14, Pedro Rafael escreveu: >> >> Pessoal como faço para colocar dentro de cada poligono o nome do >> município? Importei o shapefile do estado da paraíba e queria que >> os >> municípios (pologonos) ficassem com o nome. Existe alguma forma? >> Pedro >> >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr.br> >> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> >> >> -- >> Pedro Ribeiro de Andrade, Dr. >> Earth System Science Center (CCST) >> National Institute for Space Research (INPE) >> Sao Jose dos Campos, Brazil >> >> voice +55 12 3208 7130 >> fax +55 12 3208 6468 >> >> >> >> >> -- >> Saudações, >> Pedro Rafael Diniz Marinho - Estatístico SES-PB >> Currículo >> Lattes: >> http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 >> >> > -- > Pedro Ribeiro de Andrade, Dr. > Earth System Science Center (CCST) > National Institute for Space Research (INPE) > Sao Jose dos Campos, Brazil > > voice +55 12 3208 7130 > fax +55 12 3208 6468 > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

Re: [R-br] Criando arquivo com dados - Extracão
by Ivan Bezerra Allaman 02 Jun '11

02 Jun '11

?subset Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] [Dúvida ] Importar nome dos municípios no MAPA
by Pedro Rafael 02 Jun '11

02 Jun '11

Grande Pedro muito obrigado pela ajuda. Pedro Rafael Diniz Marinho Em 2 de junho de 2011 11:41, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Grande Pedro, muito obrigado cara. > > Forte abraço, > > Pedro Rafael Diniz Marinho > > Em 2 de junho de 2011 11:36, Pedro Andrade [via R-br] < > ml-node+3568480-1415896092-223914(a)n4.nabble.com> escreveu: > >> Existe um exemplo em http://r-spatial.sourceforge.net/gallery/ >> >> library(sp) >> library(maptools) >> nc2 <- readShapePoly(system.file("shapes/sids.shp", >> package="maptools")[1], proj4string=CRS("+proj=longlat +datum=NAD27")) >> plot(nc2) >> invisible(text(getSpPPolygonsLabptSlots(nc2), >> labels=as.character(nc2$NAME), cex=0.4)) >> >> um abraço, >> >> Pedro >> >> >> Em 02/06/2011 11:14, Pedro Rafael escreveu: >> >> > Pessoal como faço para colocar dentro de cada poligono o nome do >> > município? Importei o shapefile do estado da paraíba e queria que os >> > municípios (pologonos) ficassem com o nome. Existe alguma forma? >> > Pedro >> > >> > >> > >> > _______________________________________________ >> > R-br mailing list >> > [hidden email] <http://user/SendEmail.jtp?type=node&node=3568480&i=0> >> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> >> -- >> Pedro Ribeiro de Andrade, Dr. >> Earth System Science Center (CCST) >> National Institute for Space Research (INPE) >> Sao Jose dos Campos, Brazil >> >> voice +55 12 3208 7130 >> fax +55 12 3208 6468 >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3568480&i=1> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> >> >> ------------------------------ >> If you reply to this email, your message will be added to the discussion >> below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Importar-nome-dos-municipios-… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >> >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

[Dúvida ] Importar nome dos municípios no MAPA
by Pedro Rafael 02 Jun '11

02 Jun '11

Pessoal como faço para colocar dentro de cada poligono o nome do município? Importei o shapefile do estado da paraíba e queria que os municípios (pologonos) ficassem com o nome. Existe alguma forma? Pedro

2 1

RGoogleDocs
by Claude Thiago Arrabal 02 Jun '11

02 Jun '11

Caros, Estou utilizando o pacote RGoogleDocs para acessar uma *worksheet* do excel na minha conta do google. Entretanto me deparo com o seguinte erro: library(RGoogleDocs) sheets.con = getGoogleDocsConnection(getGoogleAuth("claude.arrabal(a)gmail.com", ********)) getWorksheets("dados",sheets.con) Erro em getDocs(con, what = "spreadsheets") : problems connecting to get the list of documents: Token invalid (401) consigo fazer a conexão, mas não estou conseguindo pegar as informações. Qualquer ajuda será bem vinda. Atenciosamente. Claude Thiago Arrabal Mestrando em Estatística Universidade Federal de São Carlos Contato: (016) 8824-8983 [image: e.png] <http://www.ppgest.ufscar.br/>[image: o.png]<http://www.orkut.com.br/Main#Profile?uid=2214128981492139780>[image: i.png] <http://www.linkedin.com/in/arrabal>[image: f.png]<http://www.facebook.com/claude.arrabal>[image: t.png] <http://www.twitter.com/arrabalclt>

1 0

Exportanto factors as numeric
by Daniel Marcelino 02 Jun '11

02 Jun '11

Caros, há algum argumento para exportar factors com labels atribuídos sem estes labels? Quando eu exporto usando o write.table a função passa os labels como strings, mas eu queria manter os valores como numéricos. Alguém sabe se dá para conseguir isso com a função write.table? Eu não achei nada. Daniel

3 7

Re: [R-br] Instalação do pacote rparallel
by Junior Beleti 01 Jun '11

01 Jun '11

Boa noite Ivan. Obrigado pela dica.. Instalei e encontrei varios resultados.. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

1 0

Re: [R-br] Instalação do pacote rparallel
by Ivan Bezerra Allaman 01 Jun '11

01 Jun '11

Boa noite Junior!! Infelizmente não é a minha área, e eu não poderei lhe ajudar. No entanto, instale o pacote "sos", talvez ache alguma coisa. Eu fiz o seguinte: library(sos) ???"parallel packages" e ???"parallel programming" Apareceu vários resultados entre eles os pacotes que o Benilton citou. Utilize uma palavra chave mais adequada para buscas, é bem provável que este pacote de busca irá lhe ajudar um pouco. Abraço. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] Instalação do pacote rparallel
by Junior Beleti 01 Jun '11

01 Jun '11

Olá Benilton, nada especifico, mas pelo que eu vi sobre ele, me parece uma boa solução para aplicações paralelas pelo fato de "automatizar" o paralelismo. Porém não encontrei nenhum exemplo de sua utilização. O snow voce utiliza em maquina multicore mesmo ou cluster? Vou dar uma olhada nos outros pacotes. Obrigado -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

RES: RES: Download de dados pelo R
by Humberto Hazin (hotmail) 01 Jun '11

01 Jun '11

Olá Victor, Obrigado pela ajuda! Mais o problema não é esse! Não sei porque ele não esta reconhecendo meu diretório, porem nesse exemplo ele vai, só que tenho que colocar arquivo por arquivo! cho<-"http://aspera.jpl.nasa.gov/download/pub/sea_surface_temperature/avhrr/ pathfinder/data_v5/monthly/night/04km/1986/bsst/198601.s04m1pfv50-bsst-16b.h df" download.file(cho, destfile = "~/Dropbox/Humberto/Analises em R/ambientais/.hdf", method="auto", quiet = FALSE, mode = "w",cacheOK = FALSE) se alguem puder souber como resolver isso agradeço abs, Humberto De: r-br-bounces(a)listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Victor Eduardo Enviada em: Wednesday, June 01, 2011 9:47 AM Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] RES: Download de dados pelo R Outra coisa, estava olhando o seu erro e repare o seguinte: C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf' No finalzinho, tem um espaço entre 2001 e .hdf. De repente pode ser isso, confere direitinho o seu paste todo e ve se está encaminhando certo para o diretório. Em 1 de junho de 2011 09:44, Victor Eduardo <victorduca08(a)gmail.com> escreveu: Humberto, eu não consegui fazer o download aqui, de repente pode ser problema pelo site mesmo. Em 1 de junho de 2011 09:04, Humberto Hazin (hotmail) <hghazin(a)hotmail.com> escreveu: Olá Eduardo Ele deve estar fora de ar nesse instante ! Tenta com esse aqui http://aspera.jpl.nasa.gov/download/pub/sea_surface_temperature/avhrr/pathfi nder/data_v5/monthly/night/04km/1986/bsst/ ou http://oceandata.sci.gsfc.nasa.gov/SeaWiFS/Mapped/Monthly/9km/chlor/ Abs, Humberto De: r-br-bounces(a)listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Victor Eduardo Enviada em: Wednesday, June 01, 2011 8:58 AM Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Download de dados pelo R Humbero, não está dando para acessar o link enviado. Em 1 de junho de 2011 08:54, Humberto Hazin (hotmail) <hghazin(a)hotmail.com> escreveu: Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

4 8

Re: [R-br] Instalação do pacote rparallel
by Junior Beleti 01 Jun '11

01 Jun '11

Olá Ivan. Este pacote não é um pacote oficial, mas encontrei artigos sobre programação paralela com ele. Estou desenvolvendo minha pesquisa sobre pacotes paralelos no R. Conhece algum pacote que possa ser utilizado? Tentei utilizar o pnmath mas não tive sucesso. Até mais. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

graficos com dois eixos y
by Gilson Sanchez 01 Jun '11

01 Jun '11

Pessoal mais uma vez procuro ajuda de vc's, estou querendo montar um gráfico com dois eixos y. Então o seguinte como eu posso fazer o seguinte: 1. Colocar nomes dos anos no eixo X (jan, fev, mar, abr...dez). 2. oC quero colocar "o" sobreescrita. 3. Como coloco o nome do segundo eixo y. plot(dados$temp, type="b", pch=15, lwd=1.5, xlab="Meses (Anos 2010)", ylab="Temperatura (oC)", ylim=c(24,30), xlim=c(1,12), main= "A", font.main=2, font.lab=2) par(new=T) plot(dados$umidade, type="b", pch=14,axes=F, frame=T, ann=F) axis(4, mtext(text="Umidade (%))", 4, line = 3)) dados meses temp umidade rad chuva 1 Jan 26.49 86.58 795.88 0.36 2 Fev 26.65 88.49 710.24 0.34 3 Mar 27.19 86.16 772.99 0.21 4 Abr 26.28 89.75 574.88 0.67 5 Mai 26.62 89.22 614.02 0.31 6 Jun 26.13 87.83 680.08 0.26 7 Jul 25.83 86.57 675.97 0.15 8 Ago 27.05 83.14 756.44 0.07 9 Set 27.60 83.02 925.57 0.14 10 Out 27.44 85.16 927.71 0.17 11 Nov 26.56 88.18 788.87 0.19 12 Dez 25.87 90.63 703.94 0.33 Obrigado se puderem ajudar -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

2 4

Download de dados pelo R
by Humberto Hazin (hotmail) 01 Jun '11

01 Jun '11

Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto

3 5

modifica dataframe
by Samuel Carvalho 01 Jun '11

01 Jun '11

Boa noite caros(a) senhores(as) membros Por gentileza alguém poderia me auxliar em um código que a princípio parece ser simples mas nao consegui uma solução Abaixo mostro o codigo e a seguir venho comentando ### Codigo R dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), ht=c(10, 13, 12.5, 18, 14)) dados dados$hi<-seq(0.1, dados$ht, by=0.5) ### O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" que corresponde aquele "idfuste", replicando também cada observação para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo que a ultima observação de cada idfuste deve ser igual ao valor de ht Em resumo queria um data.frame de saída a partir dos dados de entrada, nesta estrutura idfuste dap ht hi 1 8,5 10 0 1 8,5 10 0,5 1 8,5 10 1 1 8,5 10 1,5 1 8,5 10 2 1 8,5 10 2,5 1 8,5 10 3 1 8,5 10 3,5 1 8,5 10 4 1 8,5 10 4,5 1 8,5 10 5 1 8,5 10 5,5 1 8,5 10 6 1 8,5 10 6,5 1 8,5 10 7 1 8,5 10 7,5 1 8,5 10 8 1 8,5 10 8,5 1 8,5 10 9 1 8,5 10 9,5 1 8,5 10 10 2 11 13,2 0 2 11 13,2 0,5 2 11 13,2 1 2 11 13,2 1,5 2 11 13,2 2 2 11 13,2 2,5 2 11 13,2 3 2 11 13,2 3,5 2 11 13,2 4 2 11 13,2 4,5 2 11 13,2 5 2 11 13,2 5,5 2 11 13,2 6 2 11 13,2 6,5 2 11 13,2 7 2 11 13,2 7,5 2 11 13,2 8 2 11 13,2 8,5 2 11 13,2 9 2 11 13,2 9,5 2 11 13,2 10 2 11 13,2 10,5 2 11 13,2 11 2 11 13,2 11,5 2 11 13,2 12 2 11 13,2 12,5 2 11 13,2 13 2 11 13,2 13,2 e assim por diante para todos idfustes Desde já agradeço Para quem é da area florestal estou tentando estruturar uma base de cubagem a partir de dados de parcelas Abraço e boa noite a todos Att, Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

6 8

Dúvida em argumentos do histograma
by Victor Eduardo 01 Jun '11

01 Jun '11

Pessoal, estou com um loop e dentro do loop estou querendo utilizar o seguinte para a construção do histograma: hist(paste(var1,'_',var,sep=''),main=paste('tit.',var1,sep=''),xlim=c(0,10),ylim=c(0,85)) sendo que está dando o seguinte erro: Erro em hist.default(paste(var1, "_", var, sep = ""), main = paste("tit.", : 'x' deve ser numérico O primeiro paste me daria, por exemplo, o seguinte: obs_clcl Com isso, construiria o histograma para essa variável. Como eu poderia resolver esse problema? Abraços, Victor Eduardo

2 2

Re: [R-br] graficos com dois eixos y
by Ivan Bezerra Allaman 01 Jun '11

01 Jun '11

Gilson, Repare que no script que te mandei, antes do plot eu utilizo o parâmetro "par(mar=c(0,0,0,0))". Este parâmetro serve para alterar as margens do gráfico. Vá alterando um por um para descobrir a funcionalidade. Teste 5,0,0,0; depois 0,5,0,0 e assim por diante. Com certeza achará um valor adequado para que o nome do eixo y apareça. Pelo o código que te mandei o nome estava aparecendo, mais isso deve depender de máquina pra máquina. As vezes a margem que vc deve utilizar dentro do "mar" deve ser maior do que aqueles que estão no exemplo que te mandei. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] graficos com dois eixos y
by Ivan Bezerra Allaman 01 Jun '11

01 Jun '11

Enquanto o Walmes não responde, vai quebrando o galho com essa solução abaixo. lines <- 'meses temp umidade rad chuva Jan 26.49 86.58 795.88 0.36 Fev 26.65 88.49 710.24 0.34 Mar 27.19 86.16 772.99 0.21 Abr 26.28 89.75 574.88 0.67 Mai 26.62 89.22 614.02 0.31 Jun 26.13 87.83 680.08 0.26 Jul 25.83 86.57 675.97 0.15 Ago 27.05 83.14 756.44 0.07 Set 27.60 83.02 925.57 0.14 Out 27.44 85.16 927.71 0.17 Nov 26.56 88.18 788.87 0.19 Dez 25.87 90.63 703.94 0.33' da <- read.table(textConnection(lines), header=TRUE) par(mar=c(5,4,4,5)) plot(da$temp, type="b", pch=15, lwd=1.5, xlab="Meses (Anos 2010)", ylab="Temperatura (ºC)", ylim=c(24,30), xlim=c(1,12), axes=F, main= "A", font.main=2, font.lab=2) axis(1, at=1:12,labels=c("Jan.","Fev.","Mar.","Abr.","Maio","Jun.", "Jul.", "Ago.", "Set.", "Out.", "Nov.", "Dez.")) axis(2) par(new=T) plot(da$umidade, type="b", pch=14,axes=F, frame=T, ann=F) axis(4) mtext(text="Umidade (%)", 4, line = 3) Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

RStudio
by Diogo Borges Provete 01 Jun '11

01 Jun '11

Olá, alguém que usa o RStudio está conseguindo exibir as mensagens de erro etc em Português? Abraço Atenciosamente, Diogo Borges Provete ============================== Biólogo Mestre em Biologia Animal (UNESP) Doutorando PPG Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás - UFG Goiânia-GO CP: 131 74001-970 Brazil Tel. +55 62 3521-1732 Cel. +55 62 8231-5775 Skype: diogoprovete MSN: diogoprov(a)yahoo.com.br Personal web page Traduza conosco: <<<American Journal Experts>>> <<<D-Lang Soluções linguisticas>>> <<<Perfil no ProZ.com como tradutor freelancer>>> ============================== ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quarta-feira, 1 de Junho de 2011 9:04 Assunto: Digest R-br, volume 4, assunto 1 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: modifica dataframe (Gustavo Henrique de Carvalho) 2. Re: teste de Homogeneidade e Normalidade (Ivan Bezerra Allaman) 3. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 4. Re: teste de Homogeneidade e Normalidade (Mauro Sznelwar) 5. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 6. Re: teste de Homogeneidade e Normalidade (Mauro Sznelwar) 7. Re: teste de Homogeneidade e Normalidade (Walmes Zeviani) 8. Re: modifica dataframe (Henrique Dallazuanna) 9. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 10. Re: modifica dataframe (Samuel Carvalho) 11. Re: teste de Homogeneidade e Normalidade (Walmes Zeviani) 12. Re: teste de Homogeneidade e Normalidade (Eder David Borges da Silva) 13. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 14. Download de dados pelo R (Humberto Hazin (hotmail)) 15. Re: Download de dados pelo R (Victor Eduardo) 16. RES: Download de dados pelo R (Humberto Hazin (hotmail)) ---------------------------------------------------------------------- Message: 1 Date: Tue, 31 May 2011 21:24:23 -0300 From: Gustavo Henrique de Carvalho <gustavo.bio(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Cc: Samuel Carvalho <samukajm(a)yahoo.com.br> Subject: Re: [R-br] modifica dataframe Message-ID: <BANLkTikLCzKcVdHqgWWUSv8f+Z8UMayxKA(a)mail.gmail.com> Content-Type: text/plain; charset=UTF-8 Não sei pra que colar esse negócio grotesco aqui, mas a primeira forma que me veio à cabeça foi essa: bleh <- function(x) {unique(c(seq(0, x, by = 0.5), x))} data.frame(lapply(dados, rep, unlist(lapply(sapply(dados$ht, bleh), length))), hi = unlist(sapply(dados$ht, bleh))) 2011/5/31 Benilton Carvalho <beniltoncarvalho(a)gmail.com> > > enqto o henrique nao responde, algo como o seguinte pode lhe servir: > hi = mapply(seq, 0.1, dados[['ht']], .5) > idx = mapply(rep, 1:length(hi), sapply(hi, length)) > dados2 = cbind(dados[unlist(idx),], hi=unlist(hi)) > b > > 2011/5/31 Samuel Carvalho <samukajm(a)yahoo.com.br> >> >> Boa noite caros(a) senhores(as) membros >> Por gentileza alguém poderia me auxliar em um código que a princípio parece ser simples mas nao consegui uma solução >> Abaixo mostro o codigo e a seguir venho comentando >> ### Codigo R >> dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), >> ht=c(10, 13, 12.5, 18, 14)) >> dados >> dados$hi<-seq(0.1, dados$ht, by=0.5) >> ### >> O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" que corresponde aquele "idfuste", replicando também cada observação >> para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo que a ultima observação de cada idfuste deve ser igual ao valor de ht >> Em resumo queria um data.frame de saída a partir dos dados de entrada, nesta estrutura >> idfuste dap ht hi >> 1 8,5 10 0 >> 1 8,5 10 0,5 >> 1 8,5 10 1 >> 1 8,5 10 1,5 >> 1 8,5 10 2 >> 1 8,5 10 2,5 >> 1 8,5 10 3 >> 1 8,5 10 3,5 >> 1 8,5 10 4 >> 1 8,5 10 4,5 >> 1 8,5 10 5 >> 1 8,5 10 5,5 >> 1 8,5 10 6 >> 1 8,5 10 6,5 >> 1 8,5 10 7 >> 1 8,5 10 7,5 >> 1 8,5 10 8 >> 1 8,5 10 8,5 >> 1 8,5 10 9 >> 1 8,5 10 9,5 >> 1 8,5 10 10 >> 2 11 13,2 0 >> 2 11 13,2 0,5 >> 2 11 13,2 1 >> 2 11 13,2 1,5 >> 2 11 13,2 2 >> 2 11 13,2 2,5 >> 2 11 13,2 3 >> 2 11 13,2 3,5 >> 2 11 13,2 4 >> 2 11 13,2 4,5 >> 2 11 13,2 5 >> 2 11 13,2 5,5 >> 2 11 13,2 6 >> 2 11 13,2 6,5 >> 2 11 13,2 7 >> 2 11 13,2 7,5 >> 2 11 13,2 8 >> 2 11 13,2 8,5 >> 2 11 13,2 9 >> 2 11 13,2 9,5 >> 2 11 13,2 10 >> 2 11 13,2 10,5 >> 2 11 13,2 11 >> 2 11 13,2 11,5 >> 2 11 13,2 12 >> 2 11 13,2 12,5 >> 2 11 13,2 13 >> 2 11 13,2 13,2 >> e assim por diante para todos idfustes >> Desde já agradeço >> Para quem é da area florestal estou tentando estruturar uma base de cubagem a partir de dados de parcelas >> Abraço e boa noite a todos >> Att, >> Samuel >> >> ==================================== >> Samuel P. C. Carvalho >> Mestre em Ciências Florestais [UFLA] >> Doutorando em Recursos Florestais [ESALQ/USP] >> ============================================= >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> > > > > -- > Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ------------------------------ Message: 2 Date: Tue, 31 May 2011 17:34:14 -0700 (PDT) From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> To: R Brasil <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <947015.46346.qm(a)web161820.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Se você quer transformar as variáveis que já estão em porcentagem, deve utilizar arcsen(x%)^0,5. Agora, se as variáveis estiverem entre 0 e 1 (em decimal), aí sim vc deve utilizar arcsen(x/100)^0,5. No entanto, como a variável "germinação" é do tipo Bernouli, aconselho-o a utilizar uma análise mais rebuscada como modelos lineares generalizados por exemplo. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/ef67c401/atta…> ------------------------------ Message: 3 Date: Tue, 31 May 2011 21:46:46 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTinXDFS=T7UD7GrEKSxy-0z1-i77FA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Ok eu entendi as transformações com relação quando já estão em porcentagem.Quando as variáveis por exemplo,não estiverem entre ) e 1, vc pode me informar que transformação eu utilizo?Por que pelo que eu entendi,as transformações com arcsen(x/110)^0.5 só entram nesse intervalo,entre 0 e 1. Em 31 de maio de 2011 21:34, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br>escreveu: > Se você quer transformar as variáveis que já estão em porcentagem, deve > utilizar arcsen(x%)^0,5. Agora, se as variáveis estiverem entre 0 e 1 (em > decimal), aí sim vc deve utilizar arcsen(x/100)^0,5. > > No entanto, como a variável "germinação" é do tipo Bernouli, aconselho-o a > utilizar uma análise mais rebuscada como modelos lineares generalizados por > exemplo. > > Allaman > (S,f,P) > > *M.Sc Ivan Bezerra Allaman* > Zootecnista > Doutorando em Produção Animal/Aquicultura - UFLA > **email e msn - ivanalaman(a)yahoo.com.br > Tel: (35)3826-6608/9900-2924 > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/7c9a93de/atta…> ------------------------------ Message: 4 Date: Tue, 31 May 2011 22:11:34 -0300 From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> To: <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <00f101cc1ff8$ddcd1fe0$6700a8c0@mauro> Content-Type: text/plain; charset="iso-8859-1" Onde eu pego os arquivos para rodar? Mauro José, Sem o seu codigo fica meio complicado, porem neste link tem uma serie de analises de estatistica experimental com onde se testa homogenidade e normalidade. http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R att Em 31 de maio de 2011 20:13, Eduardo Vargas <duvargas(a)gmail.com> escreveu: Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de médias de germinação de sementes, e as varianveis são tratamentos(trat) fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em `[.data.frame`(x, complete.cases(x)) : undefined columns selected Será que alguem pode me ajudar? -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------------------------------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/ae8f9462/atta…> ------------------------------ Message: 5 Date: Tue, 31 May 2011 22:16:25 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTin=XzAn+6AVxUAnEyeNqDCZQ1Ykug(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" São meus resultados de trabalho,preciso fazer essas transformações para normalizar. Em 31 de maio de 2011 22:11, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > Onde eu pego os arquivos para rodar? > Mauro > > José, > Sem o seu codigo fica meio complicado, porem neste link tem uma serie de > analises de estatistica experimental com onde se testa homogenidade e > normalidade. > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > att > > Em 31 de maio de 2011 20:13, Eduardo Vargas <duvargas(a)gmail.com>escreveu: > >> >> Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de >> médias de germinação de sementes, e as varianveis são tratamentos(trat) >> fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de >> variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, >> porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto >> 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em >> `[.data.frame`(x, complete.cases(x)) : >> undefined columns selected >> Será que alguem pode me ajudar? >> -- >> José Eduardo Vargas Lopes de Aráujo >> Engenheiro Agrônomo >> Contato:(38)8811-6985 >> duvargas(a)gmail.com >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> >> > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/74c6398a/atta…> ------------------------------ Message: 6 Date: Tue, 31 May 2011 22:20:04 -0300 From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> To: <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <002201cc1ffa$0decce90$6700a8c0@mauro> Content-Type: text/plain; charset="iso-8859-1" Eu me referi aos arquivos da página dos Exercícios do R no leg, lá tem exemplos do código mas não diz onde estão os arquivos dos dados que foram usados! São meus resultados de trabalho,preciso fazer essas transformações para normalizar. Em 31 de maio de 2011 22:11, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: Onde eu pego os arquivos para rodar? Mauro -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/c7a2ed12/atta…> ------------------------------ Message: 7 Date: Tue, 31 May 2011 22:53:19 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTikem0jdfVDhKZ8N-n8b83QyNkgXvw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Eduardo, Salve engano da minha parte, essa transformação é estabilizadora da variância, ou seja, o foco é corrigir heterocedasticidade e não olha para normalidade dos dados. Uma vez disse em um curso que essas transformações eram coisas do passado e fui por alguns instantes odiado (por aqueles mais tradicionais obviamente). Depois passou. Acontece que essa transformação é a recomendada quando seu dado é do tipo binomial (ou seja, já se sabe que o dado não é normal), dai parte algumas álgebras, integrais, até que se obtém essa função estabilizadora da variância. Transformação estabilizadora da variãncia existe para distribuição binomial e Poisson. Não lembro das funções porque nunca usei. Considero isso metodologia do passo porque, se você sabe que seu dado é binomial (ou Poisson), hoje já existe implementado nos aplicativos estatísticos modelos capazes de considerar essa característica inerente do dado. Esse conjunto de métodos, enfim, se chama modelos lineares generalizados, do qual a distribuição normal é um caso particular (com propriedades bem interessantes). Com essa metodologia você pode obter todos os resultados experimentais que obteria com a distribuição normal, diga-se a dobradinha anova e teste de médias (com algumas adaptações). Então, não há complicações e motivos para não se usar modelos lineares generalizados. Ainda é possível, caso o número de sementes colocadas para germinar seja grande (n, normalmente é 25 ou 50), e caso a probalidade de germinação não esteja na borda (próximo de 0 ou 1), de assumir que o dado é normal e analisa-lo assim. Porém, caso p mude muito entre tratamentos pode as variâncias amostrais serem muito discrepântes. Você pode tentar uma transformação Box-Cox. Verificar os resíduos. Em termos de funções do R, você pode usar a glm(, family=binomial), fazer a análise de resíduos usuais, obter teste (sequencial) para os efeitos fixos via anova() (que não é anova mas quadro de análise de deviance), se for comparar "médias" você pode usar a multcomp::glht() e contrast::contrast(). Verifique a documentação dessas funções para saber como usa-las. Para um começo, você pode rodar o CMR disponível nas Ridículas do LEG. http://www.leg.ufpr.br/doku.php/ridiculas?&#analise_de_dados_de_proporcao_u… À disposição. 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/849b27e7/atta…> ------------------------------ Message: 8 Date: Tue, 31 May 2011 23:10:54 -0300 From: Henrique Dallazuanna <wwwhsd(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Samuel Carvalho <samukajm(a)yahoo.com.br> Subject: Re: [R-br] modifica dataframe Message-ID: <BANLkTi=77BMKw2WfJopo7xjR--he5P03VQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Samuel, Tente assim tbm: x <- sapply(dados$ht, seq, from = 0, by = 0.5) transform(dados[rep(seq(nrow(dados)), times = sapply(x, length)),], hi = unlist(x)) 2011/5/31 Samuel Carvalho <samukajm(a)yahoo.com.br> > Boa noite caros(a) senhores(as) membros > Por gentileza alguém poderia me auxliar em um código que a princípio parece > ser simples mas nao consegui uma solução > Abaixo mostro o codigo e a seguir venho comentando > ### Codigo R > dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), > ht=c(10, 13, 12.5, 18, 14)) > dados > dados$hi<-seq(0.1, dados$ht, by=0.5) > ### > O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" > que corresponde aquele "idfuste", replicando também cada observação > para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da > árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo > que a ultima observação de cada idfuste deve ser igual ao valor de ht > Em resumo queria um data.frame de saída a partir dos dados de entrada, > nesta estrutura > > idfuste dap ht hi 1 8,5 10 0 1 8,5 10 0,5 1 8,5 10 1 1 8,5 10 1,5 1 > 8,5 10 2 1 8,5 10 2,5 1 8,5 10 3 1 8,5 10 3,5 1 8,5 10 4 1 8,5 10 4,5 > 1 8,5 10 5 1 8,5 10 5,5 1 8,5 10 6 1 8,5 10 6,5 1 8,5 10 7 1 8,5 10 > 7,5 1 8,5 10 8 1 8,5 10 8,5 1 8,5 10 9 1 8,5 10 9,5 1 8,5 10 10 2 11 > 13,2 0 2 11 13,2 0,5 2 11 13,2 1 2 11 13,2 1,5 2 11 13,2 2 2 11 13,2 > 2,5 2 11 13,2 3 2 11 13,2 3,5 2 11 13,2 4 2 11 13,2 4,5 2 11 13,2 5 > 2 11 13,2 5,5 2 11 13,2 6 2 11 13,2 6,5 2 11 13,2 7 2 11 13,2 7,5 2 > 11 13,2 8 2 11 13,2 8,5 2 11 13,2 9 2 11 13,2 9,5 2 11 13,2 10 2 11 > 13,2 10,5 2 11 13,2 11 2 11 13,2 11,5 2 11 13,2 12 2 11 13,2 12,5 2 > 11 13,2 13 2 11 13,2 13,2 > > e assim por diante para todos idfustes > Desde já agradeço > Para quem é da area florestal estou tentando estruturar uma base de cubagem > a partir de dados de parcelas > Abraço e boa noite a todos > Att, > Samuel > > *====================================* > *Samuel P. C. Carvalho > *Mestre em Ciências Florestais [UFLA] > Doutorando em Recursos Florestais [ESALQ/USP] > ============================================= > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/4937eeab/atta…> ------------------------------ Message: 9 Date: Tue, 31 May 2011 23:34:58 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTiniir0GPKzWcMm_=PsYwrBfmeJTjA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Walmes,obrigado pela sua resposta,foi bastante pertinente,e acho que para eu fazer esse modelo linear generalizado é bastante complexo.Acho que meu experimento é bastante simples,fiz um experimento com 11 tratamentos + 0 o controle,com 25 sementes cada com 4 repetições.Quero comparar as médias dos tratamentos com a média do controle,Acredito eu,usando Dunnet.Meu orientador pediu que eu fizesse esse teste de normalidade(Shapiro-wilk) e homogeneidade de variâncias(bartlett).As transformações,que acho serem necessarias,eu andei olhando nos artigos que usei como referencia e todos fizeram as transformações.Mais uma vez obrigado pela ajuda,estou meio perdido. Em 31 de maio de 2011 22:53, Walmes Zeviani <walmeszeviani(a)gmail.com>escreveu: > Eduardo, > > Salve engano da minha parte, essa transformação é estabilizadora da > variância, ou seja, o foco é corrigir heterocedasticidade e não olha para > normalidade dos dados. Uma vez disse em um curso que essas transformações > eram coisas do passado e fui por alguns instantes odiado (por aqueles mais > tradicionais obviamente). Depois passou. > > Acontece que essa transformação é a recomendada quando seu dado é do tipo > binomial (ou seja, já se sabe que o dado não é normal), dai parte algumas > álgebras, integrais, até que se obtém essa função estabilizadora da > variância. Transformação estabilizadora da variãncia existe para > distribuição binomial e Poisson. Não lembro das funções porque nunca usei. > > Considero isso metodologia do passo porque, se você sabe que seu dado é > binomial (ou Poisson), hoje já existe implementado nos aplicativos > estatísticos modelos capazes de considerar essa característica inerente do > dado. Esse conjunto de métodos, enfim, se chama modelos lineares > generalizados, do qual a distribuição normal é um caso particular (com > propriedades bem interessantes). > > Com essa metodologia você pode obter todos os resultados experimentais que > obteria com a distribuição normal, diga-se a dobradinha anova e teste de > médias (com algumas adaptações). Então, não há complicações e motivos para > não se usar modelos lineares generalizados. > > Ainda é possível, caso o número de sementes colocadas para germinar seja > grande (n, normalmente é 25 ou 50), e caso a probalidade de germinação não > esteja na borda (próximo de 0 ou 1), de assumir que o dado é normal e > analisa-lo assim. Porém, caso p mude muito entre tratamentos pode as > variâncias amostrais serem muito discrepântes. Você pode tentar uma > transformação Box-Cox. Verificar os resíduos. > > Em termos de funções do R, você pode usar a glm(, family=binomial), fazer a > análise de resíduos usuais, obter teste (sequencial) para os efeitos fixos > via anova() (que não é anova mas quadro de análise de deviance), se for > comparar "médias" você pode usar a multcomp::glht() e contrast::contrast(). > Verifique a documentação dessas funções para saber como usa-las. > > Para um começo, você pode rodar o CMR disponível nas Ridículas do LEG. > > http://www.leg.ufpr.br/doku.php/ridiculas?&#analise_de_dados_de_proporcao_u… > > À disposição. > 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 > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/43f2d612/atta…> ------------------------------ Message: 10 Date: Tue, 31 May 2011 20:02:08 -0700 (PDT) From: Samuel Carvalho <samukajm(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] modifica dataframe Message-ID: <369328.15725.qm(a)web161206.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Obrigado caros pela gentileza de sempre. Vou verificar as sugestões []'s Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] ============================================= ________________________________ De: Henrique Dallazuanna <wwwhsd(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br; Samuel Carvalho <samukajm(a)yahoo.com.br> Enviadas: Terça-feira, 31 de Maio de 2011 23:10 Assunto: Re: [R-br] modifica dataframe Samuel, Tente assim tbm: x <- sapply(dados$ht, seq, from = 0, by = 0.5) transform(dados[rep(seq(nrow(dados)), times = sapply(x, length)),], hi = unlist(x)) 2011/5/31 Samuel Carvalho <samukajm(a)yahoo.com.br> Boa noite caros(a) senhores(as) membros >Por gentileza alguém poderia me auxliar em um código que a princípio parece ser simples mas nao consegui uma solução >Abaixo mostro o codigo e a seguir venho comentando >### Codigo R >dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), >ht=c(10, 13, 12.5, 18, 14)) >dados >dados$hi<-seq(0.1, dados$ht, by=0.5) >### > >O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" que corresponde aquele "idfuste", replicando também cada observação >para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo que a ultima observação de cada idfuste deve ser igual ao valor de ht > >Em resumo queria um data.frame de saída a partir dos dados de entrada, nesta estrutura > > >idfuste dap ht hi >1 8,5 10 0 >1 8,5 10 0,5 >1 8,5 10 1 >1 8,5 10 1,5 >1 8,5 10 2 >1 8,5 10 2,5 >1 8,5 10 3 >1 8,5 10 3,5 >1 8,5 10 4 >1 8,5 10 4,5 >1 8,5 10 5 >1 8,5 10 5,5 >1 8,5 10 6 >1 8,5 10 6,5 >1 8,5 10 7 >1 8,5 10 7,5 >1 8,5 10 8 >1 8,5 10 8,5 >1 8,5 10 9 >1 8,5 10 9,5 >1 8,5 10 10 >2 11 13,2 0 >2 11 13,2 0,5 >2 11 13,2 1 >2 11 13,2 1,5 >2 11 13,2 2 >2 11 13,2 2,5 >2 11 13,2 3 >2 11 13,2 3,5 >2 11 13,2 4 >2 11 13,2 4,5 >2 11 13,2 5 >2 11 13,2 5,5 >2 11 13,2 6 >2 11 13,2 6,5 >2 11 13,2 7 >2 11 13,2 7,5 >2 11 13,2 8 >2 11 13,2 8,5 >2 11 13,2 9 >2 11 13,2 9,5 >2 11 13,2 10 >2 11 13,2 10,5 >2 11 13,2 11 >2 11 13,2 11,5 >2 11 13,2 12 >2 11 13,2 12,5 >2 11 13,2 13 >2 11 13,2 13,2 > >e assim por diante para todos idfustes >Desde já agradeço >Para quem é da area florestal estou tentando estruturar uma base de cubagem a partir de dados de parcelas >Abraço e boa noite a todos >Att, >Samuel > >==================================== >Samuel P. C. Carvalho >Mestre em Ciências Florestais [UFLA] >Doutorando em Recursos Florestais [ESALQ/USP] >============================================= >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/7d0c4086/atta…> ------------------------------ Message: 11 Date: Wed, 1 Jun 2011 00:52:54 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTimC0Ux64PEz2DChLzsedsWigrmoTQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Eduardo, O código abaixo mostra como não é difícil analisar os seus dados pela abordagem que discuti #------------------------------------------------------------------------------------------ # experimento com 11 tratamentos e 1 controle, resposta binomial n=25, 4 repetições trat <- c("controle", paste("trat", 1:11, sep="")) da <- expand.grid(trat=trat, rept=1:4) da$germinadas <- rbinom(nrow(da), size=25, prob=0.8) levels(da$trat) # controle é o nível de referência # ajuste um modelo linear generalizado g0 <- glm(cbind(germinadas, 25-germinadas)~trat, data=da, family=binomial) anova(g0, test="Chisq") # testa hipótese sobre os efeitos principais summary(g0) # cada efeito é o contraste controle vs trat_i # usar correção de bonferroni nas comparações com controle #------------------------------------------------------------------------------------------ Se todo mundo continuar fazendo as coisas como todos já fizeram não vamos evoluir, certo? À disposição. 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/5331e054/atta…> ------------------------------ Message: 12 Date: Wed, 1 Jun 2011 08:13:13 -0300 From: Eder David Borges da Silva <eder(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTin6WSj1N+aP3h=J+Z9caY15_+pBKg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Mauro, a pagina completa com os dados estão em: http://www.leg.ufpr.br/~eder/cursoR/ <http://www.leg.ufpr.br/~eder/cursoR/>e o wiki com as explicações em: http://www.leg.ufpr.br/doku.php/cursos:ragronomia <http://www.leg.ufpr.br/doku.php/cursos:ragronomia>Att Em 31 de maio de 2011 22:11, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > Onde eu pego os arquivos para rodar? > Mauro > > José, > Sem o seu codigo fica meio complicado, porem neste link tem uma serie de > analises de estatistica experimental com onde se testa homogenidade e > normalidade. > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > att > > Em 31 de maio de 2011 20:13, Eduardo Vargas <duvargas(a)gmail.com>escreveu: > >> >> Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de >> médias de germinação de sementes, e as varianveis são tratamentos(trat) >> fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de >> variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, >> porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto >> 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em >> `[.data.frame`(x, complete.cases(x)) : >> undefined columns selected >> Será que alguem pode me ajudar? >> -- >> José Eduardo Vargas Lopes de Aráujo >> Engenheiro Agrônomo >> Contato:(38)8811-6985 >> duvargas(a)gmail.com >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> >> > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/8b5cb003/atta…> ------------------------------ Message: 13 Date: Wed, 1 Jun 2011 08:37:36 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTincK6wOSuWSgLyYEv_A4wmBXQ+q5A(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Obrigado mais uma vez pela informação. Em 1 de junho de 2011 00:52, Walmes Zeviani <walmeszeviani(a)gmail.com>escreveu: > Eduardo, > > O código abaixo mostra como não é difícil analisar os seus dados pela > abordagem que discuti > > > #------------------------------------------------------------------------------------------ > # experimento com 11 tratamentos e 1 controle, resposta binomial n=25, 4 > repetições > > trat <- c("controle", paste("trat", 1:11, sep="")) > da <- expand.grid(trat=trat, rept=1:4) > da$germinadas <- rbinom(nrow(da), size=25, prob=0.8) > levels(da$trat) # controle é o nível de referência > > # ajuste um modelo linear generalizado > g0 <- glm(cbind(germinadas, 25-germinadas)~trat, data=da, family=binomial) > anova(g0, test="Chisq") # testa hipótese sobre os efeitos principais > summary(g0) # cada efeito é o contraste controle vs trat_i > # usar correção de bonferroni nas comparações com controle > > > #------------------------------------------------------------------------------------------ > > Se todo mundo continuar fazendo as coisas como todos já fizeram não vamos > evoluir, certo? > > À disposição. > 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 > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/833c1853/atta…> ------------------------------ Message: 14 Date: Wed, 1 Jun 2011 08:54:31 -0300 From: "Humberto Hazin (hotmail)" <hghazin(a)hotmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Download de dados pelo R Message-ID: <COL108-DS18964BBC6881C7A035B8CEB47D0(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-1" Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/7cafde66/atta…> ------------------------------ Message: 15 Date: Wed, 1 Jun 2011 08:58:14 -0300 From: Victor Eduardo <victorduca08(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Download de dados pelo R Message-ID: <BANLkTimuvALu_kvZ3iFTzho3P4WCd+bfpQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Humbero, não está dando para acessar o link enviado. Em 1 de junho de 2011 08:54, Humberto Hazin (hotmail) <hghazin(a)hotmail.com>escreveu: > Bom dia pessoal! > > > > > > Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou > conseguindo, alguém poderia me ajudar? > > > > ano<-2001 > > cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" > > > > download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em > R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), > > "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", > method="auto") > > > > Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em > R/ambientais_", : > > cannot open destfile > 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais > ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such > file or directory' > > > > > > Abs, > > > > Humberto > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/d4a14280/atta…> ------------------------------ Message: 16 Date: Wed, 1 Jun 2011 09:04:35 -0300 From: "Humberto Hazin (hotmail)" <hghazin(a)hotmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] RES: Download de dados pelo R Message-ID: <COL108-DS15CEF73034BD68424BCDFFB47D0(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-1" Olá Eduardo Ele deve estar fora de ar nesse instante ! Tenta com esse aqui http://aspera.jpl.nasa.gov/download/pub/sea_surface_temperature/avhrr/pathfi nder/data_v5/monthly/night/04km/1986/bsst/ ou http://oceandata.sci.gsfc.nasa.gov/SeaWiFS/Mapped/Monthly/9km/chlor/ Abs, Humberto De: r-br-bounces(a)listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Victor Eduardo Enviada em: Wednesday, June 01, 2011 8:58 AM Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Download de dados pelo R Humbero, não está dando para acessar o link enviado. Em 1 de junho de 2011 08:54, Humberto Hazin (hotmail) <hghazin(a)hotmail.com> escreveu: Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/18aad303/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 4, assunto 1 ***************************************

1 0

Ordenar gráfico bipartido
by Jéssica Caroline Falcão 01 Jun '11

01 Jun '11

Olá, Eu gostaria de saber como eu ordeno o gráfico bipartido gerado pelo pacote bipartite de acordo com o número de interações das espécies na rede, em ordem decrescente. Att. Jéssica Falcão -- ------------------------------------------------------------------ Jéssica Falcão Mestranda em Ecologia e Conservação da Biodiversidade Universidade Federal de Mato Grosso IB- Departamento de Ecologia e Botânica Lab. Interação Inseto-Planta / Mirmecologia Cuiabá, Mato Grosso, BRASIL Email: jecafalcao(a)gmail.com

2 1

Re: [R-br] teste de Homogeneidade e Normalidade
by Ivan Bezerra Allaman 31 Mai '11

31 Mai '11

Se você quer transformar as variáveis que já estão em porcentagem, deve utilizar arcsen(x%)^0,5. Agora, se as variáveis estiverem entre 0 e 1 (em decimal), aí sim vc deve utilizar arcsen(x/100)^0,5. No entanto, como a variável "germinação" é do tipo Bernouli, aconselho-o a utilizar uma análise mais rebuscada como modelos lineares generalizados por exemplo. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

R e LaTex
by Cassiano 31 Mai '11

31 Mai '11

Olá a todos, Gostaria de aproveitar a lista para difulgar um trabalho sobre o Sweave - uma interface entre o R e o LaTex. O trabalho está publicado na The PracTex Journal 2010-2 e pode ser acessado no endereço: http://tug.org/pracjourn/2010-2/toc.html Att, -- ======================================= MSc Cassiano S. Rosa Universidade Federal de Viçosa -UFV Dept. Biologia Animal Pós-graduação em Entomologia Lab. Termitologia Viçosa-MG, Brasil CEP 36570-000 Tel: (31) 38994017 =======================================

2 1

Distribuição Exponencial
by Fernanda Balbino 31 Mai '11

31 Mai '11

Alexandro, Muito obrigada, os comandos dão o resultado que eu precisava sobre a distribuição exponencial. Fernanda. ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Terça-feira, 31 de Maio de 2011 12:00:02 Assunto: Digest R-br, volume 3, assunto 36 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. [Dúvida] Plotando pontos centroides no MAPA (Pedro Rafael) 2. Re: Diagrama de voronoi (Wagner bonat) 3. Re: Análise fatorial e regressão em painel (Felipe Emanoel Barletta Mendes) 4. Re: Diagrama de voronoi (Fabrício Tavares) 5. Distribuição Exponencial (Fernanda Balbino) 6. Re: Distribuição Exponencial (Fernando Mayer) 7. Re: Distribuição Exponencial (Alexandro (Yahoo)) 8. nao consigo instalar o pacote "rimage" - Ubuntu (Guilherme Moraes Ferraudo) 9. Re: nao consigo instalar o pacote "rimage" - Ubuntu (Daniel C Bezerra) 10. Re: nao consigo instalar o pacote "rimage" - Ubuntu (Benilton Carvalho) 11. Livro (Peterson Gardin) ---------------------------------------------------------------------- Message: 1 Date: Mon, 30 May 2011 12:00:20 -0300 From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] [Dúvida] Plotando pontos centroides no MAPA Message-ID: <BANLkTik2Nu2xDKh2BqDKUdYX6pa=j_W42Q(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Pessoal como faço para plotar no map as coordenadas dos centróides dos municípios. Importei o mapa com os comandos que seguem require(SpatialEpi) require(maptools) dados = read.table("C:\\Users\\GORR-PEDRO\\Desktop\\Malha Municipal PB\\dados.txt", header = TRUE) attach(dados) mapa = readShapePoly("C:\\Users\\GORR-PEDRO\\Desktop\\Malha Municipal PB\\PB.shp") geo <- coordinates(mapa) geo <- latlong2grid(geo) Gostaria de plotar no mapa o objeto "geo" Pedro -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/1eef92f9/atta…> ------------------------------ Message: 2 Date: Mon, 30 May 2011 12:18:55 -0300 From: Wagner bonat <wbonat(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Diagrama de voronoi Message-ID: <BANLkTinCo6d5NgvxV3VhbVQn+y7PMFNOBg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Veja este exemplo, fiz isto a algum tempo deve estar desatualizado require(gpclib) require(tripack) require(sp) require(spdep) ## Lendo o arquivo das bordas bordas <- read.table(" http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/points_higi.csv ",header=TRUE,sep=";",dec=".") bordas = bordas[,-1] bordas[44,] <- bordas[1,] ## Lendo o arquivo que contem as coordenadas mosq = read.table(" http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/mosqhigi.csv ",header=TRUE,sep=";",dec=".") ## Funcoes extras source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/voronoi.R") source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/grp2sp.r") coords.x <- c() coords.y <- c() cod.mosq <- unique(mosq$mosquitrap) for(i in 1:40){ coords.x[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$x) coords.y[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$y)} mapa.mosq <- voronoi(coords.x1=coords.x,coords.x2=coords.y,poligono = bordas) mapa.mosq = grp2sp(mapa.mosq,ID = cod.mosq) mapa.mosq.nb <- poly2nb(mapa.mosq) plot(mapa.mosq) -- Wagner Hugo Bonat LEG - Laboratório de Estatística e Geoinformação UFPR - Universidade Federal do Paraná -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/13fe55d8/atta…> ------------------------------ Message: 3 Date: Mon, 30 May 2011 13:13:07 -0300 From: "Felipe Emanoel Barletta Mendes" <felipe(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Análise fatorial e regressão em painel Message-ID: <500a2d6b5bb74891c63e982af2c8d472.squirrel(a)leg.ufpr.br> Content-Type: text/plain;charset=iso-8859-1 Henrique, há o pacote plm para dados em painel. ?plm Henrique Brilhante > Bom dia senhores, já acompanho a lista a algum tempo pois estou > desenvolvendo um trabalho na graduação sobre o modelo APT. > Dado meu conhecimento bem limitado do software gostaria de saber se alguém > pode me auxiliar com o problema abaixo de forma passo-a-passo, mesmo que > fora desta lista de discussão. > > VarAcao = PIBreal + Ptax + Ouro + IBOVESPA + Dow Jones + Nasdaq + > Crise2008 (dummy) > > O método seria realizar a análise fatorial dos 6 fatores e salvar os betas > para a regressão em painel considerando o efeito da crise de 2008. > > Abaixo um exemplo dos dados (3 das 40 empresas do estudo). > > Data Acao VarAcao PIBreal Ptax Ouro IBOVESPA Dow > Jones Nasdaq Crise2008 > 1 T 2005 VALE - PNA 0,0949 -0,0623 -0,0552 -0,0757 > -0,0439 -0,0838 -0,1420 0,0000 > 2 T 2005 VALE - PNA -0,1857 -0,0352 -0,1851 -0,1677 > -0,1194 -0,0739 -0,0234 0,0000 > 3 T 2005 VALE - PNA 0,3734 -0,0651 -0,1112 -0,0357 > 0,1766 -0,0411 -0,0243 0,0000 > 4 T 2005 VALE - PNA -0,0332 -0,0499 -0,0078 0,0745 > -0,0022 -0,0191 -0,0085 0,0000 > 1 T 2006 VALE - PNA 0,0960 -0,0374 -0,1294 -0,0117 > 0,0713 -0,0197 0,0036 0,0000 > 2 T 2006 VALE - PNA -0,0473 -0,0339 -0,0406 0,0280 > -0,0723 -0,0476 -0,1256 0,0000 > 3 T 2006 VALE - PNA -0,0937 -0,0238 -0,0352 -0,0752 > -0,0447 0,0111 0,0037 0,0000 > 4 T 2006 VALE - PNA 0,2958 -0,0309 -0,0593 0,0045 > 0,1564 0,0391 0,0412 0,0000 > 1 T 2007 VALE - PNA 0,1618 -0,0235 -0,0847 -0,0134 > -0,0134 -0,0566 -0,0452 0,0000 > 2 T 2007 VALE - PNA 0,1198 -0,0211 -0,0996 -0,0899 > 0,1347 0,0383 0,0287 0,0000 > 3 T 2007 VALE - PNA 0,3559 -0,0257 -0,0832 -0,0022 > 0,0690 0,0071 0,0084 0,0000 > 4 T 2007 VALE - PNA -0,0349 -0,0232 -0,0780 0,0547 > 0,0145 -0,0803 -0,0522 0,0000 > 1 T 2008 VALE - PNA -0,0143 -0,0237 -0,0537 0,0210 > -0,0878 -0,1210 -0,1941 0,0000 > 2 T 2008 VALE - PNA -0,0749 -0,0307 -0,1425 -0,1398 > 0,0159 -0,1294 -0,0460 0,0000 > 3 T 2008 VALE - PNA -0,3879 -0,0273 0,1415 0,0941 > -0,3148 -0,0771 -0,1285 1 > 4 T 2008 VALE - PNA -0,2963 -0,0874 0,1551 0,1265 > -0,3216 -0,2058 -0,2716 1 > 1 T 2009 VALE - PNA 0,1008 -0,0607 -0,0507 0,0053 > 0,0448 -0,1835 -0,0719 1 > 2 T 2009 VALE - PNA 0,1140 -0,0184 -0,2078 -0,1473 > 0,1922 0,0666 0,1449 1 > 3 T 2009 VALE - PNA 0,1976 -0,0028 -0,1213 -0,0717 > 0,1502 0,1164 0,1223 0,0000 > 4 T 2009 VALE - PNA 0,1440 -0,0073 -0,0525 0,0447 > 0,0773 0,0502 0,0460 0,0000 > 1 T 2010 VALE - PNA 0,1401 -0,0191 -0,0181 0,0517 > -0,0150 0,0122 0,0272 0,0000 > 2 T 2010 VALE - PNA -0,2695 -0,0164 -0,0208 0,0537 > -0,1762 -0,1288 -0,1521 0,0000 > 3 T 2010 VALE - PNA 0,1949 -0,0272 -0,0926 -0,0285 > 0,0993 0,0703 0,0876 0,0000 > 4 T 2010 VALE - PNA 0,0360 -0,0402 -0,0643 0,0509 > -0,0495 0,0416 0,0843 0,0000 > 1 T 2005 PETROBRAS PN 0,0422 -0,0623 -0,0552 -0,0757 > -0,0439 -0,0838 -0,1420 0,0000 > 2 T 2005 PETROBRAS PN 0,0424 -0,0352 -0,1851 -0,1677 > -0,1194 -0,0739 -0,0234 0,0000 > 3 T 2005 PETROBRAS PN 0,2955 -0,0651 -0,1112 -0,0357 > 0,1766 -0,0411 -0,0243 0,0000 > 4 T 2005 PETROBRAS PN 0,0263 -0,0499 -0,0078 0,0745 > -0,0022 -0,0191 -0,0085 0,0000 > 1 T 2006 PETROBRAS PN 0,1471 -0,0374 -0,1294 -0,0117 > 0,0713 -0,0197 0,0036 0,0000 > 2 T 2006 PETROBRAS PN 0,0141 -0,0339 -0,0406 0,0280 > -0,0723 -0,0476 -0,1256 0,0000 > 3 T 2006 PETROBRAS PN -0,0664 -0,0238 -0,0352 -0,0752 > -0,0447 0,0111 0,0037 0,0000 > 4 T 2006 PETROBRAS PN 0,2168 -0,0309 -0,0593 0,0045 > 0,1564 0,0391 0,0412 0,0000 > 1 T 2007 PETROBRAS PN -0,0797 -0,0235 -0,0847 -0,0134 > -0,0134 -0,0566 -0,0452 0,0000 > 2 T 2007 PETROBRAS PN 0,1124 -0,0211 -0,0996 -0,0899 > 0,1347 0,0383 0,0287 0,0000 > 3 T 2007 PETROBRAS PN 0,1386 -0,0257 -0,0832 -0,0022 > 0,0690 0,0071 0,0084 0,0000 > 4 T 2007 PETROBRAS PN 0,3949 -0,0232 -0,0780 0,0547 > 0,0145 -0,0803 -0,0522 0,0000 > 1 T 2008 PETROBRAS PN -0,1893 -0,0237 -0,0537 0,0210 > -0,0878 -0,1210 -0,1941 0,0000 > 2 T 2008 PETROBRAS PN 0,2043 -0,0307 -0,1425 -0,1398 > 0,0159 -0,1294 -0,0460 0,0000 > 3 T 2008 PETROBRAS PN -0,2857 -0,0273 0,1415 0,0941 > -0,3148 -0,0771 -0,1285 1 > 4 T 2008 PETROBRAS PN -0,4048 -0,0874 0,1551 0,1265 > -0,3216 -0,2058 -0,2716 1 > 1 T 2009 PETROBRAS PN 0,2109 -0,0607 -0,0507 0,0053 > 0,0448 -0,1835 -0,0719 1 > 2 T 2009 PETROBRAS PN 0,1264 -0,0184 -0,2078 -0,1473 > 0,1922 0,0666 0,1449 1 > 3 T 2009 PETROBRAS PN 0,0791 -0,0028 -0,1213 -0,0717 > 0,1502 0,1164 0,1223 0,0000 > 4 T 2009 PETROBRAS PN 0,0478 -0,0073 -0,0525 0,0447 > 0,0773 0,0502 0,0460 0,0000 > 1 T 2010 PETROBRAS PN -0,0565 -0,0191 -0,0181 0,0517 > -0,0150 0,0122 0,0272 0,0000 > 2 T 2010 PETROBRAS PN -0,2710 -0,0164 -0,0208 0,0537 > -0,1762 -0,1288 -0,1521 0,0000 > 3 T 2010 PETROBRAS PN 0,0181 -0,0272 -0,0926 -0,0285 > 0,0993 0,0703 0,0876 0,0000 > 4 T 2010 PETROBRAS PN -0,0089 -0,0402 -0,0643 0,0509 > -0,0495 0,0416 0,0843 0,0000 > 1 T 2005 ITAUUNIBANCO PN 0,0771 -0,0623 -0,0552 -0,0757 > -0,0439 -0,0838 -0,1420 0,0000 > 2 T 2005 ITAUUNIBANCO PN -0,0190 -0,0352 -0,1851 -0,1677 > -0,1194 -0,0739 -0,0234 0,0000 > 3 T 2005 ITAUUNIBANCO PN 0,2088 -0,0651 -0,1112 -0,0357 > 0,1766 -0,0411 -0,0243 0,0000 > 4 T 2005 ITAUUNIBANCO PN 0,0466 -0,0499 -0,0078 0,0745 > -0,0022 -0,0191 -0,0085 0,0000 > 1 T 2006 ITAUUNIBANCO PN 0,1268 -0,0374 -0,1294 -0,0117 > 0,0713 -0,0197 0,0036 0,0000 > 2 T 2006 ITAUUNIBANCO PN -0,0140 -0,0339 -0,0406 0,0280 > -0,0723 -0,0476 -0,1256 0,0000 > 3 T 2006 ITAUUNIBANCO PN 0,0326 -0,0238 -0,0352 -0,0752 > -0,0447 0,0111 0,0037 0,0000 > 4 T 2006 ITAUUNIBANCO PN 0,1629 -0,0309 -0,0593 0,0045 > 0,1564 0,0391 0,0412 0,0000 > 1 T 2007 ITAUUNIBANCO PN -0,0672 -0,0235 -0,0847 -0,0134 > -0,0134 -0,0566 -0,0452 0,0000 > 2 T 2007 ITAUUNIBANCO PN 0,1719 -0,0211 -0,0996 -0,0899 > 0,1347 0,0383 0,0287 0,0000 > 3 T 2007 ITAUUNIBANCO PN 0,0653 -0,0257 -0,0832 -0,0022 > 0,0690 0,0071 0,0084 0,0000 > 4 T 2007 ITAUUNIBANCO PN -0,0243 -0,0232 -0,0780 0,0547 > 0,0145 -0,0803 -0,0522 0,0000 > 1 T 2008 ITAUUNIBANCO PN -0,1238 -0,0237 -0,0537 0,0210 > -0,0878 -0,1210 -0,1941 0,0000 > 2 T 2008 ITAUUNIBANCO PN -0,0016 -0,0307 -0,1425 -0,1398 > 0,0159 -0,1294 -0,0460 0,0000 > 3 T 2008 ITAUUNIBANCO PN -0,0265 -0,0273 0,1415 0,0941 > -0,3148 -0,0771 -0,1285 1 > 4 T 2008 ITAUUNIBANCO PN -0,1913 -0,0874 0,1551 0,1265 > -0,3216 -0,2058 -0,2716 1 > 1 T 2009 ITAUUNIBANCO PN -0,0153 -0,0607 -0,0507 0,0053 > 0,0448 -0,1835 -0,0719 1 > 2 T 2009 ITAUUNIBANCO PN 0,1808 -0,0184 -0,2078 -0,1473 > 0,1922 0,0666 0,1449 1 > 3 T 2009 ITAUUNIBANCO PN 0,2339 -0,0028 -0,1213 -0,0717 > 0,1502 0,1164 0,1223 0,0000 > 4 T 2009 ITAUUNIBANCO PN 0,0709 -0,0073 -0,0525 0,0447 > 0,0773 0,0502 0,0460 0,0000 > 1 T 2010 ITAUUNIBANCO PN 0,0035 -0,0191 -0,0181 0,0517 > -0,0150 0,0122 0,0272 0,0000 > 2 T 2010 ITAUUNIBANCO PN -0,1913 -0,0164 -0,0208 0,0537 > -0,1762 -0,1288 -0,1521 0,0000 > 3 T 2010 ITAUUNIBANCO PN 0,2218 -0,0272 -0,0926 -0,0285 > 0,0993 0,0703 0,0876 0,0000 > 4 T 2010 ITAUUNIBANCO PN -0,0381 -0,0402 -0,0643 0,0509 > -0,0495 0,0416 0,0843 0,0000 > > > Agradeço a gentileza, > > Muito obrigado. > Henrique. > > brilhante.henrique(a)gmail.com > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Felipe E. Barletta Mendes (41)9189-5198 (41)3025-2150 (41)3328-7216 http://www.leg.ufpr.br/doku.php/pessoais:felipe ------------------------------ Message: 4 Date: Mon, 30 May 2011 17:28:00 -0300 From: Fabrício Tavares <fabriciotavareso(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Diagrama de voronoi Message-ID: <BANLkTin-nDpN30Z1vbkjZpW7VDhFfSk+aQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Wagner, funcionou perfeitamente! Muito obrigado! Fabrício. http://vizzual.net Em 30 de maio de 2011 12:18, Wagner bonat <wbonat(a)gmail.com> escreveu: > Veja este exemplo, fiz isto a algum tempo deve estar desatualizado > > require(gpclib) > require(tripack) > require(sp) > require(spdep) > ## Lendo o arquivo das bordas > bordas <- read.table(" > http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/points_higi.csv > ",header=TRUE,sep=";",dec=".") > bordas = bordas[,-1] > bordas[44,] <- bordas[1,] > > ## Lendo o arquivo que contem as coordenadas > mosq = read.table(" > http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/mosqhigi.csv > ",header=TRUE,sep=";",dec=".") > > ## Funcoes extras > source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/voronoi.R") > source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/grp2sp.r") > > coords.x <- c() > coords.y <- c() > cod.mosq <- unique(mosq$mosquitrap) > for(i in 1:40){ > coords.x[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$x) > coords.y[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$y)} > > mapa.mosq <- voronoi(coords.x1=coords.x,coords.x2=coords.y,poligono = > bordas) > mapa.mosq = grp2sp(mapa.mosq,ID = cod.mosq) > mapa.mosq.nb <- poly2nb(mapa.mosq) > plot(mapa.mosq) > > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/644eef83/atta…> ------------------------------ Message: 5 Date: Mon, 30 May 2011 13:51:16 -0700 (PDT) From: Fernanda Balbino <ferbalbino(a)yahoo.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Distribuição Exponencial Message-ID: <25518.42840.qm(a)web161501.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Olá, boa tarde a todos, Tenho um conjunto de dados que devido a natureza e aos gráficos de análise exploratória acredito que sigam uma distribuição exponencial. Gostaria de saber se no R é possível testar os dados para confirmar se realmente eles se comportam conforme uma distribuição exponencial. Muito obrigada. Fernanda Balbino. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/21cbd860/atta…> ------------------------------ Message: 6 Date: Mon, 30 May 2011 18:28:14 -0300 From: Fernando Mayer <fernandomayer(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Distribuição Exponencial Message-ID: <BANLkTi=ofAuy684SbjuU1doRt1z9nyY3Hw(a)mail.gmail.com> Content-Type: text/plain; charset=ISO-8859-1 Veja a função fitdistr() do pacote MASS [], --- Fernando Mayer URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com 2011/5/30 Fernanda Balbino <ferbalbino(a)yahoo.com.br>: > Olá, boa tarde a todos, > Tenho um conjunto de dados que devido a natureza e aos gráficos de análise > exploratória acredito que sigam uma distribuição exponencial. > Gostaria de saber se no R é possível testar os dados para confirmar se > realmente eles se comportam conforme uma distribuição exponencial. > Muito obrigada. > Fernanda Balbino. > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > ------------------------------ Message: 7 Date: Mon, 30 May 2011 16:11:06 -0700 (PDT) From: "Alexandro $Yahoo$" <alexandro2x(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Distribuição Exponencial Message-ID: <499551.91899.qm(a)web130108.mail.mud.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Boa noite Fernanda, seja x, a variável que tenha os seus dados, tu podes testar o seguinte: #EXEMPLO DE NÃO TER DIST. EXP #x=rnorm(1000,10,2) #ks.test(x, pexp, rate=5) #pvalor=ks.test(x, pexp, rate=5)$p.value #EXEMPLO DE DIST. EXP x <- rexp(30, rate=5) ks.test(x, pexp, rate=5) pvalor=ks.test(x, pexp, rate=5)$p.value pvalor alfa=0.05 teste<-function(pvalor) { if(pvalor<alfa) { cat("rejeita Ho >>> pvalor é ",pvalor, " NÃO TEM DIST. EXP ","\n") } else if(pvalor>=alfa) { cat("NÃO rejeita Ho >>> pvalor é ",pvalor, " TEM DIST. EXP. ","\n") } } teste(pvalor) Espero ter ajudado ! Alexandro >________________________________ >De: Fernanda Balbino <ferbalbino(a)yahoo.com.br> >Para: r-br(a)listas.c3sl.ufpr.br >Enviadas: Segunda-feira, 30 de Maio de 2011 17:51 >Assunto: [R-br] Distribuição Exponencial > > >Olá, boa tarde a todos, > > >Tenho um conjunto de dados que devido a natureza e aos gráficos de análise >exploratória acredito que sigam uma distribuição exponencial. >Gostaria de saber se no R é possível testar os dados para confirmar se realmente >eles se comportam conforme uma distribuição exponencial. > > >Muito obrigada. >Fernanda Balbino. > > > > >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/c7b1fdad/atta…> ------------------------------ Message: 8 Date: Mon, 30 May 2011 21:58:02 -0300 From: Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] nao consigo instalar o pacote "rimage" - Ubuntu Message-ID: <BANLkTim9AKvn2zsX_t=oeBffhX=0k+beKQ(a)mail.gmail.com> Content-Type: text/plain; charset=windows-1252 Boa noite pessoal, estou no Ubuntu e não consigo instalar o pacote "rimage" - Ubuntu. Olha o erro: ERROR: configuration failed for package ?rimage? * removing ?/home/gferrau/R/i486-pc-linux-gnu-library/2.10/rimage? The downloaded packages are in ?/tmp/RtmpcEgQZS/downloaded_packages? Warning message: In install.packages("rimage", dep = T) : installation of package 'rimage' had non-zero exit status O que fazer?? Abs -- Guilherme Moraes Ferraudo Jaboticabal/Campinas - SP http://lattes.cnpq.br/2096118558794430 ------------------------------ Message: 9 Date: Mon, 30 May 2011 22:10:42 -0300 From: Daniel C Bezerra <danielcbezerra(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] nao consigo instalar o pacote "rimage" - Ubuntu Message-ID: <BANLkTikZ=cNKNEQN0AH_0WbXXCQgbWd1PA(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Parece que vc está usando a versão 2.10 do R. Tente atualizar primeiro. Abs, D 2011/5/30 Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> > Boa noite pessoal, > > estou no Ubuntu e não consigo instalar o pacote "rimage" - Ubuntu. > > Olha o erro: > > ERROR: configuration failed for package ?rimage? > * removing ?/home/gferrau/R/i486-pc-linux-gnu-library/2.10/rimage? > > The downloaded packages are in > ?/tmp/RtmpcEgQZS/downloaded_packages? > Warning message: > In install.packages("rimage", dep = T) : > installation of package 'rimage' had non-zero exit status > > > O que fazer?? > > Abs > > > -- > Guilherme Moraes Ferraudo > Jaboticabal/Campinas - SP > http://lattes.cnpq.br/2096118558794430 > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/51828ce7/atta…> ------------------------------ Message: 10 Date: Tue, 31 May 2011 10:43:51 +0100 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] nao consigo instalar o pacote "rimage" - Ubuntu Message-ID: <BANLkTinyw7kdkhY-6ZSFNuJ8jk7cndNfLw(a)mail.gmail.com> Content-Type: text/plain; charset=windows-1252 Adicionalmente, a mensagem de erro parece incompleta. 2011/5/31 Daniel C Bezerra <danielcbezerra(a)gmail.com>: > Parece que vc está usando a versão 2.10 do R. Tente atualizar primeiro. > > Abs, > > D > > 2011/5/30 Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> >> >> Boa noite pessoal, >> >> estou no Ubuntu e não consigo instalar o pacote "rimage" - Ubuntu. >> >> Olha o erro: >> >> ERROR: configuration failed for package ?rimage? >> * removing ?/home/gferrau/R/i486-pc-linux-gnu-library/2.10/rimage? >> >> The downloaded packages are in >> ?/tmp/RtmpcEgQZS/downloaded_packages? >> Warning message: >> In install.packages("rimage", dep = T) : >> installation of package 'rimage' had non-zero exit status >> >> >> O que fazer?? >> >> Abs >> >> >> -- >> Guilherme Moraes Ferraudo >> Jaboticabal/Campinas - SP >> http://lattes.cnpq.br/2096118558794430 >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) ------------------------------ Message: 11 Date: Tue, 31 May 2011 07:48:21 -0700 (PDT) From: Peterson Gardin <petersongardin(a)yahoo.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Livro Message-ID: <596446.76610.qm(a)web125801.mail.ne1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Precisava deste livro, alguém tem ? "R Through Excel" Obrigado João Peterson Pereira Gardin Pesquisador Epagri - Videira - SC Doutor em Fisiologia Vegetal Engenheiro Agrônomo Fones: (49) 3566-0054 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/9f3ca333/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 3, assunto 36 ****************************************

1 0

Livro
by Peterson Gardin 31 Mai '11

31 Mai '11

Precisava deste livro, alguém tem ? "R Through Excel" Obrigado João Peterson Pereira Gardin Pesquisador Epagri - Videira - SC Doutor em Fisiologia Vegetal Engenheiro Agrônomo Fones: (49) 3566-0054

4 3

Alterar a função codebook
by Daniel Marcelino 31 Mai '11

31 Mai '11

Caros, esta função codebook tem uma saída para valores characters e factors diferente do que eu preciso para montar um livro de códigos. No código abaixo, quando a coluna é string a única informação que mostra é: "A character vector". O que eu preciso é que a saída me mostre pelo menos uns 3 exemplos de strings. Alguém sabe como alterar isso? Daniel function (dataFrame = .data) { cat("\n", attr(dataFrame, "datalabel"), "\n", "\n") x1 <- dataFrame[1, ] for (i in 1:ncol(dataFrame)) { cat(paste(names(dataFrame)[i], "\t", ":", "\t", attr(dataFrame, "var.labels")[i]), "\n") if (all(is.na(dataFrame[, i]))) { cat(paste("All elements of ", names(dataFrame)[i], " have a missing value", "\n")) } else { if (any(class(x1) == "data.frame")) { x2 <- x1[, i] } else { x2 <- x1 } if (any(class(x2) == "character") | any(class(x2) == "AsIs")) { cat("A character vector", "\n") } else { if (any(class(x2) == "difftime")) { print(summary(x2)) } else { if (is.logical(x2)) x2 <- as.factor(x2) if (any(class(x2) == "factor")) { table1 <- (t(t(table(dataFrame[, i])))) table1 <- cbind(table1, format(table1/sum(table1) * 100, digits = 3)) colnames(table1) <- c(.frequency1, .percent) if (is.null(attr(dataFrame, "val.labels")[i])) { print.noquote(table1, right = TRUE) } else { if (any(is.na(attr(dataFrame, "label.table")))) { print.noquote(table1, right = TRUE) } else { attr(dataFrame, "label.table")[which(is.na(attr(attr(dataFrame, "label.table"), "names")))] <- "" index <- attr(attr(dataFrame, "label.table"), "names") == attr(dataFrame, "val.labels")[i] index <- na.omit(index) if (suppressWarnings(!all(rownames(as.data.frame(attr(dataFrame, "label.table")[index])) == levels(x2)))) { print.noquote(table1, right = TRUE) } else { table2 <- data.frame(attr(dataFrame, "label.table")[index], table1) colnames(table2) <- c("code", colnames(table1)) cat("Label table:", attr(dataFrame, "val.labels")[i], "\n") print.noquote(table2, right = TRUE) } } } } else { print(summ(dataFrame[, i], graph = FALSE)) } } } } cat("\n", "==================", "\n") } } -- Daniel Marcelino http://danielmarcelino.zip.net Skype: d_marcelino

1 0

RCurl - Ubuntu 11.04
by Fabrício Tavares 31 Mai '11

31 Mai '11

Boa tarde pessoal. Estou com dificuldades pra instalar o pacote RCurl no Ubuntu 11.04 (uso há pouco tempo). > sessionInfo() R version 2.13.0 (2011-04-13) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=pt_BR.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tcltk_2.13.0 tools_2.13.0 > install.packages('RCurl') Installing package(s) into ‘/home/fabricio/R/x86_64-pc-linux-gnu-library/2.13’ (as ‘lib’ is unspecified) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done trying URL 'http://cran-r.c3sl.ufpr.br/src/contrib/RCurl_1.6-4.tar.gz' Content type 'application/x-gzip' length 703548 bytes (687 Kb) opened URL ================================================== downloaded 687 Kb * installing *source* package ‘RCurl’ ... checking for curl-config... no Cannot find curl-config ERROR: configuration failed for package ‘RCurl’ * removing ‘/home/fabricio/R/x86_64-pc-linux-gnu-library/2.13/RCurl’ The downloaded packages are in ‘/tmp/Rtmp8Pq7iI/downloaded_packages’ Warning message: In install.packages("RCurl") : installation of package 'RCurl' had non-zero exit status --- O que é esse "curl-config"? Estou deixando de fazer algo? (o curl 7.21 está instalado) Abs, Fabrício. http://vizzual.net

2 2

nao consigo instalar o pacote "rimage" - Ubuntu
by Guilherme Moraes Ferraudo 31 Mai '11

31 Mai '11

Boa noite pessoal, estou no Ubuntu e não consigo instalar o pacote "rimage" - Ubuntu. Olha o erro: ERROR: configuration failed for package ‘rimage’ * removing ‘/home/gferrau/R/i486-pc-linux-gnu-library/2.10/rimage’ The downloaded packages are in ‘/tmp/RtmpcEgQZS/downloaded_packages’ Warning message: In install.packages("rimage", dep = T) : installation of package 'rimage' had non-zero exit status O que fazer?? Abs -- Guilherme Moraes Ferraudo Jaboticabal/Campinas - SP http://lattes.cnpq.br/2096118558794430

3 2

Distribuição Exponencial
by Fernanda Balbino 30 Mai '11

30 Mai '11

Olá, boa tarde a todos, Tenho um conjunto de dados que devido a natureza e aos gráficos de análise exploratória acredito que sigam uma distribuição exponencial. Gostaria de saber se no R é possível testar os dados para confirmar se realmente eles se comportam conforme uma distribuição exponencial. Muito obrigada. Fernanda Balbino.

3 2