- R-br - DINF UFPR

transferir modelo da função nls para a função curve() automaticamente
by Fernando Antonio de souza 03 Mar '16

03 Mar '16

Olá a todos. Estou criando uma função para trabalhar com o pacote Rpanel no qual pretendo criar uma lista de seleção com opção para escolher vários modelos. Como vou trabalhar com vários modelos quero criar funções que automatizem todo o processo. Minha questão é a seguinte: Eu ajusto um modelo não-linear, produzo um gráfico e adiciono a curva do modelo aos gráfico, utilizando a função curve(). Gostaria de passar o modelo para a função curve() de forma que dependendo do modelo selecionado a expressão incluída na função curve() se altere. Tentei a opção modelo$call conforme código abaixo o qual me retorna o lado esquerdo da fórmula. Entretanto não funciona. A função curve retorna a seguinte mensagem de erro; "Error in curve(fm1$call[[2]][[3]], add = TRUE, col = 2) : 'expr' must be a function, or a call or an expression containing 'x' Isto se deve ao fato de a formula retornada pelo comando fm1$call[[2]][[3]] ser do typo "language" typeof(fm1$call[[2]][[3]]) [1] "language" ja tentei utilizar o comando as.expression e não obtive sucesso. Procurei na internet e nada! Como eu posso substituir o comando "fm1$call[[2]][[3]]" na função curve de forma que ela funcione? fm1 <- nls(circumference ~ A/(1+exp((B-age)/C)), Orange,start = list(A = 160, B = 700, C = 350)) plot (circumference~ age, Orange) with(as.list(coef(fm1)), curve("fm1$call[[2]][[3]]",add=TRUE, col=2)) -- ======================================================================= Fernando Souza Zootecnista, DSc. Produção Animal e-mail:nandodesouza@gmail.com https://producaoanimalcomr.wordpress.com/ ========================================================================

3 8

problema ao conveter data:hora armazenado em texto (função strptime)
by Paulo Nogueira Starzynski 03 Mar '16

03 Mar '16

Pessoal, estou tentando fazer uma conversão simples de data e hora armazena em um vetor de caracters para formato "POSIXlt" "POSIXt". Porém, está acontecendo algo que não entendo ao usar a função strptime. Se eu passo o vetor com os dados diretamete, retorna NA. Se eu copio e colo os valores na função, dá certo... ou o mesmo usando o dput, como abaixo... vejam... > class(lista.exames[,"Data de Liberação:"])[1] "character"> dt.hr <- lista.exames[,"Data de Liberação:"][1:10]> strptime(dt.hr, format = "%d/%m/%Y %H:%M:%S") [1] NA NA NA NA NA NA NA NA NA NA # não deu certo > dput(lista.exames[,"Data de Liberação:"][1:10])c("10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34")> dt.hr <- c("10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34",+ "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34",+ "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34",+ "10/06/2007 01:03:34")> strptime(dt.hr, format = "%d/%m/%Y %H:%M:%S") [1] "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" [4] "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" [7] "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" [10] "2007-06-10 01:03:34 BRT" > # deu certo!! O que pode estar ocorrendo? Não basta o meu vetor original ser do tipo character? Abraços, Paulo

1 0

Converter JSON
by Elias Carvalho 03 Mar '16

03 Mar '16

Ola Pessoal Eu preciso transformar um json em formato tabular. Alguem pode ajudar ? -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* _____________________________________________ Prof. Elias César Araújo de Carvalho CV: http://lattes.cnpq.br/4248328961021251

5 5

Merge não funciona
by Wagner Tassinari 02 Mar '16

02 Mar '16

Olá pessoal, tudo bem ? Estou com um dúvida quando executo a função merge. Tenho um banco de dados x com 1.630.328 registros. Tenho um outro banco de dados y com 503.615. Quando irei fazer o relacionamento dessas bases, x e y, com a chave primária chamada ID. w = merge(x, y, by.x=c('ID'), by.y=c('ID'), all.x = T) Neste caso, como fixei o o numero de resgistroa igual ao do banco x, o banco w deveria ter o mesmo número de x, mas isso não ocorre. O banco w fica com 1.641.990. Gostaria de saber o que está ocorrendo, e como posso "blindar" o número de registros de w sendo igual ao do banco x. Muito obrigado mais uma vez pessoal. -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

3 2

Armin - Joinville - Mudança no e-mail
by ARMIN KOENIG 29 Fev '16

29 Fev '16

Prezado(a), Estou migrando para o GMAIL, dada a lentidão de navegação no Yahoo.Peço que alterem meu e-mail para : armin.kng(a)gmail.com Muito obrigado, Armin Koenig

2 1

Artigo que não tem no periodicos capes
by Pedro Emmanuel Alvarenga Americano do Brasil 29 Fev '16

29 Fev '16

Amigos de R, Será que alguem tem na mão que possa me enviar o pdf desse artigo aqui: http://onlinelibrary.wiley.com/doi/10.1002/sim.4780141909/full Grato, Pedro Brasil

5 8

Função para remover colunas
by Elias Carvalho 29 Fev '16

29 Fev '16

Ola Pessoal Para facilitar a vida de um usuário leigo, eu criei uma função para excluir as colunas de uma tabela, conforme código abaixo: 1. remove.variables <- function(dataset.name, vars.to.remove) 2. { 3. # scan the variables to be removed 4. for (x in 1 : length(vars.to.remove)) 5. { 6. # transform command.Remove in "data$variable <- NULL" 7. command.Remove <- paste(dataset.name,"$",vars.to.remove[x], " <- NULL",sep="") 8. # execute the command 9. command.Remove <- eval(parse(text=command.Remove)) 10. } # for (x in 1 : length(vars.to.remove)) 11. } # remove.variables <- function(dataset.name, vars.to.remove) Fazendo um teste local, ou seja, executando o procedimento abaixo com as variáveisdataset.name, vars.to.remove carregadas respectivamente com o nome do meu data frame "data" e com o nome de duas variáveis "var1" e "var5" a serem removidas e executando da linha 2 a 11 tudo tunciona bem e as colunas são removidas. dataset.name <- "data" vars.to.remove<-c("var1", "var5") No entanto se carrego a função e executo os procedimentos abaixo em outro script, a função é executada sem erro, mas não remove nenhuma coluna do data frame: vars.to.remove = c("var1", "var5") dataset.name = "data" remove.variables(dataset.name, vars.to.remove) Eu acredito que tem a ver com o ambiente, acesso global algo assim mas não estou achando a solução. Agradeço se aguem puder me ajudar. -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* _____________________________________________ Prof. Elias César Araújo de Carvalho CV: http://lattes.cnpq.br/4248328961021251

6 5

Colocar identificar único em registros repetidos
by Wagner Tassinari 28 Fev '16

28 Fev '16

Olá pessoal, Tenho um banco com registros repetidos por linha e gostaria de colocar o mesmo identificador (contador) para os registros que são repetidos baseado em alguns atributos. Por exemplo, se o existirem 4 registros com o mesmo nome, nome da mãe e data de nascimento, o identificador será o mesmo. Segue um exemplo do meu banco de dados: Registro Nome Nome Mãe Data Nascimento 1 João Maria 1976-09-08 2 Manuel Selma 1988-03-10 3 João Maria 1976-09-08 4 João Maria 1976-09-08 5 Ricardo Sonia 1990-05-12 6 Marcelo Marta 1987-08-06 7 Ricardo Sonia 1990-05-12 8 João Maria 1976-09-08 9 Ricardo Sonia 1990-05-12 10 Marcelo Marta 1987-08-06 Gostaria obter esse banco com o identificador ID: Registro ID Nome Nome Mãe Data Nascimento 1 1 João Maria 1976-09-08 2 2 João Selma 1988-03-10 3 1 João Maria 1976-09-08 4 1 João Maria 1976-09-08 5 3 Ricardo Sonia 1990-05-12 6 4 Marcelo Marta 1987-08-06 7 3 Ricardo Sonia 1990-05-12 8 1 João Maria 1976-09-08 9 3 Ricardo Sonia 1990-05-12 10 4 Marcelo Marta 1987-08-06 Muito obrigado a todos, -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

3 2

Exportar RSQLite para CSV
by Paulo Dick 26 Fev '16

26 Fev '16

Boa tarde colegas listeiros, Estou usando o RSQLite para manipular uma base grande que gostaria de exportar para CSV usando o R, contudo não consegui descobrir como fazer isto. Tentei combinar write.table com dbGetQuery mas extrapola a memória. Agradeço a ajuda. *Paulo Dick* Estatístico Mestre em Epidemiologia em Saúde Pública Tel.: (55 21) 99591-2716

2 2

Gerar arquivo de saída netCDF no R
by Mateus Dias Nunes 25 Fev '16

25 Fev '16

Olá, eu preciso gerar um arquivo de saída .nc (netcdf) dos resultados da média da coluna de ozônio que carregue as informações de lat, lon e valor em cada ponto de grade da região analisada. Já consegui gerar um arquivo .txt, porém o programa que utilizo não lê arquivos modo texto. abaixo o arquivo de entrada também .nc Tentei utilizar a função create da biblioteca net https://www.dropbox.com/s/ntvmtx7wonwo839/teste.nc?dl=0 library(maps) library(ncdf4) dados <- nc_open('teste.nc') # lendo coordenadas espaço-temporal lat <- ncvar_get( dados, 'lat' ) lon <- ncvar_get( dados, 'lon' ) time <- ncvar_get( dados, 'time' ) # lendo dados coluna total de Ozônio ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' ) # dimensoes da variavel ColumnAmountO3 dims_ColumnAmountO3 <- dim(ColumnAmountO3) # tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de grade X tempo dim(ColumnAmountO3) <- c(dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] ) # calculando a média e retornado-a em 2D media_ColumnAmountO3 <- rowMeans( ColumnAmountO3) dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1],dims_ColumnAmountO3[2] ) # longitude varia de 0 a 360, convertendo para -180 a 180, essaconversão é feita para plotagem sobre o mapa for (i in 1:dim(lon)) { if (lon[i] > 180) { lon[i] <- lon[i]-360 } } desde já agradeço ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 81125154 ____________________________________________________________________________

2 1