- R-br - DINF UFPR

Cursos sobre R no Coursera (Off-topic)
by Augusto Ribas 06 Abr '14

06 Abr '14

Pessoal, o Coursera está oferecendo uma especialização em data-science focada em R. Assistir as aulas é gratuito. Você pode pagar se quiser um certificado validado, mas isso não é obrigatório. https://www.coursera.org/specialization/jhudatascience/1?utm_medium=listing… Ela é composta de nove cursos, todos sobre R, abordando vários tópicos diferentes. Amanha começam os três primeiros cursos dessa especialização. Getting and Cleaning Data https://www.coursera.org/course/getdata R Programming https://www.coursera.org/course/rprog The Data Scientist's Toolbox https://www.coursera.org/course/datascitoolbox Eu tenho feito os cursos oferecidos pela John Hopkins no coursera e tenho achado bem legal. Normalmente tem uma seleção de tópicos interessantes, exemplos bem amplos, tornando os cursos legais para diferentes áreas do conhecimento, e todos os cursos podem ser acompanhados gratuitamente. Eu acho uma oportunidade impar para aprender em casa e conhecer novas areas. Já que acredito que muita gente, deve ter dificuldade em locomoção e não pode sair pelo pais assistindo todos os cursos de R legais que temos por aqui. Achei que seria interessante comentar aqui, talvez mais pessoas se interessem. Espero que não se importem com o off-topic. Abraços a todos e continuem com as interessantes discussões na lista que também são uma excelente fonte de aprendizado. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Plotar várias curvas de densidades num único objeto gráfico
by Daniel Marcelino 05 Abr '14

05 Abr '14

Caro, estou com um problema para plotar densidades de distribuição posterior e talvez alguém aqui já tenha resolvido isso de alguma maneira. Tenho um objeto de dados de classe mcmc.list com vários paramentros, entre os quais eu gostaria de sobrepor três distribuições posteriores como abaixo. Do jeito que está, contudo, não funciona porque a cada comando plot o R mantem a curva no centro e altera as escalas. Tem alguma forma de plotar essas densidades conservando as posições delas originais? Coloquei alguns dados aqui: https://dl.dropboxusercontent.com/u/1339742/mcmc.samples.Rdata quartz(width=10,height=7) plot(mcmc.samples[, c(paste("alphaPSDB[", 2, "]", sep = ""))], trace = FALSE, density = TRUE, smooth = TRUE, lwd=3, col="blue") par(new=TRUE) plot(mcmc.samples[, c(paste("alphaPT[", 2, "]", sep = ""))], trace = FALSE, density = TRUE, smooth = TRUE, lwd=3, col="red") par(new=TRUE) plot(mcmc.samples[, c(paste("alphaPRB[", 2, "]", sep = ""))], trace = FALSE, density = TRUE, smooth = TRUE, lwd=3, col="green")

2 2

Pegar dados de uma tabela em html
by Sérgio Henrique almeida da silva ju 04 Abr '14

04 Abr '14

Prezados Como faço para pegar os dados de uma tabela que está dentro de uma Frame em HTML usando o R? O site é esse: http://portal.inep.gov.br/basica-censo-escolar-matriculanele eu preciso dizer o estado que eu quero os dados e a pesquisa me retorna uma tabela com os dados. Não existe nenhuma opção para exportar essa tabela. Sei que existe o Rcurl e o XML, mas não tenho ideia de como fazer isso. Segue um pedaço da frame: <td >Federal</td> <td align="right" >322</td> <td align="right" >334</td> <td align="right" >3541</td> <td align="right" >6855</td> <td align="right" > 16910</td> <td align="right" >10427</td> <td align="right" >48</td> <td align="right" >1301</td> <td align="right" >0</td> <td align="right" >0</td> <td align="right" >9</td> <td align="right" >70</td> <td align="right" >361</td> <td align="right" >311</td> <td align="right" >209</td> <td align="right" >83 </td> <td align="right" >0</td> <td align="right" >2</td> -- 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) 968463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 3

AIC para variograma ajustado por OLS (variofit)
by Éder Comunello 04 Abr '14

04 Abr '14

Senhores, boa tarde! Ao utilizar a geoR temos a possibilidade de ajustes de modelos para variogramas por variofit() e likfit(). Eu costumava usar Grau de Dependência Espacial (GDE) e somatório de erro como critérios para seleção de modelos ajustados por variofit() e principalmente Akaike Information Citerion (AIC) para modelos ajustados por likfit(). Todavia, recentemente estudando a função loglik.GRF(), vi que é possível calcular um valor log-likelihood para modelos ajustados por variofit(). Minha dúvida é se faz sentido empregar esse valor no cálculo de um AIC para modelos ajustados por variofit(), uma vez que não encontrei embasamento na literatura. Uma segunda dúvida é no número de parâmetros utilizados no cálculo do AIC. Observei que para obter um valor igual a o cálculado pela geoR devo considerar um parâmetro a mais do que vinha considerando ("tausq", "sigmasq" e "phi"). Poderia considerar o erro a ser estimado como um quarto parâmetro do modelo? Agradeço qualquer ajuda... Segue o CMR para apoiar a análise (modificado do help da função)... ### <code r> require(geoR) # ?loglik.GRF ## Computing the likelihood of a variogram fitted model s100.v <- variog(s100, max.dist=1); plot(s100.v) s100.vf1 <- variofit(s100.v, ini=c(1,.5), weights="equal", cov="exp"); s100.vf1 # variofit: model parameters estimated by OLS (ordinary least squares): # covariance model is: exponential # parameter estimates: # tausq sigmasq phi # 0.1540 1.1505 0.6621 # Practical Range with cor=0.05 for asymptotic range: 1.983401 # # variofit: minimised sum of squares = 0.0801 s100.vf2 <- variofit(s100.v, ini=c(1,.5), weights="equal", cov="sph"); s100.vf2 # variofit: model parameters estimated by OLS (ordinary least squares): # covariance model is: spherical # parameter estimates: # tausq sigmasq phi # 0.1841 0.8179 0.9487 # Practical Range with cor=0.05 for asymptotic range: 0.9486608 # # variofit: minimised sum of squares = 0.0727 plot(s100.v) lines.variomodel(s100.vf1, col=2) lines.variomodel(s100.vf2, col=3) lk1 <- loglik.GRF(s100, obj=s100.vf1); lk1 lk2 <- loglik.GRF(s100, obj=s100.vf2); lk2 # ?AIC # AIC = -2*log-likelihood + k*npar # log-like: maximized value of the likelihood function for the estimated model # parameters <- c("tausq", "sigmasq", "phi") + erro? # npar=3 npar=4 k=2 # for usual AIC aic1 <- -2*lk1+k*npar; aic1 aic2 <- -2*lk2+k*npar; aic2 ### </code> Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

4 9

Sobre distribuições em Estatística
by Listeiro 037 03 Abr '14

03 Abr '14

Saudações. Estou pesquisando sobre distribuições em Estatística e necessito de uma recomendação de uma página de internet ou livro em que estejam dispostas. Mais além que a Gaussiana, gostaria de observar distribuições como a elíptica, a parabólica e a hiperbólica. Comparar se há ligação com as formas quadráticas relacionadas e com as respectivas equações diferenciais parciais de nome semelhante. Não sei se algo próximo de esgotar o tema, seria o mais próximo disto. Desde já agradeço qualquer contribuição.

3 3

Outlier
by Cassiano 03 Abr '14

03 Abr '14

Prezados, Preciso "quebrar" o eixo y para mostrar um ponto que está muito fora do padrão geral. A maioria dos meus pontos variam entre 0 e 100, mas tenho um ponto em 780. Não gostaria de descartar este ponto, pois é real. Mas mostrar os dados na escala em que ele aparece, "esconde" o padrão geral dos outros pontos. Já vi em algum lugar que tem como "quebrar" (não sei o termo correto) o eixo y. Tipo, mostro de 0 a 100, faço uma "quebra" de 110 a 750, e volto a mostrar o eixo até 800, por exemplo. Vocês podem me ajudar nesta tarefa? Alguém sabe comandos para tal façanha? Estou usando o pacote "car" >library(car) >with(dados, interaction.plot(COL,Amb, ABUND, type="b", pch=15+as.numeric(Amb),col=15+as.numeric(Amb))) Obrigado -- ======================================= Cassiano S. Rosa Professor Visitante Universidade do Estado do Mato Grosso - UNEMAT Programa de Pós-graduação em Ecologia e Conservação Br 158, Km 655 - Caixa Postal 08 Antiga Sede da FAB UNEMAT - *Campus* Universitário de Nova Xavantina Nova Xavantina - MT CEP: 78690-000 Telefone: (66)3438-1224 ou 3438-2670 - FAX: (66) 3438-1224 Curriculo lattes: http://lattes.cnpq.br/3030825329622014 Scholar Google: http://scholar.google.com/citations?user=JmELzOMAAAAJ =======================================

2 1

Coeficiente de regressão parcial modelo misto linear
by Fernando Souza 03 Abr '14

03 Abr '14

Caros amigos Eu gerei um modelo de regressão utilizando o modelo misto com o comando abaixo. Eu gostaria de saber o quanto cada uma das 3 variáveis preditoras contribuiram para a redução da variação do Consumo (R2 parcial), ou seja particionar a variância total de Y em função dos preditores. Gostaria de saber como posso fazer isso? Se há ou se alguém possui um código de como fazê-lo. Li no livro do Kutner que o r2parcial é calculado como (SSE(X2)-SSE(X1,X2))/SSE(X2) para o R2 da variável x1 para um modelo com duas variáveis. No entanto no modelo mixto não consigo gerar modelos com ou sem as preditoras para obter a soma de quadrados , pois nem sempre há convergência, mantendo todos os argumentos da função lme constantes. Sendo assim há uma forma de obter esses valores? modelo0<-lme(fixed=CONSUMO~I(PV^0.75)+GEST+TEMPMAX,correlation=corAR1(form=~1+GEST|ANIMAL),weights=varIdent(form=~1|ANIMAL),data=dados,random=~1|ANIMAL,control=ctrl,na.action=na.omit,method="REML")

1 0

Parâmetros gráficos na Carta de Controle
by Edson Lira 02 Abr '14

02 Abr '14

Com a rotina abaixo consegui alterar os parâmetros do eixo x, do eixo y, só não consegui alterar os parâmetros (não encontrei nada na ajuda) do eixo da direita que contém os labels dos limites de controle (UCL, LCL), gostaria de reduzir o tamanho da fonte dos mesmos. require(qcc) old <- qcc.options() # save defaults qcc.options("cex.stats") # get a single parameter qcc.options("cex.stats"=0.8) # change parameters qcc.options("cex.axis"=0.7) # change parameters qcc.options(bg.margin="white") par(cex=0.7,cex.lab=0.8,cex.axis=0.7) qcc(M04, nsigmas=2,cex=0.7,cex.lab=0.8,cex.axis=0.7,title="",xlab="Dia",ylab="Tempo (min)") Alguma sugestão? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

Criar função utilizando "for"
by Rodrigo Muniz 25 Mar '14

25 Mar '14

Caros membros da lista, Me desculpem a ignorância, porém sou novo como usuário e estou tentando criar uma função para facilitar uma operação no meu banco de dados. considerem o exemplo: dados = as.data.frame(matrix(rnorm(200), ncol=10, nrow=50)) cor.multipla <- function(x, y){ n <- ncol(y) cor.mult <- numeric(0) + for(i in 1:n){cor.mult<- cor.test(x, y[,i]) + return(cor.mult$p.value)} } cor.multipla(dados$V1, dados[2:10]) Nesta função, eu gostaria fazer a correlação da coluna 1 com as colunas de 2 a 10, e extrair um vetor de p.value para cada correlação, no entanto, a função está me retornando apenas um valor de p.value. Alguém poderia me dizer onde eu estou errando ou mesmo se estou indo pelo caminho errado? Desde já agradeço... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 5

rodando código R na web
by Paulo Nogueira Starzynski 25 Mar '14

25 Mar '14

Caros colegas de lista, algum tempo atrás deparei-me com uma ferramenta online onde era possível executar códigos em R. Na época não guardei nos favoritos e agora não consigo voltar a encontrar. Lembra um pouco essa outra, para Python: http://www.codeskulptor.org/ Alguém sabe qual é e pode me passar a localização? Abraços, Paulo Nogueira Starzynski

4 4

Grárico de controle
by Edson Lira 25 Mar '14

25 Mar '14

Caros amigos da lista. Tenho 3507 medidas de tempo de doação distribuídas por dias conforme a tabela abaixo. 01 02 03 04 05 07 08 09 10 11 14 15 16 17 18 19 21 22 23 25 26 28 29 30 31 92 165 154 113 237 146 139 156 92 149 147 142 120 123 142 125 111 94 84 85 73 91 140 130 105 Gostaria de plotar os gráficos xbarra e o da amplitude para esses tempos. Usei a rotina abaixo como o Eder me indicou: amostra = sample(rep(1:631,5),3155) #rever esses parâmetros, pois só postei uma amostra require(qcc) tem <- qcc.groups(tempo,amostra) Plotei o gráfico e funcionou, me mostrando os tempos e os grupos de amostra. Funcioniou, só que me deparei com outro problema, essa rotina me dá o tempo e o número da amostra, como faço para ver o dia em que o tempo ou amostra de tamanho 5 foi registrada. Ou seja, necessito mostrar o gráfico com o tempo e os dias do registro do tempo em amostras de tamanho 5. No link abaixo tem um link para o arquivo que estou usando com 300 observações. http://www.datafilehost.com/d/e38f36e0 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 12

Tirar maiuscula apos o ponto
by Bernardo Rangel Tura 25 Mar '14

25 Mar '14

Boa noite a todos Primeiro de tudo quero pedir desculpas por uma duvida tao basica considerem o seguinte exemplo: texto <- 'Deus abencoou Noe e seus filhos, e lhes disse: "Sede fecundos, multiplicai-vos, enchei a terra. Sede o medo e o pavor de todos os animais da terra e de todas as aves do ceu, como de tudo o que se move na terra e de todos os peixes do mar: eles sao entregues nas vossas maos. Tudo o que se move e possui a vida vos servira de alimento, tudo isso eu vos dou, como vos dei a verdura das plantas. Mas nao comereis a carne com sua alma, isto e, o sangue. Pedirei contas porem, do sangue de cada um de vos. Pedirei contas a todos os animais e ao homem, aos homens entre si, eu pedirei contas da alma do homem. Quem derrama o sangue do homem pelo homem terra seu sangue derramado. Pois a imagem de Deus o homem foi feito. Quanto a vos, sede fecundos, multiplicai-vos, povoai a terra e dominai-a." Agora desejo transformar as maiusculas apos o ponto em minusuclas, por exemplo: ". Sede o medo" em ". sede o medo" sei que posso fazer com gsub e pensei que gsub('. [A-Z]','. [a-z]',texto) seria suficiente, mas nao foi. Onde estou errando? Antecipadamente agradeco aos R-mestres []s Tura

6 13

Re: [R-br] Criar função utilizando "for"
by Robert Iquiapaza 24 Mar '14

24 Mar '14

Caso se sinta mais confortável com código explícito: dados = data.frame(matrix(rnorm(500), ncol=10, nrow=50)) cor.multipla <- function(x, y){ x=as.numeric(x) n <- ncol(y) y=as.matrix(y,ncol=n) cor.mult <- vector() for(i in 1:n){cor.mult[i]<- cor.test(x, y[,i])$p.value} return(cor.mult) } dim(dados) cor.multipla(dados$X1, dados[2:10]) Sds Robert -----Mensagem Original----- From: Paulo Justiniano Sent: Monday, March 24, 2014 10:08 PM To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Criar função utilizando "for" O seu loop for tá sobre escrevendo resultados em cor.mult<- cor.test(x, y[,i]) e portanto somente o ultimo está sendo guardado (note tb que sua matrz de dados esá reciclando valores) Uma sugestão: x <- 1:10 y <- matrix(rnorm(200), nrow=10) apply(y, 2, function(y, x=x) cor.test(x,y)$p.value On Mon, 24 Mar 2014, Rodrigo Muniz wrote: > Caros membros da lista, Me desculpem a ignorância, porém sou novo como > usuário e estou tentando criar uma função para facilitar uma operação no > meu banco de dados. > considerem o exemplo: > > dados = as.data.frame(matrix(rnorm(200), ncol=10, nrow=50)) > cor.multipla <- function(x, y){ > n <- ncol(y) > cor.mult <- numeric(0) > + for(i in 1:n){cor.mult<- cor.test(x, y[,i]) > + return(cor.mult$p.value)} > } > > cor.multipla(dados$V1, dados[2:10]) > > Nesta função, eu gostaria fazer a correlação da coluna 1 com as colunas de > 2 a 10, e extrair um vetor de p.value para cada correlação, no entanto, > a função está me retornando apenas um valor de p.value. > Alguém poderia me dizer onde eu estou errando ou mesmo se estou indo pelo > caminho errado? > Desde já agradeço... > -- > Rodrigo A. MunizEng. Agrônomo. Ms Produção Vegetal (UENF) > Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) > Cel (19) 8300-4333 (Pessoal) > Tel (19) 3375-1196 (Residêncial) > Tel (19) 3447-8537 (Trabalho) > > _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

problemas com comparações múltiplas
by Maurício Lordêlo 24 Mar '14

24 Mar '14

Saudações a todos! Estou realizando inferência simultânea para um conjunto de dados correlacionados onde desejo comparar o efeito de quatro tratamentos. Para isto estou utilizando a função "simplesimint" do pacote "BSagri". Consigo fazer utilizando Dunnet porém quando mudo para Tukey ocorre um problema. Abaixo tenho um CMR onde tenho o vetor de coeficientes estimados (coefs), a matriz de variância e covariâncias (VarCov), o número de obervações em cada grupo (n) e a matriz de contrastes (cmat). O problema principal surge quando faço utilizando todas as linhas (são seis) da matriz de contrastes (cmat2). Já para duas, três e até quatro linhas da matriz de contrastes não há problemas. Grato caso alguém possa ajudar. Abraço library(BSagri) coefs=c("TratamentoA"=-1.0618,"TratamentoB"=1.1470,"TratamentoC"=0.3068,"TratamentoC"=0.2501) VarCov=matrix(data = NA, nrow = 4, ncol = 4, byrow = FALSE,dimnames=NULL) VarCov[1,]=c(0.3133311,-0.1618691,-0.1736806,-0.1712345) VarCov[2,]=c(-0.1618697,0.4510964,0.1595845,0.1622442) VarCov[3,]=c(-0.1736809,0.1595843,0.2729564,0.1617069) VarCov[4,]=c(-0.1712371,0.1622446,0.1617080,0.2757233) VarCov n=c("TratamentoA"=40,"TratamentoB"=20,"TratamentoC"=60,"TratamentoD"=60) cmat1<-contrMat(n=n, type="Dunnett") cmat1 test1<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat1, alternative="two.sided") plotCI(test1) cmat2<-contrMat(n=n, type="Tukey") cmat2 test2<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat2, alternative="two.sided") ##ERRO!! #Erro em uniroot(pfct, interval = interval) : # f() values at end points not of opposite sign test3<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:2,1:4], alternative="two.sided") summary(test3) plotCI(test3) test4<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:3,1:4], alternative="two.sided") summary(test4) plotCI(test4) test5<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:4,1:4], alternative="two.sided") summary(test5) plotCI(test5) test6<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:5,1:4], alternative="two.sided") ##ERRO!! #Erro em uniroot(pfct, interval = interval) : # f() values at end points not of opposite sign

1 0

função BiCopKPlot
by Natalia Martins 24 Mar '14

24 Mar '14

Prezados boa tarde, estou estudando a função BiCopKPlot utilizada para a modelagem de cópulas. Nesta função tem-se: function (u1, u2, PLOT = TRUE, ...) { if (is.null(u1) == TRUE || is.null(u2) == TRUE) stop("u1 and/or u2 are not set or have length zero.") if (length(u1) != length(u2)) stop("Lengths of 'u1' and 'u2' do not match.") if (length(u1) < 2) stop("Number of observations has to be at least 2.") if (PLOT != TRUE && PLOT != FALSE) stop("The parameter 'PLOT' has to be set to 'TRUE' or 'FALSE'.") Wi <- W(u1, u2) Hi <- H(u1, u2) Hi.sort <- sort(Hi) n <- length(u1) W.in <- rep(NA, n) for (i in 1:n) { f = function(w) { w * (-log(w)) * (w - w * log(w))^(i - 1) * (1 - w + w * log(w))^(n - i) } W.in[i] <- n * choose(n - 1, i - 1) * (integrate(f, lower = 0, upper = 1)$value) } g <- function(w) { w - w * log(w) } if (PLOT) { plot(g, xlim = c(0, 1), ylim = c(0, 1), pch = "x", xlab = expression(W[1:n]), ylab = "H", ...) points(W.in, Hi.sort, pch = "x", cex = 0.4, ...) abline(a = 0, b = 1) } else { kendall.plot.output <- list(W.in, Hi.sort) names(kendall.plot.output) <- c("W.in", "Hi.sort") return(kendall.plot.output) } } Porem, não encontrei no pdf da library(CDVine) como é definida a função H e W (em destaque). Por favor se alguem souber e puder me ajudar. Muito obrigada! -- Natália Martins Contato: (35) 91812482

3 3

Python no Coursera [offtopic]
by Paulo Nogueira Starzynski 24 Mar '14

24 Mar '14

Pessoal, me desculpem caso a mensagem esteja fora do propósito de uma lista sobre R, mas creio que Python seja interessante para muitos dos usuários do R. Começou essa semana um curso introdutório sobre a linguagem no Coursera. Pode valer a pena. Eu vou fazer. https://www.coursera.org/course/interactivepython Abraços, Paulo Nogueira Starzynski

1 0

Re: [R-br] Digest R-br, volume 39, assunto 21
by Marcos Batista 21 Mar '14

21 Mar '14

Prezados, sou biólogo e um dos objetivos da minha pesquisa é examinar os efeitos da heterogeneidade do ambiente na variabilidade genética. Para isso, eu estimei a distância genética, geodésica e ecológica e usei uma função do R - http://datadryad.org/resource/doi:10.5061/dryad.kt71r/1 - para testar a regressão múltipla dessas matrizes. O script para realizar o teste eu consegui fazer, entretanto, não encontrei uma função para fazer o gráfico da regressão linear sobre meus dados. O único gráfico que consigo fazer é entre duas variáveis é usando a função splom (matriz1, matriz2, lm=TRUE) e obtive uma resposta de erro (lm object not found). Mas, o que eu queria era fazer um gráfico de uma matriz1 x matriz2+matriz3. Alguém tem uma sugestão de função que eu possa usar? Agradeço a todos desde já. Cordialmente, Marcos Batista Marcos Roberto Dias Batista Departamento de Genética, Evolução e Bioagentes Instituto de Biologia, Unicamp Caixa Postal 6109, Campinas, São Paulo CEP 13083-970, Brasil Telefone: (19) 3521 1151 2014-03-21 12:00 GMT-03:00 <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: Calendário do pacote googleVis (Éder Comunello) > 2. Re: Grárico de controle (Edson Lira) > 3. Iteração com optim (jborssoi .) > 4. Re: Iteração com optim (Diogo Ferrari) > 5. Re: Grárico de controle (Éder Comunello) > 6. Re: Iteração com optim (walmes .) > 7. Re: Iteração com optim (Wagner Bonat) > 8. Re: Iteração com optim (jborssoi .) > 9. Re: Iteração com optim (Wagner Bonat) > 10. Re: Grárico de controle (Edson Lira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 21 Mar 2014 07:17:25 -0400 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Calendário do pacote googleVis > Message-ID: > <CABmC8gnNgPU2mWLFczTUrWYsqOf2A9OjKv5dpF7UnNM= > aTie6A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá, pessoal! Bom dia! > > Tudo indica que o problema está na API do Google e não propriamente no R. O > link que segue relata o mesmo problema numa aplicação fora do R, feita > também por um brasileiro, reforçando que o problema está ligado ao > 'Brazilian > Portuguese'. > > > https://code.google.com/p/google-visualization-api-issues/issues/detail?id=… > > Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> > Dourados, MS - [22 16.5'S, 54 49'W] > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/a91fe90a/atta… > > > > ------------------------------ > > Message: 2 > Date: Fri, 21 Mar 2014 05:47:41 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "ecomunel(a)gmail.com" <ecomunel(a)gmail.com>, > "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Grárico de controle > Message-ID: > <1395406061.56210.YahooMailNeo(a)web160103.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Eder, muito obrigado, ajudou muito, vi que tenho no mesmo dia vários > grupos. > > Como faço para ter um único resultado por dia? Ou seja, para cada dia um > único ponto no gráfico. > > > [ ]'s. > > Edson Lira > Estatístico > Manaus-Amazonas > > > > Em Terça-feira, 18 de Março de 2014 12:52, Éder Comunello < > comunello.eder(a)gmail.com> escreveu: > > Edson, boa tarde! > > Não sei se irá lhe atender, mas tentei alguma coisa com os dados de > exemplo que você havia enviado (n=300). > > ### <code r> > ### restaurando dados de exemplo! > dat <- > > data.frame(tempo=c(30.48,32.17,32.32,32.68,33.42,34.1,34.33,35.12,35.95,36.18,36.73,36.75,36.78,37.05,37.18,38.1,38.4,40.05,40.62,40.87,41.33,41.38,41.5,41.63,43.38,43.65,43.85,43.88,44.37,45.67,45.73,46.17,46.5,46.75,47.17,47.55,47.63,47.65,49.62,49.85,49.98,50.65,51.15,51.67,51.7,51.97,53.12,53.23,53.77,54.08,54.18,54.32,54.65,54.98,55.55,56.47,58.13,58.65,58.93,59.2,59.67,59.92,60.08,60.28,60.63,60.93,61.78,62.58,63.42,63.93,64.13,65.12,66.02,66.55,67.23,67.85,71.32,73.03,73.68,75.1,75.52,79.75,79.93,80.88,81.23,91.85,103.73,109.27,115.65,123.63,131.88,173.28,30.68,31.68,32.1,33.12,33.25,33.48,34.5,35.15,35.52,35.87,36.52,36.63,36.75,39.32,40.37,41.4,41.57,41.57,41.78,43.85,44.05,44.97,45.47,45.93,47.03,47.1,47.35,47.63,48.3,48.38,48.62,48.85,48.88,49,49.87,50.65,50.73,51.07,51.42,51.8,52.08,52.72,52.83,53.23,53.32,53.57,53.68,54,54.12,54.23,54.28,54.5,55.13,55.32,55.65,55.85,56.02,56.15,56.43,56.45,56.58,56.58,57.22,57.77,57.78,58.17,58.43,58.57,58 > .58,59.38,59.8,59.87,60.18,60.67,61.35,61.42,61.6,62.07,62.18,62.42,63.07,63.93,64.57,65.2,65.33,65.78,65.8,66.38,66.55,66.78,67.07,67.43,67.57,68.28,68.3,68.32,68.43,68.68,68.82,69,69.55,70.22,70.58,71.42,71.6,72.13,72.17,72.55,72.57,72.67,72.7,72.8,72.8,73.38,73.38,73.67,74.28,75.05,75.4,75.53,75.75,76.6,77.75,78.22,78.42,78.47,78.63,79.27,79.58,79.78,80.18,81.53,81.78,82.17,84.87,85.07,85.2,85.77,86.37,87.25,88.78,88.93,88.98,89.68,91,92.62,93.5,93.83,95.12,97.1,97.23,97.77,98.83,101.47,107.87,108.17,109.07,112.25,112.58,113.53,113.75,119.48,121.32,132.88,139.73,31.95,32.4,35.2,35.35,36.62,36.88,37.62,38.65,39.68,40.05,41.25,41.28,42,42.03,42.23,42.63,42.67,43.02,43.32,43.32,44.23,44.27,44.7,45.08,46.42,48.17,48.25,49.55,50.13,50.85,50.97,51.85,52.28,52.85,52.9,53.13,54.23,54.98,55.45,56.28,56.53,57.33,57.75), > dia_doa=rep(c(1:3), c(92,165,43))) > > table(dat$dia_doa) > > ### (r)andom (Pos)itions (dentro de cada dia) > {set.seed(765); rPos <- unlist(sapply(unique(dat$dia_doa), function(x) > sample(which(dat$dia_doa==x))))} > > ### set.seed() é pra controle do CMR, considere removê-lo > > rDat <- dat[rPos,]; head(rDat) > rDat$dia_doa ### observe que não muda ordem dos dias! > > dN <- unique(table(dat$dia_doa)) ### dados em cada dia! > grp <- NULL; for (n in dN) grp <- c(grp, rep(1:ceiling(n/5), each=5, > len=n)) ### grupos dentro de dias > > DAT <- cbind(rDat, grp); row.names(DAT) <- NULL; head(DAT) > DAT$id <- DAT$dia_doa*1000+DAT$grp ### cria um id/label único para grupos! > > ### Caso deseje remover dados/grupos incompletos... > grp.i <- names(table(DAT$id))[table(DAT$id)<5];grp.i ### grupos > incompletos > dat.i <- unlist(sapply(grp.i, function(x) which(DAT$id==x))); dat.i ### > dados pertencentes a grupos incompletos > > DAT[dat.i,] ### linhas dos grupos incompletos > DATCompl <- DAT[-dat.i,] ### só grupos completos > table(DATCompl$dia_doa) > table(DATCompl$id) > > ### > > require(qcc) > qcc(qcc.groups(DAT$tempo, DAT$id)) ### todos > qcc(qcc.groups(DATCompl$tempo, DATCompl$id)) ### somente completos > ### </code> > > > > > Éder Comunello <comunello.eder(a)gmail.com> > > Dourados, MS - [22 16.5'S, 54 49'W] > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/451376a0/atta… > > > > ------------------------------ > > Message: 3 > Date: Fri, 21 Mar 2014 10:26:15 -0300 > From: "jborssoi ." <joelmirb(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Iteração com optim > Message-ID: > < > CAEFzrQsPbZ_DxKqC5w3iSYP7851UZaqqCdTUUikwBTnVeSShPA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Caros, > tenho 100 amostras e preciso fazer a estimação de parâmetros usando > a função optim. Como faço para que, a cada nova estimação, o vetor com o > chute inicial ("par") seja atualizado com as estimativas obtidas no passo > anterior? > > Grato pela ajuda > Joelmir > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/85bdad8e/atta… > > > > ------------------------------ > > Message: 4 > Date: Fri, 21 Mar 2014 09:55:59 -0400 > From: Diogo Ferrari <diogoferrari(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Iteração com optim > Message-ID: > < > CAEvjS8Lwo0C5DuQUisvNfATUJ6Zy3Qcao_pa_UxM5jA5G_2eDg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Nao sei se compreendi bem sua questao, mas acredito q optim faz exatamente > isso automaticamente. > On Mar 21, 2014 9:26 AM, "jborssoi ." <joelmirb(a)gmail.com> wrote: > > > Caros, > > tenho 100 amostras e preciso fazer a estimação de parâmetros usando > > a função optim. Como faço para que, a cada nova estimação, o vetor com o > > chute inicial ("par") seja atualizado com as estimativas obtidas no passo > > anterior? > > > > Grato pela ajuda > > Joelmir > > > > _______________________________________________ > > 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/20140321/7fdd0eb1/atta… > > > > ------------------------------ > > Message: 5 > Date: Fri, 21 Mar 2014 10:09:24 -0400 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, Edson Lira > <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] Grárico de controle > Message-ID: > < > CABmC8gkK_HxJ7wGUqsUhounM84TPkQgzd6ZHcA7O0VyPWPzqAA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Edson, bom dia! > > Um ponto por dia simplifica bastante o procedimento e você pode dispensar > praticamente todo o código anteriormente enviado. > > Tendo criado (ou restaurado) o objeto de dados 'dat', utilize: > > ### qcc.groups(data, sample) > qcc.groups(dat$tempo, dat$dia_doa) > qcc(qcc.groups(dat$tempo, dat$dia_doa)) > > Reforçando que serão amostras de tamanhos distintos. > > Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> > Dourados, MS - [22 16.5'S, 54 49'W] > > > Em 21 de março de 2014 08:47, Edson Lira <edinhoestat(a)yahoo.com.br > >escreveu: > > > Eder, muito obrigado, ajudou muito, vi que tenho no mesmo dia vários > > grupos. > > > > Como faço para ter um único resultado por dia? Ou seja, para cada dia um > > único ponto no gráfico. > > > > > > [ ]'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/20140321/8124fe7c/atta… > > > > ------------------------------ > > Message: 6 > Date: Fri, 21 Mar 2014 11:15:50 -0300 > From: "walmes ." <walmeszeviani(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Iteração com optim > Message-ID: > <CAFU= > EkZMw34AABsMzb9r8Ey5KkbZi6rg1UFgpmBcc0u03wfzUA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Faça com um laço for por exemplo > > chute <- ... > for(...){ > op <- optim(chute, ...) > chute <- op$par > } > > À disposição. > Walmes. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/d078d6db/atta… > > > > ------------------------------ > > Message: 7 > Date: Fri, 21 Mar 2014 15:16:27 +0100 > From: Wagner Bonat <wbonat(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Iteração com optim > Message-ID: > <CANt= > 4MhG-eB4uOTxSCzmJrdm1aEbpmXa3LF6ax3c_OsaDAfgdw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > O que vc quer dizer com 100 amostras ?? Vc tem 100 conjuntos de dados > diferentes e quer estimar para cada um determinado modelo ... Nao consegui > entender ... Que metodo vc quer usar para inferencia ?? > > > Em 21 de março de 2014 14:55, Diogo Ferrari <diogoferrari(a)gmail.com > >escreveu: > > > Nao sei se compreendi bem sua questao, mas acredito q optim faz > exatamente > > isso automaticamente. > > On Mar 21, 2014 9:26 AM, "jborssoi ." <joelmirb(a)gmail.com> wrote: > > > >> Caros, > >> tenho 100 amostras e preciso fazer a estimação de parâmetros usando > >> a função optim. Como faço para que, a cada nova estimação, o vetor com o > >> chute inicial ("par") seja atualizado com as estimativas obtidas no > passo > >> anterior? > >> > >> Grato pela ajuda > >> Joelmir > >> > >> _______________________________________________ > >> 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. > > > > > > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/a5d6012f/atta… > > > > ------------------------------ > > Message: 8 > Date: Fri, 21 Mar 2014 11:23:04 -0300 > From: "jborssoi ." <joelmirb(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Iteração com optim > Message-ID: > < > CAEFzrQuNS+oj9FmamJJkmUtp9zp8rYhHXknMdawxnxcmRE151A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá Wagner, > é isso... tenho 100 conjuntos de dados. Eu preciso maximizar o vetor de > parâmetros para a primeira amostra. Já para a segunda (e assim por diante) > o optim deve utilizar o vetor obtido para a primeira amostra (e assim por > diante). > Imagino algo do tipo: > for(i in 1:J){ > maxL.rep <- optim(par=???, fn, gr, method="BFGS", > control=list(fnscale=-1,maxit=1000),y=W1$y[,i], u=W1$u[,i], X=d[,c(1,2)], > Z=Zi, sigma.e=1) > } > > mas não sei como fazer para que o optim atualize "par". > > abraço > > > Em 21 de março de 2014 11:16, Wagner Bonat <wbonat(a)gmail.com> escreveu: > > > O que vc quer dizer com 100 amostras ?? Vc tem 100 conjuntos de dados > > diferentes e quer estimar para cada um determinado modelo ... Nao > consegui > > entender ... Que metodo vc quer usar para inferencia ?? > > > > > > Em 21 de março de 2014 14:55, Diogo Ferrari <diogoferrari(a)gmail.com > >escreveu: > > > > Nao sei se compreendi bem sua questao, mas acredito q optim faz > exatamente > >> isso automaticamente. > >> On Mar 21, 2014 9:26 AM, "jborssoi ." <joelmirb(a)gmail.com> wrote: > >> > >>> Caros, > >>> tenho 100 amostras e preciso fazer a estimação de parâmetros usando > >>> a função optim. Como faço para que, a cada nova estimação, o vetor com > o > >>> chute inicial ("par") seja atualizado com as estimativas obtidas no > passo > >>> anterior? > >>> > >>> Grato pela ajuda > >>> Joelmir > >>> > >>> _______________________________________________ > >>> 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. > >> > > > > > > > > -- > > Wagner Hugo Bonat > > LEG - Laboratório de Estatística e Geoinformação > > UFPR - Universidade Federal do Paraná > > > > _______________________________________________ > > 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/20140321/94ab16c2/atta… > > > > ------------------------------ > > Message: 9 > Date: Fri, 21 Mar 2014 15:36:43 +0100 > From: Wagner Bonat <wbonat(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Iteração com optim > Message-ID: > <CANt= > 4Mh5gG2QAodJvBDv+3R9vZs3_wf+qRnGZ9Bo0iV-NjjYhw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Talvez isso ajude ! > > ll <- function(mu,x){-sum(dnorm(x,mean=mu, sd=1,log=TRUE))} > ll(mu = 0, x = rnorm(10)) > dados <- matrix(rnorm(1000), ncol=100, nrow=10) > saida <- list() > saida[1][[1]] <- optim(par= 0, f = ll, x = dados[,1], method="Brent", upper > = 10, lower = -10) > saida[1][[1]]$par > for(i in 2:100){ > print(saida[c(i-1)][[1]]$par) > saida[i][[1]] <- optim(par= saida[c(i-1)][[1]]$par, f = ll, x = > dados[,i], method="Brent", upper = 10, lower = -10) > } > > > > Em 21 de março de 2014 15:23, jborssoi . <joelmirb(a)gmail.com> escreveu: > > > Olá Wagner, > > é isso... tenho 100 conjuntos de dados. Eu preciso maximizar o vetor de > > parâmetros para a primeira amostra. Já para a segunda (e assim por > diante) > > o optim deve utilizar o vetor obtido para a primeira amostra (e assim por > > diante). > > Imagino algo do tipo: > > for(i in 1:J){ > > maxL.rep <- optim(par=???, fn, gr, method="BFGS", > > control=list(fnscale=-1,maxit=1000),y=W1$y[,i], u=W1$u[,i], X=d[,c(1,2)], > > Z=Zi, sigma.e=1) > > } > > > > mas não sei como fazer para que o optim atualize "par". > > > > abraço > > > > > > Em 21 de março de 2014 11:16, Wagner Bonat <wbonat(a)gmail.com> escreveu: > > > > O que vc quer dizer com 100 amostras ?? Vc tem 100 conjuntos de dados > >> diferentes e quer estimar para cada um determinado modelo ... Nao > consegui > >> entender ... Que metodo vc quer usar para inferencia ?? > >> > >> > >> Em 21 de março de 2014 14:55, Diogo Ferrari <diogoferrari(a)gmail.com > >escreveu: > >> > >> Nao sei se compreendi bem sua questao, mas acredito q optim faz > >>> exatamente isso automaticamente. > >>> On Mar 21, 2014 9:26 AM, "jborssoi ." <joelmirb(a)gmail.com> wrote: > >>> > >>>> Caros, > >>>> tenho 100 amostras e preciso fazer a estimação de parâmetros usando > >>>> a função optim. Como faço para que, a cada nova estimação, o vetor > com o > >>>> chute inicial ("par") seja atualizado com as estimativas obtidas no > passo > >>>> anterior? > >>>> > >>>> Grato pela ajuda > >>>> Joelmir > >>>> > >>>> _______________________________________________ > >>>> 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. > >>> > >> > >> > >> > >> -- > >> Wagner Hugo Bonat > >> LEG - Laboratório de Estatística e Geoinformação > >> UFPR - Universidade Federal do Paraná > >> > >> _______________________________________________ > >> 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. > > > > > > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/e0e422aa/atta… > > > > ------------------------------ > > Message: 10 > Date: Fri, 21 Mar 2014 07:57:30 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "ecomunel(a)gmail.com" <ecomunel(a)gmail.com>, > "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Grárico de controle > Message-ID: > <1395413850.17732.YahooMailNeo(a)web160101.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Era isso mesmo que precisava. > > Agora tenho as seguintes perguntas: > > Como faço para fixar os limites em 2 sigmas, usando o comando que vi no > help não funcionou, não sei por que. > Quando uso size ele fixa, mas não me permite mexer no nsigma. tem algo que > não estou sabendo fazer. > > Veja se você consegue alguma sugestão. > > Muito obrigado pela sua atenção Eder. > > > [ ]'s. > > Edson Lira > Estatístico > Manaus-Amazonas > > > > Em Sexta-feira, 21 de Março de 2014 10:09, Éder Comunello < > comunello.eder(a)gmail.com> escreveu: > > Edson, bom dia! > > Um ponto por dia simplifica bastante o procedimento e você pode dispensar > praticamente todo o código anteriormente enviado. > > Tendo criado (ou restaurado) o objeto de dados 'dat', utilize: > > ### qcc.groups(data, sample) > qcc.groups(dat$tempo, dat$dia_doa) > qcc(qcc.groups(dat$tempo, dat$dia_doa)) > > Reforçando que serão amostras de tamanhos distintos. > > > Éder Comunello <comunello.eder(a)gmail.com> > Dourados, MS - [22 16.5'S, 54 49'W] > > > > Em 21 de março de 2014 08:47, Edson Lira <edinhoestat(a)yahoo.com.br> > escreveu: > > Eder, muito obrigado, ajudou muito, vi que tenho no mesmo dia vários > grupos. > > > >Como faço para ter um único resultado por dia? Ou seja, para cada dia um > único ponto no gráfico. > > > > > > > > > >[ ]'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/20140321/3cd59f3c/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 39, assunto 21 > ***************************************** >

1 0

(sem assunto)
by Luís Fernando Pirondi Junior 21 Mar '14

21 Mar '14

Olá, escrevi o seguinte script: > A<-read.table("A.txt", header=TRUE) > A Dias Massa 1 38 370 2 121 4200 3 208 13200 4 291 17100 5 381 18000 > plot(A$Dias,A$Massa,xlab="Dias",ylab="Massa", col="blue") > modelo<-glm(Massa~Dias, data=A, family=binomial) Mensagens de aviso perdidas: In eval(expr, envir, enclos) : #sucessos não-inteiro em um glm binomial! > summary(modelo) Call: glm(formula = Massa ~ Dias, family = binomial, data = A) Deviance Residuals: 1 2 3 4 5 -0.19282 -0.03793 0.15578 0.09105 -0.11083 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -3.737847 4.166859 -0.897 0.37 Dias 0.006566 0.014028 0.468 0.64 (Dispersion parameter for binomial family taken to be 1) Null deviance: 0.339765 on 4 degrees of freedom Residual deviance: 0.083459 on 3 degrees of freedom AIC: 5.1496 Number of Fisher Scoring iterations: 6 E obtive o gráfico em anexo. Minha intenção era fazer uma regressão logística com os valores da Tabela "A", no entanto a curva não parece muito ajustada aos parâmetros. O que pode estar errado ? -- *Luís Fernando Pirondi Junior*

2 1

Iteração com optim
by jborssoi . 21 Mar '14

21 Mar '14

Caros, tenho 100 amostras e preciso fazer a estimação de parâmetros usando a função optim. Como faço para que, a cada nova estimação, o vetor com o chute inicial ("par") seja atualizado com as estimativas obtidas no passo anterior? Grato pela ajuda Joelmir

4 5

Calendário do pacote googleVis
by bruno vidigal 21 Mar '14

21 Mar '14

Olá, alguém tentou rodar o calendário do pacote googleVis? Estou utilizando o mesmo exemplo disponibilizado pelo autor mas não aparecem os valores a partir do dia 20/out para todos os 3 anos de análise. Alguém teria alguma solução? install.packages(c("devtools","RJSONIO", "knitr", "shiny", "httpuv")) library(devtools) install_github("mages/googleVis") cl1 <- gvisCalendar(Cairo, datevar="Date", numvar="Temp", options=list(width=1000)) plot(cl1) Talvez seria um problema de locale... -- Bruno C. Vidigal

5 10

Organizando os dados (limpando sinais desnecessários)
by Tropidurus Torquatus 20 Mar '14

20 Mar '14

Tenho a matriz a baixo: > str(bmatrix1) 'data.frame': 2189 obs. of 16 variables: $ numMA.ACC : num 15 15 15 15 15 14 14 8 8 8 ... $ numMA.DON : num 11 11 11 11 11 10 10 5 5 5 ... $ numMA.Mol_ID : num 2 2 2 2 2 4 4 31 31 31 ... $ numMA.ROT : num 5 5 5 5 5 4 4 6 6 6 ... $ charMA.energy : Factor w/ 2189 levels "10.745595054136636 __ ",..: 825 840 855 928 957 768 807 54 66 165 ... ..- attr(*, "names")= chr "CMP1" "CMP2" "CMP3" "CMP4" ... A coluna "$ charMA.energy” É a minha coluna de interesse, como retirar os “__” e alterá-la para coluna numérica mantendo todos as casas decimais? No mais agradeço. Thales

2 2

Teste proporção
by Andre Oliveira 19 Mar '14

19 Mar '14

Boa noite, qual função faz o teste para proporções e calcula IC além da função binom.test(12,29, p = 0.5, conf.level = 0.95)? André Oliveira Souza

4 5

Re: [R-br] IC Variancia
by Robert Iquiapaza 17 Mar '14

17 Mar '14

#http://pages.stat.wisc.edu/~yandell/st571/R/append7.pdf var.interval = function(data, conf.level = 0.95) { df = length(data) - 1 chilower = qchisq((1 - conf.level)/2, df) chiupper = qchisq((1 - conf.level)/2, df, lower.tail = FALSE) v = var(data) c(l.inf=df * v/chiupper, l.sup=df * v/chilower) } plot(x=rnorm(48,34,7)) hist(x) var(x) var.interval(x) Em 15 de março de 2014 09:08, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: > Pessoal qual função do R calcula IC para variancia? > obrigado > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código > mínimo reproduzível.

2 1

combine_factor
by Dian Andrade 17 Mar '14

17 Mar '14

Como posso combinar NA com outra categoria a <- c("4", "9","1","1", "4",NA,"9",NA) table(a, useNA = "ifany") 1 4 9 <NA> 2 2 2 2 Neste exemplo como incluir o valor do NA em ? Obrigado

3 3

Re: [R-br] Digest R-br, volume 39, assunto 17
by Gilenio Borges Fernandes 17 Mar '14

17 Mar '14

Pequeno detalhe na cálculo do coeficiente de variacao! 100*sqrt(qms)/m Gilenio Fernandes Em 17 de março de 2014 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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. Coeficiente de determinação (cv): experimentos com mais de > uma estimativa do erro (Jose Claudio Faria) > 2. Re: Coeficiente de determinação (cv): experimentos com mais > de uma estimativa do erro (walmes .) > 3. Re: Coeficiente de determinação (cv): experimentos com mais > de uma estimativa do erro (Jose Claudio Faria) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 17 Mar 2014 07:41:15 -0300 > From: Jose Claudio Faria <joseclaudio.faria(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Coeficiente de determinação (cv): experimentos com > mais de uma estimativa do erro > Message-ID: > < > CAN+Emd8zJtNHF1CKM_6oQcmVV2UX1Mmn0iEuXdQjzgZSCcOF9w(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Pessoal, > > Alguém teria um bom algoritmo em R para determinar a estimativa do > coeficiente de variação (cv - medida estatística básica e clássica > sobre a precisão) dos experimentos com mais de uma estimativa do erro > (por exemplo: parcelas subdivididas, subsubdividads, etc). > > Tenho um feito apenas para experimentos com apenas uma estimativa do erro: > > #=============================================================================== > # Name : cv > # Author : José Cláudio Faria/UESC/DCET > # Date : 26/7/2006 15:25:07 > # Version : v3 > # Aim : Calculate coefficient variation from aov objects > > #=============================================================================== > > # Arguments: > # av aov object > > cv <- function(av) > { > if(is.null(av)) > stop('Please, check the parameter!') > qmee <- with(av, sum(residuals^2) / df.residual) > cv <- 100 * sqrt(qmee) / mean(av$fitted.values) > return(round(cv, 2)) > } > > # Exemplo de uso: > > av <- aov(Sepal.Length ~ Species, data=iris) > > > cv(av) > [1] 8.81 > > Ab, > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > Jose Claudio Faria > Estatistica > UESC/DCET/Brasil > joseclaudio.faria at gmail.com > Telefones: > 55(73)3680.5545 - UESC > 55(73)9966.9100 - VIVO > 55(73)9100.7351 - TIM > 55(73)8817.6159 - OI > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > > > ------------------------------ > > Message: 2 > Date: Mon, 17 Mar 2014 10:29:17 -0300 > From: "walmes ." <walmeszeviani(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Coeficiente de determinação (cv): experimentos com > mais de uma estimativa do erro > Message-ID: > <CAFU=Eka7XYyAJZbGpXc7dP2aP7NKzXS=uyPeqcgyRj-Hg= > i9KQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Prezado Faria, > > Embora eu seja contrário à interpretação ou uso, ou melhor, abuso do CV por > certas pessoas e revistas, eu antigamente calculava assim, só para que > tenha uma idéia. Sempre em cursos eu argumento contra o CV e meus códigos > raramente o fazem aparecer. Isso porque muitas pessoas acham que um CV alto > invalidam o experimento ou a análise. Argumento que não. Enfim, segue minha > contribuição. > > > ##----------------------------------------------------------------------------- > > da <- expand.grid(a=gl(4,1), b=gl(3,1), c=gl(5,1)) > da$y <- rnorm(da$a, mean=3) > > m0 <- aov(y~a+b*c+Error(a:b), da) > summary(m0) > > a0 <- summary(m0) > str(a0) > > qms <- sapply(a0, > function(i){ > i[[1]]["Residuals",3] > }) > m <- mean(da$y) > > 100*qms/m > > > ##----------------------------------------------------------------------------- > > À disposição. > Walmes. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140317/0c3f06c6/atta… > > > > ------------------------------ > > Message: 3 > Date: Mon, 17 Mar 2014 11:39:37 -0300 > From: Jose Claudio Faria <joseclaudio.faria(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Coeficiente de determinação (cv): experimentos com > mais de uma estimativa do erro > Message-ID: > <CAN+Emd-xv=_jDu0GQKY3x8_gPx8ngw4zfdstE_2T_cEi0= > imYw(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Walmes, > > Concordo com seus argumentos, mas considero o cv uma medida útil > (comparativamente) da qualidade de um experimento frente a grupos > similares. Uso no ensino de graduação em experimentação. > > Obrigado pelo código! > (Havia encontrado alguns mas não havia gostado do código de nenhum do que > vi). > > Ab, > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > Jose Claudio Faria > Estatistica > UESC/DCET/Brasil > joseclaudio.faria at gmail.com > Telefones: > 55(73)3680.5545 - UESC > 55(73)9966.9100 - VIVO > 55(73)9100.7351 - TIM > 55(73)8817.6159 - OI > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > > > 2014-03-17 10:29 GMT-03:00 walmes . <walmeszeviani(a)gmail.com>: > > Prezado Faria, > > > > Embora eu seja contrário à interpretação ou uso, ou melhor, abuso do CV > por > > certas pessoas e revistas, eu antigamente calculava assim, só para que > tenha > > uma idéia. Sempre em cursos eu argumento contra o CV e meus códigos > > raramente o fazem aparecer. Isso porque muitas pessoas acham que um CV > alto > > invalidam o experimento ou a análise. Argumento que não. Enfim, segue > minha > > contribuição. > > > > > ##----------------------------------------------------------------------------- > > > > da <- expand.grid(a=gl(4,1), b=gl(3,1), c=gl(5,1)) > > da$y <- rnorm(da$a, mean=3) > > > > m0 <- aov(y~a+b*c+Error(a:b), da) > > summary(m0) > > > > a0 <- summary(m0) > > str(a0) > > > > qms <- sapply(a0, > > function(i){ > > i[[1]]["Residuals",3] > > }) > > m <- mean(da$y) > > > > 100*qms/m > > > > > ##----------------------------------------------------------------------------- > > > > À disposição. > > Walmes. > > > > _______________________________________________ > > 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 39, assunto 17 > ***************************************** > -- Gilenio Borges Fernandes Professor Associado IV 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

1 0

Significado de graus de liberdade
by Luiz Roberto Martins Pinto 17 Mar '14

17 Mar '14

*Pessoal,* *Perdoem-me apresentar uma pergunta que nada tem a ver com o R, mas não participo de outra lista com tantos entusiastas por análise de dados e estatística.* *Alguém tem referência (livro ou artigo) que trate do 'Significado dos graus de liberdade'?* *Estou especialmente interessado em entender os graus de liberdade em uma interação entre dois fatores A e B, isto é, para AB, por exemplo na ANOVA. * *Sei que os graus de liberdade para esta interação é igual a (a-1)*(b-1), em que 'a' é o número de níveis de A e 'b' é o número de níveis de B. Mas gostaria de encontrar os resultados entendendo o numero total de observações e o número de parâmetros que estou realmente calculando.* *Utilizando o argumento de Chong Ho Yu "df is a function of both the number of* *observations and the number of variables in one's model".* *Então estou desejoso de entender "both the number of* *observations and the number of variables" para a interação AB.* *Se não for do interesse do pessoal da lista, peço àqueles que puderem ajudar que me enviem enviem e-mail off-lista para luizroberto.uesc(a)gmail.com <luizroberto.uesc(a)gmail.com>.* *Agradeço,* *Luiz Roberto* Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

3 4

Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro
by Jose Claudio Faria 17 Mar '14

17 Mar '14

Pessoal, Alguém teria um bom algoritmo em R para determinar a estimativa do coeficiente de variação (cv - medida estatística básica e clássica sobre a precisão) dos experimentos com mais de uma estimativa do erro (por exemplo: parcelas subdivididas, subsubdividads, etc). Tenho um feito apenas para experimentos com apenas uma estimativa do erro: #=============================================================================== # Name : cv # Author : José Cláudio Faria/UESC/DCET # Date : 26/7/2006 15:25:07 # Version : v3 # Aim : Calculate coefficient variation from aov objects #=============================================================================== # Arguments: # av aov object cv <- function(av) { if(is.null(av)) stop('Please, check the parameter!') qmee <- with(av, sum(residuals^2) / df.residual) cv <- 100 * sqrt(qmee) / mean(av$fitted.values) return(round(cv, 2)) } # Exemplo de uso: > av <- aov(Sepal.Length ~ Species, data=iris) > cv(av) [1] 8.81 Ab, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9966.9100 - VIVO 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

IC Variancia
by Andre Oliveira 15 Mar '14

15 Mar '14

Pessoal qual função do R calcula IC para variancia? obrigado

1 0

Quebrar eixo/escala em gráfico de barras agrupadas no R
by Rodolfo Souza 14 Mar '14

14 Mar '14

Olá pessoal, Alguém sabe como quebrar o eixo/escala em gráfico de barras agrupadas no R? Eu preciso fazer um gráfico e tem um grupo bem maior que os demais, e nesse caso passa ter um problema de escala. Segue o CMR cp1=c(6.81, 3.78, 3.16, 0.66, 4.27, 0.75, 5.64, 2.02, 0.99, 31.00) cp2=c(4.57, 0, 0, 0, 3.81, 0, 6.6, 1.27, 0, 140) cp3=c(0.46, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 35.00) grups=c(1,2,3,4,5,6,7,8,9,10) dbh=data.frame(grups,cp1,cp2,cp3) dbhm=t(as.matrix(dbh[,-1])) colnames(dbhm)=dbh$grups barplot(dbhm, beside=T, las=1) O último grupo tem valores bem maior que os demais, por isso que fica ruim ver as variações de cada componente nos outros grupos. Eu queria fazer semelhante a esse ( http://uspace.shef.ac.uk/servlet/JiveServlet/showImage/2-12478-3980/gap_plo…) mas não sei como quando utiliza a função "beside=T". Grato pela atenção. Rodolfo Marcondes Silva Souza

3 2

Manipulação de data frame
by Rodrigo Muniz 14 Mar '14

14 Mar '14

Caros membros da lista, Desculpe por não enviar nenhum exemplo ou lista de comando. Tenho um dataframe com 1500 colunas e 700 linhas, eu preciso fazer uma razão entre todas as colunas, ou seja, C1/C2, C1/C3, C1/C4....C1499/C1500. Fazer essa divisão de forma manual apenas por indexação dos vetores é simples, porém, isso daria muito trabalho e acredito que deve haver um jeito mais fácil de fazer isso. Porém não tenho a menor ideia de como começar a fazer-lo. Por isso peço a ajuda de vocês, se alguém conhece alguma forma, algum exemplo, ou mesmo algum material que eu possa consultar, para que eu tenha esse "input" inicial. Desde já agradeço. Att... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 6

qcc- carta de controle.
by Edson Lira 13 Mar '14

13 Mar '14

Caros amigos, estou tentando fazer um gráfico xbarra no pacote qcc e estou tendo problemas. Tenho 3000 amostras de um tempo entre 20 e 180 minutos. Criei a variável amostra com a seguinte rotina: amostra<-rep(1:5, each = 5, times = 120) require(qcc) usei a rotina tem <- qcc.groups(tempo,amostra) Essa rotina está me dando problema, quando ploto o gráfico está me mostrando somente 5 grupos, quando deveria ter 600 grupos (amostras de tamanho 5). Tem um exemplo no pacote, mas não consegui ver aonde está o meu erro. ?qcc [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 3

Advanced RNA-Seq and ChIP-Seq Data Course (Campinas/Brazil/2014)
by Benilton Carvalho 12 Mar '14

12 Mar '14

Peço desculpas pelo off-topic e também a aqueles que possam receber esta mensagem duplicada. É com satisfação que anuncio a realização do "Advanced RNA-Seq and ChIP-Seq Data Course (Campinas/Brazil/2014)" na Faculdade de Ciências Médicas da Universidade Estadual de Campinas (UNICAMP) no período de 7 a 10 de abril de 2014. O curso será ministrado em Inglês e contará com a participação de instrutores experientes, responsáveis pelo curso (homônimo) no Instituto Europeu de Bioinformática, Hixton UK (EMBL-EBI-UK). Benilton Carvalho , State University of Campinas - UNICAMP Gabriella Rustici , EMBL-EBI, UK Myrto Kostadima , EMBL-EBI, UK Nicolas Delhomme , EMBL, Heidelberg Remco Loos , EMBL-EBI, UK O conteúdo programado para o curso inclui: - Introdução ao ambiente R/Bioconductor - Controle de qualidade de dados de sequenciamento - Mapeamento - Modelos estatísticos para dados de sequenciamento - Estimação de expressão gênica e uso de exons - RNA-Seq - ChIP-Seq - Determinação e anotação de picos - Análise de motifs Os interessados deverao efetuar suas candidaturas no sítio eletrônico: http://tiny.cc/rnaseqcampinas As candidaturas serão analisadas e, no máximo, 32 participantes selecionados. Estes selecionados serão contactados até o dia 31 de março e receberão as instruções para o pagamento da inscrição (R$ 300,00). Os participantes responsabilizar-se-ão por suas viagens, acomodação e alimentação. Solicito a gentileza de repassarem a informação a quem possa ter interesse. Atenciosamente, Benilton Carvalho

1 0

gráfico de dispersão em geoestatística
by Hélio Gallo Rocha 12 Mar '14

12 Mar '14

Listeiros, bom dia. Em um objeto da classe geodata (pacote geoR) , ao usar o comando plot é produzido um gráfico quadruplo, um com a posição dos pontos, dois mostrando o comportamento da variável com relação aos sentidos NS e SW e um Histograma. Como poderia separar os gráficos de comportamento? Outra dúvida. Teria como editar as fontes dos eixos x e y? aumentar a fonte ou negrito. Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

3 6

fatorial_duplo_regressao
by Odirley Campos 12 Mar '14

12 Mar '14

Srs(as), boa noite. Sou iniciante no R. Como faço para ajustar um modelo de 2º grau para os dados abaixo. Antecipademante agradeco a<-read.table("dose_fcr_ph2.txt", h=T) names(a) tra<-as.factor(a$trat) #usado apenas para verificar pressupostos da anova fcr<-as.factor(a$fcr) #4 niveis (1500, 2000, 2500 e 3000) dos<-as.factor(a$dose) #3 niveis (0, 1, 2), tratado como qualitativo ph<-as.numeric(a$ph) b<-data.frame(tra,fcr,dos,ph) b > fcr [1] 1500 1500 1500 2000 2000 2000 2500 2500 2500 3000 3000 3000 1500 1500 1500 [16] 2000 2000 2000 2500 2500 2500 3000 3000 3000 1500 1500 1500 2000 2000 2000 [31] 2500 2500 2500 3000 3000 3000 > dos [1] 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 > ph [1] 5.86 6.19 6.87 5.87 6.30 6.30 6.05 6.37 7.03 5.69 5.72 6.52 6.05 6.02 6.09 [16] 6.19 6.73 6.80 7.08 7.08 7.08 9.25 8.81 9.06 9.25 9.24 8.81 9.13 8.97 9.36 [31] 9.05 8.99 8.93 9.21 8.90 7.42 ane2<-aov(ph~dos*fcr, data=b) summary(ane2)#tenho interesse no desdobramento # Df Sum Sq Mean Sq F value Pr(>F) #dos 2 45.25 22.625 137.133 7.37e-14 *** #fcr 3 2.44 0.812 4.921 0.00836 ** #dos:fcr 6 14.08 2.347 14.223 7.24e-07 *** #Residuals 24 3.96 0.165 ane3<-aov(ph~dos/fcr, data=b)#guarda efeitos aninhados summary(ane3, split=list("dos:fcr"=list("dos0"=c(1,4,7), "dos1"=c(2,5,8), "dos2"=c(3,6,9)))) # Df Sum Sq Meam Sq F value Pr(>F) #dos 2 45.10 22.550 43.910 9.56e-09 *** #dos:fcr 9 8.69 0.966 1.881 0.10465 #dos:fcr: dos0 3 0.42 0.140 0.272 0.84507 #dos:fcr: dos1 3 7.50 2.499 4.866 0.00877 ** #dos:fcr: dos2 3 0.78 0.259 0.503 0.68359 #Residuals 24 12.33 0.514 #Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #--------------------------------------------------------- #Ajustando modelo de 2 grau > m1<-aov(ph~-1+ dos/(fcr+I(fcr^2)), data=b)# aparece a mesagem de erro abaixo #Erro em `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : #contrastes podem ser aplicados apenas a fatores com 2 ou mais níveis #Além disso: Mensagens de aviso perdidas: #In Ops.factor(fcr, 2) : ^ not meaningful for factors # está faltando alguma informação no "comando"? Odirley R. Campos Engenheiro Agrônomo UFV/MG

3 2

Teste t Student
by Luiz Roberto Martins Pinto 12 Mar '14

12 Mar '14

Olá, Estou precisando utilizar a função t.test da biblioteca {stats} para comparar respostas de tratamentos em contrastes dois-a-dois, em um experimento com 20 tratamentos. Este teste me permite comparar amostras com variancias diferentes (var.equal = FALSE), bem como utilizar o método de Satterthwaite (não preciso neste caso específico, mas é interessante em outros). As demais opções podem ser o 'default'. Preciso gerar uma matriz 20 x 20 com o p-valor referente a cada um dos pares de tratamentos, em que na diagonal principal contenha 1.00. Exemplo de matriz com p-valor T1 T2 T3 ... T 20 T1 1.00 0.85 0.03 ... 0.25 T2 0.85 1.00 0.42 ... 0.31 T3 0.03 0.42 1.00 ... 0.44 ... .... .... .... ... ... T20 0.25 0.31 0.44 ... 1.00 Ocorre que the formula interface is only applicable for the 2-sample tests, e preparar 90 arquivos com estes contrastes dois-a-dois para gerar esta matriz 20 x 20 não me parece ser a alternativa mais adequada. Mas, também não tenho conhecimento (análise e programação) suficiente de R para preparar um script para fazer isto. Alguém pode ajudar? Abaixo, endereço com os dados. http://www.datafilehost.com/d/1ec1ec2b Luiz Roberto Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

4 8

Análise de uma expressão como string
by Waldelene Moura 12 Mar '14

12 Mar '14

Sou iniciante do R. Estou precisando de uma sugestão para fazer o seguinte: Tenho numa string com uma expressão e desejo obter o resultado desta expressão. Ex: "2 * 3" me apresenta como resultado 5. "42 + 7" me apresenta como resultado 49. Existe alguma função para isto? Agradeço a quem puder me ajudar. Waldelene Moura Analista de Sistemas Universidade Federal do Maranhão - PROEN

2 2

Anava com modelo autoregressivo
by Kalill V.Páscoa 12 Mar '14

12 Mar '14

Bom dia a todos Estou analisando um experimento em DBC "split plot in the time" São 16 medições (0,68 a 10,25 anos). O experimento apresenta: - 4 blocos; - 4 tratamentos referentes aos espaçamentos de plantio; - e 2 subtratamentos dentro de cada tratamento referentes a aplicação ou não da desrama; cada subparcela tem 112 plantas. - tenho as coordenadas de cada planta; Gostaria de aplicar a ANAVA com modelo autorregressivo, contudo estou com dificuldades para construir a matriz de proximidade espacial. Estou usando a função nb2mat do pacote spdep Minha dificuldade esta no argumento neighbours. Alguém já trabalhou com essa metodologia e poderia me auxiliar? Att -- *Kalill José Viana da Páscoa* *Engenheiro Florestal do Laboratório de Estudos e Projetos em Manejo Florestal - LEMAF * *Mestrando em Engenharia FlorestalUniversidade Federal de Lavras* *Rua Geraldo dos Passos, 81 * *Água Limpa, **Lavras-MG* *CEP: 37200-000* *(35)3821-6833* *(35)9176-9332* *(35)9928-6574*

1 0

comparação de médias em experimento desbalanceado, não fatorial
by Santino Aleandro da Silva 11 Mar '14

11 Mar '14

2 1

Extração de valor-p na função hlme do pacote lcmm
by Felipe Buchbinder 11 Mar '14

11 Mar '14

Prezados, Estou utilizando a função hlme do pacote lcmm para estimar um *Growth Mixture Model.* Um exemplo de sintaxe é : gmm2 <- hlme(Y ~ time + time2, ng = 2, mixture = Y ~ time + time2, subject = 'ID',data = x.reg) Eu queria obter o valor-p do termo quadrático para cada classe, mas não estou conseguindo. (O comando *summary* me mostra esse valor-p, mas eu consigo extraí-lo.) Desde já, agradeço. Felipe Buchbinder

3 3

Loop para carregar bases
by Gerson R. Primo Jr 11 Mar '14

11 Mar '14

Oi pessoal, Recebo vários arquivos e preciso uni-los no R. Atualmente faço da seguinte forma: total = rbind(read.csv("base_1.csv", sep=";"), read.csv("base_2.csv", sep=";"), read.csv("base_3.csv", sep=";"), read.csv("base_4.csv", sep=";"), read.csv("base_5.csv", sep=";") ) Acontece que recebo esses arquivos semanalmente e a quantidade de arquivo varia (25 a 50), Alguém sabe como posso criar um 'for' que declare a quantidade de arquivo o R faz esse loop? *Att,Gerson R. Primo Jr*

5 6

Multi Layer Perceptron para Múltiplas Variáveis Independentes
by Lucas Venezian Povoa 11 Mar '14

11 Mar '14

Olá comunidade R, Eu gostaria de saber se existe algum pacote no R para gerar análises de estimativa com Redes Neurais Muti Layer Perceptron (MLP) para múltiplas variáveis independentes. O que procuro é algo parecido com o código abaixo: > data <- data.frame(a = rnorm(100, 10, 3), b = rnorm(100, 0, 1), y = runif(100)) > > lm(y ~ a + b, data) > > mlp(y ~ a + b, data, hiddenl = 1) > Obrigado pela ajuda. Lucas Venezian Povoa

3 4

Integração Google Maps
by Frederico Steinmetz Alvarez 10 Mar '14

10 Mar '14

Boa noite a todos. Será que algum dos membros da lista saberia me informar se existe algum pacote no R que faça a integração com o google maps? Estou precisando montar um mapa com informações sobre o município, algo tipo camadas. Gostaria de apresentar as informações como pontos no google maps para melhor visualização espacial dos dados. Atenciosamente -- Frederico Steinmetz Alvarez Engenheiro Eletricista - Sistemas de Produção - UFSC Mestre em Eng. Mecânica - Energia Eólica - UFPE (48) 8448-5994 (Oi) (48) 9101-5994 (Vivo)

4 5

(sem assunto)
by Fátima Lima Paula 10 Mar '14

10 Mar '14

Prezados, eu tenho um data frame assim cap=c("A1","A3","B2","B5","C6","C8","D4","F5","G8","H1","K9","A1","A3","B2","B5","C6","C8","D4","F5","G8","H1","K9") Var1=c("nao","nao","nao","nao","nao","nao","nao","nao","nao","nao","nao", "sim","sim","sim","sim","sim","sim","sim","sim","sim","sim","sim") Var2=c(2,5,7,4,1,2,4,8,9,6,7,4,5,3,2,4,5,8,7,0,1,1) df=cbind(cap,Var1,Var2) df=as.data.frame(df) df cap Var1 Var2 1 A1 nao 2 2 A3 nao 5 3 B2 nao 7 4 B5 nao 4 5 C6 nao 1 6 C8 nao 2 7 D4 nao 4 8 F5 nao 8 9 G8 nao 9 10 H1 nao 6 11 K9 nao 7 12 A1 sim 4 13 A3 sim 5 14 B2 sim 3 15 B5 sim 2 16 C6 sim 4 17 C8 sim 5 18 D4 sim 8 19 F5 sim 7 20 G8 sim 0 21 H1 sim 1 22 K9 sim 1 E gostaria que ficasse assim: cap=c("A1","A3","B2","B5","C6","C8","D4","F5","G8","H1","K9") Var1=c("nao","nao","nao","nao","nao","nao","nao","nao","nao","nao","nao") Var2=c(2,5,7,4,1,2,4,8,9,6,7) Var3=c("sim","sim","sim","sim","sim","sim","sim","sim","sim","sim","sim") Var4=c(4,5,3,2,4,5,8,7,0,1,1) df=cbind(cap,Var1,Var2,Var3,Var4) df=as.data.frame(df) df cap Var1 Var2 Var3 Var4 1 A1 nao 2 sim 4 2 A3 nao 5 sim 5 3 B2 nao 7 sim 3 4 B5 nao 4 sim 2 5 C6 nao 1 sim 4 6 C8 nao 2 sim 5 7 D4 nao 4 sim 8 8 F5 nao 8 sim 7 9 G8 nao 9 sim 0 10 H1 nao 6 sim 1 11 K9 nao 7 sim 1 Alguém poderia me ajudar? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 2

visualizar redes grandes com igraph
by Roney Fraga Souza 08 Mar '14

08 Mar '14

Caros, Estou tentando aprimorar a visualização de uma rede grande onde possa identificar a evolução de grupos. A menor rede que tenho tem 5000 vértices e 40000 edges, sendo que a maior tem 240000 vertices e mais de 2 milhões de edges (ainda não sei como trabalhar com essa rede no R, mas isso é outra história). A melhor visualização de rede que encontrei foi essa (figura a): https://dl.dropboxusercontent.com/u/61883020/rede.grande.grupos.tiff o autor utiliza o algoritmo de visalização large graph layout (LGL), disponível no pacote igraph em layout.lgl(). O problema é que não consegui combinar os argumentos da função para gerar uma rede em formato 'esférico', com baixa nível de sobreposição dos grupos. Alguém tem alguma sugestão? Atenciosamente Roney # ------------------------------ # Início código R # ------------------------------ # carregando os pacotes igraph() e plyr() library(igraph) library(plyr) # criando uma rede aleatória com quatro componentes fortemente conectados # isso vai garantir grupos homogêneos # # para realizar testes optei por criar uma rede pequena, pois, o tempo necessário para visualizar uma rede # maior aumenta significativamente g <- erdos.renyi.game(50, p=1/2) + erdos.renyi.game(50, p=1/2) + erdos.renyi.game(50, p=1/2) + erdos.renyi.game(50, p=1/2) # adicionando links entre os quatro diferentes componentes g <- add.edges( g, c(1,51, 1,52, 51,101, 51,102, 101,151, 102,152, 150,151, 1,101, 1,102, 1,151, 1,200 ) ) g # nomeando os vértices em uma sequencia de 1 a 200 (necessário para fazer o merge() nos próximos passos) V(g)$name <- 1:200 # clusterizando # A Clauset, MEJ Newman, C Moore: Finding community structure in very large networks fc <- fastgreedy.community(as.undirected(g)) # quatro grupos foram identificados, como era de ser esperar table(membership(fc)) sizes(fc) # adicionado qual grupo cada vértice pertence V(g)$grupo <- membership(fc) # pretendo colorir os edges conforme o grupo de quem esta citando, # para isso vou fazer um merge para adicionar esse atributo dos vértices nos edges # criando um data frame dos vértices e de qual grupo este pertence g.ver <- data.frame( name=V(g)$name, grupo=V(g)$grupo ) head(g.ver) # criando um data frame dos edges g.edg <- as.data.frame(get.edgelist(g)) head(g.edg) g.edg <- plyr::rename(g.edg, c('V1'='citante')) g.edg <- plyr::rename(g.edg, c('V2'='citado')) # adicionando qual grupo cada edge pertence, assumindo que se um vértice pertence ao # grupo 1 e tem 8 edges com outros vértices, todos os 8 nós pertencem ao grupo 1 g.edg2 <- merge( g.edg, g.ver, by.x = 'citante', by.y = 'name', all.x = TRUE, all.y = F ) dim(g.edg2) head(g.edg2) # adicionando cores nos edges para os diferentes grupos E(g)$grupo.citante <- g.edg2$grupo E(g)[grupo.citante==1]$color <- 'dark green' E(g)[grupo.citante==2]$color <- 'dark red' E(g)[grupo.citante==3]$color <- 'dark blue' E(g)[grupo.citante==4]$color <- 'sky blue' # adicionado cores nos vértices para os diferentes grupos V(g)[grupo==1]$color <- 'dark green' V(g)[grupo==2]$color <- 'dark red' V(g)[grupo==3]$color <- 'dark blue' V(g)[grupo==4]$color <- 'sky blue' # adicionado cores na borda dos vértices para os diferentes grupos V(g)[grupo==1]$vertex.frame.color <- 'dark green' V(g)[grupo==2]$vertex.frame.color <- 'dark red' V(g)[grupo==3]$vertex.frame.color <- 'dark blue' V(g)[grupo==4]$vertex.frame.color <- 'sky blue' # http://igraph.sourceforge.net/doc/R/layout.html # definir o layout da rede l1 <- layout.lgl( g # maxiter = 150, # The maximum number of iterations to perform (150) # maxdelta = vcount(g), # The maximum change for a vertex during an iteration (the number of vertices) # area = 30*vcount(g)^2 # The area of the surface on which the vertices are placed (square of the number of vertices) # coolexp = , # The cooling exponent of the simulated annealing (1.5) # repulserad = vcount(g) # Cancellation radius for the repulsion (the area times the number of vertices) # cellsize = # The size of the cells for the grid. When calculating the repulsion forces between vertices only vertices in the same or neighboring grid cells are taken into account (the fourth root of the number of area) # root = # The id of the vertex to place at the middle of the layout. The default value is -1 which means that a random vertex is selected ) plot( g, layout = l1, vertex.label = NA, vertex.size = 0.0000001, edge.arrow.size = FALSE, edge.arrow.width = FALSE, ) # ------------------------------ # Fim código R # ------------------------------

1 0

Re: [R-br] Teste Paralelismo no r
by Robert Iquiapaza 06 Mar '14

06 Mar '14

Geovane Assumindo que as regressões propostas nos gráficos (p8) fazem sentido e satisfazem os pressupostos do modelo de regressão linear poderia fazer o seguinte: library(car) escol='Superior/Pós graduação' model2=lm(Risco~0+(Idade*Quintis),subset(dados,Escolaridade3cat==escol)) summary(model2) # o teste para ‘Idade:Quintis5‘ na verdade está testando se os coeficientes angulares são iguais ( ou a diferença deles é zero) #para testar a diferença dos interceptos linearHypothesis(model2, "Quintis1 - Quintis5 =0") linearHypothesis(model2, "Idade:Quintis5 =0") Aqui mostrei os testes para uma escolaridade (escol), na verdade vc tem mais informação, teria que rodar um só modelo com todas as informações e fazer os testes nele: model1=lm(Risco~0+Idade*Escolaridade3cat*Quintis,data=dados) summary(model1) Boa sorte Robert From: Benilton Carvalho Sent: Thursday, March 06, 2014 2:15 PM To: r-br ; geovane barbosa Subject: Re: [R-br] Teste Paralelismo no r Geovane, Para realizar uma analise de dados, o estatistico precisa entender a (hi/e)storia dos dados. Por exemplo, sem saber do que se tratam os dados, e' dificil acreditar que a regressao linear simples que voce apresenta no seu exemplo seja adequada. A outra coisa: como/onde entra o fator "Quintis" existente no seu conjunto de dados? Enfim, supondo adequacidade em todos os sentidos (cabera' a vc verificar os pressupostos), se o interesse e' averiguar a identidade das regressoes nos tres grupos escolares (vou ignorar Quintis aqui, pois realmente nao sei o que e', nem como entra na sua analise), eu ajustaria os modelos com interacao (vis-a-vis Escolaridade) e compararia com o modelo simplificado (que contem apenas Idade): fit1 = lm(Risco ~ Idade*Escolaridade3cat, data=dados) fit0 = lm(Risco ~ Idade, data=dados) anova(fit0, fit1) benilton Em 6 de março de 2014 10:43, geovane barbosa <geovanecb(a)yahoo.com.br> escreveu: Olá pessoal , segue em anexo um exemplo reprodutível. Quero verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e banco num link. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") https://www.dropbox.com/s/ymyzh7p0xhmumda/bancoteste.csv Em Quinta-feira, 6 de Março de 2014 10:12, geovane barbosa <geovanecb(a)yahoo.com.br> escreveu: Olá Prof, segue em anexo um exemplo reprodutível. Como eu lhe disse anteriormente, queria verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e o banco de dados ficticio. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") Em Quarta-feira, 5 de Março de 2014 16:54, Benilton Carvalho <beniltoncarvalho(a)gmail.com> escreveu: A gente pede que exemplos reprodutíveis sejam dados para que possamos dar uma resposta definitiva sem a necessidade de idas e vindas. Como eu não entendi quais são os três cenários, fica a tarefa p vc: o caso a que vc se refere consiste em testar se as interações das variáveis de interesse com os referidos grupos são conjuntamente iguais a zero. b On Mar 5, 2014 11:00 AM, "geovane barbosa" <geovanecb(a)yahoo.com.br> wrote: Olá pessoal bom dia. Estou com uma questão que não consigo soluciona-la. Realize um regressão linear múltipla para comprar um grau de risco por idade e por escolaridade. Quero verificar se a inclinação proporcionada por essas retas é estatisticamente significativa, ou seja comparando a inclinação nos três cenários. bem como tambem seu intercepto. Como eu faço isso? Seria o test de Brant? Não consigo encontrar tal ferramenta no r. abraços _______________________________________________ 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. _______________________________________________ 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.

1 0

teste kappa por categorias
by branco.sk8＠ig.com.br 06 Mar '14

06 Mar '14

prezados, gentileza. como se faz para fazer o teste kappa por categorias (condicionar a outra variável?) (ex. testar concordância kappa entre dois teste em função do sexo, ou em função da idade etc) (OBS. já sei fazer o teste kappa, só não sei fazer o kappa por outras variáveis... obrigado

2 1

Re: [R-br] Teste Paralelismo no r
by geovane barbosa 06 Mar '14

06 Mar '14

Olá pessoal , segue em anexo um exemplo reprodutível. Quero verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e banco num link. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") https://www.dropbox.com/s/ymyzh7p0xhmumda/bancoteste.csv Em Quinta-feira, 6 de Março de 2014 10:12, geovane barbosa <geovanecb(a)yahoo.com.br> escreveu: Olá Prof, segue em anexo um exemplo reprodutível. Como eu lhe disse anteriormente, queria verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e o banco de dados ficticio. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") Em Quarta-feira, 5 de Março de 2014 16:54, Benilton Carvalho <beniltoncarvalho(a)gmail.com> escreveu: A gente pede que exemplos reprodutíveis sejam dados para que possamos dar uma resposta definitiva sem a necessidade de idas e vindas. Como eu não entendi quais são os três cenários, fica a tarefa p vc: o caso a que vc se refere consiste em testar se as interações das variáveis de interesse com os referidos grupos são conjuntamente iguais a zero. b On Mar 5, 2014 11:00 AM, "geovane barbosa" <geovanecb(a)yahoo.com.br> wrote: > >Olá pessoal bom dia. > > > > >Estou com uma questão que não consigo soluciona-la. > > > > >Realize um regressão linear múltipla para comprar um grau de risco por idade e por escolaridade. > > >Quero verificar se a inclinação proporcionada por essas retas é estatisticamente significativa, ou seja comparando a inclinação nos três cenários. bem como tambem seu intercepto. Como eu faço isso? Seria o test de Brant? Não consigo encontrar tal ferramenta no r. > > >abraços > > > > > >_______________________________________________ >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.

2 1

Coordenadas Geográficas no Pacote GeoR
by Helder Gramacho 06 Mar '14

06 Mar '14

Boa tarde, É possível trabalhar com coordenadas geográficas (em graus decimais) no pacote GeoR? Caso seja possível, para estimar o semivariograma utilizamos as distâncias calculadas em metros ou em graus decimais? A maioria dos trabalhos que tenho visto tem utilizado coordenadas planas UTM, e aí surgiu a dúvida: Quando a região estiver situada entre dois fusos UTM como proceder? Att, Hélder Santos *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 4

Teste Paralelismo no r
by geovane barbosa 05 Mar '14

05 Mar '14

Olá pessoal bom dia. Estou com uma questão que não consigo soluciona-la. Realize um regressão linear múltipla para comprar um grau de risco por idade e por escolaridade. Quero verificar se a inclinação proporcionada por essas retas é estatisticamente significativa, ou seja comparando a inclinação nos três cenários. bem como tambem seu intercepto. Como eu faço isso? Seria o test de Brant? Não consigo encontrar tal ferramenta no r. abraços

2 1

Colocar estatística Z na curva normal - TH
by Andre Oliveira 04 Mar '14

04 Mar '14

Pessoal, estou aqui tentando, mas sem sucesso colocar a estatística Ztab, Zcal, mas sem sucesso. media <- 170 stder <- 8 x <- seq(media-3*stder, media+3*stder, 0.01) y <- dnorm(x, media, stder) plot(x, y, type='l',axes=F,xlab="",ylab="") abline(h=0) hx <- seq(media+2*stder, media+3*stder, 0.01) hy <- dnorm(hx, media, stder) n <- length(hy) polygon(c(hx, rev(hx)), c(hy, rep(0, n)), col=2) # Tentei text(locator(n=1),"Zcal=3.53") text(locator(n=1),"Ztab=1,64") Mas, preciso aumentar a área abaixo do gráfico para plotar estes valores e não consegui! Alguém tem alguma sugestão? Outra situação é que Z_{cal} para ficar mais elegante não funcionou com locator. obrigado

3 4