Outubro 2011 - R-br

Ajuda em desenho de experimentos
by Bernardo Rangel Tura 07 Out '11

07 Out '11

Mestres em R e DOE! Preciso de uma ajuda para resolver um projeto de pesquisa. No meu experimento temos 156 níveis (ou tratamentos) que preciso analisar e só posso pesquisar 6 em cada pessoa entrevistada. Pensei em trabalhar com um modelo de blocos incompletos balanceados (BIB). Desta forma seriam necessárias 12.480 entrevistas, pois cada nível seria avaliado por 80 pessoas, mas eu precisaria agora gerar os blocos para o experimento. Como posso fazer isto no R? Uma segunda questão, esta mais difícil na minha opinião, estes 156 níveis são divididas em 3 estratos: A, B e C. Algumas pessoas sugerem que estas entrevistas teriam de ter blocos do tipo: 2 níveis A, 2 níveis B e 2 níveis C. Existe algum modelo que permita isto, qual seria? É possível implementar este plano no R como no caso anterior? Antecipadamente agradeço []s Tura

2 1

Análise de MSE, PRESS e SSE
by Marcelo Claro de Souza 07 Out '11

07 Out '11

Pessoal preciso calcular MSE, PRESS e SSE de umas regressões lineares mas não consegui encontrar os comandos no R. Tentei instalar o pacote qpcR mas está dando erro na hora de instalar o rgl. Vocês poderiam me informar os comandos ou sugerir um pacote? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

pacote RCurl não instala
by Paulo Nogueira 06 Out '11

06 Out '11

Caros, estou tentando instalar o pacote RCurl mas está dando o seguinte erro: > install.packages('RCurl') Warning in install.packages : argument 'lib' is missing: using '/home/paulo/R/i686-pc-linux-gnu-library/2.11' tentando a URL 'http://www.vps.fmvz.usp.br/CRAN/src/contrib/RCurl_1.6-10.tar.gz' Content type 'application/x-gzip' length 709176 bytes (692 Kb) URL aberta ================================================== downloaded 692 Kb * installing *source* package ‘RCurl’ ... checking for curl-config... no Cannot find curl-config ERROR: configuration failed for package ‘RCurl’ * removing ‘/home/paulo/R/i686-pc-linux-gnu-library/2.11/RCurl’ * restoring previous ‘/home/paulo/R/i686-pc-linux-gnu-library/2.11/RCurl’ Warning in install.packages : installation of package 'RCurl' had non-zero exit status The downloaded packages are in ‘/tmp/RtmpxSea5o/downloaded_packages’ > Alguém sabe o que ocorre e como posso resolver? Outros pacotes eu consigo instalar normalmente. Em tempo. Tenho alguns códigos fonte em HTML e gostaria de desmembrar a sua estrutura para popular uma base de dados com informações de um site. Penso que o RCurl possa me ajudar. Existem ainda outras recomendações? Grato, Paulo

1 1

Re: [R-br] Como manter várias janelas gráficas ativas ao mesmo tempo
by Thiago De paula protásio 05 Out '11

05 Out '11

Pessoal; Acredito que a função windows()Atenda essa necessidade tb. Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Fwd: About relevent
by Fernando Colugnati 05 Out '11

05 Out '11

Pessoal, O autor de um artigo me enviou um pacote utilizado, mas ainda não disponível no CRAN por falta de documentação. Tentei construir e instalar o pacote utilizando o Rtolls. o RCMDda um feedback como ok, mas quando tento instalar ou montar o pacote ele tenta rodar um "sh" e falha na leitura do Description File do pacote. Detalhe, isso ocorre no Windows, no Ubunto funcionou, mas preciso rodar no windows tbm. Alguma dica de como fazer? Os arquivos seguem em anexo. Abs e obrigado ---------- Forwarded message ---------- From: Carter T. Butts <buttsc(a)uci.edu> Date: 2010/2/26 Subject: Re: About relevent To: Fernando Colugnati <fernando(a)ipti.org.br> Hi, Fernando. Just a quick FYI - I sent you the working development version, which I belatedly realized is still being worked on (and may have transient issues). Attached is the latest stable version...might be best to use that until the next one is finalized. Thanks, -Carter Fernando Colugnati wrote: > Dear Carter, > Thanks a lot! For sure I'll let you know about any bug, and by the way to > solve any hard "doubt" I cannot deal with... > > All the best > > 2010/2/25 Carter T. Butts <buttsc(a)uci.edu <mailto:buttsc@uci.edu>> > > > Hi, Fernando - > > > Fernando Colugnati wrote: > > Dear Dr Butts, > > My name is Fernando Colugnati, I am a Brazilian researcher at > the Research Institutte for Technology and Innovation, working > on analysis of WEB 2.0 tools. I am statistician. > > I found your paper "A relational event framework for social > action" really interesting, and want to model Discussion Forums > using that approach. I would like to ask you if it is possible > to use the R package cited, in the article, RELEVENT, to test in > some datasets I have, having my mais interest in compare forums > related to content interest and social relationship, both from a > distance course platform, having an intersection set of authors > that played in both forums. > The package will be used for research actvities only, and if you > need I can send you more details, or share results as the > research goes forward. > > Looking forward to hear from you. > > Best regards > > > > No problem. Here is the R source package; it's not yet on CRAN > because the documentation could use some work, but the basic > routines are documented. The two main functions are rem and > rem.ego, the former being designed for dyadic models of the sort > used in the paper, and the latter being originally intended for > "egocentric" event data (although, in point of fact, it can actually > be used to fit a wider range of models, given appropriate > specification of the model statistics). The main routines have been > tested a fair amount (in their core functionality, at least), but > this is still a beta package and under active development...please > let me know if you find any bugs. > > Hope this helps! > > -Carter > > > > > -- > Fernando A.B. Colugnati > Pesquisador Associado > > Instituto de Pesquisas em Tecnologia e Inovação - IPTI > Av. São Luís, 86 - cj 192 - 19º andar > Cep 01046 000 - São Paulo - SP > Tel. 55 11 3256.2150 > Fax. 55 11 3211.2041 > www.ipti.org.br <http://www.ipti.org.br> > fernando(a)ipti.org.br <mailto:fernando@ipti.org.br> > -- Fernando A.B. Colugnati Pesquisador Associado Instituto de Pesquisas em Tecnologia e Inovação - IPTI Tel. 55 11 8704-9812 www.ipti.org.br fernando(a)ipti.org.br

2 4

Forçar breaks histograma
by Alexandre Ribeiro Leichsenring 05 Out '11

05 Out '11

Caros, Quando se especifica os breaks do histograma (através do argumento breaks), o R o leva como uma sugestão, mas às vezes se rebela e não utiliza! Há alguma forma de forçar o histograma a adotar os breaks sugeridos? Curioso é que no help, ele afirma: breaksone of: - a vector giving the breakpoints between histogram cells, - a single number giving the number of cells for the histogram, - a character string naming an algorithm to compute the number of cells (see ‘Details’), - a function to compute the number of cells. *In the last three cases the number is a suggestion only. * deixando talvez implícito que um vetor com os breakpoints não seria apenas uma sugestão... Mas, mesmo especificando os breakpoints, ele às vezes não usa. Obrigado, Alexandre

3 7

Problema com o pacote "RMySQL".
by Juliana Freitas de Mello e Silva 05 Out '11

05 Out '11

Bom dia! Eu estou tentando usar o pacote "RMySQL" com o objetivo de criar um banco de dados a partir de outros (selecionando as colunas desejadas de cada um). Fiz da seguinte forma: dados <- dbGetQuery(con,paste("SELECT OBS_inmet.datahh,FCT_10km.fdata,FCT_10km.frange,OBS_inmet.newstn,FCT_10km.newstn,OBS_inmet.TempMax,FCT_10km.ftmax FROM OBS_inmet AS OBS_inmet INNER JOIN FCT_10km AS FCT_10km ON OBS_inmet.datahh=FCT_10km.fdata AND OBS_inmet.newstn=FCT_10km.newstn")) Porém "dados" só contém 29 observações, era para ser bem mais e eu não sei o motivo. Há alguma parte errada no comando? Outra dúvida: esse banco de dados é alimentado diariamente, como eu faço para selecionar datas específicas (como o mês de Agosto inteiro, pox exemplo)? P.S.: não sei se fui clara o bastante, mas desde já, agradeço!

2 3

Criando data frame
by Carlos Mendonça 05 Out '11

05 Out '11

Amigos, tenho um data frame com seq, unid e ref. Nele tenho n observações para cada seq. Eu gostaria de saber como posso fazer um "for" para gerar um data frame contendo todas as observações do 1º desses seq, depois um outro data contendo todas as observações do 2º e assim por diante. Eu teria no final tantos data frames quantos fossem os diferentes seq desse arquivo. Um exemplo para melhor compreensão. data frame arq1 seq unid ref 001 50 1000 001 150 5000 001 220 11000 001 500 1000 002 810 21500 002 450 44000 002 990 15000 o que eu quero é um código que me gere o arq2 e o arq3 com o seguinte conteúdo. arq2 seq unid ref 001 50 1000 001 150 5000 001 220 11000 001 500 1000 arq3 seq unid ref 002 810 21500 002 450 44000 002 990 15000 Obrigado, Carlos Mendonça.

3 5

Como manter várias janelas gráficas ativas ao mesmo tempo
by Soane Mota 04 Out '11

04 Out '11

Olá, criei uma função que tem que soltar 3 gráficos diferentes, em 3 janelas diferentes. A idéia inicial era que o usuário definisse qual/quais gráficos ele queria, porém, depois de muito tentar e sem muito sucesso, se eu conseguir colocar as 3 figuras em 3 janelas gráficas diferentes ao mesmo tempo já está legal. Alguém conhece algum comando que, incorporado a minha função, libere os 3 gráficos em suas respectivas janelas? (Nem é necessário que o usuário selecione o gráfico) Abraços Soane

3 2

gerar dados
by Ricardo Pereira 04 Out '11

04 Out '11

Outra duvida. Estou tentando simular dados, tal como um exemplo em outro programa... Exemplo abaixo: proc iml; do i=1 to 1000; seed=456; var=150; w = RANNOR(seed)*SQRT(var); end; *tentei:* set.seed(456) var <- 150 for (i in 1:1000) { w [i] <- rnorm * sqrt(varg) } *ERRO:* Erro em W[j] <- 1 : objeto 'W' não encontrado

10 15