Julho 2011 - R-br - DINF UFPR

[Dúvida] Rodando Script R no terminal sem entrar no prompt.
by Pedro Rafael 22 Jul '11

22 Jul '11

Pessoal existe alguma forma de rodar o script do R no terminal Linux sem ser necessário "entar no R"? Tipo assim, referenciar o código do R no diretório que ele se encontra e mandar ele rodar sem ser necessário entrar no prompt de comando do R. Obrigado, Pedro Rafael

4 3

Problema em regressão segmentada!!
by Ivan Bezerra Allaman 22 Jul '11

22 Jul '11

Boa tarde senhores!! Estou tendo dificuldades ao ajustar uma regressão segmentada com a library ''segmented'. Não entendo qual a dificuldade que a função está encontrando para ajustar uma regressão aos meus dados. Creio que os breaks pontos chutados não estão tão fora da realidade. Segue abaixo um CMR real. Y <- c(18.335, 19.456, 19.383, 19.542, 20.718, 21.531) X <- c(0,3,6,9,12,15) plot(X,Y,type='b') mod <- lm(Y ~ X) library(segmented) mod.seg <- segmented(mod, seg.Z=~X, psi=list(X=c(3,8))) Se alguém puder ajudar, agradeço. (S,f,P) Allaman Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

alterar cores apenas na area grafico
by jose coelho 22 Jul '11

22 Jul '11

Vivam a todos Gostaria de saber com alterar a cor apenas da área do gráfico feita com plot (x,y). já tentei com plot(x,y,bg="gray) não dá erro mas também não funciona!! ao usar par(bg = "gray") fica toda a janela a cinzento e não apenas a área do gráfico!! Como posso fazer isso? Muito obrigado -- jmmcoelho

2 1

Re: [R-br] xtabs ou equivalente
by Jose Claudio Faria 22 Jul '11

22 Jul '11

Test! -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

PLOT SERIE TEMPORAL
by Enéias Façanha 22 Jul '11

22 Jul '11

PessoALL, Utilizando os dados abaixo, como eu faço para plotar essas duas séries, "MILLE_META" e "MILE_RECEITA", no mesmo gráfico, inserindo os meses no eixo x? MILLE_META MILLE_RECEITA JAN 339.155,42 275.649,07 FEV 320.427,98 267.601,66 MAR 333.703,99 279.127,36 ABR 318.293,90 184.393,12 MAI 342.795,95 262.414,92 JUN 328.074,05 268.703,78 Obrigado, -- Enéias Façanha

2 1

Minimização Condicionada
by Luciano Ramos Gonçalves 22 Jul '11

22 Jul '11

Gostaria de saber qual é a melhor função para realizar uma minimização com restrições: nlm(), optim(), optimize() ou constrOptim(), ou outra? Abaixo vai a função objeto, e as restrições descritas. Qualquer dica já ajuda. O meu problema maior foi aplicar as restrições aos métodos de otimização. Obs.: Salvo engano, tudo abaixo é plenamente reproduzível. # Tenho a função abaixo. LogLik <- function(parametros,dados){ alpha <- parametros[1] delta <- parametros[2] epsilon.c <- parametros[3] epsilon.v <- parametros[4] mu <- parametros[5] B <- dados[,2] S <- dados[,3] logl <- sum(-epsilon.c-epsilon.v+((max(c(B,S))+min(c(B,S)))/2)* (log((epsilon.v/(mu+epsilon.v)))+ log((epsilon.c/(mu+epsilon.c))))+ B*log(mu+epsilon.c)+S*log(mu+epsilon.v))+ sum( log( (alpha)*(1-delta)*exp(-mu)* ((((epsilon.v/(mu+epsilon.v))))^ (S-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.c/(mu+epsilon.c))^ (-((max(c(B,S))+min(c(B,S)))/2)))+ (alpha)*(delta)*exp(-mu)* ((((epsilon.c/(mu+epsilon.c))))^ (B-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.v/(mu+epsilon.v))^ (-((max(c(B,S))+min(c(B,S)))/2)))+ (1-alpha)* ((epsilon.c/(mu+epsilon.c))^ (B-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.v/(mu+epsilon.v))^ (S-((max(c(B,S))+min(c(B,S)))/2))) )) return(logl) } # Supondo os valores BS <- cbind(1,500,501) theta <- c(0.4,0.5,40,40,50) # Quero minimizar a função: LogLik(theta,BS) # Com as seguintes restrições, que eu não consegui implementar (as restrições são # aos valores de theta, ou "parametros" na função, em ordem): # 0 <= alpha <= 1 ("alpha" é um valor entre 0 e 1) # 0 <= delta <= 1 ("delta" é um valor entre 0 e 1) # epsilon.c > 20 ("epsilon.c" é maior que 20) # epsilon.v > 20 ("epsilon.v" é maior que 20) # mu > 30 ("mu" é maior que 30) Grato, Luciano.

1 0

[Dúvida] Salvando arquivo CSV sem aspas nas strings
by Pedro Rafael 22 Jul '11

22 Jul '11

Pessoal estou salvando um data.frame com o objeto write.table(objeto, file = "nome_do_arquivo.csv", sep = ";"). Contudo quero que as strings não contenham aspas. Como fazer para que o csv não contenham nas strings? Saudações, Pedro

3 2

xtabs ou equivalente
by Vinicius Brito Rocha 22 Jul '11

22 Jul '11

Senhores, estou usando o xtabs, no intuito de fazer alto parecido com o que se faz em uma tabela dinamica tenho um data.frame com dois fatores e uma variavel resposta continua estou fazendo algo do tipo: xtabs(y~A+B,data=base) tenho como resposta uma tabela com os níveis de A e B e uma estatística dentro de cada célula. Minha pergunta é: Como faço para mudar a estatística do xtabs? média, quantidade, mediana,.. não tem um parâmetro FUN dentro do xtabs. tem algum pacote com uma função equivalente que faça isso? 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 8

R is Hot!
by Alejandro C. Frery 21 Jul '11

21 Jul '11

Boa noite, Recomendo a leitura da matéria "R is Hot!" http://t.co/MFs3Bz9 Alejandro -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008

1 0

Incluir formula em um função
by Leo Bastos 20 Jul '11

20 Jul '11

Olá pessoal, Tenho algumas funções que eu gostaria de incluir uma formula como parametro de entrada. Da mesma forma que as funcoes lm e glm. Alguem poderia me indicar um exemplo, ou indicar um artigo/tutorial/link/whatever para que eu possa incluir uma formula com entrada de uma funcao? Por exemplo, suponha que eu tenho um data frame com um vetor de resposta y e algumas covariáveis. E suponha q eu quero uma função q retorne uma estimativa dos coeficientes de uma regressão linear. (sim, eu sei q o lm faz isso, mas eu gostaria de um exemplo de função onde eu possa usar uma formula como parametro de entrada). Uma função com um vetor y e uma matrix X é trivial, mas como fazer usando uma formula como parametro de entrada? Agradeço antecipadamente qq ajuda. Abraço, Leo Bastos

2 1