- R-br - DINF UFPR

[OFF-TOPIC] Reproduzindo mensagem de erro do R com Sweave!
by Ivan Bezerra Allaman 11 Out '11

11 Out '11

Boa noite senhores! Estou preparando um material básico para meus alunos, e gostaria de reproduzir algumas mensagens de erro do R no arquivo .tex. No entanto, algumas mensagens de erro não estou conseguindo reproduzir. Por exemplo, <<echo=TRUE,eval=FALSE>>= sqrt(-2) @ <<echo=FALSE,eval=TRUE>>= cat(try(sqrt(-2))) @ <<echo=TRUE,eval=FALSE>>= sqrt(2 ) @ <<echo=FALSE,eval=TRUE>>= cat(try(sqrt(2 ))) @ as mensagens acima, eu as reproduzo sem problema algum. No entanto, o trecho abaixo, <<echo=TRUE,eval=FALSE>>= squareroot(2) @ <<echo=FALSE,eval=TRUE>>= cat(try(squareroot(2))) @ é reproduzido, mais com um inconveniente: "Error in try(squareroot(2)) : não foi possível encontrar a função "squareroot", está aparecendo a mensagem "Error in try(squareroot(2))", não sei o porque! Já o seguinte trecho, não imprimi de forma alguma, impedindo a compilação do arquivo. <<echo=TRUE,eval=FALSE>>= sqrt 2 @ <<echo=FALSE,eval=TRUE>>= cat(try(sqrt 2)) @ Como reproduzir este erro, e impedir o inconveniente do penúltimo chuck? Desde já, agradeço a todos! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 2

Divirtam-se
by Benilton Carvalho 11 Out '11

11 Out '11

http://www.inside-r.org/howto/enter

1 0

comando weight lme
by Fernando Antonio de souza 11 Out '11

11 Out '11

Caros amigos Estou analisando dados em um modelo misto linear utilizando o pacote nlme contendo uma matrix de correlação do tipo compSymm e necessito utilizar como argumento "weights" o* inverso do quadrado do erro padrão* de cada média. Alguem sabe como faço isto. Sei que o argumento weigue utiliza algus objetos varfunction mas não compreendo muito bem seus objetos e não sei se algum deles fazem o que desejo. Alguém pode me ajudar? consumo<-lme(Fixed=Y~X, data=stpierre,random=X|Estudo) consumo<-update(consumo, random=compSymm(~X), Weight=?) Preciso usar o W2 como peso sendo: w2 =w1/ w sendo w1= o inverso do quadrado médio do erro w= média dos quadrados médios

2 1

Modelo de Efeitos Mistos: implementação
by Rodrigo Fadini 11 Out '11

11 Out '11

Prezados membros do R-br Estou auxiliando um aluno com as análises do seu mestrado, mas não tenho experiência com os Modelos de Efeitos Mistos. O aluno tenta avaliar mudanças na biomassa da floresta em função do tipo de floresta (manejada vs. não manejada) ao longo de 4 anos. Dentro de cada área (manejada e não manejada) ele possui várias parcelas de 1 hectare distribuídas sistematicamente, onde ele faz medidas nas mesmas árvores ao longo dos anos. Em cada parcela, ele divide as árvores em categorias (1-10 cm, 11 a 20 cm, e por aí vai). Eu sugeri a ele utilizar um Modelo Misto, por conta dessa estrutura hierárquica dos dados (parcelas dentro dos tipos de floresta) e também por conta das medidas repetidas às mesmas parcelas. O fato é que estou com problemas na implementação do modelo. Enfim, alguém poderia nos ajudar? Att., Rodrigo

2 2

Como salvar uma função sem precisar carregar
by Marcelo Cardoso mello 11 Out '11

11 Out '11

Como salvar a função abaixo de forma que não seja necessário carregar toda vez que quiser usar? Usar ela como se usa qualquer função do base. Obs.: No Windows. Obrigado por qualquer dica. x<-2 y<-1 f<-function(x,y){x*y+x} f(x,y)

2 1

[OFF-topic] muitos seminários sobre DOE, modelos mistos e relacionados (com vídeo)
by Walmes Zeviani 10 Out '11

10 Out '11

Pessoal, Acidentalmente tive a felicidade de encontrar seminários sobre diversos tópicos associados à delineamento experimental e modelos mistos. Quase todos possuem vídeo e slides. Vale muito a pena conferir e salvar aí na sua barra de favoritos. Aproveitem! http://www.newton.ac.uk/programmes/DAE/seminars/ À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Valores iniciais em hwe.cc - teste de verossimilhança para equilíbrio de HW
by Pedro Emmanuel Alvarenga Americano do Brasil 10 Out '11

10 Out '11

Amigos de R, Estou tentando entender valores iniciais que essa função pede. http://finzi.psych.upenn.edu/R/library/gap/html/hwe.cc.html Infelizmente o autor parece nao ter explicado bem o que significa esses valores iniciais e os exemplos tambem não ajudam muito. Eu tentei modificar progressivamente esses valores (de 0 a 100) e parece que os resultados se modificam mas nao modificaram a interpretação inicial. Mandei um email pro mesmo na semana passada que não retornou. Então estou apelando para amigos da lista que poderiam ja ter utilizado essa função e dar uma luz. Estou com receio de estar incorrendo em erro nao entendendo exatamente como esse valor modifica o resultado final, e quais valores seriam mais apropriados na minha situação. Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 1

Ajuste de regressão não-linear com bloco!
by Ivan Bezerra Allaman 09 Out '11

09 Out '11

Bom noite senhores! Estou analisando uma base de dados que se refere a crescimento de peixes avaliado em duas densidades. Em cada densidade, utilizou-se 2 tanques como bloco, sendo amostrado em torno de 50 peixes por taque. Os animais foram avaliados durante 12 meses, sendo amostrados mensalmente. Andei pesquisando na internet, mais não encontrei nenhum trabalho que tivesse ajustado um modelo não-linear em um delineamento em blocos. Alguém já viu algum trabalho, ou poderia me indicar algum caminho para tal tipo de análise? O ideal é que os parâmetros estimados para cada densidade fossem livre dos efeitos de bloco. Desde já muito obrigado pela atenção! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

Re: [R-br] [R_STAT] Teste de esfericidade
by Fernando Henrique Ferraz Pereira da Rosa 09 Out '11

09 Out '11

Michele, Acho que você estava se referindo a função sphericiy.test que eu implementei e fiz referẽncia em um material sobre ANOVA com medidas repetidas. Deixei o código disponível aqui: http://www.fernandohrosa.com.br/en/P/sphericity-test-for-covariance-matrice… Fernando, 2008/2/25 michele barbosa <mixxele(a)yahoo.com.br> > ** > > Como encontrar no help conteúdo para teste de esfericidade? > Algo referente ao comando: > *sphericity.test* > ** > obrigada > Michele > > ------------------------------ > Abra sua conta no Yahoo! Mail<http://br.rd.yahoo.com/mail/taglines/mail/*http://br.mail.yahoo.com/>, > o único sem limite de espaço para armazenamento! > > __._,_.___ > Mensagens neste tópico > <http://br.groups.yahoo.com/group/R_STAT/message/2683;_ylc=X3oDMTM1ZG1uYWw5B…>( > 1) Responder (através da web) > <http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJxdDB1bmdlBF9TAzk3N…>| Adicionar > um novo tópico > <http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmOXFwdmd2BF9TAzk3N…> > Mensagens<http://br.groups.yahoo.com/group/R_STAT/messages;_ylc=X3oDMTJmNW9qdjR2BF9TA…>| > Arquivos<http://br.groups.yahoo.com/group/R_STAT/files;_ylc=X3oDMTJnaG5zYzJzBF9TAzk3…>| > Fotos<http://br.groups.yahoo.com/group/R_STAT/photos;_ylc=X3oDMTJmbWRzZ3BjBF9TAzk…>| > Links<http://br.groups.yahoo.com/group/R_STAT/links;_ylc=X3oDMTJnamVkZWw1BF9TAzk3…>| Banco > de dados<http://br.groups.yahoo.com/group/R_STAT/database;_ylc=X3oDMTJkamg0ZWVqBF9TA…>| > Enquetes<http://br.groups.yahoo.com/group/R_STAT/polls;_ylc=X3oDMTJncjBzczVxBF9TAzk3…> > [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlZmFiNHFxBF9TAzk3NDkwNDM1BGdycElkAz…> > Alterar configurações via web<http://br.groups.yahoo.com/group/R_STAT/join;_ylc=X3oDMTJnaGRiMnBvBF9TAzk3N…>(Requer Yahoo! ID) > Alterar configurações via e-mail: Alterar recebimento para lista diária de > mensagens<R_STAT-digest@yahoogrupos.com.br?subject=+Recebimento+de+e-mail:+Lista+de+mensagens>| Alterar > formato para o tradicional<R_STAT-traditional@yahoogrupos.com.br?subject=Alterar+formato+de+distribui%C3%A7%C3%A3o:+Tradicional> > Visite seu Grupo > <http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJldjFsajNkBF9TAzk3NDkwND…>| Termos > de uso do Yahoo! Grupos <http://br.yahoo.com/info/utos.html> | Sair do > grupo <R_STAT-unsubscribe(a)yahoogrupos.com.br?subject=> > Atividade nos últimos dias > > - 1 > Novos usuários<http://br.groups.yahoo.com/group/R_STAT/members;_ylc=X3oDMTJnNzczN3RxBF9TAz…> > > Visite seu Grupo > <http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmanAxcGFzBF9TAzk3NDkwND…> > Yahoo! Mail > > Conecte-se ao mundo<http://us.lrd.yahoo.com/_ylc=X3oDMTJtcW92NGJmBF9TAzk3NDkwNDM3BF9wAzEEZ3JwSW…> > > Proteção anti-spam > > Muito mais espaço > Yahoo! Barra > > Instale grátis<http://us.lrd.yahoo.com/_ylc=X3oDMTJwczc4ZjY2BF9TAzk3NDkwNDM3BF9wAzIEZ3JwSW…> > > Buscar sites na web > > Checar seus e-mails . > Yahoo! Grupos > > Crie seu próprio grupo<http://br.groups.yahoo.com/;_ylc=X3oDMTJwcWI3bGZvBF9TAzk3NDkwNDM3BF9wAzMEZ3…> > > A melhor forma de comunicação > > . > > __,_._,___ > -- "Though this be randomness, yet there is structure in't." Fernando H Rosa - Statistician http://www.fernandohrosa.com.br / http://www.feferraz.net - Estatística, Matemática e Computação BankReview.com.br <http://www.bankreview.com.br/> - Escolha melhor seus serviços financeiros! AprendaAlemao.net <http://aprendaalemao.net/> - Seu ponto de partida para melhorar seu Alemão! @fhrosa

2 2

Função não reconhecida
by Mauro Sznelwar 09 Out '11

09 Out '11

Alguém sabe dizer o que está errado por aqui? set.seed(123) b<-5 sigma2<-2 n<-20 x<-runif(n,10,20) y<-b*x+rnorm(n,0,sqrt(sigma2)) plot(x,y) ############################################### #Función para estimar la distribución posterior f1<-function(beta0,nu,a,b,m,x,y) { A<-0 for(i in 1:n) { for(j in 1:n) { A<-A+x[j]*y[i]*(x[j]*y[i]-x[i]*y[j])/(sum(x^2)+nu) } } B0<-(nu*sum((y-beta0*x)^2))/(sum(x)+nu) med<-(sum(x*y)+nu*beta0)/(sum(x^2)+nu) dp<-sqrt((A+B0+b)/((n+a)*(sum(x^2)+nu))) beta<-rt(m,(n+a)/2)*dp+med gamma1<-(n+a)/2; delta1<-(A+B0+b)/2 sigma2<-1/rgamma(m,gamma1,delta1) w<-matrix(c(beta,sigma2),ncol=2) colnames(z)<-c("beta","sigma2") return(w) } ############################################### #Asignamos los valores para las ditribuciones # apriori y el número de simulaciones beta0<-4 nu<-0.00001 a<-1 b<-1 m<-1000 > simula<-f1(beta0,nu,a,b,m,x,y) Erro em colnames(z) <- c("beta", "sigma2") : objeto 'z' não encontrado > simula Erro: objeto 'simula' não encontrado

3 3

Solução
by Daniel Marcelino 08 Out '11

08 Out '11

Pessoal, tô com o seguinte problema: tenho eleições nos DISTRITO x, y, w, z, e vários eleitores para cada distrito desses. Na minha base tenho uma coluna que é o peso de cada voto em relação ao PARTIDO (A, B, C, D), e outras colunas com as opções de preferência de voto de cada eleitor: TURN1, TURN2, TURN3 etc. Eu preciso saber a soma total desse PESO em cada distrito, dividir isso por dois para obter meu cutting point. Depois, eu preciso calcular a soma da variável PESO para cada partido na coluna TURN1, ou seja, a soma do peso de quem votou para o PARTIDO: A, B, C, D, etc. Agora vem a parte difícil; eu preciso verificar se algum dos partidos obteve o valor do cutting point em cada DISTRITO. Se algum deles conseguiu, o distrito já é salvo no objeto final com a indicação de quem venceu. Caso nenhum partido tenha alcançado o valor do cutting point, é preciso eliminar o partido que obteve menos votos e verificar na segunda coluna TURN2 qual foi a segunda opção de voto do eleitor que votou no partido que acabo de eliminar; o o valor na linha da variável PESO do partido que acabo de eliminar é adicionado ao partido de segunda preferência. Agora, verificamos novamente se algum dos partidos restantes no distrito obteve o valor do cutting point, se sim, salvamos ele no objeto final, caso contrário, procedemos da mesma forma anterior até algum partido obter a maioria dos votos. Alguém tem algum ideia de como começar um script para rodar isso.

2 1

Re: [R-br] Análise de MSE, PRESS e SSE
by Thiago De paula protásio 07 Out '11

07 Out '11

Marcelo; Dê uma olhada na função influence.measures().Ela permite explorar mais algumas coisas da análise de regressão. Talvez isso tb te ajude.Att., 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

1 0

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 05 Out '11

05 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

Não reconhece variável
by Mauro Sznelwar 04 Out '11

04 Out '11

Alguém pode explicar por que a variável x não está sendo reconhecida neste caso? for (i in 1:30){ x[i]<--0.1+0.2*i }

7 9

Número da linha e da coluna
by Tiago Vieira 03 Out '11

03 Out '11

Boa noite a todos, Minha dúvida é a seguinte: Tenho uma matriz e quero saber qual o número da linha e da coluna que encontra-se a cela com o valor mínimo da matriz. Exemplo: > m1 <- matrix(1:12,3);m1 [,1] [,2] [,3] [,4] [1,] 1 4 7 10 [2,] 2 5 8 11 [3,] 3 6 9 12 > min(m1) [1] 1 Então quero saber a linha e a coluna onde o menor valor no "corpo" da matriz, que neste caso é 1, se encontra. É uma dúvida boba, mas desde já agradeço pela ajuda!!! Att, Tiago Vieira.

4 4

Re: [R-br] Número da linha e da coluna
by Tiago Vieira 03 Out '11

03 Out '11

Obrigado, Rodrigo! Este comando já me ajuda bastante, mas só por curiosidade, não existe outro que faça o mesmo procedimento em uma matriz e não em um vetor? ________________________________ De: Rodrigo Coster <rcoster(a)gmail.com> Para: Tiago Vieira <tiagovieira7(a)yahoo.com.br> Enviadas: Segunda-feira, 3 de Outubro de 2011 21:12 Assunto: Re: [R-br] Número da linha e da coluna ?which.min 2011/10/3 Tiago Vieira <tiagovieira7(a)yahoo.com.br> Boa noite a todos, > > >Minha dúvida é a seguinte: Tenho uma matriz e quero saber qual o número da linha e da coluna que encontra-se a cela com o valor mínimo da matriz. Exemplo: > > >> m1 <- matrix(1:12,3);m1 > [,1] [,2] [,3] [,4] >[1,] 1 4 7 10 >[2,] 2 5 8 11 >[3,] 3 6 9 12 >> min(m1) >[1] 1 > > > >Então quero saber a linha e a coluna onde o menor valor no "corpo" da matriz, que neste caso é 1, se encontra. É uma dúvida boba, mas desde já agradeço pela ajuda!!! > > >Att, > >Tiago Vieira. >_______________________________________________ >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

Foreach
by Gustavo Zeferino 03 Out '11

03 Out '11

Alguém sabe como utilizar a função *foreach*, para processamento paralelo? Estou tentando rodar código aqui, mas sempre aparece o erro Error in { : task 4 failed - "subscript out of bounds" No caso, o número da task é variável Código que estou rodando: registerDoMC(ncores) > results <- foreach(i = 1:10, .inorder=TRUE,.combine = rbind) %dopar% { > split_data <- splitData(data=data, split_rate=split_rate, > unique_users=unique_users) > training <- split_data$training > validation <- split_data$validation > error <- model(params_model, training, validation) > return(error) > } Todas as funções dentro do laço funcionam sem erros. O problema parece estar na hora de combinar o resultado final. Gustavo Zeferino

2 1

Re: [R-br] Mlp redes neurais (Fernando Neto)
by Roney Fraga 03 Out '11

03 Out '11

Fernando, talvez as dicas contidas neste endereço ajude. http://ecologia.ib.usp.br/labtrop/doku.php?id=labtrop:dicas:intalarunbutu Att -- Roney

3 3

Tempo de Processamento
by Ricardo Fernandes Campos Junior 03 Out '11

03 Out '11

Olá, pessoal. Estou usando a função pairwise.fst() do pacote adegenet. Eu a uso porque ela implementa uma forma de computar Fst independente da ploidia dos dados, diferentemente de Fst() implementado no pegas. No entanto, para meus propósitos, ela demora muito tempo. Gostaria de saber se algum de vocês conhece algum outro pacote que faça o mesmo cálculo mais rápido, alguma implementação em C para esta, ou alguma forma de diminuir o tempo de processamento significativamente. Desde já, muito obrigado. Abraços, Ricardo Fernandes Campos Junior Departamento de Genética e Biologia Evolutiva Universidade de São Paulo

3 3

Dúvida
by daiane brondani 02 Out '11

02 Out '11

Olá. Sou nova no R e estou com um problema. Estou fazendo um programa e tentando interpolar dados. Que função eu uso para fazer uma interpolação simples? Abraços

5 6

Mlp redes neurais
by Fernando Neto 02 Out '11

02 Out '11

Por favor, estou precisando usar a pacote que tenha a funcao mlp (multilayer perceptron) so q o pacote q achei como referencia esta dando Not Found. Alguem poderia me ajudar? Estou usando o R no ubuntu. Grato, Fernando

3 7

Daniele deixou uma mensagem para você...
by Badoo 01 Out '11

01 Out '11

Daniele deixou uma mensagem para você... Só você pode ler o conteúdo desta mensagem e ver quem a enviou. Delete a qualquer momento ou responda imediatamente com o sistema de troca de mensagens. Para descobrir o que diz a mensagem, siga este link. http://us1.badoo.com/01132477499/in/Ng0pKMMcVZc/?lang_id=61&m=63&mid=4e8702… Mais gente que espera pacientemente por você: Ivan (Ilhéus, Brasil) Lyly Nogueira (Ilhéus, Brasil) Wagner (Curitiba, Brasil) http://us1.badoo.com/01132477499/in/Ng0pKMMcVZc/?lang_id=61&m=63&mid=4e8702… Se o link desta mensagem não funcionar, copie e cole-o em seu navegador. Este email faz parte de nosso sistema de entrega para a mensagem enviada por Daniele. Se recebeu este email por engano, por favor ignore-o. A mensagem será deletada em breve. Divirta-se! A Equipe Badoo Você recebeu este email, porque um membro do Badoo deixou uma mensagem para você no Badoo. Esta é uma mensagem de somente envio. Respostas a esta mensagem não são monitoradas ou respondidas. Se você não quer receber mais mensagens do Badoo, por favor nos notifique: http://us1.badoo.com/impersonation.phtml?lang_id=61&mail_code=63&email=r-br…

1 0

Função não reconhece variável extraída de um modelo!
by Ivan Bezerra Allaman 30 Set '11

30 Set '11

Boa tarde senhores! Gostaria de saber como posso fazer uma função ler uma variável a partir de um modelo, nls por exemplo. Vamos supor um CMR. da <- expand.grid( A=factor( LETTERS[1: 3]), x=1:10) da$y <- 10/(1+exp(-( da$x-5)/1) )+rnorm(da$ x,0,0.1) n0 <- nls(y~As[A]/ (1+exp(-( x-x0[A])/ S[A])), data=da, start=list(As= c(10,10,10) , x0=c(5,5,5), S=c(1,1,1))) fun <- function(modelo){ mean(all.vars(summary(modelo)$formula)[1]) } Percebam que eu quero que a função calcule a média de uma variável, cuja variável foi declarada num determinado modelo. A função acima não dá certo, pois a função entende o seguinte: mean("VARIÁVEL") Para dar certo, eu acho, a variável não poderia estar entre aspas. Enfim, alguém tem alguma idéia de como eu posso fazer isso? Desde já agradeço a todos. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 4

Duvida sobre if e ifelse no R
by Augusto Ribas 30 Set '11

30 Set '11

Ola a todos, tenho uma duvida aqui. Estou lendo um livro sobre programação em R e empaquei num problema aqui. Eu tenha a função: x ≤ 0 ∈ (0, 1] > 1 f(x) −x3 x2 √x ai resolvo no intervalo de -2 a 2 com o comando if: #começo do script rm(list=ls()) # input x.values <- seq(-2, 2, by = 0.1) # for each x calculate y n <- length(x.values) y.values <- rep(0, n) for (i in 1:n) { x<-x.values[i] # your expression for y goes here if (x>1) { y<-sqrt(x) } else { if (x<=0) { y<- -x^3 } else { y<- x^2 } } #end of y expression y.values[i] <- y } # output plot(x.values, y.values, type = "l") #até aqui tudo funcionou beleza. #mas dai fui tentar fazer com ifelse a mesma coisa y.values2 <- rep(0, n) resultado<-ifelse(x.values>1,y.values2<-sqrt(x.values), ifelse(x.values<=0,y.values2<- -x.values^3, y.values2<- x.values^2) ) #o vetor resultado fica correto resultado==y.values #mas o vetor y.values2 fica errado y.values2==y.values #o resultado ficou [1] 8.000000 6.859000 5.832000 4.913000 4.096000 3.375000 2.744000 2.197000 [9] 1.728000 1.331000 1.000000 0.729000 0.512000 0.343000 0.216000 0.125000 [17] 0.064000 0.027000 0.008000 0.001000 0.000000 0.010000 0.040000 0.090000 [25] 0.160000 0.250000 0.360000 0.490000 0.640000 0.810000 1.000000 1.048809 1.303840 1.341641 1.378405 [41] 1.414214 #mas os valores que ficam salvos são no vetor y.values2 são: > y.values2 [1] 4.00 3.61 3.24 2.89 2.56 2.25 1.96 1.69 1.44 1.21 1.00 0.81 0.64 0.49 0.36 [16] 0.25 0.16 0.09 0.04 0.01 0.00 0.01 0.04 0.09 0.16 0.25 0.36 0.49 0.64 0.81 [31] 1.00 1.21 1.44 1.69 1.96 2.25 2.56 2.89 3.24 3.61 4.00 #minha duvida então é pq ele so resolve tudo elevando ao quadrado no y.values2 #o que acontece, por exemplo, com o primeiro valor de resultado que vejo 8 qd mas fica salvo 4 quando vou olhar o y.values2. #alguém pode me iluminar sobre o que esta acontecendo? -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

4 5

Distância entre strings
by Pedro Rafael 29 Set '11

29 Set '11

Pessoal qual o objeto para compara duas strings e observar a distância entre elas? Exemplo: Gostaria que Pedro Rafael fosse igual à peDro Rafael, que também fosse igual à Rafael Pedro. Existe alguma forma de fazer isto? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

7 18

DIFFTIME
by Edson Lira 29 Set '11

29 Set '11

Gente, estou reenviando este email, tenho os dados abaixo: DIA SAIDA CHEGADA 1 09/09/2011 08:11 11:20 2 09/09/2011 08:17 08:25 3 09/09/2011 08:45 09:05 4 09/09/2011 09:15 09:20 5 09/09/2011 09:30 09:35 6 09/09/2011 09:50 10:15 7 09/09/2011 10:40 10:55 8 12/09/2011 08:25 11:05 9 12/09/2011 08:30 08:45 10 12/09/2011 09:25 09:40 Todas as variáveis tem a estrutura de fatores. Vejam que as datas (DIA), estão separadas da hora da saida (SAIDA) e da hora da chegada (CHEGADA). Preciso calcular a diferença em minutos entre a hora da chegada e a hora da saida (CHEGADA-SAIDA). O Augusto me respondeu sugerindo usar o help para: ?as.POSIXct ?strptime ?difftime Verifiquei, mas não consegui, sempre tenho dificuldades em trabalhar com data. O pacote chron tem conversões, calcula diferenças, mas não consegui implementar seus comandos para o meu problema acima. []'s. Edson Lira Estatístico Manaus-Amazonas

2 1

conversão coordenadas decimais em cartesianas
by Heloíse Pavanato 29 Set '11

29 Set '11

Olá, Alguém sabe se existe um pacote para converter coordenadas geográficas em cartesianas?? Obrigada, Heloise

2 2

Problema com matrizes e vetores
by Thiago Touya 28 Set '11

28 Set '11

Senhores, preciso estimar um vetor de parâmetros (Beta) usando as saídas do svd() para comparar com a estimação pelo método dos mínimos quadrados. O grande problema é que nao sei isolar Beta(3x1) na seguinte equação: H'XBeta = H'Y = lambda'Beta, onde lambda'(1x3) Eu tenho H'(1x3), X(3x3) e Y(3x1), mas nao estou visualizando uma maneira de isolar beta, pois H'X é um vetor. X <- matrix(c(3,0,-3,0,6,3,-3,3,6),ncol=3,nrow=3) dvs <- svd(X) lambda <- (dvs$d)^2 Y <- c(5,2,3) D <- matrix(c(lambdas[1,1],0,0,0,lambdas[2,1],0,0,0,lambdas[3,1]), ncol=3, nrow=3) C <- dvs$u Vlinha <- t(dvs$v) sabe-se que pela Decomposição em valores singulares X=CDV' (neste caso em que X é simétrica e nao temos nenhum autovalor igual a zero) Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

1 0

Dúvida simples
by Victor Eduardo 27 Set '11

27 Set '11

Gente, como posso colocar o mi linha no R com o comando expression ??? Abraços!

3 4

Conexão R com Oracle
by nani.design＠bol.com.br 27 Set '11

27 Set '11

2 1

Ponto de corte em cluster
by Marcelo Cardoso mello 27 Set '11

27 Set '11

Como obter no R um ponto de corte em um cluster (UPGMA ... etc). Algum critério como pseudo test F ou t ou algum outro? Obrigado. M

2 1

Re: [R-br] Gráfico WMF
by Ivan Bezerra Allaman 27 Set '11

27 Set '11

Bom, então um caminho extremamente deselegante seria diminuir a fonte do eixo, o que acha? plot(1:40,runif(40),axes=F,xlab='Eixo X',ylab='Eixo Y') axis(1,at=c(0,seq(2,40,3),40), cex.axis=0.8) axis(2,seq(0,1,.1)) box() savePlot('c:/temp/temp.wmf') dev.off() (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

9 12

Gráficos 3D para exibir Simulações _
by Paulo Nogueira 27 Set '11

27 Set '11

Olá a todos. Acompanho a lista há bastante tempo e mesmo sem mandar mensagens aprendo muito com as discussões e soluções propostas. Acho que quase tudo que sei de R aprendi por aqui. Sou bastante grato. Minha questão é a seguinte: Estou simulando distribuições sendo que cada uma delas é gerada modificando-se de forma controlada o valor de alguns parâmetros da função geradora. Para compará-las graficamente, exibindo a evolução conforme modificam-se os parâmetros, gostaria de criar um gráfico semelhante ao que coloquei no link abaixo. http://goo.gl/17Erk Não conheço muito de gráficos 3D no R, mas gostaria que vocês indicassem funções/pacotes que eu possa pesquisar e chegar em um resultado semelhante. Toda ajuda é bem vinda. Grato, Paulo Nogueira Starzynski

2 4

Erro: pacote não encontrado
by Pedro Emmanuel Alvarenga Americano do Brasil 27 Set '11

27 Set '11

Amigos de R, O seguinte erro vem ocorrendo diversas vezes comigo com diferentes pacotes... > library('epicalc')Carregando pacotes exigidos: MASSCarregando pacotes exigidos: nnet Anexando pacote: 'epicalc' The following object(s) are masked from 'package:rms': lrtest The following object(s) are masked from 'package:Hmisc': fillin, recode > epicalc::followup.plot(c$prontuario,c$dias,c$CD,by=c$iris)Erro em get(search()[2]) : objeto 'package:epicalc' não encontrado> followup.plot(c$prontuario,c$dias,c$CD,by=c$iris)Erro em get(search()[2]) : objeto 'package:epicalc' não encontrado> search() [1] ".GlobalEnv" "package:epicalc" "package:nnet" "package:MASS" "package:foreign" [6] "package:rms" "package:Hmisc" "package:survival" "package:splines" "tools:rstudio" [11] "package:stats" "package:graphics" "package:grDevices" "package:utils" "package:datasets" [16] "package:methods" "Autoloads" "package:base" > Coincidentemente esse erro irritante que acontece com mais de um pacote vem acontecendo depois que eu comecei a usar o RStudio. Eu achei inicialmente porque eu tinha um monte de pacotes carregados, depois eu achei que era porque eu importei os dados com spss.get e não com o read.spss. Mas nada disse fez diferença nas tentativas e erros. Alguem tem uma luz pra eu tentar evitar esse erro? , Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 10

Re: [R-br] Gráfico WMF
by Ivan Bezerra Allaman 26 Set '11

26 Set '11

Leandro, Se você pediu uma sequência de 2 à 40 de 3 em 3, o último número que irá aparecer é claro que será o 38, por isso o número 40 não aparece. Caso queira que ele apareça, das diversas opções, você pode tentar, seq(2,41,3) ou seq(2,40,2). (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 5

Gráfico WMF
by Leandro Marino 26 Set '11

26 Set '11

Caros, Preciso de uma ajuda. Não sei como resolver isto. Nas configurações abaixo, necessito que o número 40 seja exibido no eixo X. Alguma sugestão? plot(1:40,runif(40),axes=F,xlab='Eixo X',ylab='Eixo Y') axis(1,at=c(0,seq(2,40,3),40)) axis(2,seq(0,1,.1)) box() savePlot('c:/temp/temp.wmf') dev.off() Aproveitando, existe outra função para salvar em WMF (obrigatoriamente) mas que não seja necessário abrir uma nova janela? Isto faz com que seja consumido muito mais memória e processamento! Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 ( (21) 2103-9600 R.:236 * <mailto:leandro@cesgranrio.org.br> leandro(a)cesgranrio.org.br P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

2 1

Problemas com vetores
by Thiago Touya 26 Set '11

26 Set '11

Senhores, em meu trabalho, criei 3 vetores, sendo dois deles formados por 0's e 1's. De posse dos três criei uma matriz. Agora estou tentando adicionar mais uma coluna à matriz. Coluna esta que servira para identificar a situação dos dados do outro vetor que foi criado. No fim das contar preciso ter uma matriz com 4 colunas sendo que as duas primeiras formadas por 0's e 1's determinam a condição da terceira coluna e essa condição será encontrada na quarta coluna. Só que a quarta coluna é formada por valores to tipo "Censurado", "Curado" e "Falha". Tentei fazer da seguinte forma: # 1º vetor VetCens <- rbinom(500,1,0.6) # 2º vetor VetCura <- rbinom(500,1,0.3) # 3º vetor VetBN <- rnbinom(500,48,mu=13.5) # Criando a matriz MatSim <- cbind(VetCens, VetCura, VetBN) # Tentativa de criação do quarto vetor com as situações VetSituac <- function(param){ VetCens <- param[1] VetCura <- param[2] situacao <- rep(0,500) for(i in 1:500){ if (VetCens[i]==1 and VetCura[i]==0) situacao[i]==3 } if (VetCens==1 and VetCura==1){ situacao="Curado" } if (VetCens==0 and VetCura==1) or (VetCens==0 and VetCura==0) situacao="Falha" } } Estou tendo problemas com essa ultima função. Pra um individuo ser considerado CURADO devemos ter para os dois primeiros vetores os valores 1,1 Pra um indivíduo ser considerado CENSURADO devemos ter para os dois primeiros vetores os valores 1,0 Pra um indivíduo ser considerado FALHA devemos ter para os dois primeiros vetores os valores 0,1 ou 0,0. Alguém tem uma sugestão? Obrigado -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

6 6

Off toppic - R Windows
by Leonard de Assis 26 Set '11

26 Set '11

Obrigado Norton 360 por detectar o r.dll como Virus, deletar da máquina sem me avisar e fazer eu perder um fim de semana tentando descobrir o que tinha de errado no meu R Só pra desabafar mesmo :( -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

2 1

Compilando e utilizando códigos em C no R
by Gustavo Marcatti 25 Set '11

25 Set '11

Pessoal esse email é para aqueles que tem dificuldades em compilar um código em C e utilizá-lo no R (principalmente usuários Windows). Tem um pacote (inline) que facilita bastante esse processo e pode ser interessante até para usuários unix. Esse pacote compila o código gerando uma .dll (windows) temporária que depois é carregada no R. O código em C pode ser escrito no próprio script do R (normalmente eu deixo em um arquivo separado .r e quando quero utilizar a função utilizo o source(".r"), então o código é compilado/carregado e a função já fica pronta para uso. Depois de encerrar o R o dll é excluída do sistema. Segue um exemplo. # Exemplo: # Algoritmo que retorna o valor mínimo de um vetor # 1º: utilizando o loop for do R min2 <- function(vetor){ n <- length(vetor) Min <- vetor[1] for(i in 1:(n-1)){ if(Min > vetor[i+1]) Min <- vetor[i+1] } Min } # 2º: utilizando um codigo escrito em C #install.packages("inline", dep = T) require(inline) #--- Compilando a função que encotra o menor valor em um vetor ------------ #--- Só é necessário compilar uma vez -------------------------------------- # Argumentos da função em C sigMin <- signature( vetor="double", n="integer", Min = "double") # Código propriamente dito codeMin <- " *Min = vetor[0]; for (int i=0; i < *n - 1; i++) { if(*Min > vetor[i+1]) *Min = vetor[i+1]; }" # Para compilar a função funcMin <- cfunction(sigMin, codeMin, convention=".C") #--- Fim # apenas para facilitar o uso da função min3 <- function(vetor){ n <- length(vetor) Min <- 0 funcMin(vetor,n,Min)$Min } # Teste: set.seed(1) vec <- runif(1000000) system.time(minimo1 <- min2(vec)) minimo1 system.time(minimo2 <- min3(vec)) minimo2 # Mais exemplos podem ser encotrados no pacote inline. Gustavo Marcatti Eng. Florestal- UFV

1 0

Re: [R-br] Problemas com vetores
by Thiago Touya 25 Set '11

25 Set '11

Senhores, muito obrigado pela ajuda. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

1 0

Offset da função lm
by Gustavo Marcatti 24 Set '11

24 Set '11

Bom dia a todos! Alguém sabe como a função lm estima os parâmetros de um modelo com um offset presente? Como ficaria o estimador dos parâmetros? Desde já obrigado. Gustavo Marcatti Engeharia Florestal - UFV

2 4