Julho 2013 - R-br - DINF UFPR

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 12 Jul '13

12 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 12 Jul '13

12 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