Dezembro 2011 - R-br

Ajuda
by Gilbert Queiroz 05 Dez '11

05 Dez '11

Pessoal, tenho a seguinte parte de um script que estou desenvolvendo: n<-c(87, 108, 80, 25) err.m<-matrix(NA, choose(nx,2),1) err.m nx for(i in 1:nx) { for(j in 2:nx) { if(i<j){ err<-sqrt((410.35/(n[i]+0.5))+(410.35/(n[j]+0.5))) print(err) } } } Informações úteis: nx = 4 (já foi definido em outra parte do script) err.m é uma matriz-coluna (6 x 1) O quero fazer é colocar os resultados de err dentro da matriz err.m Como faço isso??? O comando print me mostra que err está calculando o que quero. Mas me interessa colocar os resultados em err.m Abs.

5 6

Função com problema no Rprofile
by Emmanuel Arnhold 05 Dez '11

05 Dez '11

Criei uma função bem simples para simplificar a digitação da read.table("clipboard", h=T) que me é bastante útil. Mas a função não funciona via Rprofile. Alguém tem alguma sugestão? Testem a função do R e depois tentem inserir no Rprofile e usar novamente para ver o erro. A função é: clip=function(x=read.table("clipboard ", h=T)){x=read.table("clipboard", h=T); print(x)} clip()

2 1

Re: [R-br] teste Duncan
by Bruna Gabriela Wendpap 05 Dez '11

05 Dez '11

Há um pacote não oficial do R para análise de experimentos chamado ExpDes (Experimental Designs). A página abaixo contém link para download do pacote e tutorial em português. https://sites.google.com/site/ericbferreira/unifal/downloads-1 Att, Bruna Gabriela Wendpap > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 10, assunto 4 > To: r-br(a)listas.c3sl.ufpr.br > Date: Sat, 3 Dec 2011 12:00:01 -0200 > > 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: teste Duncan (Andre Oliveira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 2 Dec 2011 13:25:26 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, > "crabak(a)acm.org" <crabak(a)acm.org> > Subject: Re: [R-br] teste Duncan > Message-ID: > <1322861126.39856.YahooMailNeo(a)web161805.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Uma alternativa é o pacote laercio. > require(laercio) > ?LDuncan > > > > André Oliveira Souza > > > ________________________________ > De: Cesar Rabak <cesar.rabak(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br > Enviadas: Sexta-feira, 2 de Dezembro de 2011 10:42 > Assunto: Re: [R-br] teste Duncan > > Em 2/12/2011 09:31, Luiz André Zardo escreveu: > > Prezados, > > Qual pacote posso utilizar para rodar o teste Duncan para comparação de > > medias no R?? > No prompt do seu terminal R digite o seguinte: > > > RSiteSearch("Duncan test") > > Na primeira página dos resultados já vi dois pacotes. . . > > > -- Cesar Rabak > GNU/Linux User 52247. > Get counted: http://counter.li.org/ > _______________________________________________ > 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/20111202/f2833cee/atta…> > > ------------------------------ > > _______________________________________________ > 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 10, assunto 4 > ****************************************

1 0

teste Duncan
by Luiz André Zardo 05 Dez '11

05 Dez '11

Prezados, Qual pacote posso utilizar para rodar o teste Duncan para comparação de medias no R?? Obrigado.

4 3

Distribuição de Weibull
by Rodrigo Muniz 05 Dez '11

05 Dez '11

Boa noite pessoal. Estou com o seguinte script: require(MASS) getwd() setwd("C:/Users/Rodrigo/Documents/Pos_graduacao/DoutoradoUSP/Geoestatistica") Vento <- read.table("Weibull_full.txt", head=T) WVento50m <- fitdistr(Vento[[1]], densfun="weibull") WVento30m <- fitdistr(Vento[[2]], densfun="weibull") WVento15m <- fitdistr(Vento[[3]], densfun="weibull") WVento10m <- fitdistr(Vento[[4]], densfun="weibull") WVento3m <- fitdistr(Vento[[5]], densfun="weibull") O arquivo Weibull_full é um data frame com 5 coluna sendo a variável velocidade do vento em cinco diferentes níveis de altura. Estou tentanto calcular os parâmetros de forma e escala da distribuição de Weibull para cada nível de altura. No entanto, só consigo calcular para o nível de altura de 50 m (WVento50m), quando mando rodar os outros níveis dá o seguinte erro: Error in fitdistr(Vento[[2]], densfun = "weibull") : Weibull values must be > 0 Já conferi o arquivo dos dado para ver se há algo errado com os dados, mais não há. Tentei fazer a análise separada mais deu o mesmo erro. Um observação que eu faço é que embora eu não consiga rodar a função "fitdistr", eu consigo plotar os gráficos de histograma normalmente para todos os níveis de altura. Se alguém puder me ajudar eu agradeço. Abraço a todos. -- *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)

3 4

Algoritmo para simular VA GAMA
by Victor Eduardo 02 Dez '11

02 Dez '11

Gente, onde posso encontrar um algoritmo para simular VA Gama?? Eu uso o livro Simulation do Sheldon Ross, e ao que parece, lá não tem esse algoritmo. Posso simular uma gama quando o alfa for um numero inteiro por soma de exponenciais, e quando o mesmo não for, ou seja, e quando, por exemplo, o alfa for < 1, como poderia fazer? OBS: As restrições já foram adicionadas na minha função, ou seja, o alpha tem que ser positivo e o beta também. Abraços!

2 1

Estimar vários modelos não lineares de uma só vez
by Emmanuel Arnhold 02 Dez '11

02 Dez '11

Gostaria de saber se é possível fazer (e como fazer se for possível) via nls algo semelhante ao que segue abaixo via lm, ou seja, estimar vários modelos de uma só vez. No caso o modelo tem a mesma variável explicativa e os mesmos starts, mudam as variáveis respostas. x=c(0:5);y1=c(50:55);y2=c(500:505) modelo=lm(cbind(y1,y2)~x);modelo

3 3

Re: [R-br] Versão R-2.14.0 no Ubunto - RESOLVIDO
by marlene marchena 01 Dez '11

01 Dez '11

Oi! Somente para fechar este post. Depois de pedir ajuda em outros forums descubrí que tinha a versão R (2.13.01) instalada manualmente (sem usar a infraestructura apt-get) no meu diretório HOME. Por isso não conseguia desinstalá-lo. Bom, a única coisa que fiz foi deletar a pasta R 2.13.01 do diretório HOME e reinstalar a versão 2.14, isso resolviu tudo. Obrigada pela ajuda! Marlene Em 27 de novembro de 2011 20:53, Benilton Carvalho < beniltoncarvalho(a)gmail.com> escreveu: > R --vanilla > > No terminal > On Nov 27, 2011 4:45 PM, "marlene marchena" <marchenamarlene(a)gmail.com> > wrote: > >> Gente estou precisando de ajuda! >> >> Eu instalei o R 2.14 seguindo o link acima e os passos indicados pelo >> Pedro no ubuntu 11.10. Eu rodei o R no terminal e apareceu q a versão era >> a 2.14. Porém agora pouco eu percebi que após reiniciado o computador meu R >> volta para a versão que timha antes a 2.13.1 >> >> Não entendo o que fiz de errado? >> >> Alguém pode me ajudar a entender o que acontece. >> >> Grata, >> >> Marlene. >> >> >> Em 25 de novembro de 2011 13:57, Daniel C Bezerra < >> danielcbezerra(a)gmail.com> escreveu: >> >>> Olá Marlene (e quem estiver iniciando Ubuntu), >>> >>> Já que vc está usando Ubuntu 10.04 seguem mais duas sugestões que podem >>> te ajudar: >>> >>> Agora que vc instalou o repositório do R, vale à pena vc entrar em >>> Sistema>Administração>Gerenciador de Pacotes Synaptic e procurar por todo >>> os programas que tenham "r-cran" (sem aspas). Vc vai ver diversos pacotes >>> do R disponíveis. >>> >>> Vale à pena instalar todos os pacotes que vc usa por ali (pelo menos da >>> primeira vez) pois isto aumentará a probabilidade de que todas as >>> dependências que os pacotes precisam para funcionar também sejam >>> instaladas. >>> >>> Depois é só rodar no R: >>> update.packages(ask=F,dep=T,checkBuilt=T) >>> para que tudo seja atualizado. >>> >>> Por fim, mesmo assim pode acontecer que ao tentar instalar um pacote vc >>> não consiga por falta de dependências. Por exemplo imagine que vc não >>> esteja conseguindo instalar o pacote rgl. Para solucionar esta questão >>> basta rodar no terminal: >>> >>> sudo apt-get build-dep r-cran-rgl >>> >>> Assim todas as dependências serão instaladas. >>> >>> Abs, >>> >>> D >>> >>> 2011/11/25 Pedro Rafael <pedro.rafael.marinho(a)gmail.com> >>> >>>> Ok, que bom que conseguiu Marlene. Qualquer coisa, estamos as ordens. >>>> >>>> Em 24 de novembro de 2011 21:39, marlene marchena [via R-br] < >>>> ml-node+s2285057n4106046h3(a)n4.nabble.com> escreveu: >>>> >>>>> Oi Pedro, >>>>> >>>>> Obrigada pelo roteiro. Graças a teu post finalmente conseguí instalar >>>>> o R 2.14. Espero que o Marcelo tenha conseguido tb. >>>>> >>>>> Um abraço, >>>>> >>>>> Marlene. >>>>> >>>>> >>>>> >>>>> Em 22 de novembro de 2011 23:50, Pedro Rafael <[hidden email]<http://user/SendEmail.jtp?type=node&node=4106046&i=0> >>>>> > escreveu: >>>>> >>>>>> Mas o R está meio antigo, já estamos na versão 2.14. Qual a versão do >>>>>> UbuntU que você está utilizando? Lembrando que é Ubuntu e não Ubunto. Faz >>>>>> assim: >>>>>> >>>>>> sudo su >>>>>> cd /etc/apt/ >>>>>> gedit sources.list >>>>>> >>>>>> O último comando abrirá um arquivo de texto. No final do arquivo >>>>>> coloque: deb http://cran-r.c3sl.ufpr.br/bin/linux/ubuntu oneiric/ >>>>>> Após copiar salve o arquivo e feche. Posteriormente ainda no terminal >>>>>> rode: apt-get update para reconhecer o novo respositorio. >>>>>> >>>>>> Após isso rode no terminal os dois comandos abaixo: >>>>>> >>>>>> gpg --keyserver keyserver.ubuntu.com --recv-key E084DAB9 >>>>>> gpg -a --export E084DAB9 | sudo apt-key add - >>>>>> >>>>>> Novamente rode: apt-get update. Posteriormente remova a versão antiga >>>>>> do R com: apt-get remove r-base. Após a remoção do R instale novamente >>>>>> usando apt-get install r-base e você terá a versão 2.14 do R. Depois >>>>>> instale o pacote do R para ver se dar certo. >>>>>> >>>>>> >>>>>> >>>>>> Em 22 de novembro de 2011 19:40, FHRB Toledo [via R-br] <[hidden >>>>>> email] <http://user/SendEmail.jtp?type=node&node=4106046&i=1>>escreveu: >>>>>> >>>>>> Você deve adicionar o repositório no seu arquivo >>>>>>> /et/apt/sources.list ... veja a explicação em: >>>>>>> http://cran.r-project.org/bin/linux/ubuntu/README >>>>>>> >>>>>>> att, >>>>>>> FH >>>>>>> >>>>>>> 2011/11/22 Marcelo Cardoso mello <[hidden email]<http://user/SendEmail.jtp?type=node&node=4097703&i=0> >>>>>>> > >>>>>>> >>>>>>> Pessoal, executei no terminal: >>>>>>>> $ sudo apt-get install r-base >>>>>>>> $ sudo apt-get install r-recommendad >>>>>>>> $ sudo apt-get etc etc que me foi indicado por aqui >>>>>>>> >>>>>>>> Tudo ok. >>>>>>>> >>>>>>>> Mas o R instalado é da versão 2.10 >>>>>>>> >>>>>>>> Tudo bem até aí. O problema é que não instala pacotes recentes. >>>>>>>> >>>>>>>> usei: instal.packages("vegan") tudo ok ... instalou >>>>>>>> >>>>>>>> mas não instala o lme4 por exemplo. >>>>>>>> >>>>>>>> O problema deve ser da versão do R que tenho? Como faço para >>>>>>>> instalar a versão R-2.14.0 no Ubunto? De forma rápida e simples. >>>>>>>> >>>>>>>> Obrigado por qualquer dica. >>>>>>>> >>>>>>>> >>>>>>>> _______________________________________________ >>>>>>>> R-br mailing list >>>>>>>> [hidden email]<http://user/SendEmail.jtp?type=node&node=4097703&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. >>>>>>>> >>>>>>> >>>>>>> >>>>>>> _______________________________________________ >>>>>>> R-br mailing list >>>>>>> [hidden email]<http://user/SendEmail.jtp?type=node&node=4097703&i=2> >>>>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e >>>>>>> forneça código mínimo reproduzível. >>>>>>> >>>>>>> ------------------------------ >>>>>>> If you reply to this email, your message will be added to the >>>>>>> discussion below: >>>>>>> >>>>>>> http://r-br.2285057.n4.nabble.com/R-br-Versao-R-2-14-0-no-Ubunto-tp4097586p… >>>>>>> To unsubscribe from R-br, click here. >>>>>>> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >>>>>>> >>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> Saudações, >>>>>> Pedro Rafael Diniz Marinho. >>>>>> Estatístico - Secretaria de Estado da Saúde - PB. >>>>>> >>>>>> >>>>>> >>>>>> _______________________________________________ >>>>>> R-br mailing list >>>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4106046&i=2> >>>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>>> código mínimo reproduzível. >>>>>> >>>>> >>>>> >>>>> _______________________________________________ >>>>> R-br mailing list >>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4106046&i=3> >>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>> código mínimo reproduzível. >>>>> >>>>> ------------------------------ >>>>> If you reply to this email, your message will be added to the >>>>> discussion below: >>>>> >>>>> http://r-br.2285057.n4.nabble.com/R-br-Versao-R-2-14-0-no-Ubunto-tp4097586p… >>>>> 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…> >>>>> >>>> >>>> >>>> >>>> -- >>>> Saudações, >>>> Pedro Rafael Diniz Marinho. >>>> Estatístico - Secretaria de Estado da Saúde - PB. >>>> >>>> >>>> >>>> _______________________________________________ >>>> 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. >

2 1