- R-br - DINF UFPR

Ajuda em uma função
by Marcelo Zeuli 15 Jun '12

15 Jun '12

Bom dia, Algume conheçe como implementar função impulso, no R? Do tipo delta de Dirac: O sinal delta de Dirac δ(t) foi criado pelo físico inglês Paul A. Maurice Dirac (1902 - 1984) e pode ser visto como o equivalente contínuo do delta de Kronecker (matemático alemão,1823 - 1991) e como o limite de uma distribuição normal Grato, Marcelo.

1 0

Ajuda em uma função
by Sérgio Henrique almeida da silva ju 15 Jun '12

15 Jun '12

Olá Gostaria de fazer uma função e como não entendo nada, gostaria de pedir ajuda de vocês Eu uso o Alfa de Cronbach para analisar a consistência interna de uma escala, normalmente se calculo o alfa total e o alfa se uma variável for retirada. O alfa total é fácil de fazer, porém quando eu vou fazer o calculo de quando retiro uma variável tenho que fazer tudo manualmente. Ex.: a = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2) b = c(4,5,5,4,4,3,2,1,4,3,4,5,3,2,3,4,5,5,4,4,1) c = c(3,3,4,4,5,4,3,2,4,5,5,3,4,4,5,4,3,2,3,4,5) d = c(2,2,3,3,4,3,2,3,4,5,5,5,4,3,2,3,4,5,4,3,2) dados = data.frame(a,b,c,d) library(psy) # Alfa Total cronbach(dados[,1:4]) alfa total = 0.5324005 Agora preciso retirar a variável "a" e rodar novamente o calculo, depois colocar ela denovo no banco e retirar a variável "b" e assim sucessivamente, dessa forma terei o valor do alfa caso eu retire uma variável. # Alfa retirando a variável a = 0.220339 b = 0.3239247 c = 0.7669716 d = 0.2198046 Qual a melhor forma de fazer isso através de uma função? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 5

Re: [R-br] Tornar mensagens de erro "invisíveis" em função
by Ivan Bezerra Allaman 15 Jun '12

15 Jun '12

Tente isto, d<-data.frame(x=1:10, y=rep(c(FALSE,TRUE),c(4,6))) d glm(y~x, data=d, family=binomial) keepWarnings <- function(expr) { localWarnings <- list() value <- withCallingHandlers(expr, warning = function(w) { localWarnings[[length(localWarnings)+1]] <<- w invokeRestart("muffleWarning") }) list(value=value, warnings=localWarnings) } z <- keepWarnings(glm(y~x, data=d, family=binomial)) z z$value (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 2

R2 para modelos Poisson ou Binomial Negativo
by Sérgio Henrique almeida da silva ju 14 Jun '12

14 Jun '12

Olá Alguém sabe como estimar o R2 para modelos de Poisson ou Binomial negativo? Encontrei esse documento falando da possíbilidade dessa estimação: http://www.google.com.br/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0CFcQ… Porém não achei de como calcular no R. Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

7 7

Estacionariedade Regressão Espacial
by Isabela Zardo 14 Jun '12

14 Jun '12

Boa tarde, Alguém sabe como eu faço pra verificar a estacionariedade no R para a regressão espacial? Já que é preciso verificar para que o resulatdo do índice de Moran seja válido. Preciso de ajuda com urgência!!! Agradeço!! Isabela Zardo

1 0

Dúvida com quando o array tem NA
by Victor Eduardo 14 Jun '12

14 Jun '12

Galera, estava com um array de tamanho 150 onde tinham matrizes 1x1. Sendo que em alguns tinham NA. Dentro de minha função criei uma variável auxliar para identificar esses NA. Fiz do seguinte modo: aux<-NULL for(t in 1:n){ if(is.na(Y[,,t])==TRUE){ aux[t]<-0 }else{ aux[t]<-1 } } Ao identificar, o aux terá tamanho de 150, de 0 e 1. Feito essa verificação, estou pulando para a parte debaixo que é assim: for(t in 1:n){ if(aux[t]==1){ B<-matrix(Z[,,t],ncol=m) U[,,t]<-Y[,,t]-B%*%At[,,t]-D[,,t] F[,,t]<-B%*%Pt[,,t]%*%t(B)+H[,,t] K<-T[,,t]%*%Pt[,,t]%*%t(B)%*%ginv(F[,,t]) L[,,t]<-T[,,t]-K%*%B At[,,t+1]<-T[,,t]%*%At[,,t]+K%*%U[,,t]+C[,,t] Pt[,,t+1]<-T[,,t]%*%Pt[,,t]%*%t(L[,,t])+R[,,t]%*%Q[,,t]%*%t(R[,,t]) }else{ At[,,t+1]<-T[,,t]%*%At[,,t]+C[,,t] Pt[,,t+1]<-T[,,t]%*%Pt[,,t]%*%t(L[,,t])+R[,,t]%*%Q[,,t]%*%t(R[,,t]) } } return(list(At,Pt)) Só que o R está dando um erro que não estou conseguindo resolver, o erro é esse: Erro em svd(X) : infinite or missing values in 'x' Abraços!

1 1

Erro em plot()
by Ricardo Dobrovolski 13 Jun '12

13 Jun '12

Oi Pessoal, Eu quero plotar 18 gráficos de dispersão. No entanto, eu tive a seguinte mensagem de erro: "Erro em plot.new() : margens da figura muito grande" Os comandos que utilizei são os seguintes: par(mfrow=c(6,3)) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.na "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.na "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.na "]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.rn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.rn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.rn"]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.ra"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.ra"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.ra"]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) onde agro.gdp e agro.hdi são duas tabelas de dados e os termos entre aspas representam nomes das minhas colunas. Atenciosamente, RD.

3 2

Ordenar linhas em uma matriz
by Eric Ferreira 13 Jun '12

13 Jun '12

Prezados colegas, Estou com o seguinte problema: desejo ordenar uma matriz de acordo com os nomes de suas linhas. Seja o exemplo, M<-matrix(c(4,5,6,3,3),5,1) rownames(M)<-c('c','d','e','a','b') A matriz ordenada que eu desejo obter é: M1<-matrix(c(3,3,4,5,6),5,1) rownames(M1)<-c('a','b','c','d','e') Alguém conhece algum comando do tipo "sort by..."?? Obrigado e abraços, -- Dr Eric B Ferreira Exact Sciences Department Federal University of Alfenas Brazil

3 3

Re: [R-br] R2 para modelos Poisson ou Binomial Negativo
by Ivan Bezerra Allaman 13 Jun '12

13 Jun '12

Sua resposta aqui, http://r.789695.n4.nabble.com/compute-coefficient-of-determination-R-square… (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Tornar mensagens de erro "invisíveis" em função
by Emmanuel Arnhold 12 Jun '12

12 Jun '12

É possível tornar "invisíveis" as mensagens de erros de uma função? Incluindo mensagens de avisos perdidas. Emmanuel

1 0