- R-br - DINF UFPR

Ajuda com função para web-scraping no R
by Augusto Ribas 24 Jun '12

24 Jun '12

Ola pessoal. Estou com problemas com uma função aqui de web-scraping. Meu problema é o seguinte. Especies estão sempre mudando de nome. Novos generos surgem, outros somem, e assim vai. Mas tem um indexador de nome de especies que é uma mão na roda. O http://www.catalogueoflife.org/ Blz, é so colocar o nome nele que sai se o nome novo se mudou, o mesmo se nao mudou, enfim qual o nome aceito hoje em dia. Ai eu tentei fazer uma função pra fazer isso. Basicamente eu fiz ele fazer um search, criando o link, dai eu pego o primeiro resultado, entro na pagina do primeiro resultado, ai dentro da pagina do primeiro resultado eu pego o nome aceito e o autor. Eu tava tentando copiar essa função de algum tempo atras aqui da lista pra usar na minha necessidade http://r-br.2285057.n4.nabble.com/R-br-OT-Distancias-entre-municipios-td396… So que a entrada era pra ser o nome da especie, dai a saida era pra ser uma lista com o nome da especie e o autor. Ai o meu problem é o seguinte. Ps. a função ta no final da msg, pra nao bagunçar aqui. Qd eu uso ele sozinho sai blz >-sp.check(dados$Hospedeiro[1]) $sp.aceito [1] "Aetobatus narinari" $autor [1] "Euphrasen, 1790" Mas as vezes sai assim > sp.check(dados$Hospedeiro[1]) $sp.aceito ] "Aetobatus narinari" $autor [1] "Euphrasen, 1790" Não sei pq some o "[1" que deveria ter ali Mas ai qd eu uso pra pra saber so o nome da especie por exemplo > sp.check(dados$Hospedeiro[1])$sp.aceito tus narinari" > sp.check(dados$Hospedeiro[1])$sp.aceito narinari" Ele não funciona direito, come metade do nome ali, come umas letras sei la A minha ideia final é o seguinte, eu tenho um dafa.frame tipo assim: ##################################################################### > dput(exemplo) structure(list(Hospedeiro = c("Aetobatis narinari", "Arothron hispidus", "Astyanax aenaeus", "Astyanax fasciatus", "Astyanax fasciatus" ), Parasita = c("Acanthobothrium nicoyaense", "Bianium simonei", "Wallinia chavariiae", "Anacanthocotyle anacanthocotyle", "Gyrodactylus costaricensis" ), Grupo.H = c("peixes", "peixes", "peixes", "peixes", "peixes" ), Grupo.P = c("C", "D", "D", "M", "M")), .Names = c("Hospedeiro", "Parasita", "Grupo.H", "Grupo.P"), row.names = c(NA, 5L), class = "data.frame") ######################################################################## Eu gostaria de criar outra coluna, com o nome da especie aceito atualmente, usando tipo um loop com essa função que estou tentando fazer. E uma coluna com o autor que descreveu a especie tb para a coluna hospedeiros. Se alguem puder dar uma olhadinha o que estou fazendo errado. Eu quero chegar a isso: Hospedeiro sp aceita outras colunas 1 Aetobatis narinari Fulana detal ---- 2 Arothron hispidus Ciclana detal ----- 3 Astyanax aenaeus Beltrana detal ------ 4 Astyanax fasciatus Mesmo nome ----- 5 Astyanax fasciatus Segue aqui passo a passo o que eu fiz e tudo na forma de função também. ####################################### #Passo a passo ######################################## #Entrada do nome da especie especie<-c("Bufo marinus") #Separando o nome da especie gen<-strsplit(especie,"\\ ")[[1]][1] esp<-strsplit(especie,"\\ ")[[1]][2] gen esp #Criando o link link<-paste("http://www.catalogueoflife.org/col/search/all/key/",gen,"+",esp,"/match/1",sep="") link <- iconv(link, 'latin1', 'UTF-8') Encoding(link) <- 'bytes' link #lendo a pagina pagina <- readLines(url(link)) #Achando onde esta o primeiro resultado da lista, Acredito que essa linha #é o começo da tabelinha de resultado sempre, não sei o que mais da pra usar #pra achar o primeiro resultado :( n.linhas<-which(pagina%in%" <td class=\"field_header_black\">") #aqui pra testa depois na função se nao achou nada, menor que 1 ele nao achou nada no site length(n.linhas) #pegando o primeiro resultado pag.sp<-strsplit(pagina[n.linhas[1]+1],'\\"')[[1]][2] #gerando o link do primeiro resultado link2 <- paste( "http://www.catalogueoflife.org",pag.sp,sep="") link2 <- iconv(link2, 'latin1', 'UTF-8') Encoding(link2) <- 'bytes' link2 #baixando a pagina do primeiro resultado pagina2 <- readLines(url(link2)) #achando a linha onde esta o nome da especie aceito linha2<-grep('(accepted name)',pagina2) sp.final<-pagina2[linha2] #pegando o nome da especie corte1<-strsplit(sp.final,'<i>')[[1]][2] sp.aceito<-strsplit(corte1,'</i>')[[1]][1] #pegando o autor corte2<-strsplit(sp.final,'\\(')[[1]][2] autor<-strsplit(corte2,')')[[1]][1] #Resultados de interesse sp.aceito autor ####################################### #Agora uma função pra faze isso ######################################## #função sp.check sp.check<-function(especie) { especie<-as.character(especie) gen<-strsplit(especie,"\\ ")[[1]][1] esp<-strsplit(especie,"\\ ")[[1]][2] link<-paste("http://www.catalogueoflife.org/col/search/all/key/",gen,"+",esp,"/match/1",sep="") link <- iconv(link, 'latin1', 'UTF-8') Encoding(link) <- 'bytes' pagina <- readLines(url(link)) n.linhas<-which(pagina%in%" <td class=\"field_header_black\">") if(length(n.linhas)>0) { pag.sp<-strsplit(pagina[n.linhas[1]+1],'\\"')[[1]][2] link2 <- paste( "http://www.catalogueoflife.org",pag.sp,sep="") link2 <- iconv(link2, 'latin1', 'UTF-8') Encoding(link2) <- 'bytes' link2 pagina2 <- readLines(url(link2)) linha2<-grep('(accepted name)',pagina2) sp.final<-pagina2[linha2] corte1<-strsplit(sp.final,'<i>')[[1]][2] sp.aceito<-strsplit(corte1,'</i>')[[1]][1] corte2<-strsplit(sp.final,'\\(')[[1]][2] autor<-strsplit(corte2,')')[[1]][1] }else { sp.aceito<-c("Não encontrado") autor<-c("Não encontrado") } return(list(sp.aceito=sp.aceito,autor=autor)) } PS1: O as.character la no começo é que normalmente eu leio os dataframe com read.csv2 ai converte pra fator a coluna, ai quis fazer ele transforma ali pra nao da pau. PS2: Eu tb coloque um if pra dar um jeito qd ele não achar resultado, a principio tava funcionando. PS3: Tem como suprimir esse erro: In readLines(url(link2)) : linha final incompleta encontrada em 'http://www.catalogueoflife.org/col/details/species/id/7159709/synonym/71619…' PS4: Tem como o loop não parar qd a internet não funciona. tipo qd da aqueles erros 505 e a pagina nao abre. Agradeço a atenção e peço desculpas pela mensagem muito longa. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Ajuda com uso de list(). Como fazer uma lista de listas?
by Augusto Ribas 24 Jun '12

24 Jun '12

Ola pessoal. Estou aqui sofrendo para usar listas no R, resumidamente eu estou com dificuldades em criar uma lista de listas. Segue aqui um CMR, espero não esteja muito longo. #Eu tenho uma matriz que é o rede de interação entre algumas especies de hospedeiros e parasitas, por exemplo: dados<-matrix(c(1,1,1,0,0,0, 0,0,1,0,0,0, 0,0,1,0,0,0, 0,0,0,1,0,0, 0,0,0,1,1,0, 0,0,0,0,0,1),byrow=T,ncol=6,nrow=6, dimnames=list(paste("P",1:6,sep=""),paste("H",1:6,sep=""))) dados #E estava querendo ver qual o efeito das extinções dos hospedeiros nos parasitas, então eu crio uma sequencia de extinção de hospedeiros set.seed(13) extinto<-sample(colnames(dados),ncol(dados),replace=F) extinto # vejo de acordo com essas extinções, qts parasitas vão morrendo por falta de hospedeiro, então fiz uma função pra contar isso ext.cont<-function(dados,extintos) { spextintas<-sum(rowSums(dados[,!dimnames(dados)[[2]]%in%extintos,drop=F])<1) return(spextintas) } #ai salvo o resultado num vetor resultado<-rep(NA,ncol(dados)) resultado #loop faço um loop pra ver a sequencia, tipo somente a primeira especie, a primeira e a segunda, as tres primeiras até todas as especies. for(i in 1:ncol(dados)) { resultado[i]<-ext.cont(dados,extinto[1:i]) } resultado #mas ai eu posso fazer isso varias vezes pra ver como é de forma geral final<-cbind(resultado) for (j in 1:4) { extinto<-sample(colnames(dados),ncol(dados),replace=F) for(i in 1:1:ncol(dados)) { resultado[i]<-ext.cont(dados,extinto[1:i]) } final<-cbind(final,resultado) } final matplot(1:6,final,type="b",pch=1,lty=1) #so que eu gostaria de saber quem esta extinguindo também, então fiz uma função ext.list<-function(dados,extintos) { lista.extintos<-rownames(dados)[which(rowSums(dados[,!dimnames(dados)[[2]]%in%extintos,drop=F])<1)] return(lista.extintos) } # e como pode tanto não extinguir ninguem como mais de uma expecie imagino que seja interessante usar uma lista lista.extintos<-list(NA) #loop for(i in 1:ncol(dados)) { lista.extintos[[i]]<-ext.list(dados,extinto[1:i]) } lista.extintos #so que ai so para a contagem, eu conseguia salva varios resultados na forma de uma matriz. #Mas aqui eu não sei o que fazer, extiste alguma forma de eu fazer um objeto que seja uma lista de listas, pra guardar esses resultados. #Outra duvida sobre listas é no caso da lista acima, vamos supor que eu gostaria de saber em qual numero (o [[n]]) o "P2" esta a primeira vez, tem como fazer? #Por exemplo se fosse um vetor eu usaria o teste == , agora sendo uma lista como eu faço pra saber em quais elementos ([[n]]) ele esta, ou a primeira vez que aparece. #E se der pra fazer uma lista de listas, como faz? #Espero não ter sido confuso na terminologia. Muito obrigado e boa a noite a todos :) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Pre-tratamento de dados
by T Branquinho 24 Jun '12

24 Jun '12

Boa tarde, Estou montando um padrão de pre-tratamento de dados, a princípio montei para uma matrix de 39 linhas e 31 colunas, porém gostaria de fazer de uma forma que funcione para qualquer tamanho de matrix. Segue o que fiz (Em negrito estão as informações que necessito calcular) M1<-matrix(min(B1[,1]),39,1);M2<-matrix(min(B1[,2]),39,1);M3<-matrix(min(B1[,3]),39,1);M4<-matrix(min(B1[,4]),39,1);M5<-matrix(min(B1[,5]),39,1);M6<-matrix(min(B1[,6]),39,1);M7<-matrix(min(B1[,7]),39,1);M8<-matrix(min(B1[,8]),39,1);M9<-matrix(min(B1[,9]),39,1);M10<-matrix(min(B1[,10]),39,1);M11<-matrix(min(B1[,11]),39,1);M12<-matrix(min(B1[,12]),39,1);M13<-matrix(min(B1[,13]),39,1);M14<-matrix(min(B1[,14]),39,1);M15<-matrix(min(B1[,15]),39,1);M16<-matrix(min(B1[,16]),39,1);M17<-matrix(min(B1[,17]),39,1);M18<-matrix(min(B1[,18]),39,1);M19<-matrix(min(B1[,19]),39,1);M20<-matrix(min(B1[,20]),39,1);M21<-matrix(min(B1[,21]),39,1);M22<-matrix(min(B1[,22]),39,1);M23<-matrix(min(B1[,23]),39,1);M24<-matrix(min(B1[,24]),39,1);M25<-matrix(min(B1[,25]),39,1);M26<-matrix(min(B1[,26]),39,1);M27<-matrix(min(B1[,27]),39,1);M28<-matrix(min(B1[,28]),39,1);M29<-matrix(min(B1[,29]),39,1);M30<-matrix(min(B1[,30]),39,1);M31<-matrix(min(B1[,31]),39,1);Mmin<-cbind(M1,M2,M3,M4,M5,M6,M7,M8,M9,M10,M11,M12,M13,M14,M15,M16,M17,M18,M19,M20,M21,M22,M23,M24,M25,M26,M27,M28,M29,M30,M31);M1<-matrix(max(B1[,1]),39,1);M2<-matrix(max(B1[,2]),39,1);M3<-matrix(max(B1[,3]),39,1);M4<-matrix(max(B1[,4]),39,1);M5<-matrix(max(B1[,5]),39,1);M6<-matrix(max(B1[,6]),39,1);M7<-matrix(max(B1[,7]),39,1);M8<-matrix(max(B1[,8]),39,1);M9<-matrix(max(B1[,9]),39,1);M10<-matrix(max(B1[,10]),39,1);M11<-matrix(max(B1[,11]),39,1);M12<-matrix(max(B1[,12]),39,1);M13<-matrix(max(B1[,13]),39,1);M14<-matrix(max(B1[,14]),39,1);M15<-matrix(max(B1[,15]),39,1);M16<-matrix(max(B1[,16]),39,1);M17<-matrix(max(B1[,17]),39,1);M18<-matrix(max(B1[,18]),39,1);M19<-matrix(max(B1[,19]),39,1);M20<-matrix(max(B1[,20]),39,1);M21<-matrix(max(B1[,21]),39,1);M22<-matrix(max(B1[,22]),39,1);M23<-matrix(max(B1[,23]),39,1);M24<-matrix(max(B1[,24]),39,1);M25<-matrix(max(B1[,25]),39,1);M26<-matrix(max(B1[,26]),39,1);M27<-matrix(max(B1[,27]),39,1);M28<-matrix(max(B1[,28]),39,1);M29<-matrix(max(B1[,29]),39,1);M30<-matrix(max(B1[,30]),39,1);M31<-matrix(max(B1[,31]),39,1);Mmax<-cbind(M1,M2,M3,M4,M5,M6,M7,M8,M9,M10,M11,M12,M13,M14,M15,M16,M17,M18,M19,M20,M21,M22,M23,M24,M25,M26,M27,M28,M29,M30,M31);M1<-matrix(mean(B1[,1]),39,1);M2<-matrix(mean(B1[,2]),39,1);M3<-matrix(mean(B1[,3]),39,1);M4<-matrix(mean(B1[,4]),39,1);M5<-matrix(mean(B1[,5]),39,1);M6<-matrix(mean(B1[,6]),39,1);M7<-matrix(mean(B1[,7]),39,1);M8<-matrix(mean(B1[,8]),39,1);M9<-matrix(mean(B1[,9]),39,1);M10<-matrix(mean(B1[,10]),39,1);M11<-matrix(mean(B1[,11]),39,1);M12<-matrix(mean(B1[,12]),39,1);M13<-matrix(mean(B1[,13]),39,1);M14<-matrix(mean(B1[,14]),39,1);M15<-matrix(mean(B1[,15]),39,1);M16<-matrix(mean(B1[,16]),39,1);M17<-matrix(mean(B1[,17]),39,1);M18<-matrix(mean(B1[,18]),39,1);M19<-matrix(mean(B1[,19]),39,1);M20<-matrix(mean(B1[,20]),39,1);M21<-matrix(mean(B1[,21]),39,1);M22<-matrix(mean(B1[,22]),39,1);M23<-matrix(mean(B1[,23]),39,1);M24<-matrix(mean(B1[,24]),39,1);M25<-matrix(mean(B1[,25]),39,1);M26<-matrix(mean(B1[,26]),39,1);M27<-matrix(mean(B1[,27]),39,1);M28<-matrix(mean(B1[,28]),39,1);M29<-matrix(mean(B1[,29]),39,1);M30<-matrix(mean(B1[,30]),39,1);M31<-matrix(mean(B1[,31]),39,1);Mmean<-cbind(M1,M2,M3,M4,M5,M6,M7,M8,M9,M10,M11,M12,M13,M14,M15,M16,M17,M18,M19,M20,M21,M22,M23,M24,M25,M26,M27,M28,M29,M30,M31); M1<-matrix(var(B1[,1]),39,1);M2<-matrix(var(B1[,2]),39,1);M3<-matrix(var(B1[,3]),39,1);M4<-matrix(var(B1[,4]),39,1);M5<-matrix(var(B1[,5]),39,1);M6<-matrix(var(B1[,6]),39,1);M7<-matrix(var(B1[,7]),39,1);M8<-matrix(var(B1[,8]),39,1);M9<-matrix(var(B1[,9]),39,1);M10<-matrix(var(B1[,10]),39,1);M11<-matrix(var(B1[,11]),39,1);M12<-matrix(var(B1[,12]),39,1);M13<-matrix(var(B1[,13]),39,1);M14<-matrix(var(B1[,14]),39,1);M15<-matrix(var(B1[,15]),39,1);M16<-matrix(var(B1[,16]),39,1);M17<-matrix(var(B1[,17]),39,1);M18<-matrix(var(B1[,18]),39,1);M19<-matrix(var(B1[,19]),39,1);M20<-matrix(var(B1[,20]),39,1);M21<-matrix(var(B1[,21]),39,1);M22<-matrix(var(B1[,22]),39,1);M23<-matrix(var(B1[,23]),39,1);M24<-matrix(var(B1[,24]),39,1);M25<-matrix(var(B1[,25]),39,1);M26<-matrix(var(B1[,26]),39,1);M27<-matrix(var(B1[,27]),39,1);M28<-matrix(var(B1[,28]),39,1);M29<-matrix(var(B1[,29]),39,1);M30<-matrix(var(B1[,30]),39,1);M31<-matrix(var(B1[,31]),39,1);Mvar<-cbind(M1,M2,M3,M4,M5,M6,M7,M8,M9,M10,M11,M12,M13,M14,M15,M16,M17,M18,M19,M20,M21,M22,M23,M24,M25,M26,M27,M28,M29,M30,M31) N1<-(B1-Mmin)/(Mmax-Mmin);N2<-(B1-Mmean)/(Mmax-Mmean);N3<-(B1-Mmean)/(Mmax-Mmin);N4<-(B1-Mmean)/Mvar;N5<-sqrt(B1);N6<-log(B1);N7<-log(B1,10);N8<-1/sqrt(B1);N9<-1/B1;N10<-log((B1-Mmin)/(Mmax-B1)); Agradeço antecipadamente qualquer tipo de auxílio. Tales

7 10

TukeyC: Convencional Tukey Test (Perform the conventional Tukey test from aov and aov.list objects)
by Jose Claudio Faria 24 Jun '12

24 Jun '12

Olá, Gostaria de anunciar que encontra-se disponível nos repositórios CRAN um novo pacote sobre um velho tema - TukeyC: http://cran.r-project.org/web/packages/TukeyC/ dos autores: Jose Claudio Faria, Enio Jelihovschi e Ivan Bezerra Allaman. O pacote facilita bastante a aplicação do teste de Tukey para os delineamentos simples (DIC, DBC, DQL) e esquemas experimentais simples e usuais (Fatorias com no máximo 3 fatores, parcela sub-dividida e sub-sub-dividida). Para os que já conhecem o pacote ScottKnott, o TukeyC tem a mesma estrutura básica, porém aplica o teste de Tukey. A lista que armazena o resultado do teste contém vários componentes (slots) que podem ser muito bastante úteis no ensino, uso cotidiano do teste e pós-processamento. Um dos componentes armazena uma matriz (Diff_Prob) que mostra as diferença observadas entre as médias dos tratamentos e a probabilidade mínima a partir da qual a hipótese Ho poderia se rejeitada. Outro armazena o agrupamento convencional das médias dos tratamentos seguida de letras. Pode-se plotar o resultado do teste. Nesse caso serão mostrados: as médias, o máximo e mínimo de cada tratamento (na forma de barra vertical) e o agrupamento convencional dos tratamentos, ou seja, as letras que identificam os grupos de médias. Quando aplicável, realiza o teste para experimentos desbalanceados. O novo pacote passou por um longo período de desenvolvimento e testes. Vem sendo bastante usado pelos autores no ensino e em análises de dados. Com o passar do tempo e a partir das demandas, pretende-se expandir o teste para delineamentos e esquemas experimentias mais complexos . Vale a pena conferir! (sugestões são sempre bem vindas) Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Re: [R-br] Partição de soma de quadrados, experimento fatorial
by Listeiro 037 24 Jun '12

24 Jun '12

Olá. Não costumo aparecer aqui. Você quer dividir um número como soma de quadrados, isto é: n = n_1²+n_2²+n_3²+... Fiquei com esta impressão. É isto mesmo? Em Sat, 23 Jun 2012 16:47:47 -0300 Lucas Barbosa <lucasbcr(a)gmail.com> escreveu: > Boa tarde pessoal, > Estou com esses dados de um experimento com dois fatores. Quero > particionar a interação "fator.a X fator.b" em "fator.a.linear X > fator.b" e "fator.a.quadratica X fator.b" usando o argumento *split > *da *aov().* Como posso fazer isso? -- "There are three kinds of lies: lies, damn lies, and statistics." Benjamin Disraeli

3 4

Variaveis Dummy.
by Julia Maria 23 Jun '12

23 Jun '12

Pessoal, preciso montar um modelo utilizando variaveis Dummy e gostaria que se possível alguém me desse uma idéia, ou mesmo me indicasse um exemplo para verificar. O problema é o seguinte, tenho preço de aluguel em algumas cidades do Brasil. Rodei no R uma regressão linear simples, do seguinte modo: variavel resposta: y=preço variaveis explicativas:x1=renda e x2=despesas obtive pelo R toda informação necessária através de um lm(y~x1+x2) Agora preciso incorporar variaveis dummy do meu modelo do seguinte modo. Pra cada cidade, associo 0 se ela pertence a determinada região e 1 se não pertence. Exemplo: Curitiba: 1 para a regiao sul e 0 para as demais. Desse modo defini mais 5 variaveis, r1, r2, r3, r4, r5 cada uma correspondendo a uma região do pais (Norte, Sul, Centro- Oeste, Sudeste e Nordeste). De posse disso, refiz o modelo da seguinte forma: lm(y~x1+x2+r1+r2+r3+r4). Mas algo me diz que está errado e deveria ser algo do tipo: lm(y~x1*r1+x1*r2+...+x2*r4) Alguem poderia me dar uma ajuda. Sou aluna de graduação, um pouco curiosa e querendo me apaizonar pelo R.

2 1

Problema com odfWeaver
by Diogo Braga Mendes 23 Jun '12

23 Jun '12

Boa noite grupo, Comecei a usar o odfWeaver para gerar relatórios neste formato entretanto o arquivo gerado repete várias vezes o mesmo resultado. Criei um arquivo *.odt e colei este código: *<<cemNumeros, print=F, echo=T>>=* *cemNumeros<-round(rnorm(10,74,2.34))* *summary(cemNumeros)* *@* Quando uso a função, gero o arquivo porém o código acima é replicado mutas vezes e com o resultado tb, o mesmo quando uso gráficos fico com dez ou mais gráficos replicados. Tentei mudar do LibreOffice 3.5 para o BrOffice 3.3 mas não adiantou, talvez só o OpenOffice mesmo? segue minhas informações: > sessionInfo()R version 2.13.1 (2011-07-08) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 [2] LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 [4] LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 attached base packages: [1] stats graphics grDevices utils datasets methods [7] base other attached packages: [1] odfWeave_0.7.17 XML_3.4-2.2 lattice_0.19-30 loaded via a namespace (and not attached): [1] grid_2.13.1 tools_2.13.1 -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes Bacharel em Estatística, CONRE 9558 Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “O âmago da ciência, que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata, nascida da infatigável curiosidade, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

2 1

Re: [R-br] Partição de soma de quadrados, experimento fatorial
by Listeiro 037 23 Jun '12

23 Jun '12

Olá. Não costumo aparecer aqui. Você quer dividir um número como soma de quadrados, isto é: n = n_1²+n_2²+n_3²+... Fiquei com esta impressão. É isto mesmo? Em Sat, 23 Jun 2012 16:47:47 -0300 Lucas Barbosa <lucasbcr(a)gmail.com> escreveu: > Boa tarde pessoal, > Estou com esses dados de um experimento com dois fatores. Quero > particionar a interação "fator.a X fator.b" em "fator.a.linear X > fator.b" e "fator.a.quadratica X fator.b" usando o argumento *split > *da *aov().* Como posso fazer isso? -- "There are three kinds of lies: lies, damn lies, and statistics." Benjamin Disraeli

1 0

Partição de soma de quadrados, experimento fatorial
by Lucas Barbosa 23 Jun '12

23 Jun '12

Boa tarde pessoal, Estou com esses dados de um experimento com dois fatores. Quero particionar a interação "fator.a X fator.b" em "fator.a.linear X fator.b" e "fator.a.quadratica X fator.b" usando o argumento *split *da *aov().* Como posso fazer isso? dados = c(1.8, 2.1, 2, 2.1, 4.6, 5, 7.5, 7.9, 2.2, 2.4, 4.2, 4, 5.4, 5.6, 9.8, 9.2, 2.8, 3.2, 4.4, 4.8, 8.7, 8.4, 13.2, 13, 3.2, 3.6, 3.3, 3.5, 5.7, 5.8, 10.9, 11.1) dados = data.frame(resposta = dados, fator.a = factor( rep(seq(210, 225, 5), each=2, times=4) , ordered=T), fator.b = as.factor(rep(LETTERS[1:4], each=4))) a = aov(resposta ~ fator.a*fator.b, data=dados) summary(a) Valeu, gente! Lucas Cusinato

1 0

Holt-Winters e ARIMA
by kaue veras 22 Jun '12

22 Jun '12

Bom dia grupo, Estou realizando meu projeto final e preciso fazer minha previsão da serie através dos modelos de Holt-Winters e ARIMA pelo R, poderiam me auxiliar no meu trabalho? Desde já agradeço. Att, Kauê P. Veras da Cunha Estatística / 8º Período - UERJ MSN: kaueveras(a)hotmail.com E-mail: kaueveras(a)gmail.com Cel: 8254-9601

3 4

Problemas com o R
by Ana Lúcia Souza Silva Mateus 22 Jun '12

22 Jun '12

Prezados, Estou com o seguinte problema .... Rodei um script para 10 e 15 simulações, e o meu programa rodou. Quando acrescento o numero de simulações a partir deste número, o R fecha ... O que devo fazer? Grata

4 4

Treinamento - Redes Neurais
by Edmar Caldas 22 Jun '12

22 Jun '12

Bom Dia! Pessoal sou usuário do spss e não tenho muita facilidade no R. Preciso fazer um trabalho sobre TREINAMENTO INTERATIVO REDES NEURAIS e o spss não faz isso. Andei pesquisando e vi que o R pode fazer isso. Gostaria de uma ajuda. Preciso instalar com pacote? Se vocês tiverem um código “exemplo” me manem por favor. Obrigado. Edmar

2 1

Analise de associação
by Edmar Caldas 22 Jun '12

22 Jun '12

Bom Dia! Pessoal sou usuário do spss e não tenho muita facilidade no R. Preciso desenvolver uma analise de associação de cesta de produtos. Pelo que pesquisei tive que baixar o pacote ARULES e usar a função APRIORI mas não consegui, alguém poderia me ajudar por favor. Obrigado. Atenciosamente, Edmar

2 1

zip
by Vinicius Brito Rocha 21 Jun '12

21 Jun '12

Pessoal, alguém aí já usou a função zip ? tentei aqui mas não tive sucesso. -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

5 6

Ajuda em plot()
by Luciano F. Sgarbi 21 Jun '12

21 Jun '12

olá pessoal, Alguém sabe como eu faço para diminuir as margens de cada um dos gráficos dentro da janela para que cada um fique um pouco maior? segue exemplo: par(mfrow=c(2,3),pch=16) plot(sample(1:10),1:10,ylab="Y",xlab="X") plot(sample(1:10),1:10,ylab="Y",xlab="X") plot(sample(1:10),1:10,ylab="Y",xlab="X") plot(sample(1:10),1:10,ylab="Y",xlab="X") plot(sample(1:10),1:10,ylab="Y",xlab="X") plot(sample(1:10),1:10,ylab="Y",xlab="X") Abraços Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

4 3

reamostrar linhas de uma matriz
by Edson Ferreira 21 Jun '12

21 Jun '12

olá amigos do R, da matriz data abaixo, preciso de formar uma matriz com 1000 linhas onde cada linha da nova matriz é exatamente uma das linhas da matriz data, estou usando o código abaixo mas não estu conseguindo, alguém pode ajudar? data [,1] [,2] [,3] [,4] [1,] 1 0 0 1 [2,] 0 1 0 0 [3,] 1 1 1 0 [4,] 1 1 1 1 [5,] 1 1 0 1 [6,] 0 0 0 1 [7,] 0 0 0 0 simulacao<- data[sample(7, 1), ] # amostrar linha da matriz simulacao n<-1000 simulacao<-NULL padrao<-NULL for(i in 1:n){ simulacao[i]<-data[sample(7, 1), ] y<-simulacao[i] padrao<-matrix(y,nrow=n,ncol=4,byrow=FALSE) } padrao

3 2

diagnóstico de imputação multipla para fatores
by Pedro Emmanuel Alvarenga Americano do Brasil 21 Jun '12

21 Jun '12

Amigos de R, Estou imputando dados utilizando a função mice do pacote de mesmo nome. MIce é um função muito automática e facil de usar se os dados ja estão num formato facil pra função. A minha questão é explorar de forma simples como ocorreu a imputação. Para isso, há algumas funções no pacote que fzem graficos como por exemplo o bwplot.mids() então seguindo o exemplo da função: imp <- mice(boys,maxit=1) bwplot(imp) O problema é que essa função, assim como outras no mesmo pacote, só avalia as variáveis numericas contínuas. Eu imaginei que haveria alguma coisa parecida para fatores, como plot de barras, pois a maioria das variáveis que tenho são fatores, mas não encontrei. A pergunta: é se alguém tem alguma ideia ou sugestão para avaliar a imputação de variáveis formatadas fatores? Abraço, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, 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.libreoffice.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

[OFF - TOPIC] Material sobre cálculo de tamanho amostral
by Andre Oliveira 20 Jun '12

20 Jun '12

Pessaol, alguém tem algum material sobre cálculo de tamanho amostral em pesquisa de opinião, em pdf, que possa compartilhar? Ou ate mesmo sobre qualquer cálculo de tamanho amostral! andreolisouza(a)gmail.com André Oliveira Souza

3 3

Re: [R-br] Atualização do R no Linux Ubuntu 10.04
by Andre Oliveira 20 Jun '12

20 Jun '12

Nunca fiz tal compilação e sou usuário com pouco conhecimento sobre o assunto. Por favor gostaria se possível de alguma ajuda para fazer tal compilação algum tutorial ou algo semelhante. obrigado André Oliveira Souza ________________________________ De: Paulo Justiniano <paulojus(a)leg.ufpr.br> Para: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>; Andre Oliveira <andreolsouza(a)yahoo.com.br> Enviadas: Domingo, 17 de Junho de 2012 21:59 Assunto: Re: [R-br] Atualização do R no Linux Ubuntu 10.04 andré é provavel que os binários nao sejam mantidos, mas voce pode compular o R da fonte no ubuntu descontinuado sem problemas *eu faço exatamente isto no meu neste momento) On Sun, 17 Jun 2012, Andre Oliveira wrote: > Boa noite, > estou com a dúvida uso ubuntu 10.04 depois de alguns anos considero que em minha máquina está tudo configurado como > preciso criei ate uma versão persolizada em um DVD via garfio. Em abril de 2013 o 10.04 deixará de ter suporte. > Minha dúvida é a seguinte eu receberei atualização do R caso eu decida ficar com o 10.04 após abril de 2013? > > obrigado > > André Oliveira Souza > >

2 1

Erro no comando sqlFetch: formato de data.
by Carolina Pena 20 Jun '12

20 Jun '12

Olá pessoal, Estou lendo no R um arquivo ACCESS com diversas tabelas. Eu transformo cada uma das tabelas em um data.frame para trabalhar. O problema, no entanto, é que ao tentar fazer a transformação de uma das tabelas eu obtive o seguinte erro: >con <- odbcConnectAccess(file.choose()) > sqlFetch(con, banco1 ) Erro em as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format Uma das colunas do meu banco de dados (banco1) é uma data no formato brasileiro (dd/mm/YYYY). Eu li no google que o erro ocorre por causa do formato da data. Inclusive quando deleto essa coluna da base de dados o comando funciona normalmente. Eu gostaria de saber se consigo acessar no próprio R essa coluna (antes de transformar os dados em data.frame utilizando o comando sqlFetch) e utilizar o comando as.Date para tentar solucionar o problema. Obrigada! Carolina

2 1

[OFF-TOPIC] Indicação de livro
by Elton Cavalcanti 20 Jun '12

20 Jun '12

Bom dia, pessoal Vocês poderiam me indicar um bom livro de estatística que um estudante da área de saúde poderia usar, mas que não fosse tão superficial na abordagem de alguns conceitos estatísticos? Abraços, Elton Cavalcanti

2 2

Saída intergen e avaliação com pacote coda
by Dalinne Chrystian Carvalho dos Santos 20 Jun '12

20 Jun '12

Olá pessoal, tudo bem? Meu nome é Dalinne e eu tenho uma saída do programa intergen com análise Bayesiana no formato "txt" da seguinte maneira: 1 27.90086641691324 -42.05430078214862 6609.247759721497 5.547673239405648 -109.7603421548792 6665.595598994243 2 24.47414461350570 -37.16170819212519 8015.513063087009 5.783134528755414 -121.6600727394209 7681.170349762655...10000 5.741506313680731 -78.70398130578000 5932.380770110504 3.514573296713659 -109.4777360152619 6552.098731256713 Ou seja, as minhas 10.000 interações estão em apenas uma coluna.Queria saber como entrar com esse banco de dados no R mas como um objeto mcmc e depois disso como eu faço paraachar o ponto de burn-in, por meio do pacote coda do R? Obrigada pela atenção, Dalinne Chrystian Carvalho dos Santos

1 0

Pacote DiagnosisMed
by Edson Lira 20 Jun '12

20 Jun '12

Ao tentar instalar o pacote DiagnosisMed veriquei que não consta mais no CRAN. Alguém sabe me explicar o que houve? Qual pacote que pode ser utilizado para as mesmas análises que o mesmo fazia? Bom dia a todos os R-lista []'s. Edson Lira Estatístico Manaus-Amazonas

3 3

visualização cluster
by Vinicius Brito Rocha 20 Jun '12

20 Jun '12

Pessoal, estou plotando resultados de cluster kmeans em duas componentes e então colorindo os grupos, mas estava querendo fazer uma coisa diferente. Gostaria de plotar as fronteiras resultantes dos clusters. Mais ou menos como desenhar as regiões de fronteira de cada cluster. Alguém tem alguma idéia? Abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

Testes Estacionariedade Séries Temporais
by Josenilson Gomes 20 Jun '12

20 Jun '12

Olá, Gostaria de saber se alguém já usou algum teste de estacionariedade no R. Tenho interesse em usar os seguintes testes, mas não sei se já estão disponíveis no R: -Dickey-Fuller -Phillips-Perron -KPSS Todos testam se a série é estacionária... Agradeço desde já. Att., Josenílson

5 5

Simulação de distribuição bivariada
by Rodrigo Coster 19 Jun '12

19 Jun '12

Caros, o R tem alguma função para simulação de distribuições bivariadas? Seguem as funções densidades e distribuição da distribuição que tenho intesse em simular (u e v são as variaveis, tauU e tauL os parametros, todos no intervalo (0,1)) dsjc <- function(u,v,tauU,tauL) { ## Adaptado do Código para Matlab de Andrew Patton - http://publicecondukeedu/~ap172/codehtml k1 = 1/log2(2-tauU) k2 = -1/log2(tauL) CL1 = ((1 - (1 - u)^k1)^(k2 - 1)* (1 - u)^(k1 - 1)*(-1 + k1*(k2* (-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))) + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))))* (1 - (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(-k2^(-1)))^(k1^(-1))* (1 - (1 - v)^k1)^(k2 - 1)* (1 - v)^(k1 - 1)) CL2 = (((-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1)))^2) * ((1 - (1 - u)^k1)^k2 + (1 - (1 - v)^k1)^k2 - (1 - (1 - u)^k1)^k2* (1 - (1 - v)^k1)^k2)^2) CL1 = CL1/CL2 k1 = 1/log2(2-tauL) k2 = -1/log2(tauU) u = 1-u v = 1-v CL3 = ((1 - (1 - u)^k1)^(k2 - 1)* (1 - u)^(k1 - 1)*(-1 + k1*(k2* (-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))) + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))))* (1 - (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(-k2^(-1)))^(k1^(-1))* (1 - (1 - v)^k1)^(k2 - 1)* (1 - v)^(k1 - 1)) CL4 = (((-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1)))^2) * ((1 - (1 - u)^k1)^k2 + (1 - (1 - v)^k1)^k2 - (1 - (1 - u)^k1)^k2* (1 - (1 - v)^k1)^k2)^2) CL3 = CL3/CL4 CL = 0.5*(CL1+CL3) return(CL) } psjc <- function(U,V,tauU,tauL) { ## Adaptado do Código para Matlab de Andrew Patton - http://publicecondukeedu/~ap172/codehtml K = 1/log2(2-tauU); G = -1/log2(tauL); out1 = 1-((1-(((1-((1-U)^K))^(-G))+((1-((1-V)^K))^(-G))-1)^(-1/G))^(1/K)); K = 1/log2(2-tauL); # switching the upper and lower measures G = -1/log2(tauU); U = 1-U; V = 1-V; out2 = (1-U) + (1-V) - 1 + 1-((1-(((1-((1-U)^K))^(-G))+((1-((1-V)^K))^(-G))-1)^(-1/G))^(1/K)); out1 = 0.5*(out1+out2); return(out1) }

3 4

comando no R
by Eliana Silva 19 Jun '12

19 Jun '12

Pessoal, alguem me pode dizer se existe um comando no R uma t_student p_variada? e para a gama inversa tb gostaria. obrigada

2 1

Ajuda em uma função
by Sérgio Henrique almeida da silva ju 19 Jun '12

19 Jun '12

Olá Estou com mais uma dúvida no desenvolvimento de uma função Preciso analisar um conjunto de dados, porém essa análise é feita a cada duas colunas: a = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2) a1 = c(4,5,5,4,4,3,2,1,4,3,4,5,3,2,3,4,5,5,4,4,1) d = c(3,3,4,4,5,4,3,2,4,5,5,3,4,4,5,4,3,2,3,4,5) d1 = c(2,2,3,3,4,3,2,3,4,5,5,5,4,3,2,3,4,5,4,3,2) dados = data.frame(a,a1,d,d1) library(psy) x=kappa(dados[,c(1,2)]) x y=kappa(dados[,c(3,4)]) y z=kappa(dados[,c(n,n+1)]) z Porém gostaria de rodar uma função que calcule isso automaticamente, dando todos os resultados, para isso eu daria a coluna de inicio e a coluna final. Como posso fazer? -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

3 7

Duvida unzip no R
by Humberto Hazin 19 Jun '12

19 Jun '12

Olá Pessoal Tenho uma pasta que contem vários arquivos com extensão bz2! Gostaria de saber como descomprimir no R. Alguém pode me ajudar com um exemplo? Humberto

3 4

Script Matriz de Vizinhança
by Isabela Zardo 19 Jun '12

19 Jun '12

Gostaria de scripts pra matriz de vizinhança, alguém pode eme ajudar? Isabela Zardo

3 4

Extrair intervalo de uma coluna...
by andrebvs 19 Jun '12

19 Jun '12

2 2

Extrair intervalo de uma coluna...
by andrebvs 18 Jun '12

18 Jun '12

1 0

Contar atraso de um número numa matrix...
by andrebvs 18 Jun '12

18 Jun '12

3 5

Buscar por palavras no site da ABL
by Lucas Barbosa 18 Jun '12

18 Jun '12

Boa tarde, pessoal! No site da ABL ( http://www.academia.org.br/abl/cgi/cgilua.exe/sys/start.htm?sid=23) é possível buscar palavras no Vocabulário Ortográfico da Língua Portuguesa. Mas a busca não pode ser definida a partir de um endereço html, é necessário preencher o campo e clicar em "Pesquisar". Como faço para automatizar esse tipo de busca usando R? É possível? Encontrei, na stackoverflow <http://stackoverflow.com/>, problemas semelhantes resolvidos com as funções postForm ou getForm (pacote RCurl), mas não estou entendendo muito bem como elas funcionam. O html da caixa de texto e do botão "Pesquisar": <input id="palavra" size="30" name="palavra" maxlenght="30"> <input onclick="ajaxport(palavra.value,'')" value="Pesquisar" name="pesquisar" type="button"> Obrigado! Lucas

2 3

Divulgação curso: "Business Intelligence - A partir de uma perspectiva Bayesiana"
by nani.design＠bol.com.br 18 Jun '12

18 Jun '12

1 0

regressão linear
by Bernardo Rodrigues 18 Jun '12

18 Jun '12

Alguem me pode ajudar como é que se determina distribuições posteriori da regressão linear a partir de distribuição a priori nao informativa?

4 3

Atualização do R no Linux Ubuntu 10.04
by Andre Oliveira 18 Jun '12

18 Jun '12

Boa noite, estou com a dúvida uso ubuntu 10.04 depois de alguns anos considero que em minha máquina está tudo configurado como preciso criei ate uma versão persolizada em um DVD via garfio. Em abril de 2013 o 10.04 deixará de ter suporte. Minha dúvida é a seguinte eu receberei atualização do R caso eu decida ficar com o 10.04 após abril de 2013? obrigado André Oliveira Souza

3 2

Diferença entre duas abordagens de comparação usando GLM
by ASANTOS 17 Jun '12

17 Jun '12

Boa noite pessoal, Gostaria de saber se do ponto de vista da interpretação estatística do método empregado, existe diferença entre comparar 3 tratamento hipotéticos usando: 1) require(multcomp) y<- c(mapply(rpois, lambda=c(5,20,5), MoreArgs=list(n=100)))##Criação da variável resposta trat <- gl(3,100) ##Criação dos tratamentos m.1 <- glm(y~trat,family="poisson") anova.m.1 <- anova(m.1, test="Chi") anova.m.1 summary(glht(m.1, linfct=mcp(trat="Tukey"))) #Os tratamentos 1 e 3 são semelhantes e diferentes do tratamentos 2 ou 2) m.1 <-glm(y~trat,family="poisson") m.0 <- glm(y~1,family="poisson") anova(m.1,m.0,test="Chi")####Existe diferença sort(tapply(y,trat,mean))## Ordenando as médias levels(trat) ##Modelo T1 e T3 t1t3<-trat levels(t1t3) levels(t1t3)[1]<-"1+3" levels(t1t3)[3]<-"1+3" levels(t1t3) m.2<- glm(y~t1t3,family="poisson") anova(m.1,m.2,test="Chi")####Semelhantes ## Mas se faço qualquer uma das outras combinações "1+2" ou "2+3" os valores do teste Chi são significativos, podendo-se concluir que os tratamentos 1 e 3 são semelhantes e diferentes do tratamentos 2. Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 2

Trabalhando com equações
by Diego Leal 17 Jun '12

17 Jun '12

Oi pessoal, gostaria de saber como faço pra trabalhar com equações no R project (como represento uma equação, maximar função, derivar e integrar) Existe algum tutorial, ou coisa do tipo onde eu possa me informar melhor sobre o assunto, não sei nem por onde começar?

4 7

[OFF-TOPIC]Usar "for" via terminal linux?
by Ivan Bezerra Allaman 17 Jun '12

17 Jun '12

Boa noite senhores! Alguém sabe se é possível usar algum comando semelhante ao "for" via terminal no linux? Eu tenho uma sequência de imagens para recortar e estou utilizando o camando: convert crn_01.png -crop 970x474+25+60 +repage crn_01.png Algo como: for(i in 1:200){ convert crn_i.png -crop 970x474+25+60 +repage crn_i.png } já ajudaria! Alguém sabe se é possível fazer algo semelhante no terminal? Desde já muito obrigado! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 5

como juntar dois bancos em função de um intervalo
by Paula Mendes Luz 17 Jun '12

17 Jun '12

oi pessoal eu normalmente utilizo o merge() para juntar bancos em função de 1 ou mais variáveis que são iguais nos dois bancos. também normalmente tenho bancos com medidas repetidas para cada indivíduo de modo que novas linhas são criadas no ato do merge() em função das informações dos bancos originais. agora preciso juntar dois bancos em função de um intervalo de números (na verdade são datas mas isso não importa): um exemplo do que preciso fazer p o banco original com > 3000 ids seria banco1 tem o id, a entrada, a saída, e o grupo id <- c(1, 1, 1, 2, 2) ent <- c(1, 10, 20, 1, 30) sai <- c(9, 19, 40, 29, 50) grupo <- c(1, 2, 3, 1, 2) banco1 <- cbind(id, ent, sai, grupo) banco2 tem o id, o dia, tipo do teste id <- c(1, 1, 1, 1, 1, 2, 2, 2, 2) dia <- c(5, 8, 12, 18, 25, 22, 25, 29, 41) teste <- c("A", "A", "B", "B", "A", "B", "B", "C", "D") banco2 <- cbind(id, dia, teste) e o banco 3 seria o merge() dos dois em função da avaliação se o DIA do banco2 está contido no intervalo (entre ENT e SAI) do banco1 o resultado seria o banco 3 com 5 linhas para o id1 com as 5 vars da seguinte forma id ent sai grupo dia teste 1 1 9 1 5 A 1 1 9 1 8 A 1 10 19 2 12 B 1 10 19 2 18 B 1 20 40 3 25 A e 4 linhas para o id 2 alguém pode ajudar???? eu li a respeito da função findInterval() mas não consegui fazer funcionar.... agradeço desde já! abraços, paula.

2 4

simulação da normal bivariada
by Natalia Martins 17 Jun '12

17 Jun '12

Boa noite pessoal, estou tentando simular dados diários de temperatura e radiação solar usando uma normal bivariada e após a simulação eu tentei fazer uma reta de regressão para verificar se minha simulação esta de acordo, ou seja, eu esperaria uma reta pois os valores deveriam ser iguais. No entanto quando a faço o meu R^2 ajustado é ruim (sei que tambem devo fazer testes e não apenas olhar o R), alguem pode me ajudar, pois acredito que estou fazendo algo errado? a programação que estou utilizando segue a baixo: # gerando os dados normal bivariada set.seed(1234) #semente sigma <- matrix(c(19260.71398,28.92129,28.92129,2.173789), ncol=2) sim_bnorm <- rmvnorm(n=1000, mean=c(mean(RST$RS),mean(RST$Temp)), sigma=sigma) bnorm <- data.frame(sim_bnorm) simulado <- sample(bnorm$X2,31) modelo <- lm(simulado~RST$Temp) summary(modelo) plot( RST$Temp , simulado,xlab="Valores reais",ylab="Valores simulados") abline(modelo,col="red") Muito obrigada --

4 8

(sem assunto)
by Diogo Braga Mendes 17 Jun '12

17 Jun '12

O problema está provavelmente na linha: *dataset = 0 * Troque por (use uma lista): *dataset = vector("list",vetZ)* e troque: dataset*[[i]]* = data.frame(read.csv(myDir[i],sep=";",header=TRUE)); -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes Bacharel em Estatística, CONRE 9558 Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “O âmago da ciência, que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata, nascida da infatigável curiosidade, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

1 0

como juntar dois bancos em função de um intervalo
by Diogo Braga Mendes 17 Jun '12

17 Jun '12

tente isso (ao invés de cbind): banco1 <- data.frame(id, ent, sai, grupo) banco2 <- data.frame(id, dia, teste) Você não precisará converte as variáveis (dia vai como fator para o banco3). -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes Estatístico do Instituto de Previdência e Assistência do Município de Belém - IPAMB Mestre em Estatística-PPGME/UFPA Bacharel em Estatística, CONRE 9558 Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “O âmago da ciência, que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata, nascida da infatigável curiosidade, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

1 0

Dúvida sobre read.csv
by Caio Deustch 17 Jun '12

17 Jun '12

Boa noite, estou tentando fazer um contador para alocar bancos de dados em variáveis diferentes. Quando faço para 1 variável ( dataset ), esta variável aloca o banco corretamente. Quando ponho do contador, o meu vetor dataset[i] não aloca corretamente o banco. Alguma outra ideia de como fazer isso? Obrigado. myDirDefault = 'C:/Users/Noely/Documents/Caio/Bolsa/R/Caio/'; myFileConfig = paste(myDirDefault, 'Config.cfg', sep=""); z = read.table(myFileConfig, header=FALSE, sep="=") vetZ = length(z[,1]) - 2 myDir = 1:vetZ for (i in 1:(vetZ)) { myDir[i] = paste(z[1,2],z[i+1,2]) } myDir = sub("/ ","/", myDir, ignore.case=T) dataset = 0 for (i in 1:vetZ) { dataset[i] = data.frame(read.csv(myDir[i],sep=";",header=TRUE)); #quando tiro o i, funciona corretamente com a última interação do i. } [image: Imagem inline 1] -- Caio Deutsch Bacharelado em Matemática Aplicada a Negócios Faculdade de Economia Administração e Contabilidade de Ribeirão Preto - FEA -RP Universidade de São Paulo - USP Cel: 16 8812 4376

1 0

ggplot2 Browser
by Eder David Borges da Silva 16 Jun '12

16 Jun '12

Pessoal, A um tempo postei um tópico sobre gráficos animados para web, encontrei um gráfico do mesmo modo que quero fazer o meu: http://www.babynamewizard.com/voyager#ms=true&sw=f&exact=false Porem pesquisando na ggplot2 e na gridSVG não encontrei nenhum código de exemplo, se alguém tiver feito algo assim, poderia me dar um help? Att

1 0

Árvore de decisão
by Vinicius Brito Rocha 16 Jun '12

16 Jun '12

Senhores, estou trabalhando com um problema de árvore de decisão onde eu obrigatoriamente tenho que começar meu processo decisório por variáveis específicas. Por exemplo: necessariamente meu modelo precisa ser resolvido na ordem em que "inputo" as variáveis independentes model=tree(Resposta~ Var1+Var2+Var3,data=base) Primeira variável Var1, segunda variável Var2, e terceira Var3. Tenho como fazer algo parecido com isso? -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

intervalo de confiança bayesiana
by Bernardo Rodrigues 16 Jun '12

16 Jun '12

Boa noite, Sabem como posso obter intervalo de confiança bayesiano utilizando os metodos de monte carlo sendo por exemplo a distribuição a priori uniforme? valeu

2 2

Função - agrupamento de grupos
by Luciano F. Sgarbi 15 Jun '12

15 Jun '12

#olá pessoal, #Estou querendo montar uma função que me mostre pontos que estão á uma distância menor que uma pré-estabelecida. #o script e explicações seguem abaixo: ############################ ##Conjunto de dados para testar a<-rep(c(1,4),each=4) b<-c(1,2,4,5,1,2,4,5) cc<-c("p1","p2","p3","p4","p5","p6","p7","p8") mat<-data.frame(cc,a,b) rownames(mat)<-cc coord<-mat[,c(2,3)] plot(coord,pch="",xlim=c(0,6),ylim=c(0,6));text(coord[,"a"],coord[,"b"],rownames(coord))#esse é o meu conjunto de pontos. #note que se eu for juntar os pontos que estão a uma distância menor ou igual á dois, formarei dois grupos. #um formado pelos pontos: p1, p2, p3 e p4 e, outro grupo com os pontos p5, p6, p7 e p8 #################### ##iniciando os cálculos de distância library(vegan) (mat.dist<-vegdist(coord,"euclidean"))#calculo de distância entre os pontos resu.col1<-numeric(0) resu.col2<-numeric(0) (distancias<-as.numeric(mat.dist)) locais<-rownames(coord) n.linhas<-nrow(coord) ##Agora vou criar uma tabela para comparar as distâncias par-a-par for(i in 2:n.linhas){ linhas<-c(i:n.linhas) resu.col2<-c(resu.col2,locais[linhas]) resu.col1<-c(resu.col1,rep(locais[i-1],length(linhas))) } (mat.resu<-data.frame(resu.col1,resu.col2,distancias))#tabela com as distâncias par-a-par dist.corte<-2##argumentos usado para estabelecer distância para juntar pontos (min.dists<-mat.resu[which(mat.resu[,"distancias"]<=dist.corte),])##tabela com todos os pontos que estão a uma distância menor, entre si, que á pré estabelecida. #note que se for juntando os pontos forma-se dois grupos # resu.col1 resu.col2 distancias *#1 p1 p2 1#p1 com p2* *#8 p2 p3 2#o p2 que estava juntado com o p1 junta-se ao p3* *#14 p3 p4 1**#o p3 que estava juntado com o p2 junta-se ao p4* *#23 p5 p6 1#idem* *#26 p6 p7 2* *#28 p7 p8 1* ######################################## *#Alguém, por favor, sabe me dizer como de agora em diante eu faço para saber quais que são os pontos que formam cada um dos grupos?* Atenciosamente, Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

1 0

Alteração de pontos em Biplot
by Marcelo Claro de Souza 15 Jun '12

15 Jun '12

Pessoal, Normalmente eu gero meus biplots de pca segundo script abaixo: require(vegan) dados.pca=rda(dados,scale=T) simbolos=c(rep(19,X),rep(2,X),.....) plot(dados.pca,display="species",xlab="PC1 = Y%",ylab="PC2 = Y%") points(dados.pca,display="sites",pch=simbolos) legend("topright",c("nome 1","nome 2"....),pch=c(19,2,....)) Entretanto dessa vez eu preciso colocar o nome das espécies ao invés dos símbolos. Como são várias plantas da mesma espécie não da p importar o arquivo utilizando o comando row.names. Como eu faço p criar um vetor com "texto" que possa substiur o vetor simbolos? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1