- R-br - DINF UFPR

Mensagem de erro da função optim
by Rubem Kaipper Ceratti 24 Mai '12

24 Mai '12

CaRos, Ao rodar o script abaixo, a função 'optim' retorna a mensagem: 'Erro em function (par) : não foi possível encontrar a função "fn"'. Tentei, mas não consegui encontrar a razão para o erro. Alguém teria uma idéia? Grato, Rubem ## Script R: library(tweedie) # Simulação de dados period<-factor(1:4) herd<-factor(1:10) dat<-expand.grid(period=period,herd=herd) beta<-c(-.3,1.7,2.5,3.4) x<-model.matrix(~period,dat) b<-rnorm(length(levels(dat$herd)),sd=.3) z<-model.matrix(~-1+herd,dat) mu<-as.vector(exp(x%*%beta+z%*%b)) phi<-.8; p<-1.6 dat$resp<-rtweedie(nrow(dat),p,mu,phi) y<-dat$resp th<-c(rep(0,ncol(x)),1,1,1.5) beta<-th[1:ncol(x)] sig<-th[ncol(x)+1] phi<-th[ncol(x)+2] p<-th[ncol(x)+3] z.ast<-sig*z l.u<-function(u,y,z,x,beta,phi,p){ eta.u<-as.vector(x%*%beta+z%*%u) f.y<-log(dtweedie(y,p,exp(eta.u),phi)) l.yu<-sum(f.y)+sum(dnorm(u,log=T)) return(-l.yu) } uhat<-optim(rep(0,ncol(z)),l.u,method="BFGS",hessian=T, y=y,z=z.ast,x=x,beta=beta,phi=phi,p=p)

2 2

Erro em file("toto230876.tmp", "w")
by Ricardo Dobrovolski 24 Mai '12

24 Mai '12

Olá pessoal, Criei uma função que produz arquivos de saída. No entanto, ela tem sido interrompida com a seguinte mensagem de erro: Erro em file("toto230876.tmp", "w") : não é possível abrir a conexão Além disso: Mensagens de aviso perdidas: 1: In file.create("toto230876.tmp") : não foi possível criar arquivo 'toto230876.tmp', motivo 'Permission denied' 2: In file("toto230876.tmp", "w") : não foi possível abrir o arquivo 'toto230876.tmp': Permission denied Porém, isso só ocorre às vezes, em diferentes passos da função. Abraço, RD.

1 0

Atribuir nomes a códigos
by Edson Lira 24 Mai '12

24 Mai '12

Boa noite a todos. Estou com um pequeno problema. Tenho um banco aonde a variável bairro está codificada com números, e tenho um outro arquivo que consigo ler, aonde contenho os códigos e os bairros para cada código. como associar os códigos do meu banco com a relação de bairros do outro arquivo. A disposição é mostrada a seguir: IDBairro BAIRRO BZONA CodCida 1 92 Centro Presidente Figueiredo 3 2 93 Morada do Sol Presidente Figueiredo 3 3 94 Centro Balbina 3 4 95 Honório Roldão Presidente Figueiredo 3 5 96 Aida Mendonça Presidente Figueiredo 3 6 97 José Dutra Presidente Figueiredo 3 7 98 Tancredo Neves Presidente Figueiredo 3 8 99 Galo da Serra Presidente Figueiredo 3 9 100 Sol Nascente Presidente Figueiredo O banco de dados está assim: Os bairros em código: 92 92 92 92 92 92 92 92 92 92 97 93 93 93 93 93 95 93 93 93 93 95 95 95 95 95 93 95 93 95 95 95 95 95 95 95 93 93 93 95 Como associar o código do banco com o bairro do arquivo de códigos? Obrigado. Edson Lira Estatístico Manaus-Amazonas

3 2

Extração de QMerros em lme
by Marcelo Cardoso mello 23 Mai '12

23 Mai '12

Para extrair o QMerro de subparcela na lme posso usar model$sigma, mas como faço para obter o QMerro relativo as parcelas?

2 1

Pacote_gerar_mapas
by Rafael Tieppo 23 Mai '12

23 Mai '12

Necessito gerar mapas a partir de X (LONGITUDE) Y (LATITUDE) Z(teor de umidade no solo). O mapa deve ser 2D e a variável Z pretendo representar por cores através de uma legenda, ex.: 12-16% azul ; 16,01 - 18% verde ; etc Testei o GeoR e o adehabitatMA, porém nao tive os resultados pretendidos. Alguém pode me indicar um pacote? Obrigado. Atenciosamente, Eng. Agrícola Rafael Cesar Tieppo UNEMAT - Campus Tangará da Serra - MT / Departamento de Agronomia Fone:+55-65 3311-4922 site:http://www2.unemat.br/rafaeltieppo "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".

3 2

Bayesiana [OFF-TOPIC]
by ASANTOS 23 Mai '12

23 Mai '12

Bom dia Pessoal, Gostaria de saber se alguém quem trabalha com Bayesiana poderia recomendar algumas boas literaruras para estudo para que esta començando a usar esta abordagem, preferencialmente gostaria de material aplicado ao R e não ao winBugs. Obrigado, *Alexandre dos Santos *Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG

4 3

duvidas sobre teste de causalidade de granger
by Marcelo Zeuli 23 Mai '12

23 Mai '12

Bom dia, qual é o pacote mais adequado para teste de causalidade de granger? o FIAR? ********************************************* Type Package Title Functional Integration Analysis in R Version 0.4 Date 2011-04-01 Author Bjorn Roelstraete Maintainer Bjorn Roelstraete <Bjorn.Roelstraete(a)ugent.be> Description Contains DCM, (autoregressive) SEM, and multivariate Granger causality tests for analyzing fMRI connectivity data. Depends R (>= 2.11.1),lavaan(>= 0.4-7),Matrix,boot,accuracy,np ********************************************************** Grato, Marcelo Zeuli

2 2

Re: [R-br] Alteração em função de pacote
by Junior Beleti 22 Mai '12

22 Mai '12

Olá Benilton. No manual encontrei instruções de como se criar uma função em um novo pacote. No meu caso, onde apenas modifiquei a função eigen.f para eigem.f, gostaria de saber como fazer para carregá-la no momento da execução. Quando chamo a função eigen: eigen(cbind(c(1,-1),c(-1,1))) $values [1] 2 0 $vectors [,1] [,2] [1,] -0.7071068 -0.7071068 [2,] 0.7071068 -0.7071068 Já na chamada a eigem: eigem(cbind(c(1,-1),c(-1,1))) Erro: não foi possível encontrar a função "eigem" Att, Carlos.

3 2

Visualização interativa
by Eder David Borges da Silva 21 Mai '12

21 Mai '12

Pessoal gostaria de trocar uma ideia sobre visualização de dados. Ultimamente venho dando uma olhada do pacote googleVis, em especial na função gvisMotionChart, que possibilita criar gráficos onde podemos interagir com eles selecionando time, dando destaque a um nivel etc. O que estou procurando e algo como o gvisMotionChart, porem com uma especie de filtro a mais, por exemplo imaginem que tenho um conjunto de dados, emails da R-br, nos três últimos anos, onde uma coluna tenho o ano, uma segunda o participante, e uma terceira o numero de emails. Eu gostaria de fazer um gráfico selecionando por exemplo, no ano 0, apenas o user1 versus user2, onde essa seleção seria feia em uma caixa de seleção, selecionando ano 0, depois em um segunda caixa selecionando user1 e user2, e automaticamente o gráfico de barras por exemplo iria ser atualizado. para quem já ouviu falar seria um qlikview http://www.inteligencia-negocios.com/produtos/qlikdemos porem feito via R com a saida final de interface html. Oque quero compartilhar com vocês é o seguinte: alguém já viu algo assim feito em R?

2 2

Plots com dimensões diferentes
by Marcelo Claro de Souza 21 Mai '12

21 Mai '12

Preciso plotar 4 figuras no mesmo frame (2 x2) mas preciso que as figuras tenham dimensões diferentes. No caso preciso que as 2 figuras superiores ocupem apenas 25% do frame e as figuras inferiores 75%. Como posso fazer isso? Abraço. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1