- R-br - DINF UFPR

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

duvida sobre função try {base}
by Augusto Ribas 10 Mai '13

10 Mai '13

Ola, eu gostaria de saber, se tem como eu mudar de forma simples a mensagem de erro que o R envia ao tentar usar uma função e receber um erro. Por exemplo: #Imagine que eu tenho a seguinte lista, e precise aplicar uma função #Exemplo de lista de dados teste<-list(2, "a", 1 , 10 , "b") teste #Preparando um lugar para guardar os resultados saida<-list() #Se eu tentar fazer um loop com for, ele vai parar #quando chegar a letra, porque nao da para fazer o log de caracter for( i in 1:length(teste)) { saida[i] <-log(teste[i]) } #nesse caso o loop para e eu nao tenho a saida que desejo saida #Mas então eu tentei usar a função try #minha ideia era que se der para calcular, calcule #senão de uma msg de erro e siga em frente pro proximo for( i in 1:length(teste)) { saida[[i]] <-try( log(teste[[i]][1]) ) } #agora eu tenho minha saida, com erro onde não deu para calcular, isso que eu queria saida #o que eu queria saber se é possivel fazer o seguinte #se der certo, execute a função, senão imprima o erro que eu falar #nesse caso eu gostaria de uma saida desse tipo: saida.exemplo<-list(0.6931472,"Verificar",0,2.302585,"Verificar") saida.exemplo ############################################################# Esse exemplo eu peguei na função try. Mas na verdade eu estou usando o pacote taxize (http://ropensci.github.io/taxize_/) para conferir nomes de especies, e para conferir listas grandes, as vezes para quando a conexão com a internet falha, ou o nome da especie não é encontrado, resumindo eu queria saber se existe uma forma de falar pro R, fazer isso (minha função), se der erro retorne isso (minha mensagem). Agradeço desde ja a atenção. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

4 7

Estrutura de repetição - Função
by Sérgio Henrique almeida da silva ju 10 Mai '13

10 Mai '13

Prezados Quero criar uma função na qual eu digito um conjunto de notas e no final o programa me dê a média dessas notas. O programa é interrompido quando eu digito uma nota negativa. Estou fazendo dessa forma: media <- function(){ soma = 0 n=0 X1= cat("Digite uma nota: ", sep="\n") nota= scan(stdin(),n=n) while (nota >= 0){ soma=soma+nota n = n + 1 X1= cat("Digite uma nota: ", sep="\n") nota= scan(stdin(),n=n)} return (sprintf("%.1f",soma/2))} Onde posso estar errando? 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) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 1

Tinn-R news
by Jose Claudio Faria 10 Mai '13

10 Mai '13

Dears Tinn-R users, A new version of the Editor/GUI/IDE Tinn-R (2.4.1.6) was released today. News at: http://sourceforge.net/p/tinn-r/news/2013/05/tinn-r-2416-released/ The project now has its proper page: http://nbcgib.uesc.br/lec/software/des/editores/tinn-r/en Download is available in both: 1- WebPage: http://nbcgib.uesc.br/lec/software/des/editores/tinn-r/en 2- SoureForge: http://sourceforge.net/projects/tinn-r/ All the best, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

Interface GPS com R - gpsbabel e maptools
by Hélio Gallo Rocha 09 Mai '13

09 Mai '13

Testei a interface de um garmin no R com pacote maptools, mas não deu certo, com esse erro: gpsbabel not found, ou seja, parece que o gpsbabel fazia parte do maptools, mas não faz mais. Baixei o soft gpsbabel. interessante. Ele converte em vários formatos, quero ver se consegue extrair dados da portadora L1 do Garmin. Mas não fez a interface com o R. Alguém já tentou? -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0