- R-br - DINF UFPR

Nome de Função
by David Feitosa 23 Out '14

23 Out '14

Caríssimos, boa noite. Existe modo de se obter o nome de uma função em tempo de execução? Explicando melhor: eu tenho uma função da forma minhaFuncao <-function(){ #codigo } e uma variável que recebe essa função: f <- minhaFuncao Eu gostaria de obter o nome dela para ser usado em um gráfico Em Java, eu sei que poderia usar a API de Reflection para, em tempo de execução, analisar a classe, seus métodos etc e obter o nome para uso posterior. Há como fazer o equivalente em R? No meu caso, bastaria saber qual o nome da função que a variável "f" carrega. Grato. Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

4 8

Teste de estimadores de regressão
by Adriano Borges Costa 23 Out '14

23 Out '14

Caros, Estou fazendo um exercício teórico de regressão para estimar o efeito de um programa a partir da base de dados disponível no código baixo. library(foreign) lalonde <- read.dta(" http://dl.dropboxusercontent.com/u/7453595/BD_LAB/lalonde.dta") #carrega a base de dados str(lalonde) # variáveis # As variáveis re74, re75 e re78 contêm informações sobre o rendimento para os anos de 1974, 1975 e 1978, respectivamente. # A variável de tratamento t indica se o indivíduo participou do programa (t=1) ou se não participou do programa (t=0). # As demais variáveis contém características dos trabalhadores Trata-se de um experimento aleatório e quero mostrar que há aleatoriedade no experimento. Assim, estou rodando a seguinte regressão: regressao1 <- lm(re78 ~ t, data=lalonde) summary(regressao1) E em seguida incluindo os controles em uma outra regressão e gostaria de mostrar que o estimador de t (tratamento) varia pouco (cerca de U$ 100,00). regressao2 <- lm(re78 ~ t+age+educ+nodegree+black+hisp+married, data=lalonde) summary(regressao2) Assim, gostaria de fazer um teste para verificar se o estimador de t na regressão 1 é estatisticamente diferente da mesmo estimador na regressão 2. Está certo eu fazer da seguinte forma? install.packages("car") library(car) linearHypothesis(regressao2,"t=1794.3", test = "F") Desde já obrigado pelas ajudas. At. -- Adriano Borges Costa

3 3

PNAD
by Vanessa Felix 23 Out '14

23 Out '14

Pessoal, Eu tô querendo usar os dados da PNAD 2013 mas não consigo usar o pacote IBGEPesq. Fica aparecendo esse erro aqui: library(IBGEPesq) Error: package ‘IBGEPesq’ was built before R 3.0.0: please re-install it Tem alguma opção que não seja reinstalar o R? Se alguém puder ajudar... att, -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestranda em Saúde Coletiva - IMS/UERJ (32) 8854-4881 / 9142-8428 (21) 97656-3448 -------------------------------------------------------------------------------------------

4 4

Comparação entre modelo linear e não linear
by Ari Clecius 22 Out '14

22 Out '14

Caros colegas, uso alguns modelos de isotermas de adsorção, para o modelo de Langmuir existem algumas linearizações, como faço para avaliar se o modelo não linear é melhor que o modelo linearizado? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #CRM ce=c(0.66450,0.92395,1.22510,1.53040,12.89600,32.44400,70.93000,135.54800) qe=c(0.65864,0.83815,1.09135,1.42049,11.57060,15.42198,16.41500,20.42020) #########Modelo linear#################### ce_inv=1/ce qe_inv=1/qe Langmuir_Linear=lm(qe_inv~ce_inv) summary(Langmuir_Linear) names(Langmuir_Linear) Langmuir_Linear$coefficients Langmuir_Linear$coefficients[1] qmax_L=1/Langmuir_Linear$coefficients[1] qmax_L kl_L=1/qmax*Langmuir_Linear$coefficients[2] kl_L qe_est_L=qmax_L* kl_L*ce/(1+ kl_L *ce) plot(ce,qe) lines(ce,qe_est_L,col="blue") ########################################## library(nlstools) adsorption=list(ce,qe) langmuir<-nls(qe~qmax*kl*ce/(1+kl*ce),adsorption,start=list(qmax=30,kl=0.02)) summary(langmuir) qe_Est_NL=21.40602* 0.07498 *ce/(1+ 0.07498 *ce) qe_Est_NL plot(ce,qe) lines(ce,qe_Est_NL,col="red") legend("top",legend=c("LangmuirExp","Langmuir_Est_NL"), lty=c(NA,1), col=c("black","red"),pch=c(3,NA),cex=0.6) ######################################## ############################################# plot(ce,qe) lines(ce,qe_est_L,col="blue") lines(ce,qe_Est_NL,col="red") legend("topleft",legend=c("LangmuirExp","Langmuir_Est_L","Langmuir_Est_NL"), lty=c(NA,1,1), col=c("black","blue","red"),pch=c(3,NA,NA),cex=0.6) #####################################################

1 0

Dúvida sobre análise de cluster em redes neurais
by Victor Eduardo 22 Out '14

22 Out '14

Boa noite, Ainda estou aprendendo sobre o assunto e estou com uma dúvida sobre clusterização em redes neurais utilizando o R. Eu encontrei dois pacotes que utilizam essa ferramenta, que é o nnet e kohonen, só que em nenhum deles, pelo que estava vendo, passa a formação dos clusters que tenho, ou seja, qual cluster o indivíduo pertence, como é visto no algoritmo de fuzzy c-means e no kmeans. É possível visualizar os clusters que cada indivíduo pertence? Existe algum comando que forneça essa informação dentro desses dois pacotes ou tem algum outro pacote com essa informação? Atenciosamente, Victor Eduardo

1 0

OFF- Livro Multivariada
by Andre Oliveira 22 Out '14

22 Out '14

Pessoal,qual livro de multivariada para um iniciante e em português. Achei estes .... MINGOTI, S. A. Análise de dados através de métodos de estatística multivariada: uma abordagem aplicada. Belo Horizonte: Editora UFMG, 2005. FERREIRA, Daniel Furtado.Estatística multivariada. 2. ed. Lavras : Editora UFLA, 2011 O que recomendam? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

3 2

Revolution R - Paralelização
by Jefferson Ferreira-Ferreira 21 Out '14

21 Out '14

Listeiros; Alguém aí utiliza o Revolution R Enterprise? Recém o baixei. Li alguns materiais pela internet mas não estou entendendo muito bem se ele paraleliza por "default" ou se preciso usar algo como ParallelR ou doSNOW. Alguém saberia me dizer ou indicar um bom site onde posso me informar melhor sobre isso? Obrigado. Abraços -- *Jefferson Ferreira-Ferreira* Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI Jefferson.ferreira(a)mamiraua.org.br *Instituto de Desenvolvimento Sustentável Mamirauá* Ministério da Ciência, Tecnologia e Inovação Telefone: +55 97 3343-9710 *Google Maps* - Mapas deste e-mail: Exibir mapa ampliado <https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> *Contatos particulares:* *(55) 9615-0100*

5 6

Tabela de frequência com amplitudes desproporcionais
by Flávio Fagundes 21 Out '14

21 Out '14

Prezados(as), Por um acaso existe alguma função semi-pronta no R que produz uma tabela de frequência para intervalo de classes com amplitudes de tamanhos desproporcionais? Imagino ser simples a programação (vou começar uma agora), mas pode ser que já exista. Por exemplo: Tenho os seguintes dados: dados<-c(0, 7 ,7 ,8 ,8 , 9, 9, 10, 10, 10, 11, 11, 12, 13, 14, 14, 16, 17, 18, 18, 19, 19, 20, 22, 25) A função fdt me retorna: > tb<-fdt(dados, k=4) > tb Class limits f rf rf(%) cf cf(%) [0,6.31) 1 0.04 4 1 4 [6.31,12.6) 12 0.48 48 13 52 [12.6,18.9) 7 0.28 28 20 80 [18.9,25.2) 5 0.20 20 25 100 Gostaria de obter a mesma tabela acima, mas para as seguintes classes 0 |---- 10 10 |---- 13 13 |---- 20 20 |---- 25 utilizando, algo como o parâmetro breaks=c(0,10,13,20) Obrigado Flávio

3 3

Re: [R-br] skater e deldir
by Rafael Tieppo 21 Out '14

21 Out '14

Elias, Pelo o que eu entendi, vc usou a funcao deldir(), no caso Voronoi para fazer os polígono da área. Se for isso acho que passei pelo mesmo processo. Vou passar algumas etapas e uma função que improvisei e no meu caso deu certo. Vou explicar os procedimentos: 1) com a funcao deldir() vc faz seu polígonos na área sss <- deldir( .....) entao vc tem o objeto sss 2) Para gerar uma lista do objeto "sss" use a funcao tile.list() sss_list <- tile.list(sss) 3) Agora desse objeto que é uma lista, vc tem que extrair o conjunto de polígonos para isso improvisei uma função, grave ela no seu ambiente: FUN_VOR_MAT <- function (A,TT) #TT é o arquivo, N o número do POL { AA <- 1 BB <- 2 NROW <- length (A[[TT]]$x) + 1 #add primeiro ponto VORONOI_MATRIX <- matrix(0,NROW,2) VORONOI_MATRIX[1:NROW,AA] <- A[[TT]]$x[1:NROW] VORONOI_MATRIX[1:NROW,BB] <- A[[TT]]$y[1:NROW] VORONOI_MATRIX[NROW,AA] <- A[[TT]]$x[1] #+ 1 add primeiro ponto VORONOI_MATRIX[NROW,BB] <- A[[TT]]$y[1] #+ 1 add primeiro ponto if(TT == NROW){ print("17/10/2014, rafaeltieppo(a)yahoo.com.br ")} return(VORONOI_MATRIX) } 4) Agora vc vai usar a funcao: sss_areas <- lapply(1:length(sss_list), function (x) ( FUN_VOR_MAT(sss_list,x)) ) 5) #Formando Polygons sss_pols <- lapply(1:length(sss_areas), function(x) Polygons(list(Polygon(sss_areas[[x]])), ID = paste("reg", x, sep = ""))) 6) Para converter em SpatialPolygins sss_SPATIAL_POLS <- as.SpatialPolygons.PolygonsList(sss_pols, proj4string=CRS(as.character(NA))) Se tudo correu bem é só comecar a usar o poly2nb(sss_SPATIAL_POLS) Boa sorte Rafael Tieppo

2 1

Juntar colunas no R
by Giselle Davi 20 Out '14

20 Out '14

Prezados,Preciso juntar colunas no R !!! Tenho um arquivo assim: D01 D02 D03 D04 D05 D06 D07 D08 D09 D10 F01 F02 F03 F04 F05 F06 F07 F08 F091 G G G G G G . G G G G G G G G G G . .2 C C C C C T C C C C C C C C C C C C C3 T T T T T T T T T T T T T T T T T T T4 G G C G G G C G G G G G C G G C G G G5 C C C C C C C C C C C C C C C C C C C6 T T T T T T T T T T T T T T T C T T T Quando peço para juntar as colunas o arquivo ( Parte do arquivo na verdade) passa a ser : D01 D01 D01 D02 D01 D03 D01 D04 D01 D05 D01 D06 D01 D07 D01 D08 D01 D091 1 1 1 1 1 1 0 1 12 2 2 2 2 2 3 2 2 23 2 2 2 2 2 2 2 2 24 1 1 3 1 1 1 3 1 15 2 2 2 2 2 2 2 2 26 2 2 2 2 2 2 2 2 2 D01 D10 D01 F01 D01 F02 D01 F03 D01 F04 D01 F05 D01 F06 D01 F07 D01 F081 1 1 1 1 1 1 1 1 02 2 2 2 2 2 2 2 2 23 2 2 2 2 2 2 2 2 24 1 1 1 3 1 1 3 1 15 2 2 2 2 2 2 2 2 26 2 2 2 2 2 2 3 2 2 MAS, eu preciso juntar as colunas dessa forma: | D01 | com | F01 | , | D02 | com | F01 | , | D03 | com | F01 | , | D04 | com | F01 | , | D05 | com | F01 | , | D06 | com | F01 | , | D07 | com | F01 | , | D08 | com | F01 | , | D09 | com | F01 | , | D10 | com | F01 | , | | D01 | com | F02 | , | D02 | com | F02 | , | D03 | com | F02 | , | D04 | com | F02 | , | D05 | com | F02 | , | D06 | com | F02 | , | D07 | com | F02 | , | D08 | com | F02 | , | D09 | com | F02 | , | D10 | com | F02 | , | | D01 | com | F03 | , | D02 | com | F03 | , | D03 | com | F03 | , | D04 | com | F03 | , | D05 | com | F03 | , | D06 | com | F03 | , | D07 | com | F03 | , | D08 | com | F03 | , | D09 | com | F03 | , | D10 | com | F03 | , | | D01 | com | F04 | , | D02 | com | F04 | , | D03 | com | F04 | , | D04 | com | F04 | , | D05 | com | F04 | , | D06 | com | F04 | , | D07 | com | F04 | , | D08 | com | F04 | , | D09 | com | F04 | , | D10 | com | F04 | , | | D01 | com | F05 | , | D02 | com | F05 | , | D03 | com | F05 | , | D04 | com | F05 | , | D05 | com | F05 | , | D06 | com | F05 | , | D07 | com | F05 | , | D08 | com | F05 | , | D09 | com | F05 | , | D10 | com | F05 | , | | D01 | com | F06 | , | D02 | com | F06 | , | D03 | com | F06 | , | D04 | com | F06 | , | D05 | com | F06 | , | D06 | com | F06 | , | D07 | com | F06 | , | D08 | com | F06 | , | D09 | com | F06 | , | D10 | com | F06 | , | | D01 | com | F07 | , | D02 | com | F07 | , | D03 | com | F07 | , | D04 | com | F07 | , | D05 | com | F07 | , | D06 | com | F07 | , | D07 | com | F07 | , | D08 | com | F07 | , | D09 | com | F07 | , | D10 | com | F07 | , | | D01 | com | F08 | , | D02 | com | F08 | , | D03 | com | F08 | , | D04 | com | F08 | , | D05 | com | F08 | , | D06 | com | F08 | , | D07 | com | F08 | , | D08 | com | F08 | , | D09 | com | F08 | , | D10 | com | F08 | , | | D01 | com | F09 | , | D02 | com | F09 | , | D03 | com | F09 | , | D04 | com | F09 | , | D05 | com | F09 | , | D06 | com | F09 | , | D07 | com | F09 | , | D08 | com | F09 | , | D09 | com | F09 | , | D10 | com | F09 | , | Alguém poderia me ajudar ??? Desde já agradeço a atenção !!!! Att, Giselle

5 6

Krigagem: Obter valores pontuais a partir da sup. interpolada
by Helder Gramacho 20 Out '14

20 Out '14

Bom dia pessoal, Para código reproduzível desta dúvida pode-se utilizar o exemplo disponível na página dos tutoriais do geoR: http://leg.ufpr.br/geoR/tutorials/ctc.R Ao final deste tutorial chegamos ao mapa da distribuição espacial obtido por krigagem. [image: Imagem inline 1] Gostaria de saber se é possível, obter o valor da variável para coordenadas específicas por exemplo na figura acima: x=5; y=8 x=25; y=10 x=15; y=15 Se é possível, como posso implementar no R. Desde já agradeço toda colaboração, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 4

skater e deldir
by Hélio Gallo Rocha 20 Out '14

20 Out '14

Caros da lista, Estou tentando implementar o uso de SKATER do pacote spdep. Gostaria de contar com a boa vontade dos listeiros para resolver esta dificuldade. Para iniciar o skater tenho de dividir minha área em poligonos. Para tal estou usando o pacote deldir. Abaixo um exemplo da própria página CRAN. require(deldir) set.seed(42) x <- runif(20) y <- runif(20) d <- deldir(x,y) (ttt <- triang.list(d) ) plot(ttt,border="red",showrect=TRUE) (sss <- tile.list(d)) plot(sss) plot(ttt,add=TRUE,border="blue") O objeto ttt lista os conjuntos de três pontos que formam o poligono Teria uma forma automatizada para gerar estes poligonos no formato shape, pois o CRM do skater que tenho utiliza o comando "readShapePoly" do spdep. Tentei exportar para o formato XLS, mas não foi possível pois o objeto ttt não é um data.frame Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 1

Integral
by andrebvs＠bol.com.br 20 Out '14

20 Out '14

4 4

Rodar uniroot várias vezes e salvar os valores encontrados em vetor
by Rodolfo Marcondes Silva Souza 19 Out '14

19 Out '14

Olá pessoal, No meu último post vocês me ajudaram a encontrar o valor de intersecção entre duas funções (Ponto de intersecção entre duas funções no R <http://r-br.2285057.n4.nabble.com/R-br-Ponto-de-interseccao-entre-duas-func…>). Bem, agora o problema é que eu preciso que a função uniroot seja executada várias vezes os valores da intersecção seja armazenado em um vetor para que seja utilizando posteriormente. O estranho que aproveitando tentando fazer isso nas respostas do meu último post funciona, mas quando eu tento fazer o mesmo no meu código "de verdade" não funcionada e aparece o seguinte erro: Error in uniroot(f = ffpsil, interval = c(-5, -0.05), Tar = Tar[i], Tl = Tl[i], : f() values at end points not of opposite sign In addition: Warning messages: 1: In if (is.na(f.lower)) stop("f.lower = f(lower) is NA") : the condition has length > 1 and only the first element will be used 2: In if (is.na(f.upper)) stop("f.upper = f(upper) is NA") : the condition has length > 1 and only the first element will be used Se eu adiciono valores para upper, lower, f.lower e f.upper para aparecer Error in uniroot(f = ffpsil, interval = c(-5, -0.05), Tar = Tar[i], Tl = Tl[i], : f() values at end points not of opposite sign Como o meu código está um pouco grande para postar aqui, estou colocando em um arquivo no dropbox (Código-Rodolfo <https://www.dropbox.com/s/oczeupi0mc4diob/Modelagem-Rodolfo.R?dl=0>) e se vocês puderem olharem e me ajudar a encontrar o erro vai me ajudar bastante. PS.: Como eu ainda estou escrevendo o código, ele está um pouco desorganizado, mas espero que vocês entendam. Desde já obrigado pela atenção e ajuda.

1 0

Re: [R-br] Ordenamento entre colunas (Giselle Davi)
by Rafael Tieppo 18 Out '14

18 Out '14

Segue (atençao com a última vírgula antes do último colchete): MAT <- matrix(c("A2", "A2", "A1", "B3", "B2", "B1"),3,2) MAT[order(MAT[,1], MAT[,2]),] Resposta: MAT[order(MAT[,1], MAT[,2]),] [,1] [,2] [1,] "A1" "B1" [2,] "A2" "B2" [3,] "A2" "B3" Rafael Tieppo On Saturday, October 18, 2014 12:00 PM, "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> wrote: 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: Problema (Vinicius Brito Rocha) 2. Re: Ordenamento entre colunas (Giselle Davi) 3. Re: Ordenamento entre colunas (Giselle Davi) 4. Re: Ordenamento entre colunas (Tiago Souza Marçal) 5. Estatísticas do R (Andre Oliveira) 6. Re: Estatísticas do R (David Feitosa) 7. Re: Estatísticas do R (Roney Fraga Souza) 8. Re: Estatísticas do R (David Feitosa) 9. Re: Estatísticas do R (Marco Nunes) ---------------------------------------------------------------------- Message: 1 Date: Fri, 17 Oct 2014 12:17:49 -0300 From: Vinicius Brito Rocha <viniciusbritor(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Problema Message-ID: <CAJSzTedNBBJRbX2GkT3qbZjd21UfUJWWM4eb7c4vbJbNjcSzQQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Algo deve estar errado em sua planilha, pois o excel comporta 1.000.000 de registros, com essa maquina q vc tem deveria ser bem tranquilo o import. Ou vc tem varias planilhas do excel e as está importando simultaneamente, de repente assim poderia criar um problema. De qualquer forma tente converter para CSV. Leia via read.csv2 Abs Em 17/10/2014 00:21, "Tayron Juliano" <tayron.js(a)hotmail.com> escreveu: > Boa noite, pessoal! > > > > Bem, sou iniciante no uso do R e estou com um problema. Estou trabalhando > com muitos dados que estou exportando do excel e estou utilizando o pacote > XLConnect e quando coloco o programa para executar ele inicia normalmente, > quanto chega em determinado trecho aparece o seguinte erro: > > *Error: cannot allocate vector size Gb* > > > Tenho 8 de memória ram, sistema de 64 btis e o windows 7 > > > O que devo fazer??? > > > > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141017/1a87b435/atta…> ------------------------------ Message: 2 Date: Fri, 17 Oct 2014 15:19:03 +0000 (UTC) From: Giselle Davi <giselle_davi(a)yahoo.com.br> To: Paulo Justiniano <paulojus(a)leg.ufpr.br>, "R-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Ordenamento entre colunas Message-ID: <1012586073.159238.1413559143873.JavaMail.yahoo(a)jws10620.mail.bf1.yahoo.com> Content-Type: text/plain; charset="utf-8" Infelizmente não funcionou !!!Era assim: head(dados) SNPName F9 D7 F5 F4 D2 F3 F2 D5 D6 F1 D3 F6 D4 D8 D1 F7 D9 D10 F81 S1_10045 . . G G G G G G G G G G G G G G G G C2 S1_222471 C C C C C C C C T C C C C C C C C C C3 S1_222473 T T T T T T T T T T T T T T T T T T T4 S1_239225 G C G G G C G G G G C C G G G G G G G5 S1_623890 C C C C C C C C C C C C C C C C C C C6 S1_645264 T T T T T T T T T T T C T T T T T T T E ficou assim:dialelo<- dados[,sort(colnames(dados))]> head(dialelo) D1 D10 D2 D3 D4 D5 D6 D7 D8 D9 F1 F2 F3 F4 F5 F6 F7 F8 F9 SNPName1 G G G G G G G . G G G G G G G G G C . S1_100452 C C C C C C T C C C C C C C C C C C C S1_2224713 T T T T T T T T T T T T T T T T T T T S1_2224734 G G G C G G G C G G G G C G G C G G G S1_2392255 C C C C C C C C C C C C C C C C C C C S1_6238906 T T T T T T T T T T T T T T T C T T T S1_645264 Continuo aceitando sugestões !!!! Desde já agradeço !!!! Att, Giselle Em Sexta-feira, 17 de Outubro de 2014 11:23, Paulo Justiniano <paulojus(a)leg.ufpr.br> escreveu: se M é sua matrix neste caso M[,sort(colnames(M)] vai ordenar como desejado On Fri, 17 Oct 2014, Giselle Davi wrote: > Prezados, > > Como eu consigo ordenar as seguintes colunas de uma matriz no R : A5, A3, A1, A4, A2, B2, B1, B3 ? > O objetivo é que elas fiquem assim: A1, A2, A3, A4, A5, B1, B2, B3 > > Utilizei o comando order, mas o ordenamento ocorre somente dentro das colunas e não entre como eu preciso. > > > Desde já agradeço a atenção. > > > Att, > > Giselle > > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141017/55756330/atta…> ------------------------------ Message: 3 Date: Fri, 17 Oct 2014 15:24:26 +0000 (UTC) From: Giselle Davi <giselle_davi(a)yahoo.com.br> To: Paulo Justiniano <paulojus(a)leg.ufpr.br>, "R-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Ordenamento entre colunas Message-ID: <897674063.158189.1413559466523.JavaMail.yahoo(a)jws10631.mail.bf1.yahoo.com> Content-Type: text/plain; charset="utf-8" Mas, foi quase !!!! Em Sexta-feira, 17 de Outubro de 2014 12:19, Giselle Davi <giselle_davi(a)yahoo.com.br> escreveu: Infelizmente não funcionou !!!Era assim: head(dados) SNPName F9 D7 F5 F4 D2 F3 F2 D5 D6 F1 D3 F6 D4 D8 D1 F7 D9 D10 F81 S1_10045 . . G G G G G G G G G G G G G G G G C2 S1_222471 C C C C C C C C T C C C C C C C C C C3 S1_222473 T T T T T T T T T T T T T T T T T T T4 S1_239225 G C G G G C G G G G C C G G G G G G G5 S1_623890 C C C C C C C C C C C C C C C C C C C6 S1_645264 T T T T T T T T T T T C T T T T T T T E ficou assim:dialelo<- dados[,sort(colnames(dados))]> head(dialelo) D1 D10 D2 D3 D4 D5 D6 D7 D8 D9 F1 F2 F3 F4 F5 F6 F7 F8 F9 SNPName1 G G G G G G G . G G G G G G G G G C . S1_100452 C C C C C C T C C C C C C C C C C C C S1_2224713 T T T T T T T T T T T T T T T T T T T S1_2224734 G G G C G G G C G G G G C G G C G G G S1_2392255 C C C C C C C C C C C C C C C C C C C S1_6238906 T T T T T T T T T T T T T T T C T T T S1_645264 Continuo aceitando sugestões !!!! Desde já agradeço !!!! Att, Giselle Em Sexta-feira, 17 de Outubro de 2014 11:23, Paulo Justiniano <paulojus(a)leg.ufpr.br> escreveu: se M é sua matrix neste caso M[,sort(colnames(M)] vai ordenar como desejado On Fri, 17 Oct 2014, Giselle Davi wrote: > Prezados, > > Como eu consigo ordenar as seguintes colunas de uma matriz no R : A5, A3, A1, A4, A2, B2, B1, B3 ? > O objetivo é que elas fiquem assim: A1, A2, A3, A4, A5, B1, B2, B3 > > Utilizei o comando order, mas o ordenamento ocorre somente dentro das colunas e não entre como eu preciso. > > > Desde já agradeço a atenção. > > > Att, > > Giselle > > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141017/685d79ba/atta…> ------------------------------ Message: 4 Date: Fri, 17 Oct 2014 18:27:43 +0300 From: Tiago Souza Marçal <tiagosouzamarcal(a)hotmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Ordenamento entre colunas Message-ID: <BLU170-W12420DF1838548A1CEB434AD9A80(a)phx.gbl> Content-Type: text/plain; charset="utf-8" Tenta deixar o cabeçalho da seguinte forma F09 D07 F05 F04 D02 F03 F02 D05 D06 F01 D03 F06 D04 D08 D01 F07 D09 D10 F08. Att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ################################################################# Date: Fri, 17 Oct 2014 15:19:03 +0000 From: giselle_davi(a)yahoo.com.br To: paulojus(a)leg.ufpr.br; r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Ordenamento entre colunas Infelizmente não funcionou !!!Era assim: head(dados) SNPName F9 D7 F5 F4 D2 F3 F2 D5 D6 F1 D3 F6 D4 D8 D1 F7 D9 D10 F81 S1_10045 . . G G G G G G G G G G G G G G G G C2 S1_222471 C C C C C C C C T C C C C C C C C C C3 S1_222473 T T T T T T T T T T T T T T T T T T T4 S1_239225 G C G G G C G G G G C C G G G G G G G5 S1_623890 C C C C C C C C C C C C C C C C C C C6 S1_645264 T T T T T T T T T T T C T T T T T T T E ficou assim:dialelo<- dados[,sort(colnames(dados))]> head(dialelo) D1 D10 D2 D3 D4 D5 D6 D7 D8 D9 F1 F2 F3 F4 F5 F6 F7 F8 F9 SNPName1 G G G G G G G . G G G G G G G G G C . S1_100452 C C C C C C T C C C C C C C C C C C C S1_2224713 T T T T T T T T T T T T T T T T T T T S1_2224734 G G G C G G G C G G G G C G G C G G G S1_2392255 C C C C C C C C C C C C C C C C C C C S1_6238906 T T T T T T T T T T T T T T T C T T T S1_645264 Continuo aceitando sugestões !!!! Desde já agradeço !!!! Att, Giselle Em Sexta-feira, 17 de Outubro de 2014 11:23, Paulo Justiniano <paulojus(a)leg.ufpr.br> escreveu: se M é sua matrix neste casoM[,sort(colnames(M)]vai ordenar como desejadoOn Fri, 17 Oct 2014, Giselle Davi wrote:> Prezados,> > Como eu consigo ordenar as seguintes colunas de uma matriz no R : A5, A3, A1, A4, A2, B2, B1, B3 ?> O objetivo é que elas fiquem assim: A1, A2, A3, A4, A5, B1, B2, B3> > Utilizei o comando order, mas o ordenamento ocorre somente dentro das colunas e não entre como eu preciso.> > > Desde já agradeço a atenção.> > > Att,> > Giselle> > > _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne?a c?digo m?nimo reproduz?vel. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141017/23a46026/atta…> ------------------------------ Message: 5 Date: Sat, 18 Oct 2014 01:15:23 +0000 (UTC) From: Andre Oliveira <andreolsouza(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Estatísticas do R Message-ID: <1499878210.167402.1413594923888.JavaMail.yahoo(a)jws100208.mail.ne1.yahoo.com> Content-Type: text/plain; charset="utf-8" Pessoal,preciso de estatísticas do R tais como nº de bibliotecas oficiais, pacotes que o R já vem na instalação. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141018/5fbdc08a/atta…> ------------------------------ Message: 6 Date: Fri, 17 Oct 2014 22:26:32 -0300 From: David Feitosa <contato(a)davidfeitosa.com> To: r-br(a)listas.c3sl.ufpr.br, Andre Oliveira <andreolsouza(a)yahoo.com.br> Subject: Re: [R-br] Estatísticas do R Message-ID: <CAL1rLtrb=x4JEKUcdJ0TH=8tZUTYLPm0U95kuRSptLzb1Tx0nA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Deve ter nos repositorios essa info, mas... mal pergunte, qual o motivo? Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("") Em 17 de outubro de 2014 22:15, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: > Pessoal, > preciso de estatísticas do R tais como nº de bibliotecas oficiais, > pacotes que o R já vem na instalação. > > > > > André Oliveira Souza. > Graduação em Matemática, mestrado em estatística aplicada.Instituto > Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141017/c5b9fcac/atta…> ------------------------------ Message: 7 Date: Sat, 18 Oct 2014 08:16:19 +0100 From: Roney Fraga Souza <roneyfraga(a)gmail.com> To: R-br <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Estatísticas do R Message-ID: <EC356409-FBC9-4B41-B174-7506AD1768ED(a)gmail.com> Content-Type: text/plain; charset=iso-8859-1 Algumas coisas básicas tem aqui: http://www.rdocumentation.org/ Att Roney ------------------------------ Message: 8 Date: Sat, 18 Oct 2014 08:53:32 -0300 From: David Feitosa <contato(a)davidfeitosa.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Estatísticas do R Message-ID: <CAL1rLtrDOXrdMkR23u+3OU0CnjpXy-D1CN9y9caCg4HcOU+oqA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Não conhecia esse. Obg, Roney. [?] Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("") 2014-10-18 4:16 GMT-03:00 Roney Fraga Souza <roneyfraga(a)gmail.com>: > Algumas coisas básicas tem aqui: http://www.rdocumentation.org/ > > Att > Roney > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141018/77bb3e1b/atta…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: 330.gif Tipo: image/gif Tamanho: 96 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141018/77bb3e1b/atta…> ------------------------------ Message: 9 Date: Sat, 18 Oct 2014 16:09:24 +0300 From: Marco Nunes <nunes.ma(a)outlook.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Estatísticas do R Message-ID: <SNT149-W25399CDF719FA691E3DBA09DA90(a)phx.gbl> Content-Type: text/plain; charset="utf-8" Também não conhecia não.Muito bom. From: contato(a)davidfeitosa.com Date: Sat, 18 Oct 2014 08:53:32 -0300 To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Estatísticas do R Não conhecia esse. Obg, Roney. Atenciosamente, David Feitosa (\_(\(=°;°)(("")("") 2014-10-18 4:16 GMT-03:00 Roney Fraga Souza <roneyfraga(a)gmail.com>: Algumas coisas básicas tem aqui: http://www.rdocumentation.org/ Att Roney _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne?a c?digo m?nimo reproduz?vel. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141018/c9b85ace/atta…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: 330.gif Tipo: image/gif Tamanho: 96 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141018/c9b85ace/atta…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 46, assunto 18 *****************************************

1 0

Ordenamento entre colunas
by Giselle Davi 17 Out '14

17 Out '14

Prezados, Como eu consigo ordenar as seguintes colunas de uma matriz no R : A5, A3, A1, A4, A2, B2, B1, B3 ?O objetivo é que elas fiquem assim: A1, A2, A3, A4, A5, B1, B2, B3 Utilizei o comando order, mas o ordenamento ocorre somente dentro das colunas e não entre como eu preciso. Desde já agradeço a atenção. Att, Giselle

5 8

Problema
by Tayron Juliano 17 Out '14

17 Out '14

Boa noite, pessoal! Bem, sou iniciante no uso do R e estou com um problema. Estou trabalhando com muitos dados que estou exportando do excel e estou utilizando o pacote XLConnect e quando coloco o programa para executar ele inicia normalmente, quanto chega em determinado trecho aparece o seguinte erro: Error: cannot allocate vector size Gb Tenho 8 de memória ram, sistema de 64 btis e o windows 7 O que devo fazer???

4 3

Multivariada no R
by Andre Oliveira 17 Out '14

17 Out '14

Pessoal,além da vegan, bpca qual outra biblioteca vocês me recomendariam para trabalhar com multivariada. obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

3 2

Ponto de intercepção entre duas funções no R
by Rodolfo Marcondes Silva Souza 16 Out '14

16 Out '14

Olá pessoal, Gostaira de saber se é possível calcular o ponto de intercepção entre das funções/equações não lineares no R? CMR: f1=function(x,y){ x^2 - 3 * y } f2=function(x,z){ 4 * z - x^2 } x=seq(1,10,0.5) y=3 z=15 plot(x,f1(x,y), type='l', col='red', ylim=c(0,100)) lines(x,f2(x,z), col='blue') Eu preciso encontrar um função que fazendo f1=f2 possa me retornar o valor em x onde ocorre a intercepção de f1 e f2. Alguém sabe como se resolve esse problema no R? Obrigado,

3 2

Re: [R-br] Ponto de intercepção entre duas funções no R
by Robert Iquiapaza 16 Out '14

16 Out '14

Isso resolve? Sds ------------------------------ f1=function(x,y){ x^2 - 3 * y } f2=function(x,z){ 4 * z - x^2 } f3=function(x,y,z){ f1(x,y)-f2(x,z) } x=seq(1,10,0.5) y=3 z=15 plot(x,f1(x,y), type='l', col='red', ylim=c(-50,100)) lines(x,f2(x,z), col='blue') lines(x,f3(x,y,z), col='green') grid() ?uniroot (xo=uniroot(f3, c(0,10),y=y,z=z,tol= 6e-6)$root) round(c(xo,f1(xo,y),f2(xo,z),f1(xo,y)-f2(xo,z)),4) ------------------------- -----Mensagem Original----- From: Rodolfo Marcondes Silva Souza Sent: Thursday, October 16, 2014 4:51 PM To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Ponto de intercepção entre duas funções no R Olá pessoal, Gostaira de saber se é possível calcular o ponto de intercepção entre das funções/equações não lineares no R? CMR: f1=function(x,y){ x^2 - 3 * y } f2=function(x,z){ 4 * z - x^2 } x=seq(1,10,0.5) y=3 z=15 plot(x,f1(x,y), type='l', col='red', ylim=c(0,100)) lines(x,f2(x,z), col='blue') Eu preciso encontrar um função que fazendo f1=f2 possa me retornar o valor em x onde ocorre a intercepção de f1 e f2. Alguém sabe como se resolve esse problema no R? Obrigado, _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

boxplots
by Edson Lira 16 Out '14

16 Out '14

Caros amigos da lista, tenho um banco com as variáveis: certos, fase(antes/depois) e ocupa(com 5 categorias). A minha pergunta é: Tem como fazer boxplots de acertos para as fases antes e depois por ocupação. A idéia é colocar a caixa antes e depois lado a lado. O banco para análise está no link abaixo: Download novo.csv Download novo.csv Download novo.csv free. File hosted free by DataFileHost. Visualizar em www.datafileho... Visualizado por Yahoo Um grande abraço. Edson Lira Estatístico Manaus-Amazonas

2 1

processo R inicia 'sozinho'
by Elias T. Krainski 16 Out '14

16 Out '14

Olá, alguém experimentou observar um processo R iniciar sozinho? Digamos, estou usando o notebook para coisas simples (sem usar muito processamento) e então o cooler começa a trabalhar. Quando vejo (top) há um ou mais processos R (user=root) consumindo bastante CPU. Isso ocorreu após eu instalar alguns bioconductor packages. Lembro que ao instalar, apareceu a opção de atualizar vários pacotes e eu aceitei.

2 2

mensagem de erro usando função glmer
by Ana Maria Nievas 16 Out '14

16 Out '14

Oi pessoal! Tentando reproduzir um script do livro do Zuur, me deparei com a seguinte mensagem de erro: O1.glmer <-glmer(NCalls~offset(LBroodSize)+SexParent*FoodTreatment+SexParent*ArrivalTime+(1 | fNest),data = an1, family = poisson) Mensagens de aviso perdidas: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, : Model is nearly unidentifiable: very large eigenvalue - Rescale variables? Alguém sabe como solucionar? -- Ana Maria Nievas Bióloga e Mestre em Ciências, com ênfase em Ecologia Aplicada (16) 98843-7744

3 3

Adicionar linhas em xyplot multipaneis
by Fernando Souza 16 Out '14

16 Out '14

Caros amigos estou criando um gráfico utilizando o pacote látice e a função xyplot(). O gráfico é composto de dois paneis um contendo dados de temperatura e humidade em função da data do ano. Eu tenho as variáveis de temperatura (minima,maxima e média) e humidade relativa (minima, maxima e média). Eu fiz os gráficos da temperatura e humidade média em função da data e gostaria de adicionar os valores de máximo e mínimo dentro dos respectivos paneis. Eu utilizei a função (exemplo) panel= panel.lines (x,dados$TEMPMAX) associado a funçao prepanel, no entanto o linha de temperatura máxima foi colocada tanto no panel de temperatura quanto no panel de humidade relativa. Fiz uma busca na internet e não encontrei nada que pudesse solucionar meu problema. Como eu posso fazer para adicionar as linhas TEMPMAX e TEMPMIN S somente no panel de temperatura e URAMAX E URAMIN somente no panel de humidade relativa? dados<- structure(list(DATA = structure(c(14485, 14486, 14487, 14488, 14489, 14490, 14491, 14492, 14493, 14494, 14495, 14496, 14497, 14498, 14499, 14500, 14501, 14502, 14503, 14504, 14505, 14506, 14507, 14508, 14509, 14510, 14511, 14512, 14513, 14514, 14515, 14516, 14517, 14518, 14519, 14520, 14523, 14524, 14525, 14526, 14527, 14528, 14529, 14530, 14531, 14532, 14533, 14534, 14535, 14536), class = "Date"), TEMPMIN = c(11.85, 11.85, 13.65, 17.9, 18.45, 18.6, 19.6, 19.2, 19.8, 17.8, 18.7, 18.7, 15.5, 15.55, 15.65, 16.4, 18.7, 18.15, 17.55, 17.85, 18.75, 20.95, 19.3, 19.35, 19.4, 18.8, 17.65, 17.75, 16.45, 16.95, 18.4, 21.1, 18.2, 18.05, 18.05, 18.55, 19.95, 20.7, 18.6, 18.6, 16.9, 17.15, 18.6, 20.8, 18.1, 19.85, 19.85, 20.45, 20.4, 19), TEMPMAX = c(28.05, 28.05, 28.3, 30, 30.9, 30.8, 30.35, 31.2, 27.25, 27.9, 31.3, 31.3, 28.6, 28.25, 29.7, 31.05, 31.95, 32.75, 30.7, 31.65, 32.85, 32.65, 33.2, 33.05, 24.7, 24.35, 24.95, 28, 28.15, 30.9, 33.1, 34.95, 26.65, 28.55, 28.55, 32.45, 31.7, 31.4, 23.4, 23.4, 27.95, 28.05, 30.45, 29.6, 31.2, 33.45, 33.45, 30.4, 28.6, 29.95), URAMIN = c(38, 38, 37, 40, 41.5, 43, 42.5, 41, 54, 61.5, 44.5, 44.5, 41.5, 41, 19.5, 19, 18, 19.5, 26.5, 19.5, 19.5, 20, 25.5, 29, 70, 71.5, 42.5, 50, 44, 31, 23, 15.65, 46.5, 46, 46, 32.5, 27, 31.5, 76, 76, 44.5, 44, 36.5, 45.5, 34, 25.5, 25.5, 49.5, 55.5, 52.5), URAMAX = c(89, 89, 89, 94, 91, 91, 96.5, 94.5, 99, 99, 95.5, 95.5, 91, 93, 80.5, 76.5, 74.5, 74.5, 75.5, 77, 76, 80, 88, 92, 94, 93.5, 94, 91.5, 85.5, 83.5, 79, 51.75, 77.5, 81, 81, 91.5, 86.5, 88.5, 92.5, 92.5, 84, 85, 82, 86, 80, 90, 90, 95, 90.5, 97), TEMPMED = c(19.95, 19.95, 20.975, 23.95, 24.675, 24.7, 24.975, 25.2, 23.525, 22.85, 25, 25, 22.05, 21.9, 22.675, 23.725, 25.325, 25.45, 24.125, 24.75, 25.8, 26.8, 26.25, 26.2, 22.05, 21.575, 21.3, 22.875, 22.3, 23.925, 25.75, 28.025, 22.425, 23.3, 23.3, 25.5, 25.825, 26.05, 21, 21, 22.425, 22.6, 24.525, 25.2, 24.65, 26.65, 26.65, 25.425, 24.5, 24.475), URAMED = c(63.5, 63.5, 63, 67, 66.25, 67, 69.5, 67.75, 76.5, 80.25, 70, 70, 66.25, 67, 50, 47.75, 46.25, 47, 51, 48.25, 47.75, 50, 56.75, 60.5, 82, 82.5, 68.25, 70.75, 64.75, 57.25, 51, 33.7, 62, 63.5, 63.5, 62, 56.75, 60, 84.25, 84.25, 64.25, 64.5, 59.25, 65.75, 57, 57.75, 57.75, 72.25, 73, 74.75), ITH = c(65.88425, 65.88425, 67.355125, 71.62425, 72.664625, 72.7005, 73.095125, 73.418, 71.014375, 70.04575, 73.131, 73.131, 68.89775, 68.6825, 69.794625, 71.301375, 73.597375, 73.77675, 71.875375, 72.77225, 74.279, 75.714, 74.92475, 74.853, 68.89775, 68.216125, 67.8215, 70.081625, 69.2565, 71.588375, 74.20725, 77.471875, 69.435875, 70.6915, 70.6915, 73.8485, 74.314875, 74.63775, 67.391, 67.391, 69.435875, 69.687, 72.449375, 73.418, 72.62875, 75.49875, 75.49875, 73.740875, 72.4135, 72.377625)), .Names = c("DATA", "TEMPMIN", "TEMPMAX", "URAMIN", "URAMAX", "TEMPMED", "URAMED", "ITH"), row.names = c(NA, 50L), class = "data.frame") library(lattice) #-ítulo dos gráficos humidade <- "Relative Humidity (%)" temperatura <- as.expression(expression( paste("Average Air Temperature(", degree*C, ")") )) #--Personalização do cabeçario dos gráficos # (NB the colour used is the default lattice strip background colour) my.strip <- function(which.given, which.panel, ...) { strip.labels <- c(temperatura,humidade) panel.rect(0,0, 1, 1, col="gray95", border=1) panel.text(x=0.5, y=0.5, adj=c(0.5, 0.55), cex=0.95, lab=strip.labels[which.panel[which.given]]) } #--Intervalo do eixo X xlim <- range(as.Date(dados$DATA,format="%d/%m/%y",tz="UTC")) #--marcadores para as linhas grid d <- seq(from=as.Date("2009-08-29"), to=as.Date("2009-12-04"), by=365/4) #--Definição das cores das linhas col.raw <- "gray0" #colset[2] } see note above col.smo <-"gray80" #colset[1] } col.lm <- "grey20" #--Criação dos gráficos xyplot(TEMPMED +URAMED ~DATA, data=dados, scales=list(y="free", rot=0), xlim=xlim, strip=my.strip, outer=TRUE, layout=c(1, 2, 1), ylab="",xlab=NULL, panel=function(x, y, ...) { panel.fill(col ="gray100") #define a cor de fundo do panel panel.xyplot(x, y, ..., type="l", col=col.raw, lwd=1.0) # x~y panel.loess(x, y, ..., col=col.smo, span=0.14, lwd=1) # smoothed data panel.grid(h=-1,v=-1) }, key=list(text=list(c("raw data", "smoothed curve")), col=c(col.raw, col.smo), lty=c(1, 1), columns=2, cex=0.95, lines=TRUE ), )

3 3

ADAPTAR UMA MATRIZ NO R
by Tiago Souza Marçal 16 Out '14

16 Out '14

Boa noite pessoal,estou trabalhando com o desenvolvimento matricial dos efeitos dialélicos dos modelos de Griffing no R.No entanto, estou tendo dificuldade para adaptar a matriz CR (mostrada abaixo) para a matriz CR1 (mostrada abaixo) que codifica os efeitos recíprocos. Qualquer sugestão é bem vinda.diall = structure(c(1L, 2L, 3L, 4L, 2L, 5L, 6L, 7L, 3L, 6L, 8L, 9L, 4L, 7L, 9L, 10L), .Label = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10"), class = "factor")mat <- model.matrix(~-1+diall)mat[,colSums(mat)==1] = 0CR = matCR1 <- structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), .Dim = c(16L, 10L), .Dimnames = list(c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16"), c("dia1", "dia2", "dia3", "dia4", "dia5", "dia6", "dia7", "dia8", "dia9", "dia10")), assign = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), contrasts = structure(list(dia = "contr.treatment"), .Names = "dia")) Att.Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

3 3

ecospat
by juan manuel Otalora 16 Out '14

16 Out '14

Colegas boa noite, eu estou usando o pacote ecospat para identificar similaridade de nichos, no script aparece este erro: *Error in validObject(.Object) : invalid class “SpatialPoints” object: bbox should never contain infinite values...* *9 stop(msg, ": ", errors, domain = NA) * *8 validObject(.Object) * *7 initialize(value, ...) * *6 initialize(value, ...) * *5 new("SpatialPoints", coords = coords, bbox = bbox, proj4string = proj4string) * *4 SpatialPoints(spr[, 1:2]) * *3 inherits(xy, "SpatialPoints") * *2 kernelUD(SpatialPoints(spr[, 1:2]), h = "href", grid = mask, * * kern = "bivnorm") * *1 ecospat.grid.clim.dyn(scores.clim, scores.clim, scores.sp2, R)* Alguém poderia me orientar sobre possíveis soluções. Muito obrigado -- *Juan Manuel Otálora Villamil*

1 0

OT: Visto
by Leonard de Assis 14 Out '14

14 Out '14

Bom dia, Para eu participar de algum evento nos USA, só meu visto B1 B2 basta? Ou tenho que pedir outro específico

2 1

Estimar assimetria e curtose
by Andre Oliveira 13 Out '14

13 Out '14

Pessoal alguém tem alguma sugestão melhor para estimar assimetria e curtose? Este pacote não bate a definição de Ferreira,D, F. 2005. curtose=m_4/s^4 assimetria= m_3/s^3 -------------------------------------------------------------------------------------------------------------------------------------------------------- x=c(2,3,,5,6,8,9) library(fBasics) basicStats(x) André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

3 3

[OT] Curso: Variacao Genomica - Campinas
by Benilton Carvalho 13 Out '14

13 Out '14

Prezados, Escrevo para anunciar um curso, a ser realizado em Campinas/SP, no qual serao trabalhadas estrategias para analise de variantes genomicas. http://www.ebi.ac.uk/training/workshop/exploring-variation-animal-genomes-b… O curso sera' liderado por Ian Streeter e Anja Thormann, ambos membros do time Ensembl do EBI-EMBL. O curso sera' ministrado em ingles e contara' com sessoes teoricas e praticas. Peco a gentileza de que aqueles que se inscreveram, tentem novamente. Se voces precisarem de maiores informacoes, por favor, nao hesitem em entrar em contato comigo. Tambem, se possivel, encaminhem para eventuais interessados. Novamente, meus pedidos de desculpas pelo envio de mais uma mensagem a este respeito e possivel duplicacao. Espero ve-los em breve em Campinas, benilton

1 0

Referência sobre o R
by Helder Gramacho 12 Out '14

12 Out '14

Olá, boa noite, Alguém possui este trabalho e poderia disponibilizar? *FRERY, A. C.* ; PERCIANO, T. <http://lattes.cnpq.br/6877649058117092> Elementos de Processamento e Análise Estatístico de Imagens usando R. In: Marco Antônio Piteri; José Carlos Rodrigues. (Org.). Fundamentos de Visão Computacional. Santa Cruz do Rio Pardo, SP: Unesp, 2011, v. , p. 1-63. Desde já agradeço, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 1

dados climáticos históricos
by Fernando Antonio de souza 11 Out '14

11 Out '14

Caros amigos, estou necessitando de dados climáticos (temperatura max e min, e umidade relativa do ar max e min) do entorno da cidade de pedro leopoldo para o ano de 2009 ( julho a dezembro). Alguém sabe onde consigo estas informações? Att

4 4

mapa: rosa dos ventos
by Aline Lipsky 10 Out '14

10 Out '14

Boa tarde caros colegas, Estou plotando um mapa com o seguinte script: library(maptools)library(maps)library(mapdata)library(mapproj)library(sp) SHP22 <- readShapeSpatial("hidrografia_poligonos_utm22.shp") #Plotar o mapax11() par(mar=c(2,2,2,4.2)) # ajustando as margens (bottom, left, top, right)plot(SHP22, col='light blue', axes=F, ylim=c(6441000,6530000),xlim=c(380000,410000)) #adicionando eixo no mapaaxis(4, las=1, cex.axis=0.9) #adicionando latitude do lado direito do mapaaxis(1, las=1, cex.axis=0.9) #adicionando longitude abaixo do mapa Mas quando tento plotar a rosa dos ventos com a seguinte função: compassRose(430000,6525000, rot=0, cex=0.65) aparece o erro: Error: could not find function "compassRose". Alguém saberia me dizer o que posso fazer? Grata desde já pela atenção, Att., Aline

3 3

Selecionar variáveis de uma matriz de correlação
by Rodrigo Muniz 10 Out '14

10 Out '14

Caros, Tenho uma dataframe com cerca de 1500 variáveis e gostaria de saber a correlação entre essas variáveis, então gerei uma matriz de correlação. Como muitas variáveis estão altamente correlacionadas, eu gostaria de selecionar apenas as variáveis que tenham, por exemplo, no máximo uma correlação de 0,3. Para isso eu utilizei a função findCorrelation do pacote caret. Pela descrição, a função foi feita para isso, selecionar as variáveis mais correlacionadas, considerando um limite de corte. Como eu gostaria das variáveis menos correlacionadas, eu inverti minha matriz de correlação, pelo seguinte comando: as.matrix(sqrt((matrix.cor-1)^2), desta forma, a função selecionaria as variáveis menos correlacionada. No entanto, o resultado da função sempre me retorna bandas altamente correlacionada. Por isso gostaria de fazer de uma forma alternativa. Segue um pequeno exemplo para ajudar nas considerações. Se alguém souber uma forma de programar para obter esse resultado, aceito um direcionamento. Desde já agradeço pela atenção. require(caret) filtro=c(0.7) #Este parâmetro foi criado para definir o liminar de corte, como a #matriz está inversa, um valor de 0,7 que dizer que a correlação máxima entre as #variáveis deve cer de 0,3. matrix.cor <- as.data.frame(cor(bandas))# Gerar a matriz de correlação banda.selecao=bandas[,findCorrelation(as.matrix(sqrt((matrix.cor-1)^2)), cutoff=as.numeric(filtro), verbose=FALSE)] # selecionar do dataframe #original apenas as variáveis com fraca correlação. bandas=structure(list(X1417 = c(0.161042, 0.186075, 0.202966, 0.211156, 0.172813, 0.167929, 0.175015, 0.200225, 0.178143, 0.182516, 0.162152, 0.169258, 0.14955, 0.143907, 0.161655, 0.164643, 0.157703, 0.160282, 0.177145, 0.169043), X1366 = c(0.352479, 0.376561, 0.372766, 0.39279, 0.36767, 0.368855, 0.361146, 0.390109, 0.373625, 0.383793, 0.371615, 0.375047, 0.363033, 0.341382, 0.372219, 0.375438, 0.359356, 0.372921, 0.375115, 0.36681), X469 = c(0.07749, 0.098493, 0.083867, 0.101454, 0.087098, 0.088587, 0.079879, 0.083187, 0.087261, 0.091197, 0.085591, 0.091496, 0.081759, 0.093462, 0.094817, 0.089399, 0.087421, 0.08966, 0.097024, 0.084687), X1549 = c(0.233143, 0.262691, 0.266292, 0.280087, 0.247662, 0.242888, 0.237161, 0.269429, 0.25326, 0.260258, 0.239942, 0.244771, 0.224855, 0.211901, 0.237555, 0.242471, 0.232943, 0.234706, 0.248407, 0.243808), X1424 = c(0.149535, 0.176297, 0.193145, 0.197999, 0.16312, 0.156738, 0.165102, 0.190595, 0.165941, 0.171179, 0.149601, 0.158143, 0.139422, 0.131351, 0.150353, 0.153668, 0.145654, 0.149609, 0.162426, 0.158941), X1957 = c(0.045334, 0.078347, 0.055058, 0.071466, 0.04449, 0.052894, 0.066396, 0.087972, 0.053841, 0.073682, 0.068851, 0.04738, 0.044667, 0.060964, 0.061176, 0.065722, 0.052965, 0.05974, 0.07101, 0.076109), X1848 = c(0.247498, 0.280584, 0.274104, 0.278935, 0.252584, 0.256848, 0.235115, 0.274966, 0.268281, 0.267282, 0.255826, 0.23335, 0.237109, 0.230515, 0.262255, 0.258836, 0.2387, 0.265159, 0.271541, 0.265838), X1537 = c(0.218761, 0.248475, 0.254953, 0.264756, 0.233748, 0.230395, 0.224644, 0.256149, 0.236666, 0.24462, 0.223847, 0.229429, 0.210074, 0.198872, 0.220081, 0.227677, 0.217309, 0.21968, 0.231893, 0.226282), X1621 = c(0.285036, 0.31607, 0.313035, 0.331176, 0.302131, 0.300669, 0.290153, 0.323815, 0.307194, 0.318355, 0.298653, 0.303098, 0.285467, 0.269884, 0.297272, 0.301428, 0.291555, 0.300164, 0.302502, 0.299489), X1602 = c(0.273662, 0.302662, 0.304277, 0.319991, 0.29179, 0.291713, 0.281382, 0.313998, 0.296813, 0.305874, 0.286391, 0.294462, 0.27532, 0.260575, 0.287217, 0.289802, 0.279721, 0.285657, 0.291856, 0.289307), X1865 = c(0.199668, 0.229541, 0.230852, 0.245059, 0.226876, 0.215756, 0.207998, 0.254432, 0.20988, 0.226653, 0.206045, 0.194866, 0.195796, 0.183774, 0.214694, 0.231647, 0.207081, 0.208682, 0.223025, 0.199142), X1476 = c(0.144097, 0.170145, 0.187837, 0.192702, 0.158221, 0.154393, 0.160075, 0.184072, 0.162413, 0.166724, 0.147404, 0.15375, 0.133576, 0.128008, 0.146151, 0.150793, 0.143929, 0.141969, 0.158288, 0.154279), X2034 = c(0.081108, 0.110882, 0.100852, 0.109914, 0.109893, 0.06942, 0.128299, 0.139094, 0.10447, 0.125726, 0.118003, 0.075258, 0.090784, 0.055216, 0.102326, 0.101315, 0.08443, 0.123796, 0.104067, 0.102786), X1460 = c(0.136055, 0.161915, 0.177241, 0.184844, 0.149423, 0.143712, 0.150499, 0.173385, 0.150807, 0.156626, 0.138878, 0.145282, 0.125618, 0.120778, 0.138968, 0.139961, 0.133014, 0.132243, 0.150618, 0.144916), X1236 = c(0.41637, 0.440245, 0.422308, 0.453948, 0.428499, 0.436163, 0.427196, 0.454866, 0.433779, 0.451874, 0.447265, 0.44747, 0.44684, 0.422101, 0.451307, 0.452207, 0.436043, 0.454355, 0.44641, 0.437131), X1913 = c(0.049649, 0.06384, 0.080675, 0.0761, 0.062756, 0.056635, 0.074517, 0.084222, 0.054835, 0.05488, 0.06389, 0.058421, 0.054827, 0.050249, 0.063544, 0.066275, 0.057604, 0.072578, 0.070545, 0.061034), X1885 = c(0.105368, 0.136086, 0.161411, 0.146339, 0.115904, 0.121671, 0.149059, 0.145293, 0.125417, 0.145147, 0.124081, 0.112018, 0.101411, 0.096761, 0.122112, 0.114244, 0.105535, 0.122283, 0.140519, 0.121402), X492 = c(0.080821, 0.102833, 0.088743, 0.10809, 0.093344, 0.096766, 0.083556, 0.089436, 0.091601, 0.100154, 0.089925, 0.096294, 0.087353, 0.097834, 0.102614, 0.09692, 0.096883, 0.094752, 0.102917, 0.089361), X1916 = c(0.037402, 0.063136, 0.08053, 0.071535, 0.059838, 0.0555, 0.060523, 0.079581, 0.044433, 0.055137, 0.055831, 0.054455, 0.054067, 0.055332, 0.062323, 0.0651, 0.05896, 0.075877, 0.069466, 0.064714), X1249 = c(0.416616, 0.440854, 0.423568, 0.455739, 0.428821, 0.436904, 0.42868, 0.455141, 0.432762, 0.45234, 0.447037, 0.448069, 0.445743, 0.420802, 0.455476, 0.453444, 0.43817, 0.454327, 0.444599, 0.437981)), .Names = c("X1417", "X1366", "X469", "X1549", "X1424", "X1957", "X1848", "X1537", "X1621", "X1602", "X1865", "X1476", "X2034", "X1460", "X1236", "X1913", "X1885", "X492", "X1916", "X1249"), row.names = c(NA, 20L), class = "data.frame") -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) E-mail - muniz.ra(a)usp.br Cel (19) 98300-4333 (Pessoal)

3 5

Função K de Ripley
by Fernando Elias 10 Out '14

10 Out '14

Olá, algum de vocês tem o script da análise PAIR CORRELATION FUNCTION para avaliar padrões espaciais. Att. -- Fernando Elias Biólogo - CRBio 1a Região <http://www.crbiodigital.com.br/01/fernandoelias> Bolsista CAPES - Programa de Pós-graduação em Ecologia e Conservação <http://www.unemat.br/prppg/ppgec/> Br 158, Km 655 - Antiga FAB - Nova Xavantina - MT. CEP 78690-000 Caixa postal 08 CV - Currículo Lattes <http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4406350E6> CV - Google Acadêmico <http://scholar.google.com.br/citations?hl=pt-BR&user=RJwHMvQAAAAJ&view_op=l…> CV - Linkedin <http://br.linkedin.com/pub/fernando-elias/a6/36/459/> "Nada é difícil se for dividido em pequenas partes" - Henry Ford

1 0

Texto em vários gráficos
by Roberto Guevara 09 Out '14

09 Out '14

Caros usuários de R, tenho um problema quanto a acrescentar texto em gráficos que já tentei resolver pelo google e não consegui. Preciso gerar uma figura como oito gráficos de dispersão e para cada uma quero acrescentar um pequeno texto no canto contendo algumas informações. Gostaria que ficasse padronizado, com todos os textos exatamente na mesma localização correspondente em cada gráfico. Como posso fazer isso tendo em vista que as variáveis de cada gráfico são distintas e, portanto, as escalas dos eixos também o são. Há uma espécie de escala padrão que eu poderia usar na função text em todos os gráficos? Obrigado pela atenção -- Roberto Guevara Ferreira Lima Biólogo - CRBio 6ª Região: 73146/06-D Mestre em Zoologia Laboratório de Ecologia e Zoologia de Vertebrados Instituto de Ciências Biológicas Universidade Federal do Pará Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 Currículo Lattes: http://lattes.cnpq.br/7363382159007600

3 2

Re: [R-br] RES: doBy - summaryBy
by Andre Oliveira 09 Out '14

09 Out '14

Então estou no Ubuntu 14.04. Sempre deu certo esta análise não sei extamente o que esta ocorrendo me parece algum problema interno mesmo. E eu não sei resolver. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES Em Quarta-feira, 8 de Outubro de 2014 21:38, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: Eu copiei e colei seu script e rodei! Tenho um Windows XP e o R 3.1. Veja: dados=data.frame(peso,glicose) glicose, peso > require(doBy) Loading required package: doBy Loading required package: survival Loading required package: splines Loading required package: MASS > summaryBy(glicose~peso,data=dados,FUN=c(mean,sd)) peso glicose.mean glicose.sd 1 Abaixodopeso 77.57143 14.06956 2 Acimadopeso 114.33333 39.50633 3 Obeso 141.20000 62.83072 4 Pesopadrão 79.50000 14.34064 > Leonardo, obrigado mas deu erros! Erro em summaryBy(glicose ~ peso, data = dados, FUN = c(mean, sd), FUN = function(x) cbind(mean(x), : argumento formal "FUN" corresponde a múltiplos argumentos especificados Seria o motivo como entrei com meus dados? peso=scan(what= "as.character") Abaixodopeso Pesopadrão Acimadopeso Acimadopeso Obeso Acimadopeso Pesopadrão Abaixodopeso Abaixodopeso Acimadopeso Acimadopeso Obeso Pesopadrão Pesopadrão Pesopadrão Abaixodopeso Pesopadrão Obeso Pesopadrão Obeso Acimadopeso Abaixodopeso Pesopadrão Pesopadrão Abaixodopeso Pesopadrão Pesopadrão Obeso Acimadopeso Pesopadrão Acimadopeso Acimadopeso Abaixodopeso Pesopadrão Pesopadrão glicose=scan() 60 62 130 120 219 80 85 82 90 105 100 96 96 60 61 61 100 101 85 200 210 74 76 85 98 74 89 90 105 92 78 101 78 90 58 dados=data.frame(peso,glicose) dados str(dados) head(dados) attach(dados) require(doBy) summaryBy(glicose~peso,data=dados,FUN=c(mean,sd)) Erro em FUN(X[[1L]], ...) : não foi possível encontrar a função "currFUN" André Oliveira Souza. ________________________________ Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus está ativa.

3 3

doBy - summaryBy
by Andre Oliveira 08 Out '14

08 Out '14

Pessoal me socorre por favor, veja meu código que originou no erro esta abaixo. Alguém saberia dizer o motivo deste erro? Sempre usei para estratificar médias por tratamento e nunca me deparei com este erro. dados=data.frame(peso,glicose) dados str(dados) 'data.frame': 35 obs. of 2 variables: $ peso : Factor w/ 4 levels "Abaixodopeso",..: 1 4 2 2 3 2 4 1 1 2 ... $ glicose: num 60 62 130 120 219 80 85 82 90 105 ... require(doBy) summaryBy(glicose~peso,data=dados,FUN=c(mean,sd)) Erro em FUN(X[[1L]], ...) : não foi possível encontrar a função "currFUN" André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

3 4

Variáveis contínuas explicando Variáveis categóricas
by Jefferson Ferreira-Ferreira 08 Out '14

08 Out '14

Prezados; Uma dúvida básica de um iniciante em análises estatísticas. Tenho uma série de variáveis categóricas binárias e duas variáveis contínuas. Eu gostaria de saber o quando minhas variáveis categóricas explicam minhas duas variáveis contínuas. Por exemplo: c1 c2 c3 c4 x y 0 1 1 1 931 623 0 1 0 0 163 259 1 0 1 0 690 306 1 0 1 0 690 306 A pergunta é: o quanto a combinação das variáveis c explicam os valores de x e y? Ou, qual a correlação entre as variáveis c e as variáveis x e y? Podem me dar alguma ideia de análises possíveis? Obrigado. -- *Jefferson Ferreira-Ferreira* Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI Jefferson.ferreira(a)mamiraua.org.br *Instituto de Desenvolvimento Sustentável Mamirauá* Ministério da Ciência, Tecnologia e Inovação Telefone: +55 97 3343-9710 *Google Maps* - Mapas deste e-mail: Exibir mapa ampliado <https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> *Contatos particulares:* *(55) 9615-0100*

9 9

Função kruskal do pacote agricolae
by Jose Claudio Faria 08 Out '14

08 Out '14

Olá pessoal, Pacote insuficientemente documentado é um prato cheio para muitas dúvidas: basta ver o histórico de discussão em nosso grupos e no R-help. Uso muito pouco estatística livre de distribuição, mas como trabalho (também) com análise de dados, acabei caindo em um caso onde um teste de comparação múltiplas livre de distribuição está sendo necessário. Após estudar a teoria do teste Kruskal-Wallis pelo livro do Siegel (final de semana chuvoso...), fui buscar (no R) pacotes que continha sua implementação mais eleborada. Encontrei 2: 1- agricolae 2- pgirmess pgirmess não instala (no debian) devido a ausência do pacote rgdal. Adicionalmente, agricolae é mais conhecido ( o que poderia ajudar). Então fui estudar a implementação do teste no pacote agricolae. A primeira coisa que você observa é a insuficiência de documentação da função kruskal. Bom, então o jeito foi fazer um debug (ninguém merece) para ver suas entranhas. Observei que ela calcula a média dos postos de cada grupo, entre outras informações necessárias a posterior discriminação dos grupos, e chama a função "order.group" para tal tarefa. "order.group" é também insuficientemente documentada. Fora a parte funcional da função "order.group", de útil se retira que é usada nos testes: HSD, LSD, Kruskal-Wallis, Friedman ou Durbin. A função "order.group" recebe o argumento "p.adj" da função kruskal. E existem algumas opções: optei pela de bonferroni. Não consegui (pelo debug) identificar qual teste a função "order.group" usa. Vi uns posts antigos aqui na lista se referindo ao teste paramétrico "t". É isto mesmo? Bom, para socializar um pouco do que andei observando, o resultado do teste não é fornecido na mesma escala da variável em questão, mas em escala das médias dos postos de cada grupo. Se aplico a correção de bonferroni os resultados (classificação dos grupos) se assemelham bastante (em muitos casos idênticos) aos de uma ANOVA seguido do teste de Tukey. Se não aplico a correção se assemelham bastante (idem) aos do teste de Duncan. Fiz isso apenas para ter uma noção da dif. mínima significativa em termos comparativos. Escrevi esta semana para o autor e mantenedor do pacote (Felipe de Mendiburu) solicitando informações sobre qual teste a função "order.group" mas ele ainda não respondeu. Alguém conhece (com certeza) qual teste é usado na função "order.group"? É realmente o teste "t" aplicado a média dos postos? Não seria o de Mann-Whitney, ou outro? Ab, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9966.9100 - VIVO 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 1

Funções kruskal e order.group do pacote agricolae
by Jose Claudio Faria 07 Out '14

07 Out '14

Olá pessoal, Estou reenviando a mensagem devido a: 1- Erros grosseiros de linguagem detectados na primeira versão (era de madrugada e estava com sono); 2- Ter detalhado melhor alguns pontos; 3- Pela relevância do assunto, que ele fique como registro de possíveis opiniões. Pacote insuficientemente documentado é um prato cheio para muitas dúvidas e que acabam se refletem nas listas de discussão: basta ver o histórico de discussão sobre kruskal em nosso grupo e no R-help. Uso muito pouco estatística livre de distribuição, mas como trabalho (também) com análise de dados, acabei caindo em um caso onde um teste de comparação múltiplas (livre de distribuição) está sendo necessário. Após estudar a teoria do teste Kruskal-Wallis pelo livro do Siegel (final de semana chuvoso...), fui buscar (no R) pacotes que continham sua implementação mais elaborada. Encontrei 2: 1- agricolae 2- pgirmess pgirmess não instala (no debian) devido a ausência do pacote rgdal. Até seria possível armengando, mas adicionalmente, agricolae é mais conhecido (o que poderia ajudar em caso de dúvidas). Então fui estudar a implementação do teste no pacote agricolae. A primeira coisa que você observa é a insuficiência de documentação da função kruskal. Bom, então o jeito foi fazer um debug (ninguém merece) para ver suas entranhas. Observei que ela calcula a média dos postos de cada grupo (entre outras informações necessárias a posterior discriminação dos grupos) e chama a função "order.group" para tal tarefa. "order.group" é também insuficientemente documentada. Fora a parte funcional da função "order.group", de útil se retira apenas que é usada nos testes: HSD, LSD, Kruskal-Wallis, Friedman ou Durbin. A função "order.group" recebe o argumento "p.adj" da função kruskal. E existem algumas opções: optei pela de bonferroni. Não consegui (pelo debug) identificar qual teste a função "order.group" usa. Vi uns posts antigos aqui na lista se referindo ao teste paramétrico "t". Seria isto mesmo? Bom, para socializar um pouco do que andei observando, o resultado do teste não é fornecido na mesma escala da variável em questão, mas em escala das médias dos postos de cada grupo. Se aplico a correção de bonferroni os resultados (classificação dos grupos) se assemelham bastante (em muitas situações idênticos) aos da correlata paramétrica ANOVA seguido do teste de Tukey. Se não aplico a correção se assemelham bastante (idem) aos do teste de Duncan. Fiz isso apenas para ter uma noção do poder em termos comparativos. Escrevi esta semana para o autor e mantenedor do pacote (Felipe de Mendiburu) solicitando informações sobre qual teste a função "order.group" usa mas ele ainda não respondeu. Alguém conhece (com certeza) qual teste é usado na função "order.group"? É realmente o teste "t" aplicado a média dos postos? Acho estranho evocar um teste paramétrico numa metodologia não paramétrica. Adicionalmente, achei que a função faz cálculos demasiados para ser o teste "t". Não seria o de Mann-Whitney ou outro? Ab, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9966.9100 - VIVO 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

gráfico de medianas
by Diego Bilski 03 Out '14

03 Out '14

Caros, tenho um resultado da seguinte forma: structure(list(phi = c(0.15, 0.164, 0.042, 0.033, 0.161, 0.027, 0.01, 0.01, 0.208, 0.093), lambda = c(0.648, 0.654, 0.726, 0.728, 0.64, 0.696, 0.756, 0.748, 0.605, 0.678), gamma = c(0.202, 0.183, 0.232, 0.239, 0.198, 0.277, 0.234, 0.242, 0.186, 0.228), phi = c(0.08, 0.022, 0.081, 0.138, 0.139, 0.01, 0.01, 0.01, 0.135, 0.043), lambda = c(0.714, 0.79, 0.689, 0.665, 0.644, 0.75, 0.753, 0.762, 0.68, 0.714), gamma = c(0.206, 0.189, 0.23, 0.197, 0.216, 0.24, 0.237, 0.228, 0.184, 0.242)), .Names = c("phi", "lambda", "gamma", "phi", "lambda", "gamma"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10"), class = "data.frame") onde lambda, phi e gamma são parâmetros estimados para 5 variáveis, e a soma deles é um para todas as variáveis. Alguém pode me sugerir um gráfico que represente a mediana e intervalos para as 5 variáveis na mesma figura? Pensei num boxplot, mas queria ressaltar a ideia de que os 3 parâmetros somam um para cada variável, e fica bem poluído com 15 caixas. Um stacked barplot representa bem, mas fica muito feio com as barras de erro. obrigado Diego

2 1

função de interpolação
by Fernando Souza 03 Out '14

03 Out '14

Caros amigos Estou necessitando faze a interpolação de algumas pesagens tomadas em diferentes animais. Eu preciso da interpolação feita para cada animal separadamente e o intervalo entre medidas não é fixo. Eu estou utilizando a função approx() no entanto devido ao número de animais utilizados fica muito dispendioso fazer fazer esta interpolação uma a uma. Por isso gostaria de uma função onde um pudesse automatizar este procedimento. set.seed(20) dados<-data.frame(ANIMAL=factor(rep(1:5,each=4)),Peso=rnorm(20,30,4), data=sample(seq(as.Date("01/04/2009",'%d/%m/%Y'), as.Date("30/04/2009",'%d/%m/%Y'),length.out=30),20),day=1:20) Estou tentando fazer uma função que estime os pontos da interpolação agrupados por Animal . Entretanto tenho pouco conhecimento em programação para fazer isso. Tenho tentado fazer isso, sem muito sucesso. Alguém poderia me ajudar? Abraços aprendendo<-function(dados){ niveis<-levels(dados$ANIMAL) dia<-diff(dados$data) for(i in min(niveis):max(niveis)){ b<- approx(dados[dados$ANIMAL==as.numeric(i),]$Peso,dados[dados$ANIMAL==as.numeric(i),]$data),n=15) } return(b) }

3 10

Ajuda com apply
by Wagner Bonat 02 Out '14

02 Out '14

Caros o seguinte tenho um pequeno problema para manipular arrays dentro de um apply olhe o seguinte exemplo. Eu tenho este pequeno array my.array <- array(1:18, dim = c(3,3,2)) Quero aplicar esta função a cada elemento do array e obter outro array do mesmo tamanho igualzinho mas com informacao diferente.Veja so ID.median <- function(x, median,n.array){ if(n.array == 1){ saida <- ifelse(x < median, 0, 1)} if(n.array == 2){ saida <- ifelse(x < median, 10,20)} return(saida) } É claro que aqui só tem dois arrays, mas o fato relevante é que a minha depende de qual o número de array que eu estou. Agora quero usar o apply pra aplicar minha função a cada componente do meu array. apply(my.array, 1:3, ID.median, median = 10, n.arrray = ???) Como fazer isso ? alguma outra idéia ?? Valeu -- Wagner Hugo Bonat LEG - Laboratório de Estatística e Geoinformação UFPR - Universidade Federal do Paraná

2 1

Solicitação de texto
by Daniel Marcelino 02 Out '14

02 Out '14

Alguém da lista tem o texto abaixo? Não consigo encontra-lo de jeito nenhum. Deville, Jean-Claude (1991) A Theory of Quota Surveys. Suruey Methodology. Vol. 17, N° 2. pp. 163-181. Obrigado, Daniel

1 0

Parte de um número
by Vanessa Felix 02 Out '14

02 Out '14

Pessoal, Alguém pode me ajudar me informando se naquela função read.fwf tem como eu voltar e ler os primeiros dígitos de novo. Por exemplo, eu tenho um dado 20131100001500129856 20131100001500286416 20131100001500356871 20131100001500482267 20131100001500569716 20131100001500621563 20131100001500797165 Eu queria ler do 1-15 e do 5- 6, alguém tem alguma ideia pra ajudar? Eu posso ler do 1-15 e ter um vetor numérico, mas depois não sei extrair dele a quinta e sexta posição. 20131100001500797165 Desde já obrigada, -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestranda em Saúde Coletiva - IMS/UERJ (32) 8854-4881 / 9142-8428 (21) 97656-3448 -------------------------------------------------------------------------------------------

2 2

[OFF] Oportunidade de Estágio na BM&FBovespa
by Wilson Nascimento De Freitas 02 Out '14

02 Out '14

Pessoal, Estamos com uma vaga para estagiário na BM&FBovespa em modelagem de Pricing na Diretoria de Administração de Risco. Utilizamos R em diversas de nossas análises, de forma que familiaridade com R ou vontade de aprender sejam grandes diferenciais. Segue abaixo a descrição da vaga. Wilson Freitas Oportunidade de Estágio em Modelagem na BM&FBOVESPA A Gerência de Modelagem de Preços da BM&FBovespa procura jovens inteligentes e talentosos para compor o seu time por meio do programa de estágio. Descrição da área A Gerência de Modelagem de Preços desenvolve modelos de apreçamento de derivativos e metodologias para construção de dados de mercado, como curvas de juros, superfícies de volatilidade, dentre outros. Os preços e dados de mercado produzidos por estes modelos e metodologias são utilizados pela bolsa na sua gestão de risco, no cálculo dos valores de liquidação dos contratos e também nos sistemas de negociação, além de serem utilizados por diversas instituições para marcação de seus ativos. A área tem um forte viés de pesquisa quantitativa e incentiva a participação de todos nas soluções dos problemas. Responsabilidades ̶ Atuação conjunta com a equipe de modelagem de preços. ̶ Coleta de informações diárias e intra-diárias para análise de dados e simulações. ̶ Implementação e execução de relatórios de validação de modelos e validação de dados. ̶ Desenvolvimento e apresentação de estudos envolvendo técnicas utilizadas na modelagem de preços. Requisitos necessários ̶ Estar cursando graduação em áreas de ciências exatas: Física, Computação, Matemática, Estatística, Engenharia ou áreas relacionadas. ̶ Gostar de programação, análise de dados e simulação. ̶ Disponibilidade de 30 horas semanais (horário flexível). Requisitos desejáveis (não necessários) ̶ Conhecimento de linguagens de programação científica: R, Python, Matlab ̶ Conhecimento de Excel/VBA. ̶ Leitura em Inglês. Benefícios ̶ Bolsa auxílio compatível com o mercado. ̶ Vale alimentação, vale refeição, vale transporte, seguro de vida, plano de saúde, plano odontológico e férias remuneradas. Enviar currículo para: wfreitas(a)bvmf.com.br<mailto:wfreitas@bvmf.com.br>, com Assunto: “Vaga de Estagiário”. ' ________________________________ Esta mensagem pode conter informação confidencial e/ou privilegiada. Se você não for o destinatário ou a pessoa autorizada a receber esta mensagem, não deverá utilizar, copiar, alterar, divulgar a informação nela contida ou tomar qualquer ação baseada nessas informações. Se você recebeu esta mensagem por engano, por favor avise imediatamente o remetente, respondendo o e-mail e em seguida apague-o. Agradecemos sua cooperação. This message may contain confidential and/or privileged information. If you are not the addressee or authorized to receive this for the addressee, you must not use, copy, disclose, change, take any action based on this message or any information herein. If you have received this message in error, please advise the sender immediately by reply e-mail and delete this message. Thank you for your cooperation. '

1 0

Coeficiente de correlação em curvas de calibração analítica
by Ari Clecius 02 Out '14

02 Out '14

O coeficiente de correlação pode ser aplicado a dados de calibração química, já que os valores de x são escolhidos e não aleatórios? Li em um artigo que o coeficiente de correlação é aplicável apenas a variáveis aleatórias, como proceder? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #Curva de calibraçao em química analítica y=c(31527,26523,21318,15454,8902,6303,2530) x=c(50,40,30,20,10,5,2.5) lm.1=lm(y~x) library(chemCal) calplot(lm.1)

3 4

Criar função no R com if e for
by Rodolfo Souza 01 Out '14

01 Out '14

Olá pessoal, gostaria mais uma vez da ajuda de vocês. Eu estou precisando criar uma função no R que tem três condições e após testar todas as condições o programa me retorne um resultado na forma de vetor para cada valor (input) na função. Segue o código que consegui fazer até agora: x=c(-5,-4,-3,-2,-1,0,1,2,3,4,5);x fr=function(x){ r=numeric(length = length(x)); for(i in seq(along=x)){ if(x<-3){r[i]=0} else if(x>=-3 && x<=2){r[i]=3*x} else{r[i]=1} } return(r) } rs=fr(x);rs O resultado saí: 0 0 0 0 0 0 0 0 0 0 0 Se eu tenho só um valor no vetor e não utilizo o "for" consigo resolver, mas se o vetor tem mais de um valor dá erro. Desde já grato pela atenção/ajuda. -- Rodolfo Marcondes Silva Souza Engenheiro Agrônomo Mestre em Produção Vegetal - UFRPE/UAST Doutorando em Tecnologias Energéticas e Nucleares - UFPE/DEN Currículo Lattes <http://lattes.cnpq.br/4486320338778858>

5 5

Gerar vetor de dados com assimetria negativa no R
by Andre Oliveira 01 Out '14

01 Out '14

Pessoal, como posso gerar um vetor de tamanho n de dados assimétricos com assimetria negativa? Pode ser qualquer distribuição! André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 3

Função
by Andre Oliveira 01 Out '14

01 Out '14

Pessoal tenho uma matriz 10 por 10 com nome de dados. Gostaria que fazer o seguinte para i variando de 1 a 10 e guardar a média e desviopadrão. sd(dados2[,i]) mean(dados2[,i]) André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 4