Maio 2013 - R-br - DINF UFPR

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

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

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

08 Mai '13

Boa noite a todos! 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

Krigagem
by Arthur Telles Calegario 08 Mai '13

08 Mai '13

Bom dia Senhores,Estou com um pequeno problemas em relação a Krigagem e talvez os senhores tenham uma resposta.Na verdade o problema esta em preparar os dados. Estou trabalhando com evapotranspiração no Rio Doce e possuo as coordenadas das estações meteorológicas em graus. Contudo, para a construção do semivariograma, necessito que a distância entre as mesmas seja linear.Assim devo projeta-las. É um procedimento simples quando todos as estações estão em um mesmo 'datum'. Contudo, as estações que estão do 'datum' 23S são insuficientes para a construção do semivariograma. Assim necessito saber como faço para linearizar as distâncias entre estações meteorológicas cujos limites extrapolam um 'datum' apenas. Entenderam? rs. Agradeço desde já,Arthur Calegario.

11 33

Krigagem
by Tito Conte 03 Mai '13

03 Mai '13

Bom dia, Alguém poderia me ajudar com krigagem? Estou usando o fields, e tenho uma matriz com 57703, com x, y, z pontos, quando uso o comando Krig ele da vários erros. O comando é K = Krig(xy,z) xy são as posições e z é a altura Tito Conte

4 3

pp-plot
by Julia Maria 02 Mai '13

02 Mai '13

Pessoal, o R tem alguma forma pré definida para se fazer um Probability-Probability plot? Por exemplo, o pacote circular tem a função: pp.plot(x), que verifica o ajuste a uma distribuição von-mises. E se eu quiser fazer o mesmo com outras ditribuições? Grata

1 0

set.seed
by Daniel Marcelino 02 Mai '13

02 Mai '13

Caros, estou realizando um experimento e acabei constatando que a semente que eu distribuir é fundamental para eu obter as amostras binárias. Eu estou interessado em distribuir 25 unidades para tratamento e o restante para controle. Há alguma maneira de saber, sem ser por tentativa e erro, que minha distribuição terá 25 unidades em cada categoria? Por exemplo, > # 50 Bernoulli trials > set.seed(51) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 22 28 > set.seed(12345) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 23 27 > set.seed(123) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 25 25 > set.seed(2) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 25 25 -- \begin{signature} Daniel Marcelino ☁ dm.silva(a)umontreal.ca ☎ (514) 343 6111 #3799 Skype: d.marcelino ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature} "Small steps toward a much better world"

4 5

excluir linhas que tenham códigos repetidos
by Fátima Lima Paula 02 Mai '13

02 Mai '13

Pessoal, existe algum comando que exclua as linhas que apresentam códigos repetidos? Não seria deixar apenas um não. É excluir todas. Exemplo nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano","Paulo","Ana") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima","Pitanga","Amora") aih=c(6867,6798,8903,7865,8903,6745,6867,8903) ano=c(1944,1945,1923,1954,1947,1945,1946,1944) dataint=c(20121004,20120907,20120703,20120301,20120805,20121209,20120907,20120703) sexo=c("M","F","F","M","M","M","M","F") banco1=cbind(nome,sobrenome,aih,sexo,ano,dataint) banco1=data.frame(banco1) banco1 nome sobrenome aih sexo ano dataint 1 Ricardo Cunha 6867 M 1944 20121004 2 Maria Andrade 6798 F 1945 20120907 3 Tereza Silva 8903 F 1923 20120703 4 Pedro Soares 7865 M 1954 20120301 5 José Silva 8903 M 1947 20120805 6 Germano Lima 6745 M 1945 20121209 7 Paulo Pitanga 6867 M 1946 20120907 8 Ana Amora 8903 F 1944 20120703 Quero excluir as linhas que apresentam as aihs repetidas (6867 e 8903), de forma que o banco fique assim: nome sobrenome aih sexo ano dataint 1 Maria Andrade 6798 F 1945 20120907 2 Pedro Soares 7865 M 1954 20120301 3 Germano Lima 6745 M 1945 20121209 Obrigada -- "Minha felicidade depende da qualidade dos meus pensamentos"

5 5

Função maxDrawdown
by Marcos Vagner Macedo 01 Mai '13

01 Mai '13

Boa noite pessoal, Não estou conseguindo utilizar a função maxDrawdown do pacote PerformanceAnalytics. Alguém pode me ajudar? Estou utilizando dias úteis no formato zoo. Quando utilizo a função tenho o seguinte resultado: > t(round(maxDrawdown(fi_97_12[,2:4]),4)) Worst Drawdown Ibovespa NaN ALFA.DINAMICO.FIA NaN ALFA.FIA NaN Warning messages: 1: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct 2: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct 3: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct > Muito Obrigado, Marcos Macedo

1 0