Janeiro 2012 - R-br

Legenda eixo y
by Marcelo Claro de Souza 30 Jan '12

30 Jan '12

Alguém poderia me informar quais os comandos para aumentar o tamanho da fonte e deixar em negrito as legendas do eixo Y de um boxplot? Só encontrei essas informações para alterações no título. Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Re: [R-br] Importar Microdados - arquivos .txt sem separação entre as colunas
by Roney Fraga 29 Jan '12

29 Jan '12

Obrigado a todos pela ajuda. O caminho que segui foi esse http://sites.google.com/site/marcosfs2006/microdados e em seguida utilizar a função read.fwf(). Att -- Roney

1 0

Markov-switching GARCH (MS-GARCH)
by Rodrigo Coster 27 Jan '12

27 Jan '12

Povo, alguem tem, ou sabe onde encontrar, um código pra ajustar um MS GARCH? O que eu encontrei até agora foi um artigo dos autores do DEoptim ( http://www.jstatsoft.org/v40/i06/paper, aquele mesmo que eu pedi ajuda para gerar uma DLL), mas não tomei muita fé nos resultados que obtive e resolvi procurar por alternativas. []'s

1 0

arredondar soma para 100%
by Alexandro (Yahoo) 26 Jan '12

26 Jan '12

#Saudações a todos da lista R-br, fui fazer um gráfico com os dados abaixo #e a pessoa que o analisou, reclamou-me dizendo que a soma não foi 100%. #Como devo proceder para que a soma seja exatamente 100% ? Ou devo dar um print screen e mudar de 20,5% para 20,4% na categoria a? #Observe que a soma não fica 100% mesmo mudando digits para 2,3,... a=rep("a",9);b=rep("b",30);c=rep("c",5) dados=c(a,b,c);tabela=table(dados) porc=round(table(dados)/length(dados)*100,digits=1) graf=barplot(tabela,axes=FALSE,ylim=c(0,100));text(graf,porc,format(porc),pos=1);text(graf,porc," %",pos=1);text(graf,porc,format(tabela)) title(main = "Quantidade e Porcentagem das Categorias :a,b,c") #Obrigado! #Desde já agradeço! #Respeitosamente, #Alex

1 0

Re: [R-br] predict com minpack.lm
by Ivan Bezerra Allaman 26 Jan '12

26 Jan '12

Com ?nls.lm você resolveria seu problema. No entanto, segue um exemplo: f <-function(luz,a,b){ (((0.01940*luz+a+0.60375)-(((0.01940*luz+a+0.60375)^2)-4*0.01940*luz*b*(a+0.60375))^0.5)/2*b)-0.60375 } preditos <- do.call('f',c(list(luz=luz),lin1$par)) plot(luz,Pn,ylim=c(-1,5),type='l') par(new=TRUE) plot(luz,preditos,ylim=c(-1,5),col='red',xlab='',ylab='',type='l') (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}

1 0

predict com minpack.lm
by Gilson Sanchez 26 Jan '12

26 Jan '12

Pessoal mais uma vez na procura de ajuda... Estou utilizando o pacote minpack.lm para ajustar uma equação muito grande que com o nls não consegui fazer, agora estou querendo ver os valores ajustados pelo modelo, alguém pode ajudar... segue o CMR o qual estou trabalhando luz=c(0,25,50,75,100,500,1000,1500,2000) Pn=c(-0.60375, 0.04300,0.27000,0.90175,1.39225,2.70350,4.93100,4.31000, 3.80075) ps<-list(a=1,b=1) model1=function(ps,xx){(((0.01940*xx+ps$a+0.60375)-(((0.01940*xx+ps$a+0.60375)^2)-4*0.01940*xx*ps$b*(ps$a+0.60375))^0.5)/2*ps$b)-0.60375} resid1=function(observed,ps,xx){observed-model1(ps,xx)} lin1=nls.lm(ps,resid1,observed=Pn, xx=PAR) summary(lin1) predict???? Agradeço desde já! -- Bach. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 http://lattes.cnpq.br/3684553005529634 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.*** Salva un Árbol.* * *** * PIENSA**** EN EL MUNDO PARA TUS HIJOS** *

3 2

[OFF-TOPIC] Superfície de resposta
by Gustavo Marcatti 26 Jan '12

26 Jan '12

Pessoal, quais são os pressupostos para utilização de superfície de resposta? Por exemplo: as variáveis independentes podem ser aleatórias ou os valores devem ser previamente definidos pelo pesquisador, assim como num delineamento experimental. Desde já agradeço. Gustavo Marcatti Eng Florestal-UFV

2 2

Otimização de Código
by Leandro Marino 26 Jan '12

26 Jan '12

*Caros,* * * *alguém sabe me ajudar a otimizar o código do link abaixo:* * * http://pasteit.com/16646 Como referencia, coloquei o arquivo parit no link: http://pasteit.com/16647 Agradeco qualquer ajuda!!! Atenciosamente, Leandro Marino http://www.leandromarino.com.br (Fotógrafo) http://est.leandromarino.com.br/Blog (Estatístico) Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

2 3

Dúvida no uso do pacote "IBGEPesq"
by Savano Pereira 26 Jan '12

26 Jan '12

Prezados, Estou tento a seguinte mensagem de erro ao tentar carregar o pacote IBGEPesq: ------ Carregando pacotes exigidos: survey Attaching package: ‘survey’ The following object(s) are masked from ‘package:graphics’: dotchart Erro em library(pkg, character.only = TRUE) : package ‘IBGEPesq’ does not have a NAMESPACE and should be re-installed ------ Alguém já passou por algo semelhante? Como posso resolver este problema se já reinstalei o pacote? Desde já agradeço. Att, Savanno

3 7

Regressão Stepwise
by FHRB Toledo 26 Jan '12

26 Jan '12

Colegas da lista, Seguindo o livro MASS Capítulo 6, página 175 na sessão sobre escolha "automática" de modelos e recorrendo ao help da função stepAIC() fiquei em dúvida quanto ao argumento "k". Alguém tem alguma experiência nesse ponto? Não entendi direito o emprego desse argumento, no help diz que se "k" é igual a 2 efetivamente temos como critério o AIC, percebi que aumentando o "k" a resposta é um modelo mais parcimonioso, mas qual o critério para a sua escolha? No livro MASS os autores (V & R) eles fazem qualquer referência a isso... Alguém tem alguma referencia para a escolha do "k" com qualquer relação com escolha de preditores com efeito significativo pelo "F", por exemplo? att, FH

2 7