- R-br - DINF UFPR

limitando rnorm
by Hélio Gallo Rocha 01 Abr '13

01 Abr '13

Caros Teria como impor limites, máx e min para os valores gerados no uso de rnorm? Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Re: [R-br] Digest R-br, volume 25, assunto 41
by Pamela Chiroque 01 Abr '13

01 Abr '13

Disculpem alguem sabe se existe a funcao pFq (Hypergeometrica?), ajuda por favor. Obrigada 2013/3/31 <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. duvidas sobre o r (Juraci Mendes Moreira) > 2. Re: duvidas sobre o r (Daniel Marcelino) > 3. Re: duvidas sobre o r (Tiago Souza Marçal) > 4. Problemas com grid.arrange (Bernardo Rangel Tura) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 30 Mar 2013 11:06:32 -0700 (PDT) > From: Juraci Mendes Moreira <juracimendesmoreira(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] duvidas sobre o r > Message-ID: > <1364666792.97559.YahooMailClassic(a)web161005.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > como fazer analise de variancia em quadrado latino. qual o comando de anana > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130330/db11dc93/atta… > > > > ------------------------------ > > Message: 2 > Date: Sat, 30 Mar 2013 15:15:06 -0400 > From: Daniel Marcelino <dmarcelino(a)live.com> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] duvidas sobre o r > Message-ID: > <CAAjbdiK= > KN2rZWdymDiXUpHnavEfeJpDKQp0RJZ4S8i64ge88w(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Juraci, dê uma olhada nas funções ?aov, ?lme e ?lmer. > Além disso a sua pergunta não esclarece o que você precisa. O bom seria > você enviar um exemplo dos seus dados. O experimento que você está > realizando tem medidas repetidas ou não? > Por exemple, para analisar a variância em um quadrado latino simples, você > pode usar assim: > > out <- aov(time~participante+modo) > summary(out) > > > > 2013/3/30 Juraci Mendes Moreira <juracimendesmoreira(a)yahoo.com.br> > > > como fazer analise de variancia em quadrado latino. qual o comando de > anana > > _______________________________________________ > > 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. > > > > > > -- > \begin{signature} > Daniel Marcelino > ? dm.silva(a)umontreal.ca > ? (514) 343 6111 #3799 > Skype: d.marcelino > ? 3200 Jean Brillant, Office C5071 > Montreal, QC; H3T 1N8 > Canada > \end{signature} > > "Small steps toward a much better world" > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130330/1ae3897e/atta… > > > > ------------------------------ > > Message: 3 > Date: Sat, 30 Mar 2013 23:55:24 +0300 > From: Tiago Souza Marçal <tiagosouzamarcal(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] duvidas sobre o r > Message-ID: <BLU170-W124F2B19F369AEAFA04AE7ED9DC0(a)phx.gbl> > Content-Type: text/plain; charset="utf-8" > > Para fazer uma anova em dql é só usar os seguntes comantos: > trat<-factor(trat) linha<-factor(linha) coluna<-factor(coluna) > av<-aov(resp~linha+coluna+trat) # Como pode ver a linha e a coluna são > somadas aos tratamentos. Caso você queira fazer uma analise completa de um > experimento dql eu te indico a usar a função dql()::ExpDes Att. Tiago. > Date: Sat, 30 Mar 2013 11:06:32 -0700 > From: juracimendesmoreira(a)yahoo.com.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] duvidas sobre o r > > como fazer analise de variancia em quadrado latino. qual o comando de anana > _______________________________________________ > 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/20130330/41e0f417/atta… > > > > ------------------------------ > > Message: 4 > Date: Sun, 31 Mar 2013 09:02:37 -0300 > From: Bernardo Rangel Tura <tura(a)centroin.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Problemas com grid.arrange > Message-ID: <515825DD.70306(a)centroin.com.br> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > > Boa tarde, > > estou precisando fazer um gráfico chamado forest plot porem não posso > usar os já prontos no R. Já esta quase pronto utilizando este script. > > require(gridExtra) > nomes <- sample(letters,5) > logOR <- abs(rnorm(5,0,.3)) > sd_logOR <-rbeta(5,10,10) > OR <- round(exp(logOR),2) > lci_OR <- round(exp(logOR-1.96*sd_logOR),2) > uci_OR <- round(exp(logOR+1.96*sd_logOR),2) > OR <- paste(OR,' (',lci_OR,' - ',uci_OR,')',sep='') > df <- data.frame(nomes,OR,lci_OR,uci_OR) > p1 <- tableGrob(df[,1:2],cols=c('Estudos','OR (CI)'),row.just > ='left',col.just ='center', core.just ='right') > require(ggplot2) > OR <- round(exp(logOR),2) > df <- data.frame(nomes,OR,lci_OR,uci_OR) > p2 <- p2 <- ggplot(df, aes(x=nomes, y=OR, ymin=lci_OR, > ymax=uci_OR))+geom_hline(yintercept=1,lty=2)+expand_limits(y = > c(0,max(uci_OR)*1.1))+geom_pointrange()+coord_flip()#+ > geom_point(aes(size = sd_logOR)) > grid.arrange(p1,p2,ncol=2) > > Porém aqui no fim há uma desproporção entre os dados da tabela com os > dados do gráfico. > > Porventura alguém sabe como resolver isto ? > > []s > Tura > > > ------------------------------ > > _______________________________________________ > 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 25, assunto 41 > ***************************************** >

3 2

Problemas com grid.arrange
by Bernardo Rangel Tura 31 Mar '13

31 Mar '13

Boa tarde, estou precisando fazer um gráfico chamado forest plot porem não posso usar os já prontos no R. Já esta quase pronto utilizando este script. require(gridExtra) nomes <- sample(letters,5) logOR <- abs(rnorm(5,0,.3)) sd_logOR <-rbeta(5,10,10) OR <- round(exp(logOR),2) lci_OR <- round(exp(logOR-1.96*sd_logOR),2) uci_OR <- round(exp(logOR+1.96*sd_logOR),2) OR <- paste(OR,' (',lci_OR,' - ',uci_OR,')',sep='') df <- data.frame(nomes,OR,lci_OR,uci_OR) p1 <- tableGrob(df[,1:2],cols=c('Estudos','OR (CI)'),row.just ='left',col.just ='center', core.just ='right') require(ggplot2) OR <- round(exp(logOR),2) df <- data.frame(nomes,OR,lci_OR,uci_OR) p2 <- p2 <- ggplot(df, aes(x=nomes, y=OR, ymin=lci_OR, ymax=uci_OR))+geom_hline(yintercept=1,lty=2)+expand_limits(y = c(0,max(uci_OR)*1.1))+geom_pointrange()+coord_flip()#+ geom_point(aes(size = sd_logOR)) grid.arrange(p1,p2,ncol=2) Porém aqui no fim há uma desproporção entre os dados da tabela com os dados do gráfico. Porventura alguém sabe como resolver isto ? []s Tura

1 0

duvidas sobre o r
by Juraci Mendes Moreira 30 Mar '13

30 Mar '13

como fazer analise de variancia em quadrado latino. qual o comando de anana

3 2

Fwd: Redes neurais
by Bernardo Silva 29 Mar '13

29 Mar '13

Boa tarde, estou a trabalhar em redes neurais utilizando o pacote nnet. Gostaria de saber qual o comando para encontrar o melhor modelo, por exemplo qd trabalhamos em glms utilizamos o stepwise, mas em redes neurais nao sei, podem ajudar? li que tb pode ser o stepwise mas depois nao consigo implementar obrigado

1 0

nomeando pontos
by Hélio Gallo Rocha 28 Mar '13

28 Mar '13

Caros, Busquei em msg anteriores sobre plotar nomes de atributos num mapa, na msg de Alexandro Yahoo parece interessante, até consegui baixar os shapes, mas no script faz referencia a arquivos do excel que ele não anexou, assim fica dificil reproduzir. Estou usando o pacote geoR, Alguém poderia dar uma força neste sentido de plotar o nome de um atributo, tratamento por exemplo, no mapa? grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

6 11

Considerar ângulo da área para distribuir coordenadas com expand.grid ()
by ASANTOS 28 Mar '13

28 Mar '13

Boa noite pessoal, Gostaria de saber se existe uma maneira de controlar a inclinação da distribuição dos pontos em expand.grid() em função de um limite de área, como por exemplo no meu CRM abaixo, onde no exemplo 1 esta tudo Ok, mas no exemplo 2 as coordenadas não consideram a inclinação da área. Seria possível usar algum argumento para controlar isto? # Exemplo 1 x<-c(0,0,100,100,0) y<-c(0,100,100,0,0) plot(x,y, type="l", col="red") # p0<-c(10,10) # # c0 <- seq(p0[1],p0[1]+95,by=20) d0 <- seq(p0[2],p0[2]+95,by=20) d1 <- expand.grid(x=c0, y=d0) points(d1[,1],d1[,2]) # #Exemplo 2 x<-c(0,0,50,50,0) y<-c(0,100,150,50,0) plot(x,y, type="l", col="red") # p0<-c(10,10) # # c0 <- seq(p0[1],p0[1]+95,by=20) d0 <- seq(p0[2],p0[2]+95,by=20) d1 <- expand.grid(x=c0, y=d0) points(d1[,1],d1[,2]) # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 2

Unique com menos variáveis
by Fátima Lima Paula 28 Mar '13

28 Mar '13

Prezados, tenho um banco de dados com muitas variáveis e com linhas idênticas em relação a algumas variáveis. Quero retirar as linhas idênticas, mas para que isso fique melhor é necessário fazer o unique apenas com algumas variáveis. Quando faço isso são eliminadas mais linhas. Mas, por outro lado o unique me retorna um banco com menos variáveis. Como posso fazer para ficar com o banco com todas as colunas e sem as linhas idênticas. Exemplo: nome=c("A","B","C","D","E","F","F","A","C","K","L","M") idade=c("1","2","3","4","5","6","6","1","3","11","12","13") cor=c("azul","preto","verde","branco","amarelo","cinza","cinza", "marrom", "amarelo","bege","azul","azul") material=c("lápis","ignorado","caneta","papel","borracha","lápis","lápis", "lápis","ignorado","caneta","apontador","pilha") banco1=cbind(nome,idade,cor,material) # banco que tenho dim(banco1) banco2=cbind(nome,idade) # banco para fazer o unique dim(banco2) b1unique=unique(banco1) # se faço o unique com o banco todo retornam 11 linhas e 4 colunas dim(b1unique) b2unique=unique(banco2) # fazendo o unique com o banco com as 2 variáveis retornam 9 linhas e 2 colunas dim(b2unique) Para fazer o unique, de forma a eliminar o maior número de repetidas, quero utilizar apenas as variáveis nome e idade. O b1unique (com todas as variáveis) retorna 11 linhas e o b2unique retorna 9. Quero o banco com as nove linhas, mas com todas as colunas do banco1. Verifiquei que existe um argumento incomparables no unique, mas vi que não é aplicável para bancos grandes. Alguém poderia me ajudar, por favor? Obrigada Fátima

2 6

Re: [R-br] Teste para decisão do número ótimo de amostra [RESOLVIDO]
by Hélio Gallo Rocha 27 Mar '13

27 Mar '13

Caros Alexandre e Walmes. Estou usando este teste de decisão de número de amostras para ver se minha amostragem está ok. Mas fiquei na dúvida sobre o valor 10000 da linha pop.total <- rnorm(10000, mean=20, sd=15) Vou exemplificar com dados de previsão de safra em café, no caso, são relativos a altura da planta de café. Estas medidas referem a coleta em 3 pontos da lavoura, cada ponto é composto por 5 plantas, então sd e mean são de 15 plantas. Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 plantas. sd=0.2587516 mean= 3.013333 assim acho que o valor 10000 seria 1775, ok? Posso determinar o numero de amostras para para uma área de 20 ha com 37085 plantas? Grato Hélio Em 27 de março de 2013 12:30, Hélio Gallo Rocha <heliogallorocha(a)gmail.com>escreveu: > Caros Alexandre e Walmes. > > Estou usando este teste de decisão de número de amostras para ver se minha > amostragem está ok. > > Mas fiquei na dúvida sobre o valor 10000 da linha > > pop.total <- rnorm(10000, mean=20, sd=15) > > > > Vou exemplificar com dados de previsão de safra em café, no caso, são > relativos a altura da planta de café. > Estas medidas referem a coleta em 3 pontos da lavoura, cada ponto é > composto por 5 plantas, então sd e mean são de 15 plantas. > Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 > plantas. > > sd=0.2587516 > mean= 3.013333 > > assim acho que o valor 10000 seria 1775, ok? > > Posso determinar o numero de amostras para para uma área de 20 ha > com 37085 plantas? > > Grato > > Hélio > > > > Em 20 de março de 2013 18:49, ASANTOS [via R-br] < > ml-node+s2285057n4658689h27(a)n4.nabble.com> escreveu: > > Perfeito Walmes, >> >> Era isto mesmo que eu buscava, >> >> Abraço, >> >> Alexandre >> >> >> Em 20/03/2013 17:31, Walmes Zeviani escreveu: >> >> Alexandre, >> >> Segue um CMR simples: >> >> # tamanho de amostra para medida de interesse que >> # é o erro absoluto ao estimar a média da população >> # erro = abs(média_estimada-média_verdadeira) >> >> pop.total <- rnorm(10000, mean=20, sd=15) # população >> m <- mean(20) # média verdadeira >> i <- unique(floor(2^seq(1, 12, by=0.5))); i # tamanhos de amostra >> nr <- 100 # número de execuções de >> reamostragem >> fun <- function(size, pop, m){ >> # função que calcula o erro a partir de uma amostra >> x <- sample(pop, size, replace=TRUE) >> # retorna a diferença entre estimado e real >> abs(mean(x)-m) >> } >> >> results1 <- sapply(i, >> function(j){ >> replicate(nr, fun(size=j, pop=pop.total, m=m)) >> }) >> str(results1) # matriz: colunas correspondem à i >> results2 <- colMeans(results1) # tira a média para cada i >> results2 >> >> plot(results2~i) >> plot(results2~log2(i)) >> >> # qual o tamanho de amostra para ter um erro no ponto >> # médio entre o máximo (quando n=1) e o mínimo (n->infinito, erro=0) >> >> n0 <- nls(results2~A*B^(-log2(i)+1), start=list(A=1, B=1)) >> coef(n0) >> >> perc <- 0.3 # erro tolerado de 30% do máximo >> (usuário define) >> xi <- -log(perc)/log(coef(n0)["B"])+1 # tamanho de amostra correspondente >> >> plot(results2~log2(i), >> xlim=c(1, max(log2(i))), ylim=c(0,max(results2, coef(n0)["A"])), >> xlab="Tamanho da amostra (log2 unidades)", >> ylab="Erro absoluto médio ao estimar a média") >> abline(v=1, lty=2) >> with(as.list(coef(n0)), curve(A*B^(-x+1), add=TRUE)) >> with(as.list(coef(n0)), abline(h=A*c(1,perc), v=xi, lty=2, col=2)) >> >> ceiling(2^xi) # tamanho de amostra à ser considerado >> >> À 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: [hidden email]<http://user/SendEmail.jtp?type=node&node=4658689&i=0> >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> >> _______________________________________________ >> R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=1>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. >> >> >> -- >> ====================================================================== >> Alexandre dos Santos >> Proteção Florestal >> Coordenador do curso Técnico em Florestas >> Vice Coordenador do curso de Engenharia Florestal >> IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso >> Campus Cáceres >> Caixa Postal 244 >> Avenida dos Ramires, s/n >> Bairro: Distrito Industrial >> Cáceres - MT CEP: 78.200-000 >> Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)[hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=2> >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=3> >> ====================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=4> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Teste-para-decisao-do-numero-otimo-d… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >> . >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> > > > > -- > Hélio Gallo Rocha > IFSULDEMINAS - Câmpus Muzambinho > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Re: [R-br] text mining
by Vinicius Brito Rocha 27 Mar '13

27 Mar '13

tura , estou usando esse. parece bom. Obrigado. aliais se tiver material sobre o assunto poderia me mandar? abs Vinicius Enviado por Samsung Mobile -------- Mensagem original -------- De : Bernardo Rangel Tura <tura(a)centroin.com.br> Data: Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] text mining On 03/17/2013 03:49 PM, Vinicius Brito Rocha wrote: > Senhores, > > estou começando a estudar text mining. Quais pacotes, os senhores > recomendam? > > -- > /Vinicius Brito Rocha./ Vinicius, No R tem um pacote chamado tm que uso com bons resultados. Existem vários outros que você pode ver em http://cran.r-project.org/web/views/NaturalLanguageProcessing.html []s Tura _______________________________________________ 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