- 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