Março 2011 - R-br

Função para quebra de linha em título
by Fabio Mathias Corrêa 22 Mar '11

22 Mar '11

Fábio Mathias Corrêa Departamento de Estatística Universidade Estadual de Santa Cruz Tel.: 73-3680-5076 Cel.: 73-9991-8155 ----- Mensagem encaminhada ---- De: Marcio Nicolau <marcio.nicolau(a)gmail.com> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Sábado, 25 de Outubro de 2008 19:47:08 Assunto: [R_STAT] Função para quebra de linha em título Prezados, A algum tempo, questionei na lista sobre a existência de uma função para quebra de linha automática para títulos de gráficos. Após a sugestão do Benilton (função gsub), segue abaixo o código... caso necessitem. ## Pre-Processa titulos pretitulo <- function(orig, cut=55) { nparts = ceiling(nchar( orig) / cut) cuts = c(1,1:nparts * cut) spt = strsplit(orig, "\\s")[[1] ] cs = cumsum(nchar( spt)) ids <- matrix(FALSE, nparts,length( nchar(spt) )) for(i in 2:length(cuts) ) { ids[i-1,] <- (cs >= cuts[i-1] & cs < cuts[i]) } ch <- vector("character" , nparts) for(i in 1:nparts) { ch[[i]] <- paste(spt[ids[ i,]],collapse = " ") } pretitulo <- paste(ch,collapse= "\n") } Não é o melhor, computacionalmente falando, mas funciona :) Abraços Marcio Nicolau marcio.nicolau@ gmail.com __._,_.___ Mensagens neste tópico (1) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 8 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

6 6

Transformando dados com 2 fatores em uma tabela
by Augusto Ribas 22 Mar '11

22 Mar '11

Oi galera. Minha primeira postagem na nova lista que emoção :) Gostaria de ajuda com um problema simples que estou tendo aqui. #Eu queria fazer um plot assim: barplot(VADeaths, beside = TRUE, col = c("lightblue", "mistyrose", "lightcyan", "lavender", "cornsilk"), legend = rownames(VADeaths), ylim = c(0, 100)) title(main = "Death Rates in Virginia", font.main = 4) #Mas meu dados estão na forma de um data frame dessa forma... exemplo<-dados<-data.frame(cbind( Tratamento=rep(c("Controle","Ac10","Ac20","Estresse"),each=3) ,Placa=rep(1:3,4) ,Germinacao=round(rnorm(12,5,3),2) )) #Tem alguma maneira simples de fazer meu dataframe de exemplo virar uma tabela como a do "VADeaths" que é usado no help do barplot #pra fazer um plot daquele jeito, separando grupos de barrinhas tanto por tratamento como placa, no #meu caso, ou então fazer um barplot do jeito que quero(o do "VADeaths") direto com um dataframe assim do meu jeito? Obrigado -- Grato Augusto C. A. Ribas

2 2

Res: [R_STAT] Regressão de proporções com erro normal ou binomial
by Walmes Zeviani 21 Mar '11

21 Mar '11

Guilherme, Na regressão logística você observa significância porque o parâmetro de dispersão não é estimado mas sim conhecido/fixado como 1 (note que é aplicado o teste z). No modelo normal, além da estimação dos efeitos/coeficientes/parâmetros, é estimado a variância do erro aleatório e os erros padrões das estimativas estão em função disso (teste t). Se você usar uma family=quasibinomial você notará diferença. Esperamos você na [R-br]. À 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

Estudo de função!
by Elton Gean Araújo Araújo 21 Mar '11

21 Mar '11

Boa tarde! Gostaria de aplicar os conceitos de crescimento, decrescimento, máximos e mínimos de funções no R, poderiam me informar o pacote para isso? Desde já agradeço pela informação e pela nova lista! Abraço. Elton.

1 0

Especificação do modelo em lm e glm
by Diego Pujoni 21 Mar '11

21 Mar '11

#Olá pessoal, #Existem formas diferentes de se especificar um modelo na função #"lm" e "glm". Mas obtive os seguintes resultados: #Criando o data frame# X1=c(rep("A",4),rep("B",4)) X2=c(rep("C",2),rep("D",2),rep("C",2),rep("D",2)) Y=c(3,5,5,7,11,13,13,15) Table=as.data.frame(cbind(X1,X2,as.numeric(Y))) Table #Os quatro modos de se especificar o mesmo modelo fornecem a mesma deviance,# #ou seja, são iguais# glm(Y ~ X1+X2+X1:X2, data=Table, family=poisson())$deviance glm(Y ~ (X1+X2)^2, data=Table, family=poisson())$deviance glm(Y ~ X1*X2, data=Table, family=poisson())$deviance glm(Y ~ X1:X2, data=Table, family=poisson())$deviance #Mas há um deslocamento de -1.25 nos coeficientes no último modelo round(glm(Y ~ X1+X2+X1:X2, data=Table,family=poisson())$coefficients,digits=2) round(glm(Y ~ (X1+X2)^2, data=Table,family=poisson())$coefficients,digits=2) round(glm(Y ~ X1*X2, data=Table,family=poisson())$coefficients,digits=2) round(glm(Y ~ X1:X2, data=Table,family=poisson())$coefficients,digits=2) #O coeficiente de X1B:X2D sai como NA# #O quarto modo de especificar o modelo está errado? #Por que isso acontece? #Nos modelos abaixo: #Incluindo mas uma variável preditora em Table# X3=c("E","F","E","F","E","F","E","F") Table.X3=as.data.frame(cbind(X1,X2,X3,as.numeric(Y))) Table.X3 #As deviances e os coeficientes X1A:X3F e X1A:X3E também são iguais #além disso alguns coeficientes são deslocados de -1.4 (X1A:X2C,X1B:X2C,X1A:X1D), #mas os outros coeficientes não são iguais# glm(Y~X1+X2+X3+X1:X2+X1:X3,data=Table.X3,family=poisson())$deviance glm(Y~X1:X2+X1:X3,data=Table.X3,family=poisson())$deviance round(glm(Y~X1+X2+X3+X1:X2+X1:X3,data=Table.X3,family=poisson())$coefficients,digits=2) round(glm(Y~X1:X2+X1:X3,data=Table.X3,family=poisson())$coefficients,digits=2) #Especificando os coeficientes do primeiro modelo # Intercept = 1.16 # X1/X2| C D | X1/X3| E F | # A | 0 0.41| A | 0 0.41| # B |1.24 -0.25| B |1.24 -0.25| #Especificando os coeficientes do segundo modelo # Intercept = 2.56 # X1/X2| C D | X1/X3| E F | # A |-1.40 -0.99| A | 0 0.41| # B |-0.15 NA | B | 0 0.15| #Os coeficientes X1B:X3E e X1B:X3F mudam. #Alguém sabe aonde se encontra o meu erro? #Será que sou obrigado a incluir o termo de interação X2:X3? #Mas e se essas duas variáveis forem independentes? #Muito Obrigado # Diego PJ

2 1

Solução de Equação não Linear
by Carolina Pena 21 Mar '11

21 Mar '11

Boa tarde! Alguém sabe como consigo resolver uma equação não linear no R ? Por exemplo suponha que eu tenho a seguinte equação: y = 2,5*x + 3,6*x^2 + 7,8*x^4+ 10,0*x^12 Para um dado valor de y (exemplo y = 5 ) qual é o valor de x ? Obrigada, Carolina

4 4

Re: [R-br] [R_STAT] [Dúvida] Carregar dados sem ser na memória RAM
by Fernando Mayer 20 Mar '11

20 Mar '11

O pacote SOAR pode resolver. De uma olhada no manual e na vignette em http://cran-r.c3sl.ufpr.br/web/packages/SOAR/index.html []s, --- Fernando Mayer URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com 2011/3/20 pedrorafaelmarinho <pedrorafaelmarinho(a)yahoo.com.br>: > Pessoal estou com um problema ou talvez seja do R. Tenho que ler uma base de dados de 600 mb em formato DBF (DBase). O R pelo que ando lendo ele ler a base de dados e joga na memória RAM. Toda vez que carrego essa base de dados o pc trava a ponto de o mouse não funcionar. Existe alguma forma de fazer com que o R não carregue toda a base de dados diretamente na memória? Por exemplo a base de dados ficar no HD e somente carregar algumas variáveis a medida que vou precisando? Softwares simples como o Tabwin eu consigo acessar a base de dados e o pc não trava...O que está acontecendo? > O comando que usei para carregar a base foi read.dbf() usando o pacote foreign. > > > > ------------------------------------ > > Links do Yahoo! Grupos > > <*> Para visitar o site do seu grupo na web, acesse: > http://br.groups.yahoo.com/group/R_STAT/ > > <*> Para sair deste grupo, envie um e-mail para: > R_STAT-unsubscribe(a)yahoogrupos.com.br > > <*> O uso que você faz do Yahoo! Grupos está sujeito aos: > http://br.yahoo.com/info/utos.html > > >

3 4

sequencia de minutos
by Fabio Mathias Corrêa 20 Mar '11

20 Mar '11

----- Mensagem encaminhada ---- De: Henrique Dallazuanna <wwwhsd(a)gmail.com> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Terça-feira, 10 de Fevereiro de 2009 15:40:13 Assunto: Re: [R_STAT] sequencia de minutos Tente assim: 1) minutos <- format(seq(ISOdate( 2008, 2, 1, 0, 0, 0), ISOdate(2008, 2, 1, 23, 59, 59), by = "10 mins"), "%H:%M:%S") 2) use a função axis para adicionar no plot: plot(...., xaxt = 'n') axis(1, at = 1:length(minutos) , labels = minutos) 2009/2/10 Vinicius Brito Rocha <viniciusbritor@ gmail.com> preciso criar uma sequncia de tempo em intervalos de 10 em 10 minutos. > >ex: x=c(00:00, 00:10, 00:20,...,23: 50) > >quero usar isso como um eixo em um plot ou um ts.plot. > >alguem sabe como fazer 1º, para criar a sequencia e mais importante, usa-la como >eixo no gráfico. > >Att. > >-- >Vinicius Brito Rocha > > > -- Henrique Dallazuanna Curitiba-Paraná -Brasil 25° 25' 40" S 49° 16' 22" O __._,_.___ Mensagens neste tópico (2) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 4 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

2 1

Re: [R-br] Solução de Equação não Linear
by Leandro Souza 18 Mar '11

18 Mar '11

Carolina, Uma outra sugestão: fx<-function(x) (2.5*x+3.6*x^2+7.8*x^4+10*x^12) fy <- function(y) abs(fx(y)-5) optimise(fy,c(-100,100)) Abraços ________________________________ De: Daniel C Bezerra <danielcbezerra(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Sexta-feira, 18 de Março de 2011 15:44 Assunto: Re: [R-br] Solução de Equação não Linear Uma sugestão seria usar o método de Newton-Raphson: fx<-function(x) (2.5*x+3.6*x^2+7.8*x^4+10*x^12) library(animation) newton.method(fx) 2011/3/18 Carolina Pena <carolinapenaa(a)yahoo.com.br> Boa tarde! > >Alguém sabe como consigo resolver uma equação não linear no R ? > >Por exemplo suponha que eu tenho a seguinte equação: > >y = 2,5*x + 3,6*x^2 + 7,8*x^4+ 10,0*x^12 > >Para um dado valor de y (exemplo y = 5 ) qual é o valor de x ? > >Obrigada, > >Carolina > > > >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

1 0

Regressão com Group by
by Vinicius Brito Rocha 18 Mar '11

18 Mar '11

senhores, que ajustar diversos modelos via lm, usando uma espécis de Group by, pois preciso de resultados individuais para cada grupo, por isso não posso usar uma variável dummy. Preciso de um modelo para cada grupo. Alguém pode me dar uma luz? Abs. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

4 4