Junho 2015 - R-br - DINF UFPR

Curvas sobrevivência - Comparação
by Andre Oliveira 06 Jun '15

06 Jun '15

Pessoal,este é um exemplo fictício apenas para exemplificação! require(survival)tempo <- c(1, 18, 21, 21, 22, 29, 29, 35, 37, 39, 40, 50, 52, 54, 60, 80, 80, 81, 83, 84, 8)status <- c(1,0,0,0,0,0,0,1,1,0,1,0,0,0,0,0,0,0,0,0,0)grupos=c("A","F","F","F","M","A","M","A","M","F","M","A","M","F","F","F","M","F","M","F","A")km <- survfit(Surv(tempo,status) ~ grupos) # Estratificação summary(km)plot(km) survdiff(Surv(tempo, status)~grupos, rho = 1) # Teste Peto A 10% as curvas diferem, mas como comparar as curvas 2 a 2 para ver tais diferenças? A survival faz isto? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 4

Mapa / dengue
by Luiz Roberto Martins Pinto 06 Jun '15

06 Jun '15

Caros, Itabuna-BA é um dos municípios brasileiros com maior incidência de dengue. Eu desejo elaborar, periodicamente, mapas de ocorrência de dengue, para avaliar a sua evolução. Alguém da lista pode me auxiliar? Luiz Roberto. Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia-Brasil luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

6 6

Distribuir o resultado de duas condições de ifelse em duas colunas
by ASANTOS 04 Jun '15

04 Jun '15

Bom dia Pessoal, Tenho um vetor de dados artificiais chamado dados e gostaria de testar se os valores de dados são <=0.05767 e jogar o resultado (sendo que se true =1, senão 0) em uma coluna e ao mesmo tempo testar se x<=0.89000 e colocar o resultados (sendo que se true =1, senão 0) em outra coluna, para isso fiz: dados=rnorm(10) a1=0.05767 a2=0.89000 RES<-NULL for (i in 1:length(dados)){ I1 <- x[i] x[i] <- as.vector(x[i]) ifelse(x[i] <= a1, 1, 0) I2 <- x[i] x[i] <- as.vector(x[i]) ifelse (x[i] <= a2, 1, 0) RES<-rbind(RES,c(I1,I2)) } RES Mas estou fazendo algo errado, pois são estou obtendo os resultados esperados nas colunas I1 e I2, alguém poderia me ajudar? Obrigado -- ====================================================================== Alexandre dos Santos Proteção 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 Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

5 7

Shapefile sobre image.kriging
by Felinto COSTA 03 Jun '15

03 Jun '15

Prezados. Sabem alguma maneira de plotar um arquivo em formato shape sobre a imagem gerada pela função image.kriging (geoR)? A imagem que disponho já foi delimitada anteriormente na krigagem (em krige.cov, borders) mas o arquivo que desejo sobrepor possui delimitações internas (bairros de uma cidade) e o anterior apenas o perímetro externo da idade. Grato Felinto

2 1

equação no scatter3d
by Natalia Martins 03 Jun '15

03 Jun '15

Prezados boa tarde. Estou tentando colocar um modelo de regressao juntamente com o R^2 em um scatter3d, mas a equação sai como texto. Usei dois comandos "mtext3d" e "text3d" , mas não consegui inserir o modelo com nenhum. Se alguem puder me ajudar, desde já meu obrigada! scatter3d(y ~ x1 + x2) mtext3d(expression(paste(hat(y),'=7.5 + 1.59x[1] + 8.68x[2], ', R^2,'=0,9747')) text3d(expression(paste(hat(y),'=7.5 + 1.59x[1] + 8.68x[2], ', R^2,'=0,9747')) Natalia Martins

1 0

Otimizando desempenho do Servidor - dados grandes
by Pedro Emmanuel Alvarenga Americano do Brasil 02 Jun '15

02 Jun '15

Amigos de R, Eu venho com um dúvida de iniciante, porque nunca trabalhei assim. Eu não entendo muito de informática por isso estou sem luz no caminho. Bom, a minha questão é que me envolvi num projeto em que os dado são na caso de centenas de milhares de linhas, com umas 150 colunas. Estou trabalhando numa área de trabalho remota num servidor windows 2008. Eu não sei quais as especificaçõs do servidor, mas essa área esta servidor está servindo somente para esse fim. Então gostaria de uma luz ou do que ler para otimizar o tempo de processamento nas ediçoes e nas análises. Algum tutorial ou alguma dica que aproveite melhor o meu tempo de espera. O que ja fucei que é pouco... > memory.size(max = FALSE) [1] 1044.11 > memory.size(max = T) [1] 1814.94 > memory.limit(size = NA) [1] 4095 > Abraço, Pedro Brasil

3 4

Método monte carlo
by Fernando Antonio de souza 02 Jun '15

02 Jun '15

Caros amigos, Ajustei um modelo de regressão através de uma meta analise e gostaria de avaliar a sua capacidade preditiva. Infelizmente não tenho um banco de dados independente(não utilizado para ajustar modelo) para poder compara valores preditos e observados. Penso que o método de monte Carlo pode me ajudar a resolver este problema (conforme algumas leituras aleatorias pela internet e) . Gostaria de uma opinião sobre o uso desse método para esse fim. E se possível dicas ou compartilhamento de material didático( for dummies) para que eu possa introduzir ao estudo de monte Carlo. Att

2 6

Impressão do objeto Date
by David Feitosa 01 Jun '15

01 Jun '15

Boa noite! Estou deparando-me com um pequeno problema: Tenho um data.frame que desejo cortar de acordo com o valor de uma coluna. Esta coluna está em formato de data com hora: YYYY-MM-DD HH:MM:SS Minha intenção é de criar um diretório com a data e um arquivo dentro de cada diretório. Ex: <dir>/data/2015-02-02/saida.csv , <dir>/data/2015-02-03/saida.csv , <dir>/data/2015-02-04/saida.csv etc Entretanto, o que está sendo impresso, e, criado como diretório, é a representação interna do objeto Date Ex: <dir>/data/16467/saida.csv , <dir>/data/16468/saida.csv , <dir>/data/16469/saida.csv etc Já tentei usar as.character(), format(), mas estou ficando sem ideias de como "corrigir". Alguém já passou por isso ou teria alguma sugestão por favor? Segue trecho de código pertinente com um pedaço dos dados: library(R.utils) lines <- " 5928671,14,-3.831290,-38.50146,2015-02-01 00:00:01 5926098,44,-3.709390,-38.58986,2015-02-02 00:01:13 5915217,145,-3.743671,-38.48788,2015-02-03 00:01:16 5444641,174,-3.766271,-38.57452,2015-02-01 00:01:38 2250561,977,-3.738654,-38.48400,2015-02-02 00:01:47 5671754,14,-3.831290,-38.50146,2015-02-03 00:02:43 " con <- textConnection(lines) df <- read.csv(con, header = F) close(con) for (dia in unique(as.Date(df[,5]))){ dfi <- df[as.Date(df[,5]) == dia,] diretorio <- paste(getwd(), "rsc/data", format(dia, format = "%Y-%m-%d"), sep = "/", collapse = "") mkdirs(diretorio) write.csv(dfi, paste(diretorio, "saida.csv", sep = "/", collapse = "")) } Dados do ambiente: R version 3.2.0 (2015-04-16) -- "Full of Ingredients" Platform: i386-w64-mingw32/i386 (32-bit) Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

5 6