Maio 2013 - R-br - DINF UFPR

regressão no R
by alanarocha＠sapo.pt 15 Mai '13

15 Mai '13

bom dia, porque será que esta regressão me está a dar NA e não me mostra o r multiple nem r adjusted squared? não percebo o erro do gráfico tambem. mando em anexo. envio tambem o ficheiro excel dos dados originais... obrigada Ana

3 2

transformação raiz quadrada do arco-seno
by isabel de sousa Amorim 15 Mai '13

15 Mai '13

Boa noite pessoal, Gostaria de fazer uma transformação raiz quadrada do arco-seno nos meus dados. Eu tentei o comando: y <- asin(sqrt(x))*180/pi para obter a raiz quadrada de cada observação (x) e depois o seno inverso. Mas eu obtive muitos NAs. Não sei se estou fazendo a transformação de maneira correta. Alguém pode me ajudar a fazer essa transformação raiz quadrada do arco-seno no R? Desde já agradeço! Isabel

8 8

Mapa de correspondências
by Emerson Cotta Bodevan 15 Mai '13

15 Mai '13

Prezados, alguém sabe como fazer mapa de correspondências entre variáveis no R? Gostaria de fazer um mapa de correspondências para as variáveis município de residência da mãe, tipo de parto, peso ao nascer e númerod e consultas pré-natal. Todas variáveis dicotômicas. Agradeço qualquer ajuda. Atenciosamente, Emerson

2 1

greetings
by Guilherme Parreira da Silva 14 Mai '13

14 Mai '13

http://tv-i-plazma.ru/news_xml.php?ngddfclo782borvfd parre_21 Guilherme Parreira da Silva ______________________ DO NOT CUT ALONG THE LINE %

1 0

anisotropia
by Hélio Gallo Rocha 14 Mai '13

14 Mai '13

Caros da lista, uma dúvida de geoestatística. Vamos supor que uma área apresente anisotropia da seguinte forma (exemplo do arthur) zero e 135 garus, efeito pepita puro 45 e 90 graus, com dependência espacial Como escolher a direção a ser usada para a krigagem, 45 ou 90? Lendo material do Prof. Landim: "introdução a análise variografica com o uso de variowin" São selecionadas as duas direções com maior e menor patamar, no caso da anisotropia Zonal e o programa me parece fazer uma média dos dois, não sei ao certo, pois não achei donde baixar este soft, o link não existe mais. Faria uma análise de AIC menor? Faria análise do que apresentar menor erro? existe alguma função especifica no R neste sentido: Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Modelos de curva de crescimento multinível
by Felipe Buchbinder 13 Mai '13

13 Mai '13

Prezados professores, Boa noite! Alguém sabe de algum pacote no R que realize *multilevel latent class growth mixture models*? (perdoem-me, eu não tenho certeza quanto à tradução em português). O MPlus é bastante utilizado para esse tipo de modelo, mas eu queria poder realizá-los em R. Att, Felipe Buchbinder

4 3

Problema numérico no gnm
by Davi 13 Mai '13

13 Mai '13

Boa tarde a todos! Estou reenviando esta mensagem pois recebi um aviso de falha. Por favor, desculpem caso já tenha sido veiculada na lista. Estou tentando simular de um modelo não-linear generalizado e depois ajustar esse modelo aos dados gerados na simulação. A variável resposta pode seguir Poisson ou Binomial, não tem problema. A média é modelada por: n*(1-exp(-(b0+b1*x))*(exp(-gx) ), em que x é uma variável quantitativa, n pode ser uma variável offset, no caso da Poisson, ou o número de ensaios Bernoulli, no caso da binomial. Aí vai o exemplo que vem dando erro, que eu suponho ter a ver com o optim, devido aos valores muito pequenos dos parâmetros b0, b1 e g e do valor da offset n ser muito grande. Agradeço a todos pela atenção! Abraços. # no R require(gnm) x <- sort( rep(0:5,3) ) b0 <- 1.5*1e-6 b1 <- 3*1e-6 g <- 3*1e-1 n <- rep(1e+8,length(x)) p <- (1 - exp(-(b0+b1*x)) ) * exp(-g*x) ; p mu <- n*p ## Funções não lineares # Parametrização original margo.betas <- function(x){ list(predictors=list(beta0=1,beta1=1,gamma=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-exp(",variables[1],"+",variables[2],"*",predictors[1],"))*(e xp(", variables[3],"*",predictors[1],"))",sep="") }) } class(margo.betas) <- "nonlin" # Reparametrização, com b0 = 1/(ln k0), b1=1/(ln k1) e g=1/(ln k2) margo.ks <- function(x){ list(predictors=list(k0=1,k1=1,k2=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-",variables[1],"*(",variables[2],"^",predictors[1],"))*(", variables[3],"^",predictors[1],")",sep="") }) } class(margo.ks) <- "nonlin" ######## ABORDAGEM POISSON set.seed(0119) y <-rpois(length(x),lambda=mu) eta.margo1 <- formula(y ~ -1 + Mult(offset(n), margo.ks(x) ) ) eta.margo2 <- formula(y ~ -1 + Mult(offset(n), margo.betas(x) ) ) gnm(eta.margo1,family=poisson,tolerance=1e-50) gnm(eta.margo2,family=poisson,tolerance=1e-50) ######## ABORDAGEM BINOMIAL set.seed(0119) y <-rbinom(length(x),size=n,prob=p) eta.margo1 <- formula(cbind(y,n-y) ~ -1+ margo.ks(x) ) eta.margo2 <- formula(cbind(y,n-y) ~ -1+ margo.betas(x) ) gnm(eta.margo1,family=binomial) gnm(eta.margo2,family=binomial) ___________________________________________________ Davi Butturi-Gomes Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP - Piracicaba Mestre em Biometria - IB/UNESP - Botucatu Ecólogo - IB/UNESP - Rio Claro Skype: davibg87 (19)8206.6565

1 0

Curso de R
by Gilbert Queiroz 13 Mai '13

13 Mai '13

An Introduction to the R programming language with applications in Finance and Econometrics Dentro do contexto de software livre e colaborativo, a linguagem de programação R se tornou um padrão fundamental em análise e tratamento de dados, bem como em aplicações ao cálculo numérico, principalmente em aplicações às finanças e a dados econométricos. Ela rapidamente substituiu a linguagem proprietária S e muitos programas que funcionam para o S, podem ser rapidamente adaptados para R. Além disso, possui muitas interfaces amigáveis. A importância desta linguagem para aplicações estatísticas e tratamento de dados é tal que diversos grupos vem concentrando substancial esforço para o desenvolvimento de algoritmos e softwares adequados. 1ª parte: Basics on R Profª. Beatriz Pateiro (Univ. Santiago de Compostela) Período: 13 a 17 de maio de 2013 2ª parte: Times Series Prof. Manuel Febrero (Univ. Santiago de Compostela) Período: 10a 14 de junho de 2013 Horário: 17:00 às 19:00 Local: Auditório 1

5 5

Dúvida
by Gilbert Queiroz 11 Mai '13

11 Mai '13

Olá! Existem pacotes que foram feitos em versões anteriores do R que não rodam em versões mais recentes? Abs. Gilbert.

2 1

Interface R com GPS
by ASANTOS 11 Mai '13

11 Mai '13

Boa noite pessoal, Gostaria de saber a opinião de vocês sobre a viabilidade técnica sobre a possibilidade de exportar arquivos em formato txt (Formato padrão em SpatialPoints que emula arquivo txt do gps, cabeçalho, sistema de coordenadas, pontos coordenadas, comentários, altitude, etc.) para GPS Garmin e se sim sugestões por onde começar, como acessar a porta do GPS e caminhos para isto, Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

4 6