- R-br - DINF UFPR

função matlines()
by Fernando Antonio de souza 12 Jul '13

12 Jul '13

Boa tarde Estou com uma dúvida! por que ao utilizar a função matlines() a linha adcionda ao gráfico é quebrada e não contínua (parece um conjunto de retas ligando os pontos). Como eu posso solucionar isso? grato

3 7

erro de objecto
by alanarocha＠sapo.pt 12 Jul '13

12 Jul '13

boa tarde porque me dá este erro? all.equal(integrate(f=fdp_expweibull,par=c(1,1,1),0,Inf)$value, 1) Error in match.fun(f) : object 'fdp_expweibull' not found Ana Quoting r-br-request(a)listas.c3sl.ufpr.br: > 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: Simulação estocástica (Walmes Zeviani) > 2. Re: Coordenadas e distâncias (Tito Conte) > 3. Re: Soma de Quadrado Tipo III no R para analisar um fatorial > desbalaceado (Ivan Bezerra Allaman) > 4. Re: Soma de Quadrado Tipo III no R para analisar um fatorial > desbalaceado (Cesar Rabak) > 5. Gráfico de intervalo de confiança de modelo não linear > (Fernando Antonio de souza) > 6. Ordenação de boxplots pela mediana no ggplot2 (Nicolay Cunha) > 7. Re: Gráfico de intervalo de confiança de modelo não linear > (walmes .) > 8. Re: Ordenação de boxplots pela mediana no ggplot2 (walmes .) > 9. Re: Ordenação de boxplots pela mediana no ggplot2 (Nicolay Cunha) > 10. Re: Gráfico de intervalo de confiança de modelo não linear > (Fernando Antonio de souza) > 11. Dados semanais (ARMIN KOENIG) > 12. Re: Dados semanais (Edson Lira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 11 Jul 2013 12:14:10 -0300 > From: Walmes Zeviani <walmes(a)ufpr.br> > To: Hildemar Calixto <hildemarcalixto(a)yahoo.com>, > r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Simulação estocástica > Message-ID: > <CAFU=Eka_ka3W7uF_i3rzRMmdkgNDK+SDDMv_mxGLEqTpGdTSng(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Hildemar, > > Envie sua mensagem para o endereço da lista (r-br(a)listas.c3sl.ufpr.br) e > não para os administradores dela. Envite enviar arquivos anexos. A lista > recomenda que você hospede em algum lugar (na public do DropBox por > exemplo) e envie o link para acesso. > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > > 2013/7/11 Hildemar Calixto <hildemarcalixto(a)yahoo.com> > >> Prezados, gostaria de sua ajuda para resolver um problema de simulação, no >> qual utilizo o método da transformação inversa para gerar valores >> de uma v.a. com distribuição hipergeométrica. Conforme material anexo >> (pág. 80) fiz uma função para tal fim, porém está fornece alguns valores >> "NA" e não sei como solucionar. >> >> Agradeço desde já a ajuda. >> >> -------------------------------- >> Hildemar Calixto >> Graduando em Estatística >> Universidade Federal do Ceará >> -------------------------------- >> > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/0e9d581e/atta…> > > ------------------------------ > > Message: 2 > Date: Thu, 11 Jul 2013 12:17:13 -0300 > From: Tito Conte <tito.conte(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Coordenadas e distâncias > Message-ID: > <CACqq46xQ0Hm8PnPE0jHiChvqu8vhRtiRpocdX9G3HqKEAAH=9g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > utilize um laço, > o comando for ou while resolvem, dependendo do que você deseje fazer > > Tito Conte > > > > 2013/7/8 Roberto Guevara <robertoguevara.bio(a)gmail.com> > >> Olá pessoal do R-br. >> >> Meus dados consistem basicamente em coordenadas e o horário em que >> cada ponto foi obtido, como exemplificado a seguir: >> >> obs dia min x y >> 1 1 400 425 475 >> 2 1 410 425 425 >> 3 1 420 475 325 >> 4 1 430 525 325 >> 5 1 435 525 275 >> 6 1 443 525 275 >> >> obs se refere a cada unidade de amostra, dia a cada dia de coleta (um >> total de 73), min é o horário em unidades de minutos e x e y cada >> elemento do ponto cartesiano. >> >> O que preciso obter desses dados é a distribuição de distâncias >> euclidianas entre um ponto escolhido aleatoriamente e todos os pontos >> seguintes de um mesmo dia. Ou seja, preciso escolher aleatoriamente um >> ponto por dia, calcular, para todos os dias, a euclidiana entre este >> ponto e todos os pontos posteriores do mesmo dia e relacionar a >> distância espacial com a diferença de tempo entre os pontos e o dia em >> que cada distância está associada. >> >> Eu já elaborei um código pra isso. Caso tenham sugestões para melhorar >> o código, ficarei agradecido. Porém, não consigo passar para uma >> próxima etapa. Como poderia fazer para repetir esse procedimento o >> número de vezes que fosse necessário? Por exemplo, existem milhares de >> combinações de pontos escolhidos aleatoriamente ao longo dos 73 dias. >> Eu gostaria de repetir o procedimento 1000 vezes para ter uma ideia da >> variação de distribuição de distâncias que posso obter em decorrência >> da escolha aleatória de um ponto de origem por dia. >> >> Código já elaborado: >> >> arquivo<-read.table("arquivo.txt",head=T) >> >> fixo<-numeric() #vetor com origens de cada dia >> ult<-numeric() #vetor com últimos pontos de cada dia >> xy<-list() #lista para armazenar a "parte dos pontos de cada dia" >> delimitadas por fixo e ult >> hor<-list() #lista para armazenar a "parte dos horarios de cada dia" >> delimitadas por fixo e ult >> distancias<-list()#distancias euclidianas >> tempo<-list()#diferença de tempo entre os pontos >> dia<-list()#dia de coleta >> for (i in 1:max(unique(arquivo$dia))){ #repetindo até o total de dias >> fixo[i]<-sample(arquivo$obs[arquivo$dia==i],1)#sorteio ponto de >> origem >> ult[i]<-max(arquivo$obs[arquivo$dia==i])#armazenando último ponto >> de cada dia >> xy[[i]]<-arquivo[fixo[i]:ult[i],4:5]#armazenando parte dos pontos >> de cada dia >> hor[[i]]<-arquivo[fixo[i]:ult[i],6]#armazenando parte dos horários >> de cada dia >> distancias[[i]]<-dist(xy[[i]])[1:nrow(xy[[i]])-1]#somente as >> distâncias em relação ao primeiro ponto interessam >> tempo[[i]]<-dist(horarios[[i]])[1:length(horarios[[i]])-1]#somente >> diferenças ao primeiro ponto >> dia[[i]]<-rep(i,length(arquivo$dia[fixo[i]:ult[i]])-1)$cada dia de >> coleta >> } >> >> distancias<-unlist(distancias) >> tempo<-unlist(tempo) >> dia<-unlist(dia) >> >> dtd<-data.frame(distancias,tempo,dia) >> >> Este data.frame seria como uma unidade de amostra das milhares de >> distribuições que são possíveis de obter. Como otimizar para que >> qualquer número desejado não seja dispendioso de obter? >> >> 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 >> _______________________________________________ >> 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/20130711/69054e0a/atta…> > > ------------------------------ > > Message: 3 > Date: Thu, 11 Jul 2013 08:23:33 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Soma de Quadrado Tipo III no R para analisar um > fatorial desbalaceado > Message-ID: > <1373556213.31227.YahooMailNeo(a)web162301.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Esta dúvida já foi sanada a "n" tópicos atrás. Para obter a mesma > soma de quadrados do tipo III do SAS, você precisa usar a mesma > reparametrização do "DEMO", ou seja, SOMA ZERO. > > (s,f,p) > Allaman > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/31c1bdcd/atta…> > > ------------------------------ > > Message: 4 > Date: Thu, 11 Jul 2013 15:11:17 -0200 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Soma de Quadrado Tipo III no R para analisar um > fatorial desbalaceado > Message-ID: > <CAKrF98kGO-9eeDX00x9Q59MrHfTqj+LXqGQwqCE-CVcq-+7Tww(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Ivan, > > Estou entendendo que você acha que já vale a pena termos um FAQ desta lista? > > []s > -- > Cesar Rabak > > > > 2013/7/11 Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > >> Esta dúvida já foi sanada a "n" tópicos atrás. Para obter a mesma soma de >> quadrados do tipo III do SAS, você precisa usar a mesma reparametrização do >> "DEMO", ou seja, SOMA ZERO. >> >> (s,f,p) >> Allaman >> * >> * >> \begin{signature} >> <<>>= >> Prof. Dr. Ivan Bezerra Allaman >> Universidade Estadual de Santa Cruz >> Departamento de Ciências Exatas e Tecnológicas >> Ilhéus/BA - Brasil >> Fone: +55 73 3680-5596 >> E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com >> @ >> \end{signature} >> >> _______________________________________________ >> 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/20130711/03299d22/atta…> > > ------------------------------ > > Message: 5 > Date: Thu, 11 Jul 2013 17:17:52 -0300 > From: Fernando Antonio de souza <nandodesouza(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Gráfico de intervalo de confiança de modelo não linear > Message-ID: > <CAFrWFs=1ozgQT9Y7vfJaczLPwkX8LNsj+taEiFFuj+1unxR0Kg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Caros amigos > > Estou tentando criar um gráfico de um modelo não linear utilizando a função > gnls (pacote nlme). Utilizei essa função por que meus dados não possuem > variâncias homogêneas e a função gnls me permite lidar melhor com isso. > > Pois bem. Dei uma lida no post do Walmes em seu blogo (Ridículas) sobre > como fazer ( > http://ridiculas.wordpress.com/2011/05/19/bandas-de-confianca-para-modelo-d…) > no entanto não consigo obter a matriz gradiente e Hessiana neste > modelo > e por isso não consigo dar continuidade. > > Abaixo segue o CMR, bem como os comando fornecido pelo Walmes até onde > consegui evoluir. Este é meu gargalo! > > Desde já agradeço a atenção > > GM<-structure(list(Animal = c(1L, 60L, 62L, 68L, 72L, 7L, 73L, 98L, > 99L, 100L, 91L, 92L, 93L, 51L, 83L, 89L, 86L, 87L, 94L, 95L, > 17L, 82L, 84L, 97L, 3L, 76L, 78L, 48L, 52L, 55L, 69L, 13L, 22L, > 33L, 63L, 12L, 67L, 77L, 80L, 41L, 42L, 46L, 50L, 25L, 53L, 57L, > 64L, 28L, 35L, 58L, 61L, 21L, 24L, 54L, 59L, 20L, 30L, 38L, 45L, > 10L, 11L, 19L, 9L, 14L, 18L, 32L, 56L, 2L, 4L, 5L, 23L, 27L, > 29L, 40L, 47L, 6L, 8L, 31L, 34L, 36L, 39L, 44L, 65L), Manejo = c(1L, > 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, > 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, > 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, > 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, > 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > 2L, 2L), Gest = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 90L, > 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, > 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, > 110L, 110L, 110L, 110L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, > 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 140L, 140L, > 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, > 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, > 140L, 140L, 140L, 140L), Fetos = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), NaGLAMA = c(0.515523054, > 0.274319898, 0.260406616, 0.345393579, 0.941416811, 3.242298603, > 0.543581754, NA, NA, 0.310913983, 0.103119122, 0.242179757, 1.511982572, > 0.248561812, 1.377624119, 1.239927094, NA, 0.347712662, NA, 0.907352204, > NA, 1.629736193, 0.707560235, 1.260423046, 0.6423459, 1.314054397, > NA, 0.728530524, 0.398189115, 0.749234643, 0.544891837, NA, 5.434432736, > 1.273166074, NA, 1.055026755, 0.937103858, 1.297767845, 0.691943258, > 1.336532109, NA, 1.69233929, 2.762457777, 1.639361362, 1.220678864, > 0.696944641, 0.555373873, 2.694735252, 1.593673698, 4.8664716, > 1.14279477, 0.952757874, 1.538877833, 3.417410969, NA, 3.827188165, > 4.90094219, NA, 1.254578538, 3.997852624, 1.562503764, 2.689082873, > 8.459237956, 8.787103701, 2.535628503, 4.715973548, 2.187817632, > 1.860875669, 2.656679208, 4.929889812, 2.548334487, 4.177756252, > 4.293566946, 5.717494612, 1.55775636, 4.922141642, 4.650077264, > 6.112568346, 5.235183032, 1.907048841, 2.814221997, 5.477504341, > 5.465705692)), .Names = c("Animal", "Manejo", "Gest", "Fetos", > "NaGLAMA"), row.names = c(NA, -83L), class = "data.frame") > > library(nlme) > NAGLM1<-gnls(NaGLAMA~a*exp(c*Gest),start=list(a=0.04,c=0.03),weight=varIdent(form=~1|Gest),na.action=na.omit,data=GM[-c(59,57,61),],subset=c(Gest>0 > & Fetos==1)) > > > #Códigos fornecidos pelo Walmes Adaptado > > criando função que representa o modelo, retorna gradidente e hessiano > > expo.der <- deriv3(~a*exp(c*Gest), > c("a", "c"), > function(Gest, a,c) NULL) > str(expo.der) > > #----------------------------------------------------------------------------- > # diagnose gráfica e primeiro chute > > plot(NaGLAMA~Gest, data=GM, xlab="idade Gestacional (dias)", > ylab="Conteudo de sodio(g)") > a <- 0.04; c <-0.03; > curve(expo.der(x, a, c), add=TRUE, col=2) > start <- list(a=a, c=c) > > #----------------------------------------------------------------------------- > # ajustando o modelo aos dados a partir dos valores iniciais via gráfico > > n0 <- gnls(NaGLAMA~expo.der(Gest, > a,c),weight=varIdent(form=~1|Gest),na.action=na.omit, > data=GM[-c(59,57,61),], start=start) > summary(n0) > confint(n0) > #---------------------------------------------------------------------------------------------------- > #Não consigo acessar estas matrizes. Não consegui progredi. > str(n0) > str(n0$m$fitted()) > c(n0$m$fitted()) > attr(n0$m$fitted(), "gradient") > attr(n0$m$fitted(), "hessian") > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/daa99be6/atta…> > > ------------------------------ > > Message: 6 > Date: Thu, 11 Jul 2013 18:20:10 -0400 > From: Nicolay Cunha <nicolaycunha(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Ordenação de boxplots pela mediana no ggplot2 > Message-ID: > <CAAPK02Akk8Yzv6bprcRXRf7ZyMig0AjeXC2TgZ9odhDr=LjFTg(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > # Prezados, > # tenho uma dúvida em relação a ordenação dos boxplots em relação à > mediana em um gráfico de três fatores. > # quero que os boxplots fiquem ordenados da maior para a menor mediana > (escada), dentro de cada nível do fator que confere as cores no > gráfico > # abaixo o CMR > > > ################################################################################################################################ > #install.packages("ggplot2") > library(ggplot2) > cmr <- read.table ('http://pastebin.com/raw.php?i=61J9W4Gt', header = T) > str(cmr) > > #ordenando fatores > cmr$height <- factor(cmr$height, levels = c('h','m','l')) > cmr$Pulse <- factor (cmr$Pulse, levels = c("1LF_Low_Short", > "2PF_Low_Short", "2PF_Medium_Medium", > "3GF_Medium_Short","3GF_Medium_Long")) > > #plot > ggplot(cmr, aes( height , variable, > colour = factor(pop, levels = pop[order(factor(Pulse))]), > fill = factor(Pulse)))+ > xlab ('Anther levels')+ > ylab ('Length (mm)')+ > scale_colour_manual( values = rep("black", 24), legend = FALSE) + > scale_fill_manual(name = 'Flood Pulse intensity', values = > c("#E5E5E5", "#ACACAC", "#747474", "#3B3B3B", "#030303"))+ > geom_boxplot(mapping = NULL, data = NULL, stat = "boxplot", position = > "dodge", outlier.colour = "black", > outlier.shape = 8, outlier.size = 1.5, notch = FALSE, notchwidth = > 0.5, show_guide = TRUE)+ > theme_bw() > ################################################################################################################################ > > > # a dúvida é, como ordenar os boxplots dentro de cada nível de Flood > em ordem descrescente da mediana (varias escadinhas)? > # [] > > > -- > Nicolay Leme da Cunha > > Biólogo, Mestre, Doutorando em Ecologia e Conservação > Universidade Federal de Mato Grosso do Sul, 79070-900 > Campo Grande, MS, Brasil > E-mail: nicolaycunha(a)gmail.com > lattes.cnpq.br/5916316648872099 > > > ------------------------------ > > Message: 7 > Date: Thu, 11 Jul 2013 19:32:00 -0300 > From: "walmes ." <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Gráfico de intervalo de confiança de modelo não > linear > Message-ID: > <CAFU=EkZZ7TPMxfvYfgh8d2GohxdFeVBSZibMRqDs1zgnU1-SvA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Para você extrair n0 tem que ser de classe 'nls' e você tá usando a 'gnls'. > Então as matrizes por fora assim > > x <- seq(0,140,20) > args(expo.der) > m <- expo.der(Gest=x, a=coef(n0)["a"], c=coef(n0)["c"]) > > c(m) > attr(m, "gradient") > attr(m, "hessian") > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/acbdb490/atta…> > > ------------------------------ > > Message: 8 > Date: Thu, 11 Jul 2013 20:31:54 -0300 > From: "walmes ." <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Ordenação de boxplots pela mediana no ggplot2 > Message-ID: > <CAFU=Eka7H0-6v2=zivwm0B_3G-r-R9msY+P++WsPA+JQczhbzQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Eu sei fazer com a lattice. Está nesse script > > http://www.leg.ufpr.br/~walmes/cursoR/cnj/script2.R > > Vá executando desde o começo por causa das manipuções com a base de dados > que é feita em sequência e tals. A parte que interessa é em torno da linha > 405. > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/fcfcbeba/atta…> > > ------------------------------ > > Message: 9 > Date: Thu, 11 Jul 2013 20:20:04 -0400 > From: Nicolay Cunha <nicolaycunha(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Ordenação de boxplots pela mediana no ggplot2 > Message-ID: > <CAAPK02Cd18FgjkUwaFPCz7dOomB13=GtTVBGfQKDNK_kYE_aTw(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Obrigado Walmes, se não houver como fazer isso no ggplot2 partirei pro > lattice como você sugere. > Nicolay > > > > 2013/7/11 walmes . <walmeszeviani(a)gmail.com>: >> Eu sei fazer com a lattice. Está nesse script >> >> http://www.leg.ufpr.br/~walmes/cursoR/cnj/script2.R >> >> Vá executando desde o começo por causa das manipuções com a base de dados >> que é feita em sequência e tals. A parte que interessa é em torno da linha >> 405. >> >> À disposição. >> Walmes. >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: walmes(a)ufpr.br >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código >> mínimo reproduzível. > > > > -- > Nicolay Leme da Cunha > > Biólogo, Mestre, Doutorando em Ecologia e Conservação > Universidade Federal de Mato Grosso do Sul, 79070-900 > Campo Grande, MS, Brasil > E-mail: nicolaycunha(a)gmail.com > lattes.cnpq.br/5916316648872099 > > > ------------------------------ > > Message: 10 > Date: Thu, 11 Jul 2013 21:56:57 -0300 > From: Fernando Antonio de souza <nandodesouza(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Gráfico de intervalo de confiança de modelo não > linear > Message-ID: > <CAFrWFsmgR4r0YkFBbHw-7jvL3NoavXXnohOBJT60peUR5m4u-A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Mais uma vez, muito obrigado Walmes! Acabei de rodar aqui e pelo jeito é > muito mais simples que pelo nls. > > Muito obrigado > > > > > Em 11 de julho de 2013 19:32, walmes . <walmeszeviani(a)gmail.com> escreveu: > >> Para você extrair n0 tem que ser de classe 'nls' e você tá usando a >> 'gnls'. Então as matrizes por fora assim >> >> x <- seq(0,140,20) >> args(expo.der) >> m <- expo.der(Gest=x, a=coef(n0)["a"], c=coef(n0)["c"]) >> >> c(m) >> attr(m, "gradient") >> attr(m, "hessian") >> >> À disposição. >> Walmes. >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: walmes(a)ufpr.br >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> 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/20130711/c9902af6/atta…> > > ------------------------------ > > Message: 11 > Date: Thu, 11 Jul 2013 18:16:53 -0700 (PDT) > From: ARMIN KOENIG <arminkoenig19(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Dados semanais > Message-ID: > <1373591813.54894.YahooMailNeo(a)web162004.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Prezados, > Estou lidando com uma série temporal, que tem 55 semanas. > Para realizar os cálculos de previsão, utilizei: > frequency=7, start=c(....), end=c(...) e funcionou. > Mas, quando pedi plot (gerar o gráfico), o R devolveu-me um eixo X > muito estranho. > Minha dúvida é: como faço para que o "começo" (semana) e o "final" > (também semana), seja plotado adequadamente? > Obrigado pela ajuda, > Armin > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/61ca3c78/atta…> > > ------------------------------ > > Message: 12 > Date: Fri, 12 Jul 2013 05:51:41 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, ARMIN > KOENIG <arminkoenig19(a)yahoo.com.br> > Subject: Re: [R-br] Dados semanais > Message-ID: > <1373633501.42969.YahooMailNeo(a)web161206.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Poste um CMR. Umas linhas de seu banco de dados, por exemplo no > DatafileHost, DropBox e poste o link para nossos amigos possam > ajudá-lo, inclusive o seu script também pode ser postado. > > > [ ]'s. > > > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: ARMIN KOENIG <arminkoenig19(a)yahoo.com.br> > Para: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Quinta-feira, 11 de Julho de 2013 21:16 > Assunto: [R-br] Dados semanais > > > > Prezados, > Estou lidando com uma série temporal, que tem 55 semanas. > Para realizar os cálculos de previsão, utilizei: > frequency=7, start=c(....), end=c(...) e funcionou. > Mas, quando pedi plot (gerar o gráfico), o R devolveu-me um eixo X > muito estranho. > Minha dúvida é: como faço para que o "começo" (semana) e o "final" > (também semana), seja plotado adequadamente? > Obrigado pela ajuda, > Armin > _______________________________________________ > 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/20130712/41966ff1/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 31, assunto 15 > *****************************************

2 1

Dados semanais
by ARMIN KOENIG 12 Jul '13

12 Jul '13

Prezados, Estou lidando com uma série temporal, que tem 55 semanas. Para realizar os cálculos de previsão, utilizei: frequency=7, start=c(....), end=c(...) e funcionou. Mas, quando pedi plot (gerar o gráfico), o R devolveu-me um eixo X muito estranho. Minha dúvida é: como faço para que o "começo" (semana) e o "final" (também semana), seja plotado adequadamente? Obrigado pela ajuda, Armin

2 1

Gráfico de intervalo de confiança de modelo não linear
by Fernando Antonio de souza 11 Jul '13

11 Jul '13

Caros amigos Estou tentando criar um gráfico de um modelo não linear utilizando a função gnls (pacote nlme). Utilizei essa função por que meus dados não possuem variâncias homogêneas e a função gnls me permite lidar melhor com isso. Pois bem. Dei uma lida no post do Walmes em seu blogo (Ridículas) sobre como fazer ( http://ridiculas.wordpress.com/2011/05/19/bandas-de-confianca-para-modelo-d…) no entanto não consigo obter a matriz gradiente e Hessiana neste modelo e por isso não consigo dar continuidade. Abaixo segue o CMR, bem como os comando fornecido pelo Walmes até onde consegui evoluir. Este é meu gargalo! Desde já agradeço a atenção GM<-structure(list(Animal = c(1L, 60L, 62L, 68L, 72L, 7L, 73L, 98L, 99L, 100L, 91L, 92L, 93L, 51L, 83L, 89L, 86L, 87L, 94L, 95L, 17L, 82L, 84L, 97L, 3L, 76L, 78L, 48L, 52L, 55L, 69L, 13L, 22L, 33L, 63L, 12L, 67L, 77L, 80L, 41L, 42L, 46L, 50L, 25L, 53L, 57L, 64L, 28L, 35L, 58L, 61L, 21L, 24L, 54L, 59L, 20L, 30L, 38L, 45L, 10L, 11L, 19L, 9L, 14L, 18L, 32L, 56L, 2L, 4L, 5L, 23L, 27L, 29L, 40L, 47L, 6L, 8L, 31L, 34L, 36L, 39L, 44L, 65L), Manejo = c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), Gest = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L), Fetos = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), NaGLAMA = c(0.515523054, 0.274319898, 0.260406616, 0.345393579, 0.941416811, 3.242298603, 0.543581754, NA, NA, 0.310913983, 0.103119122, 0.242179757, 1.511982572, 0.248561812, 1.377624119, 1.239927094, NA, 0.347712662, NA, 0.907352204, NA, 1.629736193, 0.707560235, 1.260423046, 0.6423459, 1.314054397, NA, 0.728530524, 0.398189115, 0.749234643, 0.544891837, NA, 5.434432736, 1.273166074, NA, 1.055026755, 0.937103858, 1.297767845, 0.691943258, 1.336532109, NA, 1.69233929, 2.762457777, 1.639361362, 1.220678864, 0.696944641, 0.555373873, 2.694735252, 1.593673698, 4.8664716, 1.14279477, 0.952757874, 1.538877833, 3.417410969, NA, 3.827188165, 4.90094219, NA, 1.254578538, 3.997852624, 1.562503764, 2.689082873, 8.459237956, 8.787103701, 2.535628503, 4.715973548, 2.187817632, 1.860875669, 2.656679208, 4.929889812, 2.548334487, 4.177756252, 4.293566946, 5.717494612, 1.55775636, 4.922141642, 4.650077264, 6.112568346, 5.235183032, 1.907048841, 2.814221997, 5.477504341, 5.465705692)), .Names = c("Animal", "Manejo", "Gest", "Fetos", "NaGLAMA"), row.names = c(NA, -83L), class = "data.frame") library(nlme) NAGLM1<-gnls(NaGLAMA~a*exp(c*Gest),start=list(a=0.04,c=0.03),weight=varIdent(form=~1|Gest),na.action=na.omit,data=GM[-c(59,57,61),],subset=c(Gest>0 & Fetos==1)) #Códigos fornecidos pelo Walmes Adaptado criando função que representa o modelo, retorna gradidente e hessiano expo.der <- deriv3(~a*exp(c*Gest), c("a", "c"), function(Gest, a,c) NULL) str(expo.der) #----------------------------------------------------------------------------- # diagnose gráfica e primeiro chute plot(NaGLAMA~Gest, data=GM, xlab="idade Gestacional (dias)", ylab="Conteudo de sodio(g)") a <- 0.04; c <-0.03; curve(expo.der(x, a, c), add=TRUE, col=2) start <- list(a=a, c=c) #----------------------------------------------------------------------------- # ajustando o modelo aos dados a partir dos valores iniciais via gráfico n0 <- gnls(NaGLAMA~expo.der(Gest, a,c),weight=varIdent(form=~1|Gest),na.action=na.omit, data=GM[-c(59,57,61),], start=start) summary(n0) confint(n0) #---------------------------------------------------------------------------------------------------- #Não consigo acessar estas matrizes. Não consegui progredi. str(n0) str(n0$m$fitted()) c(n0$m$fitted()) attr(n0$m$fitted(), "gradient") attr(n0$m$fitted(), "hessian")

2 2

Ordenação de boxplots pela mediana no ggplot2
by Nicolay Cunha 11 Jul '13

11 Jul '13

# Prezados, # tenho uma dúvida em relação a ordenação dos boxplots em relação à mediana em um gráfico de três fatores. # quero que os boxplots fiquem ordenados da maior para a menor mediana (escada), dentro de cada nível do fator que confere as cores no gráfico # abaixo o CMR ################################################################################################################################ #install.packages("ggplot2") library(ggplot2) cmr <- read.table ('http://pastebin.com/raw.php?i=61J9W4Gt', header = T) str(cmr) #ordenando fatores cmr$height <- factor(cmr$height, levels = c('h','m','l')) cmr$Pulse <- factor (cmr$Pulse, levels = c("1LF_Low_Short", "2PF_Low_Short", "2PF_Medium_Medium", "3GF_Medium_Short","3GF_Medium_Long")) #plot ggplot(cmr, aes( height , variable, colour = factor(pop, levels = pop[order(factor(Pulse))]), fill = factor(Pulse)))+ xlab ('Anther levels')+ ylab ('Length (mm)')+ scale_colour_manual( values = rep("black", 24), legend = FALSE) + scale_fill_manual(name = 'Flood Pulse intensity', values = c("#E5E5E5", "#ACACAC", "#747474", "#3B3B3B", "#030303"))+ geom_boxplot(mapping = NULL, data = NULL, stat = "boxplot", position = "dodge", outlier.colour = "black", outlier.shape = 8, outlier.size = 1.5, notch = FALSE, notchwidth = 0.5, show_guide = TRUE)+ theme_bw() ################################################################################################################################ # a dúvida é, como ordenar os boxplots dentro de cada nível de Flood em ordem descrescente da mediana (varias escadinhas)? # [] -- Nicolay Leme da Cunha Biólogo, Mestre, Doutorando em Ecologia e Conservação Universidade Federal de Mato Grosso do Sul, 79070-900 Campo Grande, MS, Brasil E-mail: nicolaycunha(a)gmail.com lattes.cnpq.br/5916316648872099

2 2

Re: [R-br] Soma de Quadrado Tipo III no R para analisar um fatorial desbalaceado
by Ivan Bezerra Allaman 11 Jul '13

11 Jul '13

Esta dúvida já foi sanada a "n" tópicos atrás. Para obter a mesma soma de quadrados do tipo III do SAS, você precisa usar a mesma reparametrização do "DEMO", ou seja, SOMA ZERO. (s,f,p) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

Coordenadas e distâncias
by Roberto Guevara 11 Jul '13

11 Jul '13

Olá pessoal do R-br. Meus dados consistem basicamente em coordenadas e o horário em que cada ponto foi obtido, como exemplificado a seguir: obs dia min x y 1 1 400 425 475 2 1 410 425 425 3 1 420 475 325 4 1 430 525 325 5 1 435 525 275 6 1 443 525 275 obs se refere a cada unidade de amostra, dia a cada dia de coleta (um total de 73), min é o horário em unidades de minutos e x e y cada elemento do ponto cartesiano. O que preciso obter desses dados é a distribuição de distâncias euclidianas entre um ponto escolhido aleatoriamente e todos os pontos seguintes de um mesmo dia. Ou seja, preciso escolher aleatoriamente um ponto por dia, calcular, para todos os dias, a euclidiana entre este ponto e todos os pontos posteriores do mesmo dia e relacionar a distância espacial com a diferença de tempo entre os pontos e o dia em que cada distância está associada. Eu já elaborei um código pra isso. Caso tenham sugestões para melhorar o código, ficarei agradecido. Porém, não consigo passar para uma próxima etapa. Como poderia fazer para repetir esse procedimento o número de vezes que fosse necessário? Por exemplo, existem milhares de combinações de pontos escolhidos aleatoriamente ao longo dos 73 dias. Eu gostaria de repetir o procedimento 1000 vezes para ter uma ideia da variação de distribuição de distâncias que posso obter em decorrência da escolha aleatória de um ponto de origem por dia. Código já elaborado: arquivo<-read.table("arquivo.txt",head=T) fixo<-numeric() #vetor com origens de cada dia ult<-numeric() #vetor com últimos pontos de cada dia xy<-list() #lista para armazenar a "parte dos pontos de cada dia" delimitadas por fixo e ult hor<-list() #lista para armazenar a "parte dos horarios de cada dia" delimitadas por fixo e ult distancias<-list()#distancias euclidianas tempo<-list()#diferença de tempo entre os pontos dia<-list()#dia de coleta for (i in 1:max(unique(arquivo$dia))){ #repetindo até o total de dias fixo[i]<-sample(arquivo$obs[arquivo$dia==i],1)#sorteio ponto de origem ult[i]<-max(arquivo$obs[arquivo$dia==i])#armazenando último ponto de cada dia xy[[i]]<-arquivo[fixo[i]:ult[i],4:5]#armazenando parte dos pontos de cada dia hor[[i]]<-arquivo[fixo[i]:ult[i],6]#armazenando parte dos horários de cada dia distancias[[i]]<-dist(xy[[i]])[1:nrow(xy[[i]])-1]#somente as distâncias em relação ao primeiro ponto interessam tempo[[i]]<-dist(horarios[[i]])[1:length(horarios[[i]])-1]#somente diferenças ao primeiro ponto dia[[i]]<-rep(i,length(arquivo$dia[fixo[i]:ult[i]])-1)$cada dia de coleta } distancias<-unlist(distancias) tempo<-unlist(tempo) dia<-unlist(dia) dtd<-data.frame(distancias,tempo,dia) Este data.frame seria como uma unidade de amostra das milhares de distribuições que são possíveis de obter. Como otimizar para que qualquer número desejado não seja dispendioso de obter? 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

2 1

Re: [R-br] Simulação estocástica
by Walmes Zeviani 11 Jul '13

11 Jul '13

Hildemar, Envie sua mensagem para o endereço da lista (r-br(a)listas.c3sl.ufpr.br) e não para os administradores dela. Envite enviar arquivos anexos. A lista recomenda que você hospede em algum lugar (na public do DropBox por exemplo) e envie o link para acesso. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2013/7/11 Hildemar Calixto <hildemarcalixto(a)yahoo.com> > Prezados, gostaria de sua ajuda para resolver um problema de simulação, no > qual utilizo o método da transformação inversa para gerar valores > de uma v.a. com distribuição hipergeométrica. Conforme material anexo > (pág. 80) fiz uma função para tal fim, porém está fornece alguns valores > "NA" e não sei como solucionar. > > Agradeço desde já a ajuda. > > -------------------------------- > Hildemar Calixto > Graduando em Estatística > Universidade Federal do Ceará > -------------------------------- >

1 0

Soma de Quadrado Tipo III no R para analisar um fatorial desbalaceado
by Tales Fernandes 11 Jul '13

11 Jul '13

Pessoal, mais uma vez recorro a experiência de vocês. Estou analisando um conjunto de dados de um experimento em esquema fatorial com dois fatores e em DBC. Os dados estão desbalanceados, porém sem caselas vazias. Quero testar a seguinte hipótese: efeito de determinado parâmetro corrigido para todos os outros. Que é a famosa soma de quadrados tipo III (estou certo ?). Mas como testar esta hipótese no R? segue meu CMR: #################### library(car) dados <- read.table("material.txt",header = TRUE) names(dados) > names(dados) [1] "Lamina" "Dose" "Bloco" "V1" M1<- lm(V1~ Bloco + Lamina*Dose, data=dados) # soma de quadrados tipo I anova(M1) # soma de quadrados tipo II Anova(M1, type="II") # soma de quadrados tipo III Anova(M1, type="III") ################### O meu problema é que: a SQ tipo I bate com a tipo I do SAS; a SQ tipo II bate com a tipo II do SAS, mas a SQ tipo III não bate com a tipo III do SAS. E eu sei que a soma de quadrado correta para testar a minha hipótese é a SQ tipo III, mas não quero usar o SAS. Aquele tipo III da função Anova da biblioteca car corresponde ao mesmo tipo III do SAS ? Posso usar seus resultados como se fossem a SQ tipo III? Então porque os resultados são diferentes? Outra questão: Neste mesmo experimento foram avaliadas outras variáveis, e em uma delas não foi possível medir todas as combinações, gerando 2 caselas vazias. Daí quando faço a Anova tipo III dá o seguinte erro: > Anova(M2, type="III") Erro em Anova.III.lm(mod, error, singular.ok = singular.ok, ...) : there are aliased coefficients in the model Alguém sabe o que significa este erro ? Muito obrigado! Tales Jesus Fernandes Doutorando em Estatística UFLA Universidade Federal de Lavras

2 3

Micro-dados do Censo 2012
by Rogério Santos 11 Jul '13

11 Jul '13

Bom dia, Alguém sabe informar se os micro-dados do IBGE já estão disponíveis para venda? Obrigado -- ------------------------------------------------------------ *Rogério Santos* *Estatístico* ------------------------------------------------------------

2 1

Autocorrelação de uma variável em Séries Temporais
by Jodavid Ferreira 10 Jul '13

10 Jul '13

Como verificar autocorrelação de uma unica variáveis? ou seja, de todos as observações entre si. Conheço o método gráfico "acf" mas gostaria outro método, com valores. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

1 0

Pesos Censo 2010
by Raphael Saldanha 08 Jul '13

08 Jul '13

Prezados, Tenho uma dúvida mais conceitual sobre os microdados do censo 2010. Como devo trabalhar com os pesos dos registros? Mais exatamente no meu caso, preciso fazer algumas tabulações de variáveis categóricas, por exemplo, sexo e cor, com os microdados. Como devo proceder? ---Atenciosamente,Raphael Saldanha rfsaldanha(a)outlook.com

4 4

Sugestao de pacote
by Marcelo Laia 08 Jul '13

08 Jul '13

Preciso realizar uma análise sobre a resistência de dado genótipo a uma praga. O experimento foi o seguinte: genotipo T0 T1 T2 1 10 8 3 2 10 10 10 3 10 2 0 4 10 1 0 5 10 10 10 (...) n 10 8 2 Em cada genótipo foram adicionados 10 insetos no tempo zero e passou-se a contar o número de sobreviventes a cada dia. Genótipos onde os insetos morreram de início são mais resistentes que aqueles em que os insetos sobreviveram mais tempo. Há repetição desse mesmo experimento ao longo do tempo. Foi observado que em determinada época do ano a mortalidade varia. Poderia sugerir um pacote que auxilia nesse tipo de análise? Ou sugerir material bibliográfico? Funções, etc...? Obrigado -- Marcelo Brasil (Brazil, for English Speakers) Linux user number 487797

2 3

Tamanho de uma Variável!
by Sérgio Henrique almeida da silva ju 08 Jul '13

08 Jul '13

Qual comando do R que me retorna o tamanho de uma variável? Sei que o length me dá o tamanho do banco! ex.: var <- c(123, 123, 456,456") tamanho da variável = 3 Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

5 8

Como obter coeficientes padronizados em modelos mistos Poisson?
by Juliana Schietti 08 Jul '13

08 Jul '13

Olá, Gostaria de saber se é possivel obter os coeficientes beta padronizados para as variáveis preditoras (efeitos fixos) em um modelo misto para distribuição Poisson. Minha pergunta é mais sobre estatística do que sobre uso do R, mas espero que alguém possa ajudar. Obrigada, Juliana

2 1

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 08 Jul '13

08 Jul '13

Jônatan, a função rasterToshapefile é interessante, mas não atribuí o índices i e j ao shapefile criado. Encontrei outra solução no script abaixo e espero auxilar outras pessoas. Apenas uma ressalva, ao final gero um shapefile de pontos, com os atributos i, j e o valor da célula, mas depois associo ao shapefile de polígonos no Quantum GIS através da função Join de Tabelas. Outra solução possível é juntar a sua solução de rasterToPolygon e depois realizar o Join de tabela com os shapefiles de pontos. Realmente seria interessante que a função rasterToPolygon levasse as informações de linha e coluna para o shapefile criado. Mais uma vez obrigado Jôntan e Elias, vocês me orientaram bastante na solução deste problema. require(raster) require(rgdal) require(shapefiles) require(maptools) ### ENTRADA DE DADOS ### ## criando um raster r <- raster(ncol=319,nrow=319) atribuindo valores r[] <- (1:ncell(r)) ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) matriz dim(matriz) tail(matriz) summary (matriz) names(matriz) # a matriz não possui nome de campos # Data Frame possui nome dos atributos class(matriz) require(sp) pontos<-data.frame(matriz) head(pontos) class(pontos) dim(pontos) names(pontos) coordinates(pontos) <- c("row","col") # c("Norte_Sul","Leste_Oeste") class(pontos) dim(pontos) getSlots("SpatialPointsDataFrame") slot(pontos, "bbox") require(maptools) writePointsShape(pontos,"matriz319pt", factor2char = TRUE, max_nchar=254) readShapePoints("matriz319pt", proj4string = CRS(as.character(NA)),verbose = FALSE,repair=FALSE)

3 2

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 07 Jul '13

07 Jul '13

Mauro, falta o carregar o pacote GDAL. Por favor, insira o comando no script require(rgdal). Att, Alex Santos Por que não rodaram os últimos comandos? > require(raster) > Loading required package: raster > Loading required package: sp > > r <- raster(nc=20,nr=10) > r[] <- rep(1:5,each=40) > plot(r) > res <- data.frame(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) > names(res) <- c("row","col","cell","v") > ## criando stack com 4 layers > s <- stack(r,r,r,r) > ## atribuindo valores para cada layer baseado no res > s[] <- as.matrix(res) > plot(s) > #p<-rasterToPoints(s, spatial=T) > pp<- rasterToPolygons(s,dissolve=F) > #spplot(pp) > # writeOGR(p,getwd(),"p2",driver="ESRI Shapefile") > writeOGR(pp,getwd(),"matriz_pl",driver="ESRI Shapefile") > Error: could not find function "writeOGR" > > pescrito <- readOGR(getwd(), "matriz_pl") > Error: could not find function "readOGR" > pescrito@data > Error: object 'pescrito' not found >

1 0

Table com zeros
by Fátima Lima Paula 07 Jul '13

07 Jul '13

Olá pessoal. Tenho uma variável cnes.1 (número do hospital) num banco. Quando faço table(banco$cnes.1), retornam vários com quantidade zero. Acho que é porque existe uma variável cnes.2 que apresenta quantidades para esses hospitais que no cnes.1 estão com zeros. Existe uma forma de se evitar isso? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 2

Dúvidas sobre a amostra
by edmar Caldas 07 Jul '13

07 Jul '13

Como eu sei que a minha amostra representa a população. Ex: tenho 10 mil observações minha população. Retirei aleatoriamente 1000 observações que é minha amostra Faço uma inferência. Um teste T. A partir do resultado do p-valor <0,05 posso dizer que minha amostra pode representar minha população com um IC 95%?

3 2

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 07 Jul '13

07 Jul '13

Jônatan, excelente! no seu script apenas substitui rasterToPoints para rasterToPolygons no final. O objetivo do raster é apenas para extrair as informações do índice (i, j). Os valores do pixel funcionam como ID's. Agora, o resultado final já é georreferenciado e as quadrículas estão prontas para uma posterior análise geoestatística. Att, Alex Santos ### ENTRADA DE DADOS ### r <- raster(nc=20,nr=10) r[] <- rep(1:5,each=40) plot(r) res <- data.frame(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) names(res) <- c("row","col","cell","v") ## criando stack com 4 layers s <- stack(r,r,r,r) ## atribuindo valores para cada layer baseado no res s[] <- as.matrix(res) plot(s) #p<-rasterToPoints(s, spatial=T) pp<- rasterToPolygons(s,dissolve=F) #spplot(pp) # writeOGR(p,getwd(),"p2",driver="ESRI Shapefile") writeOGR(pp,getwd(),"matriz_pl",driver="ESRI Shapefile") pescrito <- readOGR(getwd(), "matriz_pl") pescrito@data

1 0

Re: [R-br] Fwd: Re: Fwd: Re: [Dúvida] Criando pacotes para o R.
by Pedro Rafael 06 Jul '13

06 Jul '13

Obrigado pela opinião Leonard Mendonça. [ ], Pedro Rafael Diniz Marinho. Em 17 de novembro de 2012 14:45, Pedro Rafael < pedro.rafael.marinho(a)gmail.com> escreveu: > Obrigado pela opinião Leonard Mendonça. > > [ ], > Pedro Rafael Diniz Marinho. > > > Em 17 de novembro de 2012 12:13, Leonard Mendonça de Assis [via R-br] < > ml-node+s2285057n4657102h39(a)n4.nabble.com> escreveu: > > No meu entendimento é uma nova versão >> >> []s >> Leonard de Assishttp://about.me/ldeassis >> >> Em 17/11/2012 11:32, Pedro Rafael escreveu: >> >> Benilton voce poderia me tirar uma dúvida? Já li a documentacao do R mas >> nao ficou claro para mim. Caso eu queira fazer uma nova versao do pacote >> pelo que eu entendi eu devo fazer upload do novo arquivo tar.gz informando >> a nova versao. Contudo, caso eu precise fazer pequenos ajustes, por >> exemplo, ajustes na documentacao. Como fazer isso sem mudar a versao do >> pacote? Isso pode ser feito? >> ---------- Mensagem encaminhada ---------- >> >> De: "beniltoncarvalho [via R-br]" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657102&i=0> >> > >> Data: 16/11/2012 17:50 >> Assunto: Re: [R-br] Fwd: Re: [Dúvida] Criando pacotes para o R. >> Para: "Pedro Rafael Diniz Marinho" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657102&i=1> >> > >> >> >> Aguarde 24-48h, é o tempo do ciclo dos "builds". >> On 16 Nov 2012 20:47, "Pedro Rafael" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657087&i=0>> >> wrote: >> >>> Pessoal os binarios para windows e Mac sao gerados pelo CRAN? Vi em >>> alguns pacotes que entraram nos espelhos hoje nao tem esses binarios. >>> Existe uma mensagem do lado dizendo que os pacotes sao novos. Os binarios >>> nao sao gerados logo quando o pacofe entra no CRAN? >>> ---------- Mensagem encaminhada ---------- >>> De: "beniltoncarvalho [via R-br]" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657087&i=1>> >>> >>> >>> Data: 16/11/2012 17:01 >>> Assunto: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> Para: "Pedro Rafael Diniz Marinho" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657087&i=2>> >>> >>> >>> >>> Se vc rodar >>> >>> R CMD check pacote >>> >>> ou seu equivalente para o windows... o que sao os erros e o conteudo >>> dos logs? >>> >>> 2012/11/16 tiago souza marçal <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=0> >>> > >>> >>> O meu path esta configurado da seguinte forma: >>> >>> C:\Perl\site\bin;C:\Perl\bin;C:\Rtools\bin;C:\MiKTeX2.9\miktex\bin;C:\R\R-2.15.2\bin\i386;C:\MinGW\bin;C:\HTMLHelpWorkshop;C:\Program >>> Files\Common Files\Microsoft >>> Shared\WindowsLive;%LINDOAPI_HOME%bin\win32;C:\Windows\system32;C:\Windows;C:\Windows\System32\Wbem;C:\Windows\System32\WindowsPowerShell\v1.0\;C:\Program >>> Files\Windows Live\Shared;C:\Program Files\Internet Explorer >>> >>> Você vê algum problema? >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Fri, 16 Nov 2012 13:58:58 -0200 >>> >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=1> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=2> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> na minha vm windows, eu instalei primeiro o miktex, depois instalei o >>> rtools. Deve ser alguma loucura do windows (Path faltando, por exemplo) >>> >>> []s >>> Leonard de Assishttp://about.me/ldeassis >>> >>> Em 16/11/2012 13:11, tiago souza marçal escreveu: >>> >>> O Rtools eu tenho o problema e quando eu peço para gerar o PDF ( R CMD >>> Rd2pdf ) ele me retorna uma mensagem de erro falando do MikTex 1.9, >>> entretanto, eu estou utilizando o MikTex 2.9. >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Thu, 15 Nov 2012 20:33:21 -0200 >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=3> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=4> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> RTools, errei o nome >>> >>> voce acha no CRAN >>> >>> []s >>> Leonard de Assishttp://about.me/ldeassis >>> >>> Em 15/11/2012 09:42, tiago souza marçal escreveu: >>> >>> Leonard onde eu consigo baixar esse rutil? >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Thu, 15 Nov 2012 09:28:46 -0200 >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=5> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=6> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> No Windows tem que instalar o rutil recomendado para a versão do R que >>> estamos usando. >>> >>> Tenho uma VM windows 64 bits comigo e todos os problemas que encontrei >>> quando fuçei com pacotes foram referentes ao R não localizar todos os >>> 'bichinhos' do Rutil os quais precisava. >>> >>> []s >>> Leonard de Assishttp://about.me/ldeassis >>> >>> Em 15/11/2012 08:06, tiago souza marçal escreveu: >>> >>> Estou com um problema parecido com o seu, porém, no Windows . No meu >>> caso o check do PDF está ok mas quando eu peço para ele ser gerado (R >>> CMD Rd2pdf) uma mensagem de erro e retornada informando uma versão >>> antiga do MikTex que eu não estou usando. >>> >>> Mas no seu caso você observou mensagem de erro no processo de >>> instalação do pacote (R CMD INSTALL) ? Caso isso tenha ocorrido este >>> pode ter sido o problema, em função de algum erro na edição do arquivo.Rd. >>> Observe que foi retornado pra você um WARNING no check dos arquivos.Rd >>> com algumas mensagens de erro. >>> >>> Espero ter te ajudado. >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Wed, 14 Nov 2012 23:27:37 -0300 >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=7> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=8> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> Pessoal estou estudando como criar pacotes no R e tive o seguinte erro >>> na hora de dar o check. Consigo gerar o arquivo tar.gz mas o pdf da >>> documentação não é criado. Alguém saberia me informar onde estaria o erro? >>> >>> OBS: Estou em uma máquina Linux. >>> >>> * using log directory >>> ‘/home/pedro/Documentos/Gama_Log_Logistica/FUNCAO_CODIGO/PRDM.Rcheck’ >>> * using R version 2.15.1 (2012-06-22) >>> * using platform: x86_64-pc-linux-gnu (64-bit) >>> * using session charset: UTF-8 >>> * checking for file ‘PRDM/DESCRIPTION’ ... OK >>> * checking extension type ... Package >>> * this is package ‘PRDM’ version ‘1.0’ >>> * checking package namespace information ... OK >>> * checking package dependencies ... OK >>> * checking if this is a source package ... OK >>> * checking if there is a namespace ... OK >>> * checking for executable files ... OK >>> * checking whether package ‘PRDM’ can be installed ... OK >>> * checking installed package size ... OK >>> * checking package directory ... OK >>> * checking for portable file names ... OK >>> * checking for sufficient/correct file permissions ... OK >>> * checking DESCRIPTION meta-information ... OK >>> * checking top-level files ... OK >>> * checking index information ... OK >>> * checking package subdirectories ... OK >>> * checking R files for non-ASCII characters ... OK >>> * checking R files for syntax errors ... OK >>> * checking whether the package can be loaded ... OK >>> * checking whether the package can be loaded with stated dependencies >>> ... OK >>> * checking whether the package can be unloaded cleanly ... OK >>> * checking whether the namespace can be loaded with stated dependencies >>> ... OK >>> * checking whether the namespace can be unloaded cleanly ... OK >>> * checking for unstated dependencies in R code ... OK >>> * checking S3 generic/method consistency ... OK >>> * checking replacement functions ... OK >>> * checking foreign function calls ... OK >>> * checking R code for possible problems ... OK >>> * checking Rd files ... WARNING >>> prepare_Rd: PRDM-package.Rd:30: All text must be in a section >>> prepare_Rd: PRDM-package.Rd:31: All text must be in a section >>> prepare_Rd: prdm.Rd:36-42: Dropping empty section \value >>> * checking Rd metadata ... OK >>> * checking Rd cross-references ... OK >>> * checking for missing documentation entries ... OK >>> * checking for code/documentation mismatches ... OK >>> * checking Rd \usage sections ... OK >>> * checking Rd contents ... OK >>> * checking for unstated dependencies in examples ... OK >>> * checking examples ... OK >>> * checking PDF version of manual ... ERROR >>> Rd conversion errors: >>> Warning: unknown option '-' >>> Warning in file(con, "r") : >>> cannot open file '/home/pedro/Documentos/MESTRADO': No such file or >>> directory >>> Error in file(con, "r") : cannot open the connection >>> >>> >>> [ ], >>> Pedro Rafael Diniz Marinho. >>> >>> >>> Em 14 de novembro de 2012 20:43, Pedro Rafael <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=9> >>> > escreveu: >>> >>> Grande Paulo Justiniano, como vai? Muito obrigado pelas suas >>> informações. Uso o Linux faz uns 6 anos mas nunca contribui ou publiquei >>> meus códigos. Estava vendo que a linguagem Julia vem com a licença MIT que >>> eu não acho muito interessante porque quem quiser pode vender o código. >>> Minha dúvida sobre a licença GPL ou GPL2 é se uma pessoa poderá desmanchar >>> algum código que fiz ou mudar completamente o foco do pacote, retirar nome >>> dos autores, etc. Existe um pacote hoje muito usado no R que foi feito por >>> uma pessoa e depois de um tempo parou de ser mantido. Chegou um outro cara >>> de fora e reimplementou o pacote e tirou o nome dos autores antigos e >>> depois muita luta o nome dessas pessoas permaneceu como "autores antigos". >>> Isso ocorre com certa frequência? Existe alguma periodicidade de manutenção >>> do pacote que deve ser respeitada? >>> >>> [ ], >>> Pedro Rafael Diniz Marinho. >>> >>> >>> Em 14 de novembro de 2012 20:28, Paulo Justiniano [via R-br] <[hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=10>>escreveu: >>> >>> >>> >>> On Wed, 14 Nov 2012, Pedro Rafael wrote: >>> >>> > Pessoal, digamos que faço um pacote para o R. Como funciona o professo >>> de submissão? >>> >>> 1. Leia atentamente o manual "Writing R extensions" que descreve isto >>> em >>> detalhes. >>> >>> 2 certifique-se que voce testou o seu pacote : >>> >>> R CMD check SEU_PACOTE >>> >>> e corrija todos os eventuais problemas >>> >>> > Digamos que eu submeta um pacote "meu_pacote" e depois queira >>> > fazer modificações no mesmo. Para acrescentar essas modificações eu >>> tenho que submeter novamente o pacote com o mesmo nome para que as >>> atualizações >>> > fiquem disponíveis? >>> >>> Sim. Para isto voce alter a numero da versao no arquivo DESCRIPTION >>> >>> > >>> > Outra dúvida que tenho é sobre a segurança que tenho que meu pacote >>> não seja modificado. Qualquer pessoa pode pegar esse pacote e alterar ele >>> > completamente sem nenhuma permissão minha ou qualquer modificação tem >>> que ser acompanhada por intermédio da pessoa que submeteu o pacote pela >>> primeira >>> > vez? >>> > >>> >>> Bem vindo ao Universo do Software livre >>> >>> Note que se voce quiser voce pode baixar os fontes do R e mudar o que >>> quiser >>> >>> Falando mais seriamente no DESCRIPTION voce define a licenca. >>> Leia sobre licxancas GPL e análogas >>> >>> Minha sugestao: internalize completamente o conceito de software livre e >>> livre-se destas preocupacoes >>> >>> > >>> > São dúvidas bestas mas que nunca entendi direito. >>> >>> Nao sao somente suas. >>> Não é por acaso que há uma diversidade de tipos de licensas, não sem >>> controvérsias! >>> >>> > >>> > Obrigado à todos. >>> > [ ], >>> > Pedro Rafael Diniz Marinho. >>> > >>> > >>> > >>> _______________________________________________ >>> R-br mailing list >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657009&i=0> >>> 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. >>> >>> ------------------------------ >>> If you reply to this email, your message will be added to the >>> discussion below: >>> >>> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Criando-pacotes-para-o-R-tp46… >>> To unsubscribe from R-br, click here. >>> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >>> >>> >>> >>> >>> _______________________________________________ R-br mailing list [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=11> >>> 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[hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=12>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 [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=13> >>> 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[hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=14>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 [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=15> >>> 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[hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=16>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 [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=17> >>> 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 >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=18> >>> 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 >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=19> >> 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. >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Criando-pacotes-para-o-R-tp46… >> To unsubscribe from R-br, click here. >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657087&i=3> >> 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 >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657087&i=4> >> 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. >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Criando-pacotes-para-o-R-tp46… >> To unsubscribe from R-br, click here. >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> >> >> _______________________________________________ >> R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4657102&i=2>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 >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657102&i=3> >> 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. >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Criando-pacotes-para-o-R-tp46… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >> . >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> > >

5 18

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 06 Jul '13

06 Jul '13

Jônatan, respondendo as suas perguntas: ## estás concluindo isso da saída da linha de comando que há um limite de linhas para print. Sim, realmente estou concluindo em função do limite de linhas para print. ## jogue em uma variavel o resultado do cbind e verifique res <- cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) dim(res) tail(res) Abaixo está os resultados do código que você solicitou. [24999,] 84 99 24999 9895 [ atingido getOption("max.print") -- últimas 65001 linhas omitidas ]] > matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) > dim(matriz) [1] 90000 4 > tail(matriz) row col [89995,] 300 295 89995 0 [89996,] 300 296 89996 0 [89997,] 300 297 89997 0 [89998,] 300 298 89998 0 [89999,] 300 299 89999 0 [90000,] 300 300 90000 0 Em relação a dica do shapefile, tentarei utilizar a função writeOGR() sugerida pelo autor, apenas não sei como criar os atributos. Conforme você sugeriu o link: http://listas.inf.ufpr.br/pipermail/r-br/2013-April/010226.html A questão agora é transferir este resultado "matriz" para um shapefile. Caso eu tenha êxito ao gerar o shapefile, a coluna referente aos valores do pixel do meus raster contem valores que eu posso utilizar como join. Mais uma vez obrigado! Alex Santos

2 1

[Dúvida] Operador Lógico (==).
by Pedro Rafael 06 Jul '13

06 Jul '13

Pessoal eu tenho uma dúvida. Considere o código abaixo: function(par,x){ beta = par[1] c = par[2] a = par[3] a * beta * c * exp(-(beta*x)^c) * (beta*x)^(c-1) * (1 - exp(-(beta*x)^c))^(a-1) } integrate(f=fdp_expweibull,par=c(1,1,1),0,Inf)$value == 1 Porque o resultado é false? [ ], Pedro Rafael Diniz Marinho.

2 3

Re: [R-br] Rcmdr
by Darren Norris 05 Jul '13

05 Jul '13

Marcelo, Acho que os dados nao estao no formato correto – por exemplo para fazer um “ANOVA” tem que ter pelo menos um factor no dataframe ou Rcmdr nao deixar rodar. Veja este pequeno exemplo (https://skydrive.live.com/redir?resid=F0A9FA3480208398!396) mostrando de como carregar dados da pacote "car" e como os menus ficam ativas. Se isto da certo entao o problema eh com os dados e nao com Rcmdr. Google “Introdução ao R Commander" e tambem veja este guias: http://paginapessoal.utfpr.edu.br/jonas/disciplinas/estatistica-experimenta… e Fox, J. (2012) Getting started with R Commander http://socserv.mcmaster.ca/jfox/Misc/Rcmdr/Getting-Started-with-the-Rcmdr.p… E proximo vez por favor resultados de “ sessionInfo() ” Abs, Darren From: Marcelo Lacerda Rezende Sent: Friday, July 05, 2013 4:50 PM To: Darren Norris Subject: Re: [R-br] Rcmdr Fica cinza. Entro no data set, seleciono um conjunto de dados, mas continua cinza... O que pode ser? Obrigado pela resposta! Marcelo Em 5 de julho de 2013 10:09, Darren Norris <doon75(a)hotmail.com> escreveu: atraves os menus???: ANOVA: estatisticas ->Medias estatisticas ->Ajuste de modelos modelos->testes de hipotese Kruskal-Wallis : estatisticas ->testes nao-parametricos Cinza?? Entao tem que carregar um conjunto de dados....... Abs, Darren Message: 1 Date: Thu, 4 Jul 2013 11:55:12 -0500 From: Marcelo Lacerda Rezende <marcelo.l.rezende(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Rcmdr Message-ID: <CAEYb-yvWcoHceTNHjaf6VEA2ij-=xTBg6-vnFxvenCbV7gPDPQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Obrigado Alexandre! Farei isso. Att, Marcelo Em 2 de julho de 2013 20:33, Alexandre Serpa <serpa.alexandre(a)gmail.com>escreveu: Olá Marcelo, Você pode tentar instalar algum dos plugins que tem no site do rcmdr http://www.rcommander.com/ Mas o teste está lá, mesmo que não tenha um menu para ele. basta você chamar via console. abs Alexandre Serpa Psicólogo Especialista em Métodos Computacionais Estatísticos Mestre em Educação Doutorando em Psicologia email: serpa.alexandre(a)gmail.com Em 2 de julho de 2013 19:52, Marcelo Lacerda Rezende < marcelo.l.rezende(a)gmail.com> escreveu: Prezados, instalei o Rcmdr mas alguns testes não ficaram disponíveis como o Kruskal-Wallis, Anova ... Por favor, será que devo instalar algum outro pacote ou é outro tipo de problema? Att, Marcelo _______________________________________________ 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. -- Prof Dr Marcelo Lacerda Rezende Visiting Researcher Scholar Department of Agricultural Economics -Texas A&M University - USA Professor Federal University of Alfenas - Brazil Institute of Exact Sciences -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130704/75ca532a/atta…> ------------------------------ Message: 2 Date: Thu, 4 Jul 2013 10:11:04 -0700 (PDT) From: Edson Lira <edinhoestat(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Ordenar meses Message-ID: <1372957864.57177.YahooMailNeo(a)web161206.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Esta aí uma rotina que uso em um banco de dados para doadores de sangue tip_doa$mês<-factor(tip_doa$mes, levels=c("jan","fev","mar","abr","mai","jun" ,"jul","ago","set","out","nov","dez")) tip_doa$mês esta variável recebe ordenada a variável tip_doa$mes [ ]'s. Edson Lira Estatístico Manaus-Amazonas ________________________________ De: Augusto Ribas <ribas.aca(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quinta-feira, 4 de Julho de 2013 10:30 Assunto: Re: [R-br] Ordenar meses Ola, se são apenas meses, talvez fosse mais fácil usar eles como inteiros, número inteiros de 1 a 12, ai é so colocar na ordem numérica crescente que vai ficar organizado. Se os nomes estiverem em português, eu não tenho certeza como fazer, mas se você passar para o inglês, uma possibilidade talvez seria usar o as.Date do R. Mas vc precisaria de datas completas, como por exemplo: datas <- as.Date(c("01/03/2013","01/02/2013","01/05/2013","01/04/2013","01/06/2013","01/01/2013"),"%d/%m/%Y") class(datas) Notem que os dados são de uma classe Data agora. Outra coisa, é que o formato ai eu descrevi no segundo argumento la, esse aqui "%d/%m/%Y", olhe o help da função, da para os nomes dos meses, nomes abreviados etc. Aqui tem um guia legal http://www.statmethods.net/input/dates.html Umas vez que vc converteu para data, da para usar comandos como rank ou sort nos dados. sort(datas) [1] "2013-01-01" "2013-02-01" "2013-03-01" "2013-04-01" "2013-05-01" "2013-06-01" order(datas) [1] 6 2 1 4 3 5 E a partir daqui é so usar essa informação para ordenar sua matriz. # matriz<-matrix(rnorm(6*4),ncol=4,nrow=6) matriz matriz[order(datas),] # Bem é isso, existem outros pacotes para lidar com datas também, como por exemplo o lubridate, talvez vale a pena dar uma olhada, se for precisar muito desse tipo de coisa Espero ter ajudado. Em 2 de julho de 2013 14:28, Aline Lipsky <line_lips(a)hotmail.com> escreveu: Boa tarde, Eu tenho uma matriz 12 x 43, e os nomes das linhas são o meses em ordem alfabética. Como faço para por em ordem cronológica? Grata, Aline F. Lipsky _______________________________________________ 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. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056 _______________________________________________ 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/20130704/6bc109f6/atta…> ------------------------------ Message: 3 Date: Thu, 4 Jul 2013 16:16:56 -0300 From: Arthur Telles Calegario <tcalegario(a)hotmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Exportar Mapa de Krigagem Message-ID: <BAY160-W58B5B8CC9C244450273FDCC27C0(a)phx.gbl> Content-Type: text/plain; charset="utf-8" Bom dia Pessoal,Realizei a krigagem de uns dados, seguindo http://www.leg.ufpr.br/geoR/geoRdoc/geoRintro.html, e preciso exportar pro ArGis. Alguem sabe como?Tentei exportar como tiff e referenciar no GIS. Mas o Gis não leu a figura, não sei pq.Agradeço,Arthur. Date: Thu, 4 Jul 2013 10:11:04 -0700 From: edinhoestat(a)yahoo.com.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Ordenar meses Esta aí uma rotina que uso em um banco de dados para doadores de sangue tip_doa$mês<-factor(tip_doa$mes, levels=c("jan","fev","mar","abr","mai","jun" ,"jul","ago","set","out","nov","dez")) tip_doa$mês esta variável recebe ordenada a variável tip_doa$mes [ ]'s. Edson Lira Estatístico Manaus-Amazonas De: Augusto Ribas <ribas.aca(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quinta-feira, 4 de Julho de 2013 10:30 Assunto: Re: [R-br] Ordenar meses Ola, se são apenas meses, talvez fosse mais fácil usar eles como inteiros, número inteiros de 1 a 12, ai é so colocar na ordem numérica crescente que vai ficar organizado. Se os nomes estiverem em português, eu não tenho certeza como fazer, mas se você passar para o inglês, uma possibilidade talvez seria usar o as.Date do R. Mas vc precisaria de datas completas, como por exemplo: datas <- as.Date(c("01/03/2013","01/02/2013","01/05/2013","01/04/2013","01/06/2013","01/01/2013"),"%d/%m/%Y")class(datas) Notem que os dados são de uma classe Data agora.Outra coisa, é que o formato ai eu descrevi no segundo argumento la, esse aqui "%d/%m/%Y", olhe o help da função, da para os nomes dos meses, nomes abreviados etc. Aqui tem um guia legal http://www.statmethods.net/input/dates.html Umas vez que vc converteu para data, da para usar comandos como rank ou sort nos dados. sort(datas)[1] "2013-01-01" "2013-02-01" "2013-03-01" "2013-04-01" "2013-05-01" "2013-06-01" order(datas)[1] 6 2 1 4 3 5 E a partir daqui é so usar essa informação para ordenar sua matriz. #matriz<-matrix(rnorm(6*4),ncol=4,nrow=6) matrizmatriz[order(datas),]# Bem é isso, existem outros pacotes para lidar com datas também, como por exemplo o lubridate, talvez vale a pena dar uma olhada, se for precisar muito desse tipo de coisa Espero ter ajudado. Em 2 de julho de 2013 14:28, Aline Lipsky <line_lips(a)hotmail.com> escreveu: Boa tarde, Eu tenho uma matriz 12 x 43, e os nomes das linhas são o meses em ordem alfabética. Como faço para por em ordem cronológica? Grata, Aline F. Lipsky _______________________________________________ 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. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056 _______________________________________________ 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/20130704/220bf983/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 31, assunto 6 **************************************** -- Prof Dr Marcelo Lacerda Rezende Visiting Researcher Scholar Department of Agricultural Economics -Texas A&M University - USA Professor Federal University of Alfenas - Brazil Institute of Exact Sciences

1 0

Re: [R-br] Rcmdr
by Darren Norris 05 Jul '13

05 Jul '13

atraves os menus???: ANOVA: estatisticas ->Medias estatisticas ->Ajuste de modelos modelos->testes de hipotese Kruskal-Wallis : estatisticas ->testes nao-parametricos Cinza?? Entao tem que carregar um conjunto de dados....... Abs, Darren Message: 1 Date: Thu, 4 Jul 2013 11:55:12 -0500 From: Marcelo Lacerda Rezende <marcelo.l.rezende(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Rcmdr Message-ID: <CAEYb-yvWcoHceTNHjaf6VEA2ij-=xTBg6-vnFxvenCbV7gPDPQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Obrigado Alexandre! Farei isso. Att, Marcelo Em 2 de julho de 2013 20:33, Alexandre Serpa <serpa.alexandre(a)gmail.com>escreveu: > Olá Marcelo, > > Você pode tentar instalar algum dos plugins que tem no site do rcmdr > http://www.rcommander.com/ > > Mas o teste está lá, mesmo que não tenha um menu para ele. basta você > chamar via console. > > abs > > Alexandre Serpa > > Psicólogo > Especialista em Métodos Computacionais Estatísticos > Mestre em Educação > Doutorando em Psicologia > email: serpa.alexandre(a)gmail.com > > > Em 2 de julho de 2013 19:52, Marcelo Lacerda Rezende < > marcelo.l.rezende(a)gmail.com> escreveu: > >> >> Prezados, >> instalei o Rcmdr mas alguns testes não ficaram disponíveis como o >> Kruskal-Wallis, Anova ... >> Por favor, será que devo instalar algum outro pacote ou é outro tipo de >> problema? >> Att, >> Marcelo >> >> _______________________________________________ >> 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. > -- Prof Dr Marcelo Lacerda Rezende Visiting Researcher Scholar Department of Agricultural Economics -Texas A&M University - USA Professor Federal University of Alfenas - Brazil Institute of Exact Sciences -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130704/75ca532a/atta…> ------------------------------ Message: 2 Date: Thu, 4 Jul 2013 10:11:04 -0700 (PDT) From: Edson Lira <edinhoestat(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Ordenar meses Message-ID: <1372957864.57177.YahooMailNeo(a)web161206.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Esta aí uma rotina que uso em um banco de dados para doadores de sangue tip_doa$mês<-factor(tip_doa$mes, levels=c("jan","fev","mar","abr","mai","jun" ,"jul","ago","set","out","nov","dez")) tip_doa$mês esta variável recebe ordenada a variável tip_doa$mes [ ]'s. Edson Lira Estatístico Manaus-Amazonas ________________________________ De: Augusto Ribas <ribas.aca(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quinta-feira, 4 de Julho de 2013 10:30 Assunto: Re: [R-br] Ordenar meses Ola, se são apenas meses, talvez fosse mais fácil usar eles como inteiros, número inteiros de 1 a 12, ai é so colocar na ordem numérica crescente que vai ficar organizado. Se os nomes estiverem em português, eu não tenho certeza como fazer, mas se você passar para o inglês, uma possibilidade talvez seria usar o as.Date do R. Mas vc precisaria de datas completas, como por exemplo: datas <- as.Date(c("01/03/2013","01/02/2013","01/05/2013","01/04/2013","01/06/2013","01/01/2013"),"%d/%m/%Y") class(datas) Notem que os dados são de uma classe Data agora. Outra coisa, é que o formato ai eu descrevi no segundo argumento la, esse aqui "%d/%m/%Y", olhe o help da função, da para os nomes dos meses, nomes abreviados etc. Aqui tem um guia legal http://www.statmethods.net/input/dates.html Umas vez que vc converteu para data, da para usar comandos como rank ou sort nos dados. > sort(datas) [1] "2013-01-01" "2013-02-01" "2013-03-01" "2013-04-01" "2013-05-01" "2013-06-01" > order(datas) [1] 6 2 1 4 3 5 E a partir daqui é so usar essa informação para ordenar sua matriz. # matriz<-matrix(rnorm(6*4),ncol=4,nrow=6) matriz matriz[order(datas),] # Bem é isso, existem outros pacotes para lidar com datas também, como por exemplo o lubridate, talvez vale a pena dar uma olhada, se for precisar muito desse tipo de coisa Espero ter ajudado. Em 2 de julho de 2013 14:28, Aline Lipsky <line_lips(a)hotmail.com> escreveu: > >Boa tarde, > > >Eu tenho uma matriz 12 x 43, e os nomes das linhas são o meses em ordem >alfabética. >Como faço para por em ordem cronológica? > > >Grata, > > >Aline F. Lipsky >_______________________________________________ >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. > -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056 _______________________________________________ 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/20130704/6bc109f6/atta…> ------------------------------ Message: 3 Date: Thu, 4 Jul 2013 16:16:56 -0300 From: Arthur Telles Calegario <tcalegario(a)hotmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Exportar Mapa de Krigagem Message-ID: <BAY160-W58B5B8CC9C244450273FDCC27C0(a)phx.gbl> Content-Type: text/plain; charset="utf-8" Bom dia Pessoal,Realizei a krigagem de uns dados, seguindo http://www.leg.ufpr.br/geoR/geoRdoc/geoRintro.html, e preciso exportar pro ArGis. Alguem sabe como?Tentei exportar como tiff e referenciar no GIS. Mas o Gis não leu a figura, não sei pq.Agradeço,Arthur. Date: Thu, 4 Jul 2013 10:11:04 -0700 From: edinhoestat(a)yahoo.com.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Ordenar meses Esta aí uma rotina que uso em um banco de dados para doadores de sangue tip_doa$mês<-factor(tip_doa$mes, levels=c("jan","fev","mar","abr","mai","jun" ,"jul","ago","set","out","nov","dez")) tip_doa$mês esta variável recebe ordenada a variável tip_doa$mes [ ]'s. Edson Lira Estatístico Manaus-Amazonas De: Augusto Ribas <ribas.aca(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quinta-feira, 4 de Julho de 2013 10:30 Assunto: Re: [R-br] Ordenar meses Ola, se são apenas meses, talvez fosse mais fácil usar eles como inteiros, número inteiros de 1 a 12, ai é so colocar na ordem numérica crescente que vai ficar organizado. Se os nomes estiverem em português, eu não tenho certeza como fazer, mas se você passar para o inglês, uma possibilidade talvez seria usar o as.Date do R. Mas vc precisaria de datas completas, como por exemplo: datas <- as.Date(c("01/03/2013","01/02/2013","01/05/2013","01/04/2013","01/06/2013","01/01/2013"),"%d/%m/%Y")class(datas) Notem que os dados são de uma classe Data agora.Outra coisa, é que o formato ai eu descrevi no segundo argumento la, esse aqui "%d/%m/%Y", olhe o help da função, da para os nomes dos meses, nomes abreviados etc. Aqui tem um guia legal http://www.statmethods.net/input/dates.html Umas vez que vc converteu para data, da para usar comandos como rank ou sort nos dados. > sort(datas)[1] "2013-01-01" "2013-02-01" "2013-03-01" "2013-04-01" > "2013-05-01" "2013-06-01" > order(datas)[1] 6 2 1 4 3 5 E a partir daqui é so usar essa informação para ordenar sua matriz. #matriz<-matrix(rnorm(6*4),ncol=4,nrow=6) matrizmatriz[order(datas),]# Bem é isso, existem outros pacotes para lidar com datas também, como por exemplo o lubridate, talvez vale a pena dar uma olhada, se for precisar muito desse tipo de coisa Espero ter ajudado. Em 2 de julho de 2013 14:28, Aline Lipsky <line_lips(a)hotmail.com> escreveu: Boa tarde, Eu tenho uma matriz 12 x 43, e os nomes das linhas são o meses em ordem alfabética. Como faço para por em ordem cronológica? Grata, Aline F. Lipsky _______________________________________________ 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. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056 _______________________________________________ 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/20130704/220bf983/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 31, assunto 6 ****************************************

1 0

Ordenar meses
by Aline Lipsky 04 Jul '13

04 Jul '13

Boa tarde, Eu tenho uma matriz 12 x 43, e os nomes das linhas são o meses em ordem alfabética. Como faço para por em ordem cronológica? Grata, Aline F. Lipsky

4 3

Rcmdr
by Marcelo Lacerda Rezende 04 Jul '13

04 Jul '13

Prezados, instalei o Rcmdr mas alguns testes não ficaram disponíveis como o Kruskal-Wallis, Anova ... Por favor, será que devo instalar algum outro pacote ou é outro tipo de problema? Att, Marcelo

2 2

Exemplo de fatorial duplo com um fator qualitativo e um fator quantitativo
by Alisson Lucrecio 04 Jul '13

04 Jul '13

Caro Colegas r-br, Boa tarde. Eu estou precisando de um exemplo de análise de um experimento em fatorial duplo com um fator qualitativo e um quantitativo, alguém teria algum script para me passar. Eu estou com duvida na aov, devo considerar os dois como quantitativo ou um como quantitativo e outro como qualitativo? obrigado. Alisson Lucrécio da Costa

2 1

uma pergunta sobre médias
by alanarocha＠sapo.pt 04 Jul '13

04 Jul '13

Olá Edson e todos, aqui envio o link dos dados que preciso para responder a esta pergunta. a coluna A é a coluna com os dados originais mas foi-me pedido para organizar em 4 grupos como esta na coluna B a E a pergunta para os dados é esta 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. que codigo eu tenho de fazer no R para saber a evolução e que gráfico tenho de fazer para ver essa mesma evolução? link dos dados: https://dl.dropboxusercontent.com/u/9337305/grupos_pthi.xlsx obrigada a todos Ana

1 0

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 03 Jul '13

03 Jul '13

Primeiramente obrigado a todos! Jônatan, por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo no shapefile. matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, assim, em relação a conversão para shapefile, qual seria a melhor função para criar o shapefile (preferencialmente tipo polígono) com os atributos da matriz acima. Mais uma vez agradeço a lista. Att, Alex Santos Jônatam escreveu (alterei apenas a quantidade de linhas e colunas): require(raster) ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) >

2 1

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 03 Jul '13

03 Jul '13

Primeiramente obrigado a todos! Por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo num shapefile (tipo polígono ou ponto). matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, em relação a conversão para shapefile, qual seria a melhor função no R para criar o shapefile com os atributos da matriz acima. Abaixo está o código elaborado por Jônatam, com alterações apenas na quantidade de linhas e colunas, agora resta apenas converter o resultado "matriz" para shapefile. ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r)) ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) # resultado matriz matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) matriz Mais uma vez agradeço a lista. Att, Alex Santos

1 0

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 03 Jul '13

03 Jul '13

Primeiramente obrigado a todos! Jônatan, por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo no shapefile. matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, assim, em relação a conversão para shapefile, qual seria a melhor função para criar o shapefile (preferencialmente tipo polígono) com os atributos da matriz acima. Mais uma vez agradeço a lista. Att, Alex Santos Jônatam escreveu (alterei apenas a quantidade de linhas e colunas): require(raster) ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))

1 0

Como extrair índices linha e coluna de um raster
by Alex 03 Jul '13

03 Jul '13

Primeiramente obrigado a todos! Jônatan, por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo no shapefile. matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, assim, em relação a conversão para shapefile, qual seria a melhor função para criar o shapefile (preferencialmente tipo polígono) com os atributos da matriz acima. Mais uma vez agradeço a lista. Att, Alex Santos Jônatam escreveu (alterei apenas a quantidade de linhas e colunas): require(raster) ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))

1 0

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 03 Jul '13

03 Jul '13

Primeiramente obrigado a todos! Jônatan, por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo no shapefile. matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, assim, em relação a conversão para shapefile, qual seria a melhor função para criar o shapefile (preferencialmente tipo polígono) com os atributos da matriz acima. Jônatam escreveu (alterei apenas a quantidade de linhas e colunas): require(raster) ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) Mais uma vez agradeço a lista. Att, Alex Santos

1 0

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 03 Jul '13

03 Jul '13

Primeiramente obrigado a todos! Jônatan, por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo no shapefile. matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, assim, em relação a conversão para shapefile, qual seria a melhor função para criar o shapefile (preferencialmente tipo polígono) com os atributos da matriz acima. Mais uma vez agradeço a lista. Att, Alex Santos Jônatam escreveu (alterei apenas a quantidade de linhas e colunas): require(raster) ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))

1 0

Substituir caracter
by Sérgio Henrique almeida da silva ju 03 Jul '13

03 Jul '13

Prezados No meu banco eu tenho uma variável (CPF) que tem pontos e traços e eu gostaria de retirar esses caracteres. Como faço? Ex.: > CPF = c("092.321.678-00", "089.567.543-00","456.657.543-00") > CPF 092.321.678-00" "089.567.543-00" "456.657.543-00" Quero que fique assim: "09232167800" "08956754300" "45665754300" Vi que tem o comando GSUB, mas não entendi muito bem como usá-lo. Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

3 3

Lendo R dentro do Perl
by Deisy Gysi 02 Jul '13

02 Jul '13

Boa tarde! Estou precisando rodar um programa em R dentro de um programa em Perl. Para isso, estou usando a função *system(CMD R BACH --slave arquivo.R)* do perl, porém, eu precisava que na shell aparecesse as saídas em R, não que apenas ele rodasse sem que eu pudesse ve-las, uma vez que dentro desse programa em R tenho vários avisos importantes para o usuário (se os arquivos foram devidamente importados, e análises realizadas)... Embora eu tenha criado uma log (.Rout) que é aberta no final do R rodar, eu preciso que os avisos tenham saída em tempo real na shell... Alguém tem alguma ideia de como solucionar esse problema? Obrigada, Deisy -- Atenciosamente, Deisy Morselli Gysi Bacharel em Biotecnologia Bacharel em Estatística

2 2

trabalhar com médias
by alanarocha＠sapo.pt 02 Jul '13

02 Jul '13

vou enviar o ficheiro de dados, para esta pergunta: 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. Aqui adicionar a coluna contendo a pthi T+12 https://dl.dropboxusercontent.com/u/9337305/grupos_pthi.xlsx Ana Quoting r-br-request(a)listas.c3sl.ufpr.br: > 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: interaction plot (Natalia Martins) > 2. diferenca em horas (Edson Lira) > 3. Re: diferenca em horas (andrebvs(a)bol.com.br) > 4. Instalação rmysql ubuntu (Leandro Marino) > 5. Re: Instalação rmysql ubuntu (Daniel C Bezerra) > 6. Re: Instalação rmysql ubuntu (Leandro Marino) > 7. Re: diferenca em horas (Edson Lira) > 8. trabalhar com médias (alanarocha(a)sapo.pt) > 9. trabalhar com médias (alanarocha(a)sapo.pt) > 10. Re: trabalhar com médias (Edson Lira) > 11. Re: Dúvida sobre stepwise (Victor Eduardo) > 12. [Dúvida] Função dentro de função. (Pedro Rafael) > 13. Re: [Dúvida] Função dentro de função. (Benilton Carvalho) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 1 Jul 2013 12:57:09 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: Lista R <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] interaction plot > Message-ID: > <CAGPGL2E0xhe7SPTBB8rXkH4dBVrsjoon6B1YPDLBSkdFbe1myQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Ordenei, mas mesmo assim o R ainda me dá em ordem alfabetica, > mudei o nome mas gostaria de saber se eu consigo editar o grafico, pois nas > opções do interaction nao consegui. > > Att > > > Em 1 de julho de 2013 11:02, Pedro Emmanuel Alvarenga Americano do Brasil < > emmanuel.brasil(a)gmail.com> escreveu: > >> Sugestão 1: formate como fator e ordene os fatores como conveniente. >> >> Sugestão 2: epicalc::followup.plot >> >> Sugestão 3: >> https://stat.ethz.ch/pipermail/r-sig-epi/attachments/20130401/7fcbc976/atta… >> >> >> Dr. Pedro Emmanuel A. A. do Brasil >> http://blog.ipec.fiocruz.br/lapclin-chagas/ >> Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 >> ResearchGate.net: https://www.researchgate.net/profile/Pedro_Brasil2/ >> Instituto Nacional de Infectologia/Instituto de Pesquisa Clínica Evandro >> Chagas >> Fundação Oswaldo Cruz >> Rio de Janeiro - Brasil >> Av. Brasil 4365, >> CEP 21040-360, >> Tel 55 21 3865-9648 >> e-mail: pedro.brasil(a)ipec.fiocruz.br >> e-mail: emmanuel.brasil(a)gmail.com >> >> ---Apoio aos softwares livres >> www.zotero.org - gerenciamento de referências bibliográficas. >> www.broffice.org ou www.libreoffice.org - textos, planilhas ou >> apresentações. >> www.epidata.dk - entrada de dados. >> www.r-project.org - análise de dados. >> www.ubuntu.com - sistema operacional >> >> >> Em 1 de julho de 2013 10:48, Natalia Martins <nsmbarreto(a)gmail.com>escreveu: >> >>> Prezados, bom dia! >>> Estou fazendo uns graficos utilizando a função interaction. plot. >>> O grafico sai em ordem alfabética no eixo x, como eu poderia proceder >>> para que saísse em outra ordem? >>> >>> interaction.plot(CP,Clone,Altura, xlab = 'Corpo de prova', ylab = >>> 'Altura', col=c("blue","red"),lty = 1) >>> >>> >>> O >>> brigada, >>> >>> >>> Natália da Silva Martins >>> >>> >>> _______________________________________________ >>> 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. >> > > > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestra em Ciências: Area de concentração: Estatística e Experimentação > Agronômica - ESALQ/ USP > Doutoranda em Ciências: Area de concentração: Estatística e Experimentação > Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/bbbcb1a6/atta…> > > ------------------------------ > > Message: 2 > Date: Mon, 1 Jul 2013 09:16:45 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] diferenca em horas > Message-ID: > <1372695405.42748.YahooMailNeo(a)web161205.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Caros amigos da Lista, tenho um banco com as variáveis abaixo: > > > Dia_da_Semana Dia_do_Mês > Entrada Saída Semana > > No link abaixo está um banco teste para vocês analisarem. > > > > > http://www1.datafilehost.com/d/320c3ab5 > > > Quero calcular a diferença de tempo entre a entrada e a saída. Para > poder analisar tempo médio por semana. > > Não estou conseguindo implementar nos pacotes chron e lubridate. > > Alguém poderia me ajudar? > > > [ ]'s. > > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/40e14eae/atta…> > > ------------------------------ > > Message: 3 > Date: Mon, 1 Jul 2013 13:36:42 -0300 > From: andrebvs(a)bol.com.br > To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] diferenca em horas > Message-ID: <51d1b01aabf49_a12961f478131(a)a4-winter19.tmail> > Content-Type: text/plain; charset="iso-8859-1" > > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/25fe940c/atta…> > > ------------------------------ > > Message: 4 > Date: Mon, 1 Jul 2013 15:26:50 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAKSaaFkxo81fXiDZ1No8yq0RErKD6grXobjNmMbhO7WMaNSbsQ(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Caros, > > executei a instalação do Rmysql no ubuntu sem sucesso. já executei o > primeiro passo das instruções mas nao sei como fazer o segundo passo... > alguma sugestão? > > Configuration error: > could not find the MySQL installation include and/or library > directories. Manually specify the location of the MySQL > libraries and the header files and re-run R CMD INSTALL. > > INSTRUCTIONS: > > 1. Define and export the 2 shell variables PKG_CPPFLAGS and > PKG_LIBS to include the directory for header files (*.h) > and libraries, for example (using Bourne shell syntax): > > export PKG_CPPFLAGS="-I<MySQL-include-dir>" > export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" > > Re-run the R INSTALL command: > > R CMD INSTALL RMySQL_<version>.tar.gz > > 2. Alternatively, you may pass the configure arguments > --with-mysql-dir=<base-dir> (distribution directory) > or > --with-mysql-inc=<base-inc> (where MySQL header files reside) > --with-mysql-lib=<base-lib> (where MySQL libraries reside) > in the call to R INSTALL --configure-args='...' > > R CMD INSTALL --configure-args='--with-mysql-dir=DIR' > RMySQL_<version>.tar.gz > > ERROR: configuration failed for package ?RMySQL? > * removing ?/usr/local/lib/R/site-library/RMySQL? > > The downloaded source packages are in > ?/tmp/RtmpIdWQKu/downloaded_packages? > Warning message: > In install.packages("RMySQL") : > installation of package ?RMySQL? had non-zero exit status > > a instalação do my sql se deu com o comando: > *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 > php5-mysql phpmyadmin* > * > * > * > * > * > * > <http://www.showmetech.com.br> > > *LEANDRO MARINO* | Showmetech > Estatístico | Fotógrafo > Mobile: 55 21 9845-7707 > Email: contato(a)leandromarino.com.br > Profiles SMT: facebook <http://www.facebook.com/showmetech> | > twitter<http://www.twitter.com/showmetech> > | Google+ <https://plus.google.com/115440786884851522659> > > <http://www.leandromarino.com.br> > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/1fd91564/atta…> > > ------------------------------ > > Message: 5 > Date: Mon, 1 Jul 2013 16:11:49 -0300 > From: Daniel C Bezerra <danielcbezerra(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAJGh7uzcsYu1B9_DuhUDpE2+nODVwd1UobWsBVOZoHsjM3thGw(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Rode no terminal : > > sudo apt-get build-dep r-cran-rmysql > > E depois tente instalar o pacote novamente. > > Abs, > > Daniel > On Jul 1, 2013 3:27 PM, "Leandro Marino" <leandromarino(a)leandromarino.com.br> > wrote: > >> Caros, >> >> executei a instalação do Rmysql no ubuntu sem sucesso. já executei o >> primeiro passo das instruções mas nao sei como fazer o segundo passo... >> alguma sugestão? >> >> Configuration error: >> could not find the MySQL installation include and/or library >> directories. Manually specify the location of the MySQL >> libraries and the header files and re-run R CMD INSTALL. >> >> INSTRUCTIONS: >> >> 1. Define and export the 2 shell variables PKG_CPPFLAGS and >> PKG_LIBS to include the directory for header files (*.h) >> and libraries, for example (using Bourne shell syntax): >> >> export PKG_CPPFLAGS="-I<MySQL-include-dir>" >> export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" >> >> Re-run the R INSTALL command: >> >> R CMD INSTALL RMySQL_<version>.tar.gz >> >> 2. Alternatively, you may pass the configure arguments >> --with-mysql-dir=<base-dir> (distribution directory) >> or >> --with-mysql-inc=<base-inc> (where MySQL header files reside) >> --with-mysql-lib=<base-lib> (where MySQL libraries reside) >> in the call to R INSTALL --configure-args='...' >> >> R CMD INSTALL --configure-args='--with-mysql-dir=DIR' >> RMySQL_<version>.tar.gz >> >> ERROR: configuration failed for package ?RMySQL? >> * removing ?/usr/local/lib/R/site-library/RMySQL? >> >> The downloaded source packages are in >> ?/tmp/RtmpIdWQKu/downloaded_packages? >> Warning message: >> In install.packages("RMySQL") : >> installation of package ?RMySQL? had non-zero exit status >> >> a instalação do my sql se deu com o comando: >> *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 >> php5-mysql phpmyadmin* >> * >> * >> * >> * >> * >> * >> <http://www.showmetech.com.br> >> >> *LEANDRO MARINO* | Showmetech >> Estatístico | Fotógrafo >> Mobile: 55 21 9845-7707 >> Email: contato(a)leandromarino.com.br >> Profiles SMT: facebook <http://www.facebook.com/showmetech> | >> twitter<http://www.twitter.com/showmetech> >> | Google+ <https://plus.google.com/115440786884851522659> >> >> <http://www.leandromarino.com.br> >> >> >> _______________________________________________ >> 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/20130701/37ca7970/atta…> > > ------------------------------ > > Message: 6 > Date: Mon, 1 Jul 2013 17:22:17 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAKSaaF=8D=FQ9vebGcsDZEgn4842zTQ1Z7UfT_eyDw44DKGLtg(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > valeu daniel > funcionou! :) > > vou guardar esta dica! :) > > <http://www.showmetech.com.br> > > *LEANDRO MARINO* | Showmetech > Estatístico | Fotógrafo > Mobile: 55 21 9845-7707 > Email: contato(a)leandromarino.com.br > Profiles SMT: facebook <http://www.facebook.com/showmetech> | > twitter<http://www.twitter.com/showmetech> > | Google+ <https://plus.google.com/115440786884851522659> > > <http://www.leandromarino.com.br> > > > > 2013/7/1 Daniel C Bezerra <danielcbezerra(a)gmail.com> > >> Rode no terminal : >> >> sudo apt-get build-dep r-cran-rmysql >> >> E depois tente instalar o pacote novamente. >> >> Abs, >> >> Daniel >> On Jul 1, 2013 3:27 PM, "Leandro Marino" < >> leandromarino(a)leandromarino.com.br> wrote: >> >>> Caros, >>> >>> executei a instalação do Rmysql no ubuntu sem sucesso. já executei o >>> primeiro passo das instruções mas nao sei como fazer o segundo passo... >>> alguma sugestão? >>> >>> Configuration error: >>> could not find the MySQL installation include and/or library >>> directories. Manually specify the location of the MySQL >>> libraries and the header files and re-run R CMD INSTALL. >>> >>> INSTRUCTIONS: >>> >>> 1. Define and export the 2 shell variables PKG_CPPFLAGS and >>> PKG_LIBS to include the directory for header files (*.h) >>> and libraries, for example (using Bourne shell syntax): >>> >>> export PKG_CPPFLAGS="-I<MySQL-include-dir>" >>> export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" >>> >>> Re-run the R INSTALL command: >>> >>> R CMD INSTALL RMySQL_<version>.tar.gz >>> >>> 2. Alternatively, you may pass the configure arguments >>> --with-mysql-dir=<base-dir> (distribution directory) >>> or >>> --with-mysql-inc=<base-inc> (where MySQL header files reside) >>> --with-mysql-lib=<base-lib> (where MySQL libraries reside) >>> in the call to R INSTALL --configure-args='...' >>> >>> R CMD INSTALL --configure-args='--with-mysql-dir=DIR' >>> RMySQL_<version>.tar.gz >>> >>> ERROR: configuration failed for package ?RMySQL? >>> * removing ?/usr/local/lib/R/site-library/RMySQL? >>> >>> The downloaded source packages are in >>> ?/tmp/RtmpIdWQKu/downloaded_packages? >>> Warning message: >>> In install.packages("RMySQL") : >>> installation of package ?RMySQL? had non-zero exit status >>> >>> a instalação do my sql se deu com o comando: >>> *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 >>> php5-mysql phpmyadmin* >>> * >>> * >>> * >>> * >>> * >>> * >>> <http://www.showmetech.com.br> >>> >>> *LEANDRO MARINO* | Showmetech >>> Estatístico | Fotógrafo >>> Mobile: 55 21 9845-7707 >>> Email: contato(a)leandromarino.com.br >>> Profiles SMT: facebook <http://www.facebook.com/showmetech> | >>> twitter<http://www.twitter.com/showmetech> >>> | Google+ <https://plus.google.com/115440786884851522659> >>> >>> <http://www.leandromarino.com.br> >>> >>> >>> _______________________________________________ >>> 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/20130701/5bace641/atta…> > > ------------------------------ > > Message: 7 > Date: Mon, 1 Jul 2013 21:32:56 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "andrebvs(a)bol.com.br" <andrebvs(a)bol.com.br>, > "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] diferenca em horas > Message-ID: > <1372739576.48970.YahooMailNeo(a)web161203.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > André, valeu! > > O problema é que quando trabalhamos com arquivos texto, csv, as > datas e/ou tempos são vistos como fatores. Não consegui implementar > sua idéia, fiz o seguinte: > > Transformei o arquivo em excell (xlsx) > > Nesta forma o R entende como POSIXct e atribui uma data para as horas > Usei a rotina difftime > > m$hs<-format(as.Date(m$Saída,"H:M")) > m$he<-format(as.Date(m$Entrada,"H:M")) > > m$ht<-difftime(m$Saída,m$Entrada) > m$htrab<-as.numeric(m$ht)/60 > > Muito obrigado!!! > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: "andrebvs(a)bol.com.br" <andrebvs(a)bol.com.br> > Para: r-br(a)listas.c3sl.ufpr.br; Edson Lira <edinhoestat(a)yahoo.com.br> > Enviadas: Segunda-feira, 1 de Julho de 2013 12:36 > Assunto: Re: [R-br] diferenca em horas > > > > Não precisa de pacote, acho que é isso abaixo: > > # Seus dados salvo em "EnterExit.txt": > > Entrada Saída > 13.00 17.00 > 15.01 16.38 > . > . > . > 15.12 17.40 > > dados <- read.table("EnterExit.txt",h=T) > attach(dados) > > dif <- mean(Saída-Entrada) > > [1] 1.767143 > > Att. > André BVS > > > ________________________________ > > Em 01/07/2013 13:16, Edson Lira < edinhoestat(a)yahoo.com.br > escreveu: > Caros amigos da Lista, tenho um banco com as variáveis abaixo: > > Dia_da_Semana Dia_do_Mês Entrada Saída Semana > > No link abaixo está um banco teste para vocês analisarem. > > > http://www1.datafilehost.com/d/320c3ab5 > > > Quero calcular a diferença de tempo entre a entrada e a saída. Para > poder analisar tempo médio por semana. > > Não estou conseguindo implementar nos pacotes chron e lubridate. > > Alguém poderia me ajudar? > > > [ ]'s. > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/8664afa5/atta…> > > ------------------------------ > > Message: 8 > Date: Tue, 02 Jul 2013 12:41:40 +0100 > From: alanarocha(a)sapo.pt > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] trabalhar com médias > Message-ID: > <20130702124140.Horde.yQjHRL5SZhxR0rx0GXnDs6A(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao > longo do tempo. > Ana > > > > ------------------------------ > > Message: 9 > Date: Tue, 02 Jul 2013 12:41:52 +0100 > From: alanarocha(a)sapo.pt > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] trabalhar com médias > Message-ID: > <20130702124152.Horde.fDPba5IzaPJR0ryACYDm7LA(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao > longo do tempo. > Ana > > > > ------------------------------ > > Message: 10 > Date: Tue, 2 Jul 2013 05:33:24 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] trabalhar com médias > Message-ID: > <1372768404.42539.YahooMailNeo(a)web161206.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Allan, um CMR é o que nossos amigos precisam para te ajudar. > > Poderia postar o teu banco em um site tipo dropbox, datafile host e > enviar o link no email para consulta. > > De qualquer forma, talvez você possa usar a rotina abaixo:Examples > > #AGGREGATE > No R > ?aggregate > > ## Compute the averages for the variables in 'state.x77', grouped > ## according to the region (Northeast, South, North Central, West) that > ## each state belongs to. > aggregate(state.x77, list(Region = state.region), mean) ## Compute > the averages according to region and the occurrence of more > ## than 130 days of frost. > aggregate(state.x77, list(Region = state.region, Cold = > state.x77[,"Frost"] > 130), mean) > ## (Note that no state in 'South' is THAT cold.) ## example with > character variables and NAs > testDF <- data.frame(v1 = c(1,3,5,7,8,3,5,NA,4,5,7,9), v2 = > c(11,33,55,77,88,33,55,NA,44,55,77,99) ) > by1 <- c("red","blue",1,2,NA,"big",1,2,"red",1,NA,12) > by2 <- c("wet","dry",99,95,NA,"damp",95,99,"red",99,NA,NA) > aggregate(x = testDF, by = list(by1, by2), FUN = "mean") # and if > you want to treat NAs as a group > fby1 <- factor(by1, exclude = "") > fby2 <- factor(by2, exclude = "") > aggregate(x = testDF, by = list(fby1, fby2), FUN = "mean") ## > Formulas, one ~ one, one ~ many, many ~ one, and many ~ many: > aggregate(weight ~ feed, data = chickwts, mean) > aggregate(breaks ~ wool + tension, data = warpbreaks, mean) > aggregate(cbind(Ozone, Temp) ~ Month, data = airquality, mean) > aggregate(cbind(ncases, ncontrols) ~ alcgp + tobgp, data = esoph, > sum) ## Dot notation: > aggregate(. ~ Species, data = iris, mean) > aggregate(len ~ ., data = ToothGrowth, mean) ## Often followed by xtabs(): > ag <- aggregate(len ~ ., data = ToothGrowth, mean) > xtabs(len ~ ., data = ag) ## Compute the average annual approval > ratings for American presidents. > aggregate(presidents, nfrequency = 1, FUN = mean) > ## Give the summer less weight. > aggregate(presidents, nfrequency = 1, FUN = weighted.mean, w = c(1, > 1, 0.5, 1)) Outra rotina: > No R > ?apply > #APPLY > > ## Compute row and column sums for a matrix: > x <- cbind(x1 = 3, x2 = c(4:1, 2:5)) > dimnames(x)[[1]] <- letters[1:8] > apply(x, 2, mean, trim = .2) > col.sums <- apply(x, 2, sum) > row.sums <- apply(x, 1, sum) > rbind(cbind(x, Rtot = row.sums), Ctot = c(col.sums, sum(col.sums))) > stopifnot( apply(x, 2, is.vector)) ## Sort the columns of a matrix > apply(x, 2, sort) ##- function with extra args: > cave <- function(x, c1, c2) c(mean(x[c1]), mean(x[c2])) > apply(x,1, cave, c1="x1", c2=c("x1","x2")) ma <- matrix(c(1:4, 1, > 6:8), nrow = 2) > ma > apply(ma, 1, table) #--> a list of length 2 > apply(ma, 1, stats::quantile)# 5 x n matrix with rownames > stopifnot(dim(ma) == dim(apply(ma, 1:2, sum))) ## Example with > different lengths for each call > z <- array(1:24, dim=2:4) > zseq <- apply(z, 1:2, function(x) seq_len(max(x))) > zseq ## a 2 x 3 matrix > typeof(zseq) ## list > dim(zseq) ## 2 3 > zseq[1,] > apply(z, 3, function(x) seq_len(max(x))) > # a list without a dim attribute > > > > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: "alanarocha(a)sapo.pt" <alanarocha(a)sapo.pt> > Para: r-br <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Terça-feira, 2 de Julho de 2013 7:41 > Assunto: [R-br] trabalhar com médias > > > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi > 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução > ao longo do tempo. > Ana > > _______________________________________________ > 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/20130702/97771208/atta…> > > ------------------------------ > > Message: 11 > Date: Tue, 2 Jul 2013 09:40:56 -0300 > From: Victor Eduardo <victorduca08(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Dúvida sobre stepwise > Message-ID: > <CAA=hdgaQ+nsWV+N0tVyFFwhBtp-B2ryPzCFPEoQsU9-Ks9Rm8g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá! Desculpe a demora! Muito obrigado pela ajuda! Vou dar uma olhada no > link aqui, qualquer dúvida eu falo. > > > Abraços e boa semana, > > > Victor Eduardo > > > Em 26 de junho de 2013 09:42, walmes . <walmeszeviani(a)gmail.com> escreveu: > >> Dê uma googlada com variações do termo de busca "R stepwise p values >> selection". Vários resultados são retornados. Eu gostei desse. Veja se é >> útil. >> >> >> http://stackoverflow.com/questions/3701170/stepwise-regression-using-p-valu… >> >> À disposição. >> Walmes. >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: walmes(a)ufpr.br >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> 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/20130702/ba5c8e6b/atta…> > > ------------------------------ > > Message: 12 > Date: Tue, 2 Jul 2013 11:47:18 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] [Dúvida] Função dentro de função. > Message-ID: > <CAKwavrkJmG5YZu2u=WzGJGZHOJiX8Ann11Q+69msic0b9LxKOg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Não tenho nenhum exemplo para dar. Minha dúvida é teórica. Gostaria de > saber se é uma boa ou má prática de programação definir função dentro de > função em R. Suponhamos que tenho duas funções f1 e f2. Para f2 rodar ela > precisa de f1. Não preciso de f1 no Environment Global. f1 é apenas uma > função necessária para f2. O "correto" é definir f1 fora ou dentro da > função f2? > [ ], > Pedro Rafael Diniz Marinho. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130702/9060cf66/atta…> > > ------------------------------ > > Message: 13 > Date: Tue, 2 Jul 2013 11:51:22 -0300 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] [Dúvida] Função dentro de função. > Message-ID: > <CAO-arWN4K=AJa-dAYobFandQXF_Bszu9pYwzhprV5NToH2P5ZA(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > e' questao de gosto... e gosto nao se discute. > > eu, particularmente, prefiro nao faze-lo por nao gostar de depender do > escopo lexical... mas ha' quem goste exatamente por isso... > > b > > Em 2 de julho de 2013 11:47, Pedro Rafael > <pedro.rafael.marinho(a)gmail.com> escreveu: >> Não tenho nenhum exemplo para dar. Minha dúvida é teórica. Gostaria de saber >> se é uma boa ou má prática de programação definir função dentro de função em >> R. Suponhamos que tenho duas funções f1 e f2. Para f2 rodar ela precisa de >> f1. Não preciso de f1 no Environment Global. f1 é apenas uma função >> necessária para f2. O "correto" é definir f1 fora ou dentro da função f2? >> [ ], >> Pedro Rafael Diniz Marinho. >> >> _______________________________________________ >> 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. > > > ------------------------------ > > 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 31, assunto 2 > ****************************************

1 0

trabalhar com médias
by alanarocha＠sapo.pt 02 Jul '13

02 Jul '13

vou enviar o ficheiro de dados, para esta pergunta: 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. Aqui adicionar a coluna contendo a pthi T+12 https://dl.dropboxusercontent.com/u/9337305/grupos_pthi.xlsx Ana Quoting r-br-request(a)listas.c3sl.ufpr.br: > 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: interaction plot (Natalia Martins) > 2. diferenca em horas (Edson Lira) > 3. Re: diferenca em horas (andrebvs(a)bol.com.br) > 4. Instalação rmysql ubuntu (Leandro Marino) > 5. Re: Instalação rmysql ubuntu (Daniel C Bezerra) > 6. Re: Instalação rmysql ubuntu (Leandro Marino) > 7. Re: diferenca em horas (Edson Lira) > 8. trabalhar com médias (alanarocha(a)sapo.pt) > 9. trabalhar com médias (alanarocha(a)sapo.pt) > 10. Re: trabalhar com médias (Edson Lira) > 11. Re: Dúvida sobre stepwise (Victor Eduardo) > 12. [Dúvida] Função dentro de função. (Pedro Rafael) > 13. Re: [Dúvida] Função dentro de função. (Benilton Carvalho) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 1 Jul 2013 12:57:09 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: Lista R <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] interaction plot > Message-ID: > <CAGPGL2E0xhe7SPTBB8rXkH4dBVrsjoon6B1YPDLBSkdFbe1myQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Ordenei, mas mesmo assim o R ainda me dá em ordem alfabetica, > mudei o nome mas gostaria de saber se eu consigo editar o grafico, pois nas > opções do interaction nao consegui. > > Att > > > Em 1 de julho de 2013 11:02, Pedro Emmanuel Alvarenga Americano do Brasil < > emmanuel.brasil(a)gmail.com> escreveu: > >> Sugestão 1: formate como fator e ordene os fatores como conveniente. >> >> Sugestão 2: epicalc::followup.plot >> >> Sugestão 3: >> https://stat.ethz.ch/pipermail/r-sig-epi/attachments/20130401/7fcbc976/atta… >> >> >> Dr. Pedro Emmanuel A. A. do Brasil >> http://blog.ipec.fiocruz.br/lapclin-chagas/ >> Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 >> ResearchGate.net: https://www.researchgate.net/profile/Pedro_Brasil2/ >> Instituto Nacional de Infectologia/Instituto de Pesquisa Clínica Evandro >> Chagas >> Fundação Oswaldo Cruz >> Rio de Janeiro - Brasil >> Av. Brasil 4365, >> CEP 21040-360, >> Tel 55 21 3865-9648 >> e-mail: pedro.brasil(a)ipec.fiocruz.br >> e-mail: emmanuel.brasil(a)gmail.com >> >> ---Apoio aos softwares livres >> www.zotero.org - gerenciamento de referências bibliográficas. >> www.broffice.org ou www.libreoffice.org - textos, planilhas ou >> apresentações. >> www.epidata.dk - entrada de dados. >> www.r-project.org - análise de dados. >> www.ubuntu.com - sistema operacional >> >> >> Em 1 de julho de 2013 10:48, Natalia Martins <nsmbarreto(a)gmail.com>escreveu: >> >>> Prezados, bom dia! >>> Estou fazendo uns graficos utilizando a função interaction. plot. >>> O grafico sai em ordem alfabética no eixo x, como eu poderia proceder >>> para que saísse em outra ordem? >>> >>> interaction.plot(CP,Clone,Altura, xlab = 'Corpo de prova', ylab = >>> 'Altura', col=c("blue","red"),lty = 1) >>> >>> >>> O >>> brigada, >>> >>> >>> Natália da Silva Martins >>> >>> >>> _______________________________________________ >>> 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. >> > > > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestra em Ciências: Area de concentração: Estatística e Experimentação > Agronômica - ESALQ/ USP > Doutoranda em Ciências: Area de concentração: Estatística e Experimentação > Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/bbbcb1a6/atta…> > > ------------------------------ > > Message: 2 > Date: Mon, 1 Jul 2013 09:16:45 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] diferenca em horas > Message-ID: > <1372695405.42748.YahooMailNeo(a)web161205.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Caros amigos da Lista, tenho um banco com as variáveis abaixo: > > > Dia_da_Semana Dia_do_Mês > Entrada Saída Semana > > No link abaixo está um banco teste para vocês analisarem. > > > > > http://www1.datafilehost.com/d/320c3ab5 > > > Quero calcular a diferença de tempo entre a entrada e a saída. Para > poder analisar tempo médio por semana. > > Não estou conseguindo implementar nos pacotes chron e lubridate. > > Alguém poderia me ajudar? > > > [ ]'s. > > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/40e14eae/atta…> > > ------------------------------ > > Message: 3 > Date: Mon, 1 Jul 2013 13:36:42 -0300 > From: andrebvs(a)bol.com.br > To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] diferenca em horas > Message-ID: <51d1b01aabf49_a12961f478131(a)a4-winter19.tmail> > Content-Type: text/plain; charset="iso-8859-1" > > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/25fe940c/atta…> > > ------------------------------ > > Message: 4 > Date: Mon, 1 Jul 2013 15:26:50 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAKSaaFkxo81fXiDZ1No8yq0RErKD6grXobjNmMbhO7WMaNSbsQ(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Caros, > > executei a instalação do Rmysql no ubuntu sem sucesso. já executei o > primeiro passo das instruções mas nao sei como fazer o segundo passo... > alguma sugestão? > > Configuration error: > could not find the MySQL installation include and/or library > directories. Manually specify the location of the MySQL > libraries and the header files and re-run R CMD INSTALL. > > INSTRUCTIONS: > > 1. Define and export the 2 shell variables PKG_CPPFLAGS and > PKG_LIBS to include the directory for header files (*.h) > and libraries, for example (using Bourne shell syntax): > > export PKG_CPPFLAGS="-I<MySQL-include-dir>" > export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" > > Re-run the R INSTALL command: > > R CMD INSTALL RMySQL_<version>.tar.gz > > 2. Alternatively, you may pass the configure arguments > --with-mysql-dir=<base-dir> (distribution directory) > or > --with-mysql-inc=<base-inc> (where MySQL header files reside) > --with-mysql-lib=<base-lib> (where MySQL libraries reside) > in the call to R INSTALL --configure-args='...' > > R CMD INSTALL --configure-args='--with-mysql-dir=DIR' > RMySQL_<version>.tar.gz > > ERROR: configuration failed for package ?RMySQL? > * removing ?/usr/local/lib/R/site-library/RMySQL? > > The downloaded source packages are in > ?/tmp/RtmpIdWQKu/downloaded_packages? > Warning message: > In install.packages("RMySQL") : > installation of package ?RMySQL? had non-zero exit status > > a instalação do my sql se deu com o comando: > *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 > php5-mysql phpmyadmin* > * > * > * > * > * > * > <http://www.showmetech.com.br> > > *LEANDRO MARINO* | Showmetech > Estatístico | Fotógrafo > Mobile: 55 21 9845-7707 > Email: contato(a)leandromarino.com.br > Profiles SMT: facebook <http://www.facebook.com/showmetech> | > twitter<http://www.twitter.com/showmetech> > | Google+ <https://plus.google.com/115440786884851522659> > > <http://www.leandromarino.com.br> > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/1fd91564/atta…> > > ------------------------------ > > Message: 5 > Date: Mon, 1 Jul 2013 16:11:49 -0300 > From: Daniel C Bezerra <danielcbezerra(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAJGh7uzcsYu1B9_DuhUDpE2+nODVwd1UobWsBVOZoHsjM3thGw(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Rode no terminal : > > sudo apt-get build-dep r-cran-rmysql > > E depois tente instalar o pacote novamente. > > Abs, > > Daniel > On Jul 1, 2013 3:27 PM, "Leandro Marino" <leandromarino(a)leandromarino.com.br> > wrote: > >> Caros, >> >> executei a instalação do Rmysql no ubuntu sem sucesso. já executei o >> primeiro passo das instruções mas nao sei como fazer o segundo passo... >> alguma sugestão? >> >> Configuration error: >> could not find the MySQL installation include and/or library >> directories. Manually specify the location of the MySQL >> libraries and the header files and re-run R CMD INSTALL. >> >> INSTRUCTIONS: >> >> 1. Define and export the 2 shell variables PKG_CPPFLAGS and >> PKG_LIBS to include the directory for header files (*.h) >> and libraries, for example (using Bourne shell syntax): >> >> export PKG_CPPFLAGS="-I<MySQL-include-dir>" >> export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" >> >> Re-run the R INSTALL command: >> >> R CMD INSTALL RMySQL_<version>.tar.gz >> >> 2. Alternatively, you may pass the configure arguments >> --with-mysql-dir=<base-dir> (distribution directory) >> or >> --with-mysql-inc=<base-inc> (where MySQL header files reside) >> --with-mysql-lib=<base-lib> (where MySQL libraries reside) >> in the call to R INSTALL --configure-args='...' >> >> R CMD INSTALL --configure-args='--with-mysql-dir=DIR' >> RMySQL_<version>.tar.gz >> >> ERROR: configuration failed for package ?RMySQL? >> * removing ?/usr/local/lib/R/site-library/RMySQL? >> >> The downloaded source packages are in >> ?/tmp/RtmpIdWQKu/downloaded_packages? >> Warning message: >> In install.packages("RMySQL") : >> installation of package ?RMySQL? had non-zero exit status >> >> a instalação do my sql se deu com o comando: >> *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 >> php5-mysql phpmyadmin* >> * >> * >> * >> * >> * >> * >> <http://www.showmetech.com.br> >> >> *LEANDRO MARINO* | Showmetech >> Estatístico | Fotógrafo >> Mobile: 55 21 9845-7707 >> Email: contato(a)leandromarino.com.br >> Profiles SMT: facebook <http://www.facebook.com/showmetech> | >> twitter<http://www.twitter.com/showmetech> >> | Google+ <https://plus.google.com/115440786884851522659> >> >> <http://www.leandromarino.com.br> >> >> >> _______________________________________________ >> 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/20130701/37ca7970/atta…> > > ------------------------------ > > Message: 6 > Date: Mon, 1 Jul 2013 17:22:17 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAKSaaF=8D=FQ9vebGcsDZEgn4842zTQ1Z7UfT_eyDw44DKGLtg(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > valeu daniel > funcionou! :) > > vou guardar esta dica! :) > > <http://www.showmetech.com.br> > > *LEANDRO MARINO* | Showmetech > Estatístico | Fotógrafo > Mobile: 55 21 9845-7707 > Email: contato(a)leandromarino.com.br > Profiles SMT: facebook <http://www.facebook.com/showmetech> | > twitter<http://www.twitter.com/showmetech> > | Google+ <https://plus.google.com/115440786884851522659> > > <http://www.leandromarino.com.br> > > > > 2013/7/1 Daniel C Bezerra <danielcbezerra(a)gmail.com> > >> Rode no terminal : >> >> sudo apt-get build-dep r-cran-rmysql >> >> E depois tente instalar o pacote novamente. >> >> Abs, >> >> Daniel >> On Jul 1, 2013 3:27 PM, "Leandro Marino" < >> leandromarino(a)leandromarino.com.br> wrote: >> >>> Caros, >>> >>> executei a instalação do Rmysql no ubuntu sem sucesso. já executei o >>> primeiro passo das instruções mas nao sei como fazer o segundo passo... >>> alguma sugestão? >>> >>> Configuration error: >>> could not find the MySQL installation include and/or library >>> directories. Manually specify the location of the MySQL >>> libraries and the header files and re-run R CMD INSTALL. >>> >>> INSTRUCTIONS: >>> >>> 1. Define and export the 2 shell variables PKG_CPPFLAGS and >>> PKG_LIBS to include the directory for header files (*.h) >>> and libraries, for example (using Bourne shell syntax): >>> >>> export PKG_CPPFLAGS="-I<MySQL-include-dir>" >>> export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" >>> >>> Re-run the R INSTALL command: >>> >>> R CMD INSTALL RMySQL_<version>.tar.gz >>> >>> 2. Alternatively, you may pass the configure arguments >>> --with-mysql-dir=<base-dir> (distribution directory) >>> or >>> --with-mysql-inc=<base-inc> (where MySQL header files reside) >>> --with-mysql-lib=<base-lib> (where MySQL libraries reside) >>> in the call to R INSTALL --configure-args='...' >>> >>> R CMD INSTALL --configure-args='--with-mysql-dir=DIR' >>> RMySQL_<version>.tar.gz >>> >>> ERROR: configuration failed for package ?RMySQL? >>> * removing ?/usr/local/lib/R/site-library/RMySQL? >>> >>> The downloaded source packages are in >>> ?/tmp/RtmpIdWQKu/downloaded_packages? >>> Warning message: >>> In install.packages("RMySQL") : >>> installation of package ?RMySQL? had non-zero exit status >>> >>> a instalação do my sql se deu com o comando: >>> *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 >>> php5-mysql phpmyadmin* >>> * >>> * >>> * >>> * >>> * >>> * >>> <http://www.showmetech.com.br> >>> >>> *LEANDRO MARINO* | Showmetech >>> Estatístico | Fotógrafo >>> Mobile: 55 21 9845-7707 >>> Email: contato(a)leandromarino.com.br >>> Profiles SMT: facebook <http://www.facebook.com/showmetech> | >>> twitter<http://www.twitter.com/showmetech> >>> | Google+ <https://plus.google.com/115440786884851522659> >>> >>> <http://www.leandromarino.com.br> >>> >>> >>> _______________________________________________ >>> 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/20130701/5bace641/atta…> > > ------------------------------ > > Message: 7 > Date: Mon, 1 Jul 2013 21:32:56 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "andrebvs(a)bol.com.br" <andrebvs(a)bol.com.br>, > "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] diferenca em horas > Message-ID: > <1372739576.48970.YahooMailNeo(a)web161203.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > André, valeu! > > O problema é que quando trabalhamos com arquivos texto, csv, as > datas e/ou tempos são vistos como fatores. Não consegui implementar > sua idéia, fiz o seguinte: > > Transformei o arquivo em excell (xlsx) > > Nesta forma o R entende como POSIXct e atribui uma data para as horas > Usei a rotina difftime > > m$hs<-format(as.Date(m$Saída,"H:M")) > m$he<-format(as.Date(m$Entrada,"H:M")) > > m$ht<-difftime(m$Saída,m$Entrada) > m$htrab<-as.numeric(m$ht)/60 > > Muito obrigado!!! > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: "andrebvs(a)bol.com.br" <andrebvs(a)bol.com.br> > Para: r-br(a)listas.c3sl.ufpr.br; Edson Lira <edinhoestat(a)yahoo.com.br> > Enviadas: Segunda-feira, 1 de Julho de 2013 12:36 > Assunto: Re: [R-br] diferenca em horas > > > > Não precisa de pacote, acho que é isso abaixo: > > # Seus dados salvo em "EnterExit.txt": > > Entrada Saída > 13.00 17.00 > 15.01 16.38 > . > . > . > 15.12 17.40 > > dados <- read.table("EnterExit.txt",h=T) > attach(dados) > > dif <- mean(Saída-Entrada) > > [1] 1.767143 > > Att. > André BVS > > > ________________________________ > > Em 01/07/2013 13:16, Edson Lira < edinhoestat(a)yahoo.com.br > escreveu: > Caros amigos da Lista, tenho um banco com as variáveis abaixo: > > Dia_da_Semana Dia_do_Mês Entrada Saída Semana > > No link abaixo está um banco teste para vocês analisarem. > > > http://www1.datafilehost.com/d/320c3ab5 > > > Quero calcular a diferença de tempo entre a entrada e a saída. Para > poder analisar tempo médio por semana. > > Não estou conseguindo implementar nos pacotes chron e lubridate. > > Alguém poderia me ajudar? > > > [ ]'s. > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/8664afa5/atta…> > > ------------------------------ > > Message: 8 > Date: Tue, 02 Jul 2013 12:41:40 +0100 > From: alanarocha(a)sapo.pt > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] trabalhar com médias > Message-ID: > <20130702124140.Horde.yQjHRL5SZhxR0rx0GXnDs6A(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao > longo do tempo. > Ana > > > > ------------------------------ > > Message: 9 > Date: Tue, 02 Jul 2013 12:41:52 +0100 > From: alanarocha(a)sapo.pt > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] trabalhar com médias > Message-ID: > <20130702124152.Horde.fDPba5IzaPJR0ryACYDm7LA(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao > longo do tempo. > Ana > > > > ------------------------------ > > Message: 10 > Date: Tue, 2 Jul 2013 05:33:24 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] trabalhar com médias > Message-ID: > <1372768404.42539.YahooMailNeo(a)web161206.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Allan, um CMR é o que nossos amigos precisam para te ajudar. > > Poderia postar o teu banco em um site tipo dropbox, datafile host e > enviar o link no email para consulta. > > De qualquer forma, talvez você possa usar a rotina abaixo:Examples > > #AGGREGATE > No R > ?aggregate > > ## Compute the averages for the variables in 'state.x77', grouped > ## according to the region (Northeast, South, North Central, West) that > ## each state belongs to. > aggregate(state.x77, list(Region = state.region), mean) ## Compute > the averages according to region and the occurrence of more > ## than 130 days of frost. > aggregate(state.x77, list(Region = state.region, Cold = > state.x77[,"Frost"] > 130), mean) > ## (Note that no state in 'South' is THAT cold.) ## example with > character variables and NAs > testDF <- data.frame(v1 = c(1,3,5,7,8,3,5,NA,4,5,7,9), v2 = > c(11,33,55,77,88,33,55,NA,44,55,77,99) ) > by1 <- c("red","blue",1,2,NA,"big",1,2,"red",1,NA,12) > by2 <- c("wet","dry",99,95,NA,"damp",95,99,"red",99,NA,NA) > aggregate(x = testDF, by = list(by1, by2), FUN = "mean") # and if > you want to treat NAs as a group > fby1 <- factor(by1, exclude = "") > fby2 <- factor(by2, exclude = "") > aggregate(x = testDF, by = list(fby1, fby2), FUN = "mean") ## > Formulas, one ~ one, one ~ many, many ~ one, and many ~ many: > aggregate(weight ~ feed, data = chickwts, mean) > aggregate(breaks ~ wool + tension, data = warpbreaks, mean) > aggregate(cbind(Ozone, Temp) ~ Month, data = airquality, mean) > aggregate(cbind(ncases, ncontrols) ~ alcgp + tobgp, data = esoph, > sum) ## Dot notation: > aggregate(. ~ Species, data = iris, mean) > aggregate(len ~ ., data = ToothGrowth, mean) ## Often followed by xtabs(): > ag <- aggregate(len ~ ., data = ToothGrowth, mean) > xtabs(len ~ ., data = ag) ## Compute the average annual approval > ratings for American presidents. > aggregate(presidents, nfrequency = 1, FUN = mean) > ## Give the summer less weight. > aggregate(presidents, nfrequency = 1, FUN = weighted.mean, w = c(1, > 1, 0.5, 1)) Outra rotina: > No R > ?apply > #APPLY > > ## Compute row and column sums for a matrix: > x <- cbind(x1 = 3, x2 = c(4:1, 2:5)) > dimnames(x)[[1]] <- letters[1:8] > apply(x, 2, mean, trim = .2) > col.sums <- apply(x, 2, sum) > row.sums <- apply(x, 1, sum) > rbind(cbind(x, Rtot = row.sums), Ctot = c(col.sums, sum(col.sums))) > stopifnot( apply(x, 2, is.vector)) ## Sort the columns of a matrix > apply(x, 2, sort) ##- function with extra args: > cave <- function(x, c1, c2) c(mean(x[c1]), mean(x[c2])) > apply(x,1, cave, c1="x1", c2=c("x1","x2")) ma <- matrix(c(1:4, 1, > 6:8), nrow = 2) > ma > apply(ma, 1, table) #--> a list of length 2 > apply(ma, 1, stats::quantile)# 5 x n matrix with rownames > stopifnot(dim(ma) == dim(apply(ma, 1:2, sum))) ## Example with > different lengths for each call > z <- array(1:24, dim=2:4) > zseq <- apply(z, 1:2, function(x) seq_len(max(x))) > zseq ## a 2 x 3 matrix > typeof(zseq) ## list > dim(zseq) ## 2 3 > zseq[1,] > apply(z, 3, function(x) seq_len(max(x))) > # a list without a dim attribute > > > > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: "alanarocha(a)sapo.pt" <alanarocha(a)sapo.pt> > Para: r-br <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Terça-feira, 2 de Julho de 2013 7:41 > Assunto: [R-br] trabalhar com médias > > > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi > 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução > ao longo do tempo. > Ana > > _______________________________________________ > 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/20130702/97771208/atta…> > > ------------------------------ > > Message: 11 > Date: Tue, 2 Jul 2013 09:40:56 -0300 > From: Victor Eduardo <victorduca08(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Dúvida sobre stepwise > Message-ID: > <CAA=hdgaQ+nsWV+N0tVyFFwhBtp-B2ryPzCFPEoQsU9-Ks9Rm8g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá! Desculpe a demora! Muito obrigado pela ajuda! Vou dar uma olhada no > link aqui, qualquer dúvida eu falo. > > > Abraços e boa semana, > > > Victor Eduardo > > > Em 26 de junho de 2013 09:42, walmes . <walmeszeviani(a)gmail.com> escreveu: > >> Dê uma googlada com variações do termo de busca "R stepwise p values >> selection". Vários resultados são retornados. Eu gostei desse. Veja se é >> útil. >> >> >> http://stackoverflow.com/questions/3701170/stepwise-regression-using-p-valu… >> >> À disposição. >> Walmes. >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: walmes(a)ufpr.br >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> 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/20130702/ba5c8e6b/atta…> > > ------------------------------ > > Message: 12 > Date: Tue, 2 Jul 2013 11:47:18 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] [Dúvida] Função dentro de função. > Message-ID: > <CAKwavrkJmG5YZu2u=WzGJGZHOJiX8Ann11Q+69msic0b9LxKOg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Não tenho nenhum exemplo para dar. Minha dúvida é teórica. Gostaria de > saber se é uma boa ou má prática de programação definir função dentro de > função em R. Suponhamos que tenho duas funções f1 e f2. Para f2 rodar ela > precisa de f1. Não preciso de f1 no Environment Global. f1 é apenas uma > função necessária para f2. O "correto" é definir f1 fora ou dentro da > função f2? > [ ], > Pedro Rafael Diniz Marinho. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130702/9060cf66/atta…> > > ------------------------------ > > Message: 13 > Date: Tue, 2 Jul 2013 11:51:22 -0300 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] [Dúvida] Função dentro de função. > Message-ID: > <CAO-arWN4K=AJa-dAYobFandQXF_Bszu9pYwzhprV5NToH2P5ZA(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > e' questao de gosto... e gosto nao se discute. > > eu, particularmente, prefiro nao faze-lo por nao gostar de depender do > escopo lexical... mas ha' quem goste exatamente por isso... > > b > > Em 2 de julho de 2013 11:47, Pedro Rafael > <pedro.rafael.marinho(a)gmail.com> escreveu: >> Não tenho nenhum exemplo para dar. Minha dúvida é teórica. Gostaria de saber >> se é uma boa ou má prática de programação definir função dentro de função em >> R. Suponhamos que tenho duas funções f1 e f2. Para f2 rodar ela precisa de >> f1. Não preciso de f1 no Environment Global. f1 é apenas uma função >> necessária para f2. O "correto" é definir f1 fora ou dentro da função f2? >> [ ], >> Pedro Rafael Diniz Marinho. >> >> _______________________________________________ >> 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. > > > ------------------------------ > > 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 31, assunto 2 > ****************************************

1 0

trabalhar com médias
by alanarocha＠sapo.pt 02 Jul '13

02 Jul '13

vou enviar o ficheiro de dados, para esta pergunta: 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. Aqui adicionar a coluna contendo a pthi T+12 https://dl.dropboxusercontent.com/u/9337305/grupos_pthi.xlsx Ana Quoting r-br-request(a)listas.c3sl.ufpr.br: > 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: interaction plot (Natalia Martins) > 2. diferenca em horas (Edson Lira) > 3. Re: diferenca em horas (andrebvs(a)bol.com.br) > 4. Instalação rmysql ubuntu (Leandro Marino) > 5. Re: Instalação rmysql ubuntu (Daniel C Bezerra) > 6. Re: Instalação rmysql ubuntu (Leandro Marino) > 7. Re: diferenca em horas (Edson Lira) > 8. trabalhar com médias (alanarocha(a)sapo.pt) > 9. trabalhar com médias (alanarocha(a)sapo.pt) > 10. Re: trabalhar com médias (Edson Lira) > 11. Re: Dúvida sobre stepwise (Victor Eduardo) > 12. [Dúvida] Função dentro de função. (Pedro Rafael) > 13. Re: [Dúvida] Função dentro de função. (Benilton Carvalho) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 1 Jul 2013 12:57:09 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: Lista R <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] interaction plot > Message-ID: > <CAGPGL2E0xhe7SPTBB8rXkH4dBVrsjoon6B1YPDLBSkdFbe1myQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Ordenei, mas mesmo assim o R ainda me dá em ordem alfabetica, > mudei o nome mas gostaria de saber se eu consigo editar o grafico, pois nas > opções do interaction nao consegui. > > Att > > > Em 1 de julho de 2013 11:02, Pedro Emmanuel Alvarenga Americano do Brasil < > emmanuel.brasil(a)gmail.com> escreveu: > >> Sugestão 1: formate como fator e ordene os fatores como conveniente. >> >> Sugestão 2: epicalc::followup.plot >> >> Sugestão 3: >> https://stat.ethz.ch/pipermail/r-sig-epi/attachments/20130401/7fcbc976/atta… >> >> >> Dr. Pedro Emmanuel A. A. do Brasil >> http://blog.ipec.fiocruz.br/lapclin-chagas/ >> Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 >> ResearchGate.net: https://www.researchgate.net/profile/Pedro_Brasil2/ >> Instituto Nacional de Infectologia/Instituto de Pesquisa Clínica Evandro >> Chagas >> Fundação Oswaldo Cruz >> Rio de Janeiro - Brasil >> Av. Brasil 4365, >> CEP 21040-360, >> Tel 55 21 3865-9648 >> e-mail: pedro.brasil(a)ipec.fiocruz.br >> e-mail: emmanuel.brasil(a)gmail.com >> >> ---Apoio aos softwares livres >> www.zotero.org - gerenciamento de referências bibliográficas. >> www.broffice.org ou www.libreoffice.org - textos, planilhas ou >> apresentações. >> www.epidata.dk - entrada de dados. >> www.r-project.org - análise de dados. >> www.ubuntu.com - sistema operacional >> >> >> Em 1 de julho de 2013 10:48, Natalia Martins <nsmbarreto(a)gmail.com>escreveu: >> >>> Prezados, bom dia! >>> Estou fazendo uns graficos utilizando a função interaction. plot. >>> O grafico sai em ordem alfabética no eixo x, como eu poderia proceder >>> para que saísse em outra ordem? >>> >>> interaction.plot(CP,Clone,Altura, xlab = 'Corpo de prova', ylab = >>> 'Altura', col=c("blue","red"),lty = 1) >>> >>> >>> O >>> brigada, >>> >>> >>> Natália da Silva Martins >>> >>> >>> _______________________________________________ >>> 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. >> > > > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestra em Ciências: Area de concentração: Estatística e Experimentação > Agronômica - ESALQ/ USP > Doutoranda em Ciências: Area de concentração: Estatística e Experimentação > Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/bbbcb1a6/atta…> > > ------------------------------ > > Message: 2 > Date: Mon, 1 Jul 2013 09:16:45 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] diferenca em horas > Message-ID: > <1372695405.42748.YahooMailNeo(a)web161205.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Caros amigos da Lista, tenho um banco com as variáveis abaixo: > > > Dia_da_Semana Dia_do_Mês > Entrada Saída Semana > > No link abaixo está um banco teste para vocês analisarem. > > > > > http://www1.datafilehost.com/d/320c3ab5 > > > Quero calcular a diferença de tempo entre a entrada e a saída. Para > poder analisar tempo médio por semana. > > Não estou conseguindo implementar nos pacotes chron e lubridate. > > Alguém poderia me ajudar? > > > [ ]'s. > > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/40e14eae/atta…> > > ------------------------------ > > Message: 3 > Date: Mon, 1 Jul 2013 13:36:42 -0300 > From: andrebvs(a)bol.com.br > To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] diferenca em horas > Message-ID: <51d1b01aabf49_a12961f478131(a)a4-winter19.tmail> > Content-Type: text/plain; charset="iso-8859-1" > > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/25fe940c/atta…> > > ------------------------------ > > Message: 4 > Date: Mon, 1 Jul 2013 15:26:50 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAKSaaFkxo81fXiDZ1No8yq0RErKD6grXobjNmMbhO7WMaNSbsQ(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Caros, > > executei a instalação do Rmysql no ubuntu sem sucesso. já executei o > primeiro passo das instruções mas nao sei como fazer o segundo passo... > alguma sugestão? > > Configuration error: > could not find the MySQL installation include and/or library > directories. Manually specify the location of the MySQL > libraries and the header files and re-run R CMD INSTALL. > > INSTRUCTIONS: > > 1. Define and export the 2 shell variables PKG_CPPFLAGS and > PKG_LIBS to include the directory for header files (*.h) > and libraries, for example (using Bourne shell syntax): > > export PKG_CPPFLAGS="-I<MySQL-include-dir>" > export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" > > Re-run the R INSTALL command: > > R CMD INSTALL RMySQL_<version>.tar.gz > > 2. Alternatively, you may pass the configure arguments > --with-mysql-dir=<base-dir> (distribution directory) > or > --with-mysql-inc=<base-inc> (where MySQL header files reside) > --with-mysql-lib=<base-lib> (where MySQL libraries reside) > in the call to R INSTALL --configure-args='...' > > R CMD INSTALL --configure-args='--with-mysql-dir=DIR' > RMySQL_<version>.tar.gz > > ERROR: configuration failed for package ?RMySQL? > * removing ?/usr/local/lib/R/site-library/RMySQL? > > The downloaded source packages are in > ?/tmp/RtmpIdWQKu/downloaded_packages? > Warning message: > In install.packages("RMySQL") : > installation of package ?RMySQL? had non-zero exit status > > a instalação do my sql se deu com o comando: > *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 > php5-mysql phpmyadmin* > * > * > * > * > * > * > <http://www.showmetech.com.br> > > *LEANDRO MARINO* | Showmetech > Estatístico | Fotógrafo > Mobile: 55 21 9845-7707 > Email: contato(a)leandromarino.com.br > Profiles SMT: facebook <http://www.facebook.com/showmetech> | > twitter<http://www.twitter.com/showmetech> > | Google+ <https://plus.google.com/115440786884851522659> > > <http://www.leandromarino.com.br> > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/1fd91564/atta…> > > ------------------------------ > > Message: 5 > Date: Mon, 1 Jul 2013 16:11:49 -0300 > From: Daniel C Bezerra <danielcbezerra(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAJGh7uzcsYu1B9_DuhUDpE2+nODVwd1UobWsBVOZoHsjM3thGw(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Rode no terminal : > > sudo apt-get build-dep r-cran-rmysql > > E depois tente instalar o pacote novamente. > > Abs, > > Daniel > On Jul 1, 2013 3:27 PM, "Leandro Marino" <leandromarino(a)leandromarino.com.br> > wrote: > >> Caros, >> >> executei a instalação do Rmysql no ubuntu sem sucesso. já executei o >> primeiro passo das instruções mas nao sei como fazer o segundo passo... >> alguma sugestão? >> >> Configuration error: >> could not find the MySQL installation include and/or library >> directories. Manually specify the location of the MySQL >> libraries and the header files and re-run R CMD INSTALL. >> >> INSTRUCTIONS: >> >> 1. Define and export the 2 shell variables PKG_CPPFLAGS and >> PKG_LIBS to include the directory for header files (*.h) >> and libraries, for example (using Bourne shell syntax): >> >> export PKG_CPPFLAGS="-I<MySQL-include-dir>" >> export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" >> >> Re-run the R INSTALL command: >> >> R CMD INSTALL RMySQL_<version>.tar.gz >> >> 2. Alternatively, you may pass the configure arguments >> --with-mysql-dir=<base-dir> (distribution directory) >> or >> --with-mysql-inc=<base-inc> (where MySQL header files reside) >> --with-mysql-lib=<base-lib> (where MySQL libraries reside) >> in the call to R INSTALL --configure-args='...' >> >> R CMD INSTALL --configure-args='--with-mysql-dir=DIR' >> RMySQL_<version>.tar.gz >> >> ERROR: configuration failed for package ?RMySQL? >> * removing ?/usr/local/lib/R/site-library/RMySQL? >> >> The downloaded source packages are in >> ?/tmp/RtmpIdWQKu/downloaded_packages? >> Warning message: >> In install.packages("RMySQL") : >> installation of package ?RMySQL? had non-zero exit status >> >> a instalação do my sql se deu com o comando: >> *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 >> php5-mysql phpmyadmin* >> * >> * >> * >> * >> * >> * >> <http://www.showmetech.com.br> >> >> *LEANDRO MARINO* | Showmetech >> Estatístico | Fotógrafo >> Mobile: 55 21 9845-7707 >> Email: contato(a)leandromarino.com.br >> Profiles SMT: facebook <http://www.facebook.com/showmetech> | >> twitter<http://www.twitter.com/showmetech> >> | Google+ <https://plus.google.com/115440786884851522659> >> >> <http://www.leandromarino.com.br> >> >> >> _______________________________________________ >> 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/20130701/37ca7970/atta…> > > ------------------------------ > > Message: 6 > Date: Mon, 1 Jul 2013 17:22:17 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Instalação rmysql ubuntu > Message-ID: > <CAKSaaF=8D=FQ9vebGcsDZEgn4842zTQ1Z7UfT_eyDw44DKGLtg(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > valeu daniel > funcionou! :) > > vou guardar esta dica! :) > > <http://www.showmetech.com.br> > > *LEANDRO MARINO* | Showmetech > Estatístico | Fotógrafo > Mobile: 55 21 9845-7707 > Email: contato(a)leandromarino.com.br > Profiles SMT: facebook <http://www.facebook.com/showmetech> | > twitter<http://www.twitter.com/showmetech> > | Google+ <https://plus.google.com/115440786884851522659> > > <http://www.leandromarino.com.br> > > > > 2013/7/1 Daniel C Bezerra <danielcbezerra(a)gmail.com> > >> Rode no terminal : >> >> sudo apt-get build-dep r-cran-rmysql >> >> E depois tente instalar o pacote novamente. >> >> Abs, >> >> Daniel >> On Jul 1, 2013 3:27 PM, "Leandro Marino" < >> leandromarino(a)leandromarino.com.br> wrote: >> >>> Caros, >>> >>> executei a instalação do Rmysql no ubuntu sem sucesso. já executei o >>> primeiro passo das instruções mas nao sei como fazer o segundo passo... >>> alguma sugestão? >>> >>> Configuration error: >>> could not find the MySQL installation include and/or library >>> directories. Manually specify the location of the MySQL >>> libraries and the header files and re-run R CMD INSTALL. >>> >>> INSTRUCTIONS: >>> >>> 1. Define and export the 2 shell variables PKG_CPPFLAGS and >>> PKG_LIBS to include the directory for header files (*.h) >>> and libraries, for example (using Bourne shell syntax): >>> >>> export PKG_CPPFLAGS="-I<MySQL-include-dir>" >>> export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" >>> >>> Re-run the R INSTALL command: >>> >>> R CMD INSTALL RMySQL_<version>.tar.gz >>> >>> 2. Alternatively, you may pass the configure arguments >>> --with-mysql-dir=<base-dir> (distribution directory) >>> or >>> --with-mysql-inc=<base-inc> (where MySQL header files reside) >>> --with-mysql-lib=<base-lib> (where MySQL libraries reside) >>> in the call to R INSTALL --configure-args='...' >>> >>> R CMD INSTALL --configure-args='--with-mysql-dir=DIR' >>> RMySQL_<version>.tar.gz >>> >>> ERROR: configuration failed for package ?RMySQL? >>> * removing ?/usr/local/lib/R/site-library/RMySQL? >>> >>> The downloaded source packages are in >>> ?/tmp/RtmpIdWQKu/downloaded_packages? >>> Warning message: >>> In install.packages("RMySQL") : >>> installation of package ?RMySQL? had non-zero exit status >>> >>> a instalação do my sql se deu com o comando: >>> *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 >>> php5-mysql phpmyadmin* >>> * >>> * >>> * >>> * >>> * >>> * >>> <http://www.showmetech.com.br> >>> >>> *LEANDRO MARINO* | Showmetech >>> Estatístico | Fotógrafo >>> Mobile: 55 21 9845-7707 >>> Email: contato(a)leandromarino.com.br >>> Profiles SMT: facebook <http://www.facebook.com/showmetech> | >>> twitter<http://www.twitter.com/showmetech> >>> | Google+ <https://plus.google.com/115440786884851522659> >>> >>> <http://www.leandromarino.com.br> >>> >>> >>> _______________________________________________ >>> 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/20130701/5bace641/atta…> > > ------------------------------ > > Message: 7 > Date: Mon, 1 Jul 2013 21:32:56 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "andrebvs(a)bol.com.br" <andrebvs(a)bol.com.br>, > "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] diferenca em horas > Message-ID: > <1372739576.48970.YahooMailNeo(a)web161203.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > André, valeu! > > O problema é que quando trabalhamos com arquivos texto, csv, as > datas e/ou tempos são vistos como fatores. Não consegui implementar > sua idéia, fiz o seguinte: > > Transformei o arquivo em excell (xlsx) > > Nesta forma o R entende como POSIXct e atribui uma data para as horas > Usei a rotina difftime > > m$hs<-format(as.Date(m$Saída,"H:M")) > m$he<-format(as.Date(m$Entrada,"H:M")) > > m$ht<-difftime(m$Saída,m$Entrada) > m$htrab<-as.numeric(m$ht)/60 > > Muito obrigado!!! > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: "andrebvs(a)bol.com.br" <andrebvs(a)bol.com.br> > Para: r-br(a)listas.c3sl.ufpr.br; Edson Lira <edinhoestat(a)yahoo.com.br> > Enviadas: Segunda-feira, 1 de Julho de 2013 12:36 > Assunto: Re: [R-br] diferenca em horas > > > > Não precisa de pacote, acho que é isso abaixo: > > # Seus dados salvo em "EnterExit.txt": > > Entrada Saída > 13.00 17.00 > 15.01 16.38 > . > . > . > 15.12 17.40 > > dados <- read.table("EnterExit.txt",h=T) > attach(dados) > > dif <- mean(Saída-Entrada) > > [1] 1.767143 > > Att. > André BVS > > > ________________________________ > > Em 01/07/2013 13:16, Edson Lira < edinhoestat(a)yahoo.com.br > escreveu: > Caros amigos da Lista, tenho um banco com as variáveis abaixo: > > Dia_da_Semana Dia_do_Mês Entrada Saída Semana > > No link abaixo está um banco teste para vocês analisarem. > > > http://www1.datafilehost.com/d/320c3ab5 > > > Quero calcular a diferença de tempo entre a entrada e a saída. Para > poder analisar tempo médio por semana. > > Não estou conseguindo implementar nos pacotes chron e lubridate. > > Alguém poderia me ajudar? > > > [ ]'s. > Edson Lira > Estatístico > Manaus-Amazonas > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130701/8664afa5/atta…> > > ------------------------------ > > Message: 8 > Date: Tue, 02 Jul 2013 12:41:40 +0100 > From: alanarocha(a)sapo.pt > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] trabalhar com médias > Message-ID: > <20130702124140.Horde.yQjHRL5SZhxR0rx0GXnDs6A(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao > longo do tempo. > Ana > > > > ------------------------------ > > Message: 9 > Date: Tue, 02 Jul 2013 12:41:52 +0100 > From: alanarocha(a)sapo.pt > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] trabalhar com médias > Message-ID: > <20130702124152.Horde.fDPba5IzaPJR0ryACYDm7LA(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao > longo do tempo. > Ana > > > > ------------------------------ > > Message: 10 > Date: Tue, 2 Jul 2013 05:33:24 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] trabalhar com médias > Message-ID: > <1372768404.42539.YahooMailNeo(a)web161206.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Allan, um CMR é o que nossos amigos precisam para te ajudar. > > Poderia postar o teu banco em um site tipo dropbox, datafile host e > enviar o link no email para consulta. > > De qualquer forma, talvez você possa usar a rotina abaixo:Examples > > #AGGREGATE > No R > ?aggregate > > ## Compute the averages for the variables in 'state.x77', grouped > ## according to the region (Northeast, South, North Central, West) that > ## each state belongs to. > aggregate(state.x77, list(Region = state.region), mean) ## Compute > the averages according to region and the occurrence of more > ## than 130 days of frost. > aggregate(state.x77, list(Region = state.region, Cold = > state.x77[,"Frost"] > 130), mean) > ## (Note that no state in 'South' is THAT cold.) ## example with > character variables and NAs > testDF <- data.frame(v1 = c(1,3,5,7,8,3,5,NA,4,5,7,9), v2 = > c(11,33,55,77,88,33,55,NA,44,55,77,99) ) > by1 <- c("red","blue",1,2,NA,"big",1,2,"red",1,NA,12) > by2 <- c("wet","dry",99,95,NA,"damp",95,99,"red",99,NA,NA) > aggregate(x = testDF, by = list(by1, by2), FUN = "mean") # and if > you want to treat NAs as a group > fby1 <- factor(by1, exclude = "") > fby2 <- factor(by2, exclude = "") > aggregate(x = testDF, by = list(fby1, fby2), FUN = "mean") ## > Formulas, one ~ one, one ~ many, many ~ one, and many ~ many: > aggregate(weight ~ feed, data = chickwts, mean) > aggregate(breaks ~ wool + tension, data = warpbreaks, mean) > aggregate(cbind(Ozone, Temp) ~ Month, data = airquality, mean) > aggregate(cbind(ncases, ncontrols) ~ alcgp + tobgp, data = esoph, > sum) ## Dot notation: > aggregate(. ~ Species, data = iris, mean) > aggregate(len ~ ., data = ToothGrowth, mean) ## Often followed by xtabs(): > ag <- aggregate(len ~ ., data = ToothGrowth, mean) > xtabs(len ~ ., data = ag) ## Compute the average annual approval > ratings for American presidents. > aggregate(presidents, nfrequency = 1, FUN = mean) > ## Give the summer less weight. > aggregate(presidents, nfrequency = 1, FUN = weighted.mean, w = c(1, > 1, 0.5, 1)) Outra rotina: > No R > ?apply > #APPLY > > ## Compute row and column sums for a matrix: > x <- cbind(x1 = 3, x2 = c(4:1, 2:5)) > dimnames(x)[[1]] <- letters[1:8] > apply(x, 2, mean, trim = .2) > col.sums <- apply(x, 2, sum) > row.sums <- apply(x, 1, sum) > rbind(cbind(x, Rtot = row.sums), Ctot = c(col.sums, sum(col.sums))) > stopifnot( apply(x, 2, is.vector)) ## Sort the columns of a matrix > apply(x, 2, sort) ##- function with extra args: > cave <- function(x, c1, c2) c(mean(x[c1]), mean(x[c2])) > apply(x,1, cave, c1="x1", c2=c("x1","x2")) ma <- matrix(c(1:4, 1, > 6:8), nrow = 2) > ma > apply(ma, 1, table) #--> a list of length 2 > apply(ma, 1, stats::quantile)# 5 x n matrix with rownames > stopifnot(dim(ma) == dim(apply(ma, 1:2, sum))) ## Example with > different lengths for each call > z <- array(1:24, dim=2:4) > zseq <- apply(z, 1:2, function(x) seq_len(max(x))) > zseq ## a 2 x 3 matrix > typeof(zseq) ## list > dim(zseq) ## 2 3 > zseq[1,] > apply(z, 3, function(x) seq_len(max(x))) > # a list without a dim attribute > > > > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: "alanarocha(a)sapo.pt" <alanarocha(a)sapo.pt> > Para: r-br <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Terça-feira, 2 de Julho de 2013 7:41 > Assunto: [R-br] trabalhar com médias > > > > Olá, > eu tenho estes dados de 4 grupos: > A. pthi <150 B. pthi 150-300 C. pthi > 300-600 D. >600 > 112 277 317 1153 > 126 208 454 1199 > 210 475 2794 > 257 444 2081 > 246 493 1576 > 264 599 901 > 589 871 > 583 722 > 599 1110 > 548 1412 > 454 732 > 382 601 > 533 1130 > 335 642 > 598 682 > 448 1035 > 355 1203 > 462 1136 > 500 1109 > 354 805 > 369 912 > 547 1454 > 842 > 2384 > 760 > 675 > 810 > 1215 > 1565 > 1585 > 942 > 827 > 1015 > 631 > 647 > 709 > 2192 > 1392 > 1028 > 1201 > 892 > 1287 > 1579 > 1071 > 669 > 1103 > 964 > 926 > 1501 > 979 > 1178 > 625 > 678 > 1054 > 855 > 983 > 1076 > 876 > > > que codigos eu posso fazer no R para responder a estas perguntas: > 7. Obter a Dose média em cada grupo e o gráfico da sua evolução > ao longo do tempo. > Ana > > _______________________________________________ > 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/20130702/97771208/atta…> > > ------------------------------ > > Message: 11 > Date: Tue, 2 Jul 2013 09:40:56 -0300 > From: Victor Eduardo <victorduca08(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Dúvida sobre stepwise > Message-ID: > <CAA=hdgaQ+nsWV+N0tVyFFwhBtp-B2ryPzCFPEoQsU9-Ks9Rm8g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá! Desculpe a demora! Muito obrigado pela ajuda! Vou dar uma olhada no > link aqui, qualquer dúvida eu falo. > > > Abraços e boa semana, > > > Victor Eduardo > > > Em 26 de junho de 2013 09:42, walmes . <walmeszeviani(a)gmail.com> escreveu: > >> Dê uma googlada com variações do termo de busca "R stepwise p values >> selection". Vários resultados são retornados. Eu gostei desse. Veja se é >> útil. >> >> >> http://stackoverflow.com/questions/3701170/stepwise-regression-using-p-valu… >> >> À disposição. >> Walmes. >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: walmes(a)ufpr.br >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> 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/20130702/ba5c8e6b/atta…> > > ------------------------------ > > Message: 12 > Date: Tue, 2 Jul 2013 11:47:18 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] [Dúvida] Função dentro de função. > Message-ID: > <CAKwavrkJmG5YZu2u=WzGJGZHOJiX8Ann11Q+69msic0b9LxKOg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Não tenho nenhum exemplo para dar. Minha dúvida é teórica. Gostaria de > saber se é uma boa ou má prática de programação definir função dentro de > função em R. Suponhamos que tenho duas funções f1 e f2. Para f2 rodar ela > precisa de f1. Não preciso de f1 no Environment Global. f1 é apenas uma > função necessária para f2. O "correto" é definir f1 fora ou dentro da > função f2? > [ ], > Pedro Rafael Diniz Marinho. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130702/9060cf66/atta…> > > ------------------------------ > > Message: 13 > Date: Tue, 2 Jul 2013 11:51:22 -0300 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] [Dúvida] Função dentro de função. > Message-ID: > <CAO-arWN4K=AJa-dAYobFandQXF_Bszu9pYwzhprV5NToH2P5ZA(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > e' questao de gosto... e gosto nao se discute. > > eu, particularmente, prefiro nao faze-lo por nao gostar de depender do > escopo lexical... mas ha' quem goste exatamente por isso... > > b > > Em 2 de julho de 2013 11:47, Pedro Rafael > <pedro.rafael.marinho(a)gmail.com> escreveu: >> Não tenho nenhum exemplo para dar. Minha dúvida é teórica. Gostaria de saber >> se é uma boa ou má prática de programação definir função dentro de função em >> R. Suponhamos que tenho duas funções f1 e f2. Para f2 rodar ela precisa de >> f1. Não preciso de f1 no Environment Global. f1 é apenas uma função >> necessária para f2. O "correto" é definir f1 fora ou dentro da função f2? >> [ ], >> Pedro Rafael Diniz Marinho. >> >> _______________________________________________ >> 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. > > > ------------------------------ > > 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 31, assunto 2 > ****************************************

1 0

[Dúvida] Função dentro de função.
by Pedro Rafael 02 Jul '13

02 Jul '13

Não tenho nenhum exemplo para dar. Minha dúvida é teórica. Gostaria de saber se é uma boa ou má prática de programação definir função dentro de função em R. Suponhamos que tenho duas funções f1 e f2. Para f2 rodar ela precisa de f1. Não preciso de f1 no Environment Global. f1 é apenas uma função necessária para f2. O "correto" é definir f1 fora ou dentro da função f2? [ ], Pedro Rafael Diniz Marinho.

3 4

comando ftable
by geovane barbosa 02 Jul '13

02 Jul '13

Olá pessoa estou executando o seguinte comando. ftable(mytable), como eu faço para ele me retornar os percentuais? abraços

6 8

Comparações múltiplas em modelos lineares generalizados
by Gilenio Borges Fernandes 02 Jul '13

02 Jul '13

Prezados(as) Comparações múltiplas (inferência simultânea) feitas através de funções paramétricas com o pacote multcomp do R são válidas para modelos lineares ou para o preditor linear dos modelos lineares generalizados. No caso de modelo de Poisson (ligacao log), modelo logit, etc., os resultados (conclusões) nos preditores lineares podem ser estendidos para a resposta média? Veja o exemplo seguinte do modelo logit com 4 tratamentos e uma covariável onde fiz os intervalos simultâneos na resposta média: Trat = rep(gl(4, 4, label = c('T1','T2', 'T3', 'T4')),3) y=c(1,0,1,1,0,1,0,0,1,1,0,1,1,0,0,0,0,0,0,1,1,1,0,0,0,1,1,1,1,0,0,1,0,0,0,1,1,1,0,0,1,1,1, 1,0,0,1,0) y_ant=c(0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,1,1,0,1,0,0,1,1,1,1,0,0,0,0,0,0,1,1,1,0,0,1, 1,1,1,0,0,1,0) d1=data.frame(Trat,y_ant,y) attach(d1) gmod1 <- glm(y ~ Trat+y_ant, data = d1,family = binomial()) summary(gmod1) # # Comparacao de pares de tratamentos (no preditor linear) library(multcomp) comp1=glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE)) summary(comp1, test = adjusted("Westfall")) # Wastfall (1997) # # Teste e IC's para diferencas de pares de tratamentos no preditor linear gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) summary(gmod1_ci) gmod1_ci$confint # # Comparacao de pares de tratamentos (nas prob condicionais ajustadas) gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) gmod1_ci$confint <- apply(gmod1_ci$confint, 1, binomial()$linkinv) gmod1_ci$confint # gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) gmod1_ci$confint <- apply(gmod1_ci$confint, 2, binomial()$linkinv) gmod1_ci$confint Desde já Grato, -- Gilenio Borges Fernandes Professor Associado III Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

2 1

Pacote(mclust)
by Carlos Pereira 02 Jul '13

02 Jul '13

estou tendo problema com a função densityMclust(data, ...) do pacote(mclust)o objeto data pode ser vetor ou matriz, porem quando uso um vetor diferente do exemplo do pacote dar erro. Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

2 1

Dúvida sobre stepwise
by Victor Eduardo 02 Jul '13

02 Jul '13

Tem como fazer stepwise com teste de significância? Parece que no R, o stepwise faz pelo critério do AIC

2 2

trabalhar com médias
by alanarocha＠sapo.pt 02 Jul '13

02 Jul '13

Olá, eu tenho estes dados de 4 grupos: A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 112 277 317 1153 126 208 454 1199 210 475 2794 257 444 2081 246 493 1576 264 599 901 589 871 583 722 599 1110 548 1412 454 732 382 601 533 1130 335 642 598 682 448 1035 355 1203 462 1136 500 1109 354 805 369 912 547 1454 842 2384 760 675 810 1215 1565 1585 942 827 1015 631 647 709 2192 1392 1028 1201 892 1287 1579 1071 669 1103 964 926 1501 979 1178 625 678 1054 855 983 1076 876 que codigos eu posso fazer no R para responder a estas perguntas: 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. Ana

2 1

trabalhar com médias
by alanarocha＠sapo.pt 02 Jul '13

02 Jul '13

Olá, eu tenho estes dados de 4 grupos: A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 112 277 317 1153 126 208 454 1199 210 475 2794 257 444 2081 246 493 1576 264 599 901 589 871 583 722 599 1110 548 1412 454 732 382 601 533 1130 335 642 598 682 448 1035 355 1203 462 1136 500 1109 354 805 369 912 547 1454 842 2384 760 675 810 1215 1565 1585 942 827 1015 631 647 709 2192 1392 1028 1201 892 1287 1579 1071 669 1103 964 926 1501 979 1178 625 678 1054 855 983 1076 876 que codigos eu posso fazer no R para responder a estas perguntas: 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. Ana

1 0

diferenca em horas
by Edson Lira 02 Jul '13

02 Jul '13

Caros amigos da Lista, tenho um banco com as variáveis abaixo: Dia_da_Semana Dia_do_Mês Entrada Saída Semana No link abaixo está um banco teste para vocês analisarem. http://www1.datafilehost.com/d/320c3ab5 Quero calcular a diferença de tempo entre a entrada e a saída. Para poder analisar tempo médio por semana. Não estou conseguindo implementar nos pacotes chron e lubridate. Alguém poderia me ajudar? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

Instalação rmysql ubuntu
by Leandro Marino 01 Jul '13

01 Jul '13

Caros, executei a instalação do Rmysql no ubuntu sem sucesso. já executei o primeiro passo das instruções mas nao sei como fazer o segundo passo... alguma sugestão? Configuration error: could not find the MySQL installation include and/or library directories. Manually specify the location of the MySQL libraries and the header files and re-run R CMD INSTALL. INSTRUCTIONS: 1. Define and export the 2 shell variables PKG_CPPFLAGS and PKG_LIBS to include the directory for header files (*.h) and libraries, for example (using Bourne shell syntax): export PKG_CPPFLAGS="-I<MySQL-include-dir>" export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" Re-run the R INSTALL command: R CMD INSTALL RMySQL_<version>.tar.gz 2. Alternatively, you may pass the configure arguments --with-mysql-dir=<base-dir> (distribution directory) or --with-mysql-inc=<base-inc> (where MySQL header files reside) --with-mysql-lib=<base-lib> (where MySQL libraries reside) in the call to R INSTALL --configure-args='...' R CMD INSTALL --configure-args='--with-mysql-dir=DIR' RMySQL_<version>.tar.gz ERROR: configuration failed for package ‘RMySQL’ * removing ‘/usr/local/lib/R/site-library/RMySQL’ The downloaded source packages are in ‘/tmp/RtmpIdWQKu/downloaded_packages’ Warning message: In install.packages("RMySQL") : installation of package ‘RMySQL’ had non-zero exit status a instalação do my sql se deu com o comando: *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 php5-mysql phpmyadmin* * * * * * * <http://www.showmetech.com.br> *LEANDRO MARINO* | Showmetech Estatístico | Fotógrafo Mobile: 55 21 9845-7707 Email: contato(a)leandromarino.com.br Profiles SMT: facebook <http://www.facebook.com/showmetech> | twitter<http://www.twitter.com/showmetech> | Google+ <https://plus.google.com/115440786884851522659> <http://www.leandromarino.com.br>

2 2

interaction plot
by Natalia Martins 01 Jul '13

01 Jul '13

Prezados, bom dia! Estou fazendo uns graficos utilizando a função interaction. plot. O grafico sai em ordem alfabética no eixo x, como eu poderia proceder para que saísse em outra ordem? interaction.plot(CP,Clone,Altura, xlab = 'Corpo de prova', ylab = 'Altura', col=c("blue","red"),lty = 1) O brigada, Natália da Silva Martins

2 2