Abril 2012 - R-br - DINF UFPR

Soma de quadrados de regressão..
by Julia Maria 19 Abr '12

19 Abr '12

Pessoal, como sempre estou com algumas dúvidas básicas. Eu ainda estou na graduação e estou me adaptando ao R. Será que alguem poderia me tirar uma dúvida? Quando faço uma analise de um modelo de regressão no R, por exemplo, do modelo: m1<-lm(Y~-1+x1+x2+x3+I(x1^2) +I(x2^2)+I(x3^2) (as variáveis estao definidas corretamente em um arquivo a parte) e peço anova ele solta as SQs de cada coeficiente. Gostaria de saber se há alguma forma de obter um quadro de análise de variância da regressão, com a soma de quadrados da regressao linear e da regressao quadratica, ao inves de cada coeficiente.

3 2

Fwd: [paleontropica] Macroevolution in R Short Course: Call for Applicants
by D'Apolito 19 Abr '12

19 Abr '12

We are pleased to announce an intensive short course on using R to perform comparative methods to be held in Santa Barbara on June 11th to June 15th. This course is funded by the National Science Foundation, and a number of stipends to cover or defray travel, room, and board are available to qualified students and post-docs. Topics covered will include an introduction to the R programming language, tree manipulation, independent contrasts and phylogenetic generalized least squares, ancestral state reconstruction, models of character evolution, diversification analyses, and community phylogenetic analysis. If you are interested please submit your CV along with a short (maximum 1 page) description of your research interests, background, and reasons for taking the course. Admission is competitive, and the best applications come from students with data sets to analyze. International applicants are welcome. To apply visit this URL: tinyurl.com/macro-in-R Application deadline: April 30th. Individuals from cultural, racial, linguistic, geographic and socioeconomic backgrounds that are currently underrepresented in science are especially encouraged to apply. Please contact the co-organizers, Michael Alfaro (michaelalfaro(a)ucla.edu) and Luke Harmon (lukeh(a)uidaho.edu) with any questions. -- David Bapst Dept of Geophysical Sciences University of Chicago 5734 S. Ellis Chicago, IL 60637 http://home.uchicago.edu/~dwbapst/ http://cran.r-project.org/web/packages/paleotree/index.html _______________________________________________ Paleonet mailing list Paleonet(a)nhm.ac.uk http://mailman.nhm.ac.uk/mailman/listinfo/paleonet __._,_.___ Responder a remitente<alc1006@hotmail.com?subject=RE%3A%20Macroevolution%20in%20R%20Short%20Course%3A%20Call%20for%20Applicants>| Responder a grupo<paleontropica@gruposyahoo.com?subject=RE%3A%20Macroevolution%20in%20R%20Short%20Course%3A%20Call%20for%20Applicants>| Responder mediante la Web<http://espanol.groups.yahoo.com/group/paleontropica/post;_ylc=X3oDMTJxbmVxa…>| Crear un tema nuevo<http://espanol.groups.yahoo.com/group/paleontropica/post;_ylc=X3oDMTJmNWo4a…> Mensajes con este tema<http://espanol.groups.yahoo.com/group/paleontropica/message/1633;_ylc=X3oDM…>( 1) Actividad reciente: Visita tu grupo<http://espanol.groups.yahoo.com/group/paleontropica;_ylc=X3oDMTJmMmg2N3E1BF…> [image: Yahoo! Grupos]<http://espanol.groups.yahoo.com/;_ylc=X3oDMTJldGQxdTA0BF9TAzk3NDkwNDU3BGdyc…> Cambiar a: Solo Texto<paleontropica-traditional@gruposyahoo.com?subject=Change+Delivery+Format:+Traditional>, ReseÃ±a Diaria<paleontropica-digest@gruposyahoo.com?subject=Email+Delivery:+Digest>• Cancelar suscripción<paleontropica-unsubscribe(a)gruposyahoo.xn--com?subject=cancelar+suscripcin-99c>• Condiciones de uso <http://e1.docs.yahoo.com/info/utos.html> . __,_._,___ -- _____________________________________ *Carlos D'Apolito Júnior - Intern* Center for Tropical Paleoecology and Archeology-CTPA Smithsonian Tropical Research Institute-STRI Panama City (507) 212-8019 CV: http://lattes.cnpq.br/4557754746424026

1 0

duvida comando read.fwf
by Alexandre Heinemann 18 Abr '12

18 Abr '12

Caros Colegas, Gostaria de saber se alguem poderia me ajudar a entender por que o comando abaixo não funciona? Abaixo segue exemplo com 1 arquivo, mas tenho mais de 100 arquivos para ler. http://www.datafilehost.com/download-7e562cad.html (link para download do arquivo utilizado como exemplo (somente 1 arquivo)) #codigo nasaAirTempFiles <- file.path("NASA", "Files", paste("temperature", 1:1, ".txt",sep="")) head(nasaAirTempFiles) [1] "NASA/Files/temperature1.txt" tempDF <- read.fwf(nasaAirTempFiles[1],skip=7,widths=c(-12, rep(7, 24))) #esse comando não funciona??? Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") : cannot open file 'NASA/Files/temperature1.txt': No such file or directory tempDF <- read.fwf("cloudhigh1.txt",skip=7,widths=c(-12, rep(7, 24))) # qdo substituo o nome do arquivo "nasaAirTempFiles[1]" para "cloudhigh1.txt" o nome original do arquivo (link para download) funciona head(tempDF) V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V20 V21 V22 V23 V24 1 26.0 23 23.0 17.0 19.5 17.0 16.0 16.0 16.0 19.0 18.0 19.0 19.5 19.5 18.0 16.0 14.0 14.0 14.0 18.0 20.5 22.0 25.5 25.5 2 20.0 20 18.5 16.5 18.0 15.0 15.0 15.0 16.0 15.0 15.0 17.5 17.0 17.0 17.0 17.5 17.5 18.0 19.0 19.0 20.0 22.0 23.5 23.5 3 16.0 16 14.0 12.5 13.5 14.5 13.5 13.5 13.5 13.0 13.0 16.5 18.5 20.0 21.0 21.0 21.0 20.5 20.5 19.5 19.0 19.0 19.0 20.5 4 13.0 11 11.0 11.0 14.0 13.5 11.5 12.0 15.0 18.0 18.0 18.5 20.5 20.5 18.5 18.0 17.0 15.0 15.0 14.0 13.0 13.0 13.5 12.5 5 7.5 8 10.5 10.5 15.0 16.5 19.5 20.0 20.5 19.5 17.0 16.0 16.0 16.0 12.5 9.5 8.5 7.5 5.5 4.5 4.5 6.5 6.5 10.0 6 8.0 11 13.0 19.5 19.5 26.0 27.5 25.5 20.0 14.0 9.5 8.5 10.0 8.0 4.5 5.0 5.0 4.5 3.5 2.5 2.5 4.0 8.0 12.5 Att., Alexandre

4 4

Dúvida acessar base de dados
by Marco Pollo 18 Abr '12

18 Abr '12

Caros, Algum direcionamento para eu conseguir acessar um arquivo com extensão .GDB (FireBird Database) direto de uma base de dados. Marco Pollo Almeida

2 1

Re: [R-br] Duvida String
by Sérgio Henrique almeida da silva ju 18 Abr '12

18 Abr '12

---------- Mensagem encaminhada ---------- De: Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> Data: 18 de abril de 2012 11:22 Assunto: Re: [R-br] Duvida String Para: r-br(a)listas.c3sl.ufpr.br Uma alternativa que encontrei foi f = function(v, dat, v1=sapply(gregexpr(dat,v),"[",i=1)){ n = ifelse(v1 > 1, nchar(v), NA) substr(v, v1,n) } dados1$causaa = f(dados1$linhaa, "A39") dados1$causab = f(dados1$linhab, "A39") dados1$causac = f(dados1$linhac, "A39") dados1$causad = f(dados1$linhad, "A39") dados1$causaii = f(dados1$linhaii, "A39") dados1$causa<- with(dados1, gsub('NA', '', paste(dados1$causaa,dados1$causab,dados1$causac,dados1$causad,dados1$causaii, sep=''))) Como faço para buscar mais de um código? Tem que fazer um monte de comando como esses? Existe uma forma mais simples de fazer isso? Segue o banco anexo Desculpe enviar o arquivo anexo, mas é uma tentativa para vocês me ajudarem, já que o dput não é uma boa saida! 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) 94429486/78101651 id: 123*20942 -- 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) 94429486/78101651 id: 123*20942

2 3

Agregar colunas
by Sérgio Henrique almeida da silva ju 18 Abr '12

18 Abr '12

Olá Estou tentando agregar três colunas em uma, porém quando faço isso os NAs também se agregam, como faço para descosiderá-los Ex: # Dados tmp = data.frame(x=c('A1', NA, 'A2', NA, 'A3', NA, NA), y=c(NA, NA, NA, 'A4', NA, 'A5', NA), z=c(NA, 'A6', NA, NA, NA, NA, 'A7')) # Agregando com o PASTE tmp$W <-paste(tmp$x,tmp$y,tmp$z, sep = "") # Como os dados ficam dput(tmp) structure(list(x = structure(c(1L, NA, 2L, NA, 3L, NA, NA), .Label = c("A1", "A2", "A3"), class = "factor"), y = structure(c(NA, NA, NA, 1L, NA, 2L, NA), .Label = c("A4", "A5"), class = "factor"), z = structure(c(NA, 1L, NA, NA, NA, NA, 2L), .Label = c("A6", "A7"), class = "factor"), W = c("A1NANA", "NANAA6", "A2NANA", "NAA4NA", "A3NANA", "NAA5NA", "NANAA7")), .Names = c("x", "y", "z", "W"), row.names = c(NA, -7L), class = "data.frame") # Gostaria que minha coluna W ficasse assim e não como acima: W = c("A1", "A6", "A2", "A4", "A3", "A5", "A7") Como faço 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) 94429486/78101651 id: 123*20942

2 3

Configuração de e-mails - Lista
by Danilo Cuzzuol Pedrini 17 Abr '12

17 Abr '12

Prezados, Gostaria de saber se é possível escolhermos de receber apenas um relatório diário sobre o que foi enviado à Lista? Pelo menos para mim, a caixa de e-mails fica cheia de assuntos os quais não posso ajudar e/ou tenho interesse. A coisa fica pior quando me ausento uns 4 dias do e-mail. Atenciosamente, Danilo Cuzzuol Pedrini www.danilocp.blogspot.com

2 1

Comparar duas regressões no mesmo plot
by Marcelo Claro de Souza 17 Abr '12

17 Abr '12

Plotei a seguinte regressão hipotética e gostaria de saber se existe alguma forma se verificar se as reg 1 e 2 são significativamente diferentes, talvez comparando ângulo de inclinação das retas? Tb gostaria de saber como faço p inserir na legenda o que significa cada reta, pois só consegui colocar o que significam os pontos. x1=c(1,2,3,4,5,6,7,8,9,10) x2=c(13,14,15,16,17,18,19) y1=c(30,40,50,60,70,80,90,100,110,120) y2=c(20,22,23,24,25,26,27) x=c(1,2,3,4,5,6,7,8,9,10,13,14,15,16,17,18,19) y=c(30,40,50,60,70,80,90,100,110,120,20,22,23,24,25,26,27) reg1=lm(x1~y1) reg2=lm(x2~y2) simbolos=c(rep(19,10),rep(2,7)) plot(x~y,pch=simbolos) abline(reg1,lty=2) abline(reg2,lty=1) legend("bottomright",c("D","W"),pch=c(19,2)) Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

ajuda com sample
by DANIELA CUNHA DA SE 17 Abr '12

17 Abr '12

Olá Pessoal, Gostaria de saber se alguém pode me ajudar: Quero sortear uma amostra dentro de um data.frame com varias colunas.....porém quero resgatar da amostra selecionada as informações das outras colunas. Exemplo a1<-sample(dados$parc,10) [1] 21 51 59 3 30 65 10 2 13 47 exemplo a primeira unidade selecionada é a parcela 21 (da coluna dados$parc do data.frame), como relacionar a essas parcelas os valores das outras colunas? Obrigada pela atenção.

2 1

Ajuda com mapas
by Renato Braga 17 Abr '12

17 Abr '12

Prezados, Sou iniciante com mapas: vocês podem me encaminhar/indicar material de estudo? Como plotar as bordas de um mapa, definir limites (municípios, por exemplo)? Obrigado,

2 1