- R-br - DINF UFPR

SoundexBR
by Fátima Lima Paula 15 Dez '12

15 Dez '12

Prezados, recebi do Daniel Marcelino o soundex para ser utilizado no RecordLinkage. Gostaria de saber se esse soundex é validado e pode ser utilizado sem problemas. Obrigada Fátima

3 4

Mallows - estimador robusto
by Natália Bordin Barbieri 15 Dez '12

15 Dez '12

Olá, preciso ajustar o seguinte modelo de estimação robusta, através do método de "Mallows' para o modelo de regressão logística que segue: require(robust) DM.glmRob3<-glmRob(DM ~ INFLAMACAO + IDADE + SEXOM + RACACOR + HIPERT + IMC + RCQ + TRIGT, family=binomial, data=DM, weights=NULL, method ="mallows", model = TRUE, control = glmRob.control) Porém, ao rodar, o que retorna é o seguinte erro: 'The covariance matrix has become singular during the iterations of the MCD algorithm. There are 284 observations (in the entire dataset of 400 obs.) lying on the hyperplane with equation a_1*(x_i1 - m_1) + ... + a_p*(x_ip - m_p) = 0 with (m_1,...,m_p) the mean of these observations and coefficients a_i from the vector a <- c(1, 0, 0, 0, 0, 0, 0, 0) Erro em solve.default(V) rotina Lapack dgesv: sistema é exatamente singular' Alguém já trabalhou com isso?Sabe me ajudar? Obrigada, Natália

2 7

Analise dados proporção e contagem
by José Francisco Arruda e Silva 14 Dez '12

14 Dez '12

Ola pessoal.Alguem poderia me dar uma dicas de como posso analisar dados de proporcao e contagem, uma vez que necessito saber se existe diferenca entre os tratamentos (dietas) no que se refere ao numero de individuos que nasceram e da razao sexual. segue abaixo uma pequena amostra. dados <- expand.grid(rept=1:15, trat=LETTERS[1:4]) dados$razao_sexual <- c(1,0.5,0.8,1,0.5,0.8,1,1,0.8,0.6,0.8,1,0.5,1,0.8,1,0.4,1,1,0.4,0.8,0.8,1,1,0.6,0.6,0.8,0.8,0.8,1,0.8,1,1,0.6,1,1,0.8,1,0.8,1,0.8,1,0.8, 1,0.8,1,0.6,0.7,0.7,1,1,1,0.8,1,1,1,1,1,0.8,1) dados$nascidos<-c(9,10,9,11,13,8,10,11,10,9,10,10,9,10,9,10,14,10,9,10,10,7,10,9,10,8,10,10,12,10,14,10,9,10,10,11,10,10,9,10,10,12,11,10,15,10,10,10,11, 9,10,10,9,10,9,9,10,10,10,9) Obrigado. Chicao.

3 4

AJUSTE DOS BETAS DE REGRESSÃO
by tiago souza marçal 14 Dez '12

14 Dez '12

Caros colegas, não estou conseguindo estimar os betas dos modelos quadrático e cúbico da regressão através da função reg.poly() do pacote ExpDes. y<-c(13.337250,24.96667,22.34600,24.31983) x<-c(1000000, 10000000, 50000000, 100000000) Glresíduo<- 72 Sqresíduo<- 1159.9 Gltratamento<- 3 Sqtratamento<- 167.6 reg.poly(y,x,72,1159.9,3,167.6) Agradeço pela atenção Att. Tiago.

2 1

Pacote RNetLogo
by Felipe Buchbinder 14 Dez '12

14 Dez '12

Olá, Estou tentando utilizar o pacote RNetLogo. Eu consigo abrir o NetLogo a partir do R mas, quando tento qualquer outro comando, recebo a seguinte mensagem: Error in .jcall(nl.obj, "V", "loadModel", .jnew("java/lang/String", model.path)) : RcallMethod: invalid object parameter Alguém sabe o que está acontecendo? Obrigado, Felipe

1 0

Seleção de modelo misto
by Diego Pujoni 14 Dez '12

14 Dez '12

Olá pessoal, tenho um experimento de medidas repetidas que estou analisando através de um modelo misto. Nove unidades amostrais independentes (frascos com meio de cultura com algas) foram alocadas aleatoriamente em 3 grupos ("c","t1","t2"). Não há a necessidade de inclusão de efeito aleatório de intercepto, pois as nove unidades amostrais são homogêneas entre si (retiradas da mesma amostra mãe). A concentração de algas foi mensurada a cada dois dias durante 10 dias. O objetivo é testar diferenças entre os grupos (i.e. tratamentos). Estimei um modelo somente com o intercepto e com a interação Group*Day, para testar qual é o melhor library(nlme) library(lme4) Day = rep(c(0,2,4,6,8,10),each=9) Group = as.factor(rep(c("c","c","c","t1","t1","t1","t2","t2","t2"),6)) Individual = rep(1:9,6) X = c(0.71,0.72,0.71,0.72,0.72,0.72,0.70,0.69,0.70,0.72,0.72, 0.71,0.72,0.72,0.71,0.71,0.70,0.71,0.73,0.73,0.69,0.74, 0.69,0.73,0.67,0.71,0.69,0.71,0.71,0.72,0.70,0.71,0.70, 0.52,0.64,0.60,0.70,0.73,0.73,0.67,0.66,0.71,0.47,0.56, 0.54,0.65,0.73,0.73,0.67,0.71,0.58,0.44,0.52,0.58) xyplot(X~Day, groups=Group) LME = lme(X ~ 1, random = ~Day|Individual) Erro em lme.formula(X ~ 1, random = ~Day | Individual) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) LME1 = lme(X ~ Group*Day, random = ~Day|Individual) Erro em lme.formula(X ~ Group * Day, random = ~Day | Individual) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) LMER = lmer(X ~ 1 + (Day|Individual)) LMER1 = lmer(X ~ Group*Day + (Day|Individual)) AIC(LMER) [1] -179.0302 AIC(LMER1) [1] -151.1938 anova(LMER,LMER1) Data: Models: LMER: X ~ 1 + (Day | Individual) LMER1: X ~ Group * Day + (Day | Individual) Df AIC BIC logLik Chisq Chi Df Pr(>Chisq) LMER 5 -187.2 -177.26 98.602 LMER1 10 -203.2 -183.31 111.600 25.996 5 8.939e-05 *** xyplot(fitted(LMER)~Day, groups=Group) xyplot(fitted(LMER1)~Day, groups=Group) 1a pergunta: Por que a função nlme:lme não converge, mas a lme4:lmer converge? A segunda é melhor que a primeira? 2a pergunta: Por que a função "anova" dá valores distintos de AIC para os dois modelos. Se olharmos o valor do AIC de cada um, o modelo LMER é melhor que o LMER1, mas a "anova" diz o contrário. 3a pergunta: Por que os valores ajustados do modelo somente com o intercepto (LMER) variam ao longo do tempo? Muito obrigado pela ajuda Diego PJ

1 0

mínimos quadrados bivariado
by Luiz Henrique da Conceicao Leal 14 Dez '12

14 Dez '12

Prezados Existe algum pacote no R que realize análise de regressão através dos mínimos quadrados bivariado - Bivariate Least Squares (BLS)?

3 3

[Dúvida] Tratamento de erro.
by Pedro Rafael 13 Dez '12

13 Dez '12

Pessoal como realizar tratamento de erros em R? Estou com uma função em que é executado um loop. Neste laço é executado o comando arima que em determinada interação poderá não convergir dando assim uma mensagem de erro que ocasionará o termino do laço. Gostaria de ignorar esses erros. O código segue: nrep = 50 matriz_result = matrix(,nrow=nrep,ncol=7) for(i in 1:nrep){ par = trunc(runif(5,0,4)) resultado = arima(retirar(serie,3)$Serie,order=c(par[1],2,par[2]), seasonal = list(order = c(par[3], par[4], par[5]))) aic = resultado$aic matriz_result[i,] = c(par[1],1,par[2:5],aic) } Gostaria de que no caso de erro a minha do objeto matriz_result seja completada por NA. Grato [ ], Pedro Rafael Diniz Marinho.

2 2

Qual pacote que tem o comando scripchart
by jairan dionizio 13 Dez '12

13 Dez '12

Bom dia todos e todas! gostaria de uma ajudinha. Qual pacote do contém comando scripchart? Jairan Dionizio Skype:jairan.dionizio Facebook: Jairan Dionizio Twitter: @jairandionizio Tel.:(84) 9423-5286 "A força não provém da capacidade física e sim de uma vontade indomável." (Mahatma Gandhi) Cuidados ao enviar e-mail Ao reenviar apague o meu e-mail e o meu nome, apague também os endereços dos amigos antes de reenviar qualquer e-mail, encaminhe como cópia oculta (Cco ou Bcc) aos SEUS destinatários.Agindo sempre assim dificultaremos a disseminação de vírus, spams, banners e outros.

3 2

Problemas Git-Hub
by Leandro Marino 13 Dez '12

13 Dez '12

Caros, alguém saberia me dizer o motivo do código abaixo não estar funcionando?! classcol = c(rep('character',5),rep('numeric',2),rep('character',2)) url <- 'https://raw.github.com/leandromarino/TBAmoMunic/master/Munic.dat' munic <- read.table(url,colClasses=classcol,sep='\t',header=T) Obrigado. Leandro

5 11