Agosto 2013 - R-br

Re: [R-br] Selecionar variáveis conforme linhas
by Éder Comunello 27 Ago '13

27 Ago '13

Roney, bom dia! Não deve ser a melhor forma, mas o código abaixo pode ser um começo... # <BEGIN> #install.packages('repmis') library(repmis) FinURL <- paste0(" https://dl.dropboxusercontent.com/u/61883020/2013-08-25-r-br.csv") data <- repmis::source_data(FinURL, sep = ",", header = TRUE) dim(data); head(data) data$new = 0; data$prop = 0; data$medlnk= 0 str(data) for (i in 1:length(data[,1])){ # para cada linha col = data[i,2]+2 # coluna do grupo dado dado = data[i,col] # dado para a coluna calculada nlink = length(data[1,3:6]>0) # núm de grupos com valor (ligações) data[i,8] = dado # aloja o dado data[i,9] = dado/(data[i,7]-dado) # total - dado do grupo data[i,10] = data[i,7]/nlink # total/núm ligações } head(data) #<END> > head(data) pessoa pessoa.grupo g.1 g.2 g.3 g.4 tc new prop medlnk 1 a 3 360 913 1407 674 3354 1407 0.7226502 838.50 2 b 2 235 693 588 384 1900 693 0.5741508 475.00 3 c 4 355 435 285 524 1599 524 0.4874419 399.75 4 e 3 37 129 258 80 504 258 1.0487805 126.00 5 f 1 317 40 14 95 466 317 2.1275168 116.50 6 g 2 17 216 88 27 348 216 1.6363636 87.00 -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

1 0

path analysis - pacote agricolae
by Heloíse Pavanato 27 Ago '13

27 Ago '13

Boa noite, Estou com uma dúvida sobre path analysis no pacote agricolae. No exemplo do manual do pacote se usa a correlação de Pearson para depois computar os coeficientes da análise. Minha dúvida é se posso utliizar correlação de Spearman ao invés, já que meus dados não são normais, gostaria de saber se isso bastaria. Se alguém tiver alguma outra sugestão será muito bem-vinda. O CRM segue abaixo: # dados: y <- c(77, 92, 94, 95, 79, 87, 91, 70, 65, 69, 68, 56, 71, 67, 62, 95, 91, 78, 90, 85, 83, 78, 65, 77, 87, 43, 55, 63, 66, 87, 67, 70, 69, 30, 28, 44, 67, 89, 82, 36, 60, 52, 14, 19, 51, 22, 37, 15, 86, 79, 31, 35, 15, 96, 89, 88, 92, 90, 18, 15, 13, 93, 95, 83, 95, 100, 100, 84, 90, 93, 74, 78, 73, 69, 75, 67, 68, 100, 95, 94, 94, 88, 93, 83, 73, 88, 95, 71, 69, 71, 74, 100, 71, 77, 73, 56, 52, 62, 75, 93, 88, 44, 72, 60, 27, 33, 69, 35, 44, 27, 93, 87, 56, 52, 27, 98, 94, 90, 96, 95, 26, 20, 16, 94, 97, 100) temp <- c(26.26, 26.26, 26.26, 26.26, 26.26, 26.26, 26.26, 25.83, 25.83, 25.83, 25.83, 25.83, 25.83, 25.83, 25.83, 26.40, 26.40, 26.40, 26.40, 26.40, 26.40, 26.40, 26.40, 26.13, 26.13, 26.13, 26.13, 26.13, 26.13, 26.13, 26.13, 25.60, 25.60, 25.60, 25.60, 25.60, 25.60, 25.60, 25.60, 24.96, 24.96, 24.96, 24.96, 24.96, 24.96, 24.96, 24.96, 25.30, 25.30, 25.30, 25.30, 25.30, 25.30, 25.30, 25.30, 26.76, 26.76, 26.76, 26.76, 26.76, 26.76, 26.76, 26.76, 24.36, 24.36, 24.36, 24.36, 24.36, 24.36, 24.36, 23.80, 23.80, 23.80, 23.80, 23.80, 23.80, 23.80, 24.26, 24.26, 24.26, 24.26, 24.26, 24.26, 24.26, 24.26, 23.70, 23.70, 23.70, 23.70, 23.70, 23.70, 23.70, 23.70, 23.06, 23.06, 23.06, 23.06, 23.06, 23.06, 23.06, 23.06, 21.93, 21.93, 21.93, 21.93, 21.93, 21.93, 21.93, 21.93, 22.30, 22.30, 22.30, 22.30, 22.30, 22.30, 22.30, 22.30, 24.00, 24.00, 24.00, 24.00, 24.00, 24.00, 24.00, 24.00, 24.00) pluv <- c(136.6, 136.6, 136.6, 136.6, 136.6, 136.6, 136.6, 136.0, 136.0, 136.0, 136.0, 136.0, 136.0, 136.0, 136.0, 112.3, 112.3, 112.3, 112.3, 112.3, 112.3, 112.3, 112.3, 65.3, 65.3, 65.3, 65.3, 65.3, 65.3, 65.3, 65.3, 66.0, 66.0, 66.0, 66.0, 66.0, 66.0, 66.0, 66.0, 65.0, 65.0, 65.0, 65.0, 65.0, 65.0, 65.0, 65.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 53.0, 254.6, 254.6, 254.6, 254.6, 254.6, 254.6, 254.6, 256.6, 256.6, 256.6, 256.6, 256.6, 256.6, 256.6, 231.3, 231.3, 231.3, 231.3, 231.3, 231.3, 231.3, 231.3, 142.6, 142.6, 142.6, 142.6, 142.6, 142.6, 142.6, 142.6, 126.6, 126.6, 126.6, 126.6, 126.6, 126.6, 126.6, 126.6, 118.0, 118.0, 118.0, 118.0, 118.0, 118.0, 118.0, 118.0, 110.6, 110.6, 110.6, 110.6, 110.6, 110.6, 110.6, 110.6, 57.6, 57.6, 57.6, 57.6, 57.6, 57.6, 57.6, 57.6, 57.6) # shapiro-wilk shapiro.test(y) shapiro.test(x$temperatura) # não é normal! shapiro.test(x$pluviosidade) x <- data.frame(temp,pluv) # path analysis require(agricolae) cor.y <- correlation(y,x)$correlation cor.x <- correlation(x)$correlation path.analysis(cor.x,cor.y) Grata, Heloise

1 0

Fábio, let's connect on LinkedIn
by Iago Carvalho Cunha via LinkedIn 26 Ago '13

26 Ago '13

LinkedIn ------------ Iago Carvalho Cunha requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Iago Carvalho Cunha http://www.linkedin.com/e/-33g73l-hku2i88j-z/vhj1xG42iBdUFK5d_HMEGL4iFMC511… View profile of Iago Carvalho Cunha http://www.linkedin.com/e/-33g73l-hku2i88j-z/rso/282741063/Q3Ds/name/827774… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-hku2i88j-z/plh/http%3A%2F%2Fhelp%2Elinked… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Mapa com as Regiões de Minas Gerais
by Alisson Lucrecio 24 Ago '13

24 Ago '13

Caro colegas r-br, Boa tarde. Eu estou precisando de uma mapa com as coordenas da regiões de Minas Gerais e de preferencia com as principais cidades, Algum de você saberiam me informar? Obrigado. Alisson Lucrécio da Costa

3 3

Extrair string usando padrão de caracteres
by Thiago V. dos Santos 23 Ago '13

23 Ago '13

Olá pessoal, Digamos que eu tenha uma lista com os nomes de vários arquivos: i <- c("./maize-yield.csv", "./rice-yield.csv", "./soybean-yield.csv", "./sugarcane-yield.csv", "./wheat-yield.csv") e eu precise extrair os caracteres após a barra ("/") e antes do hífen ("-") , ou seja, o nome das culturas. Qual comando ou combinação de comandos eu devo usar para chegar ao seguinte resultado? > ii [1] "maize" "rice" "soybean" "sugarcane" "wheat" Agradece desde já, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

3 4

get coordinates of patches (clumps) in a raster map
by Ludmila Rattis 21 Ago '13

21 Ago '13

Prezados, Preciso extrair as coordenadas dos centroides de diferentes clumps de um raster. Para isso, usei o seguinte código: p <- data.frame(rasterToPoints(meu_raster)) a<- p[p$layer > 0,] c<-sapply(split(p[, c("x", "y")], p$layer), colMeans) Contudo, esse código não está funcionando. Alguém saberia me dizer como posso obter tais coordenadas? []'s *Ludmila Rattis* Programa de Pós-graduação em Ecologia/UNICAMP Conservation Biogeography Lab http://www.wix.com/rdloyola/lab <http://sites.ffclrp.usp.br/ficus>

2 6

Livros ou materiais sobre VarSelRF
by Alisson Lucrecio 21 Ago '13

21 Ago '13

Caro Colegas, Eu estou precisando estudar sobre a função varselRF para seleções de variáveis, assim estou precisando de algum livro ou material. Você teriam algo que poderiam me recomendar? Obrigado. Alisson Lucrécio da Costa

1 0

comparação codigo R
by alanarocha＠sapo.pt 21 Ago '13

21 Ago '13

Bom dia, eu tenho o seguinte excerto de dados: Clínica 1s09 1s10 1s11 1s12 1s13 1 Arcos de Valdevez 70,5% 73,9% 75,7% 81,0% 77,9% 2 Braga 72,5% 79,9% 77,3% 73,4% 81,8% 3 Fafe 74,7% 74,7% 71,1% 74,4% 80,0% 4 Maia 76,7% 72,7% 77,3% 78,1% 80,5% 5 VNGaia 84,6% 84,3% 83,4% 81,8% 73,3% 6 SMFeira 73,0% 73,8% 76,7% 79,4% 68,5% pretendo comparar para saber qual o melhor em média, comparar tambem seus desvios padroes, variancias, medianas e amplitudes. Qual a melhor função para fazer isto e se posso comparar os 5 semestres juntos ou se tem de ser 2 a 2. obrigada. cumprimentos Ana Rocha

1 0

GIS - matriz de distância entre polígonos
by Ludmila Rattis 20 Ago '13

20 Ago '13

Oi pessoal, Estou tentando criar uma matriz de distância entre todos os polígonos de um shape. Alguém tem ideia de qual pacote ou função posso usar? [ ]'s *Ludmila Rattis* Programa de Pós-graduação em Ecologia/UNICAMP Conservation Biogeography Lab http://www.wix.com/rdloyola/lab <http://sites.ffclrp.usp.br/ficus>

2 1

ExpDes
by Natalia Martins 20 Ago '13

20 Ago '13

Boa tarde prezados, estou fazendo uma anova com 3 fatores e bloco fat3.rbd(hibrido, populacao, zona, bloco, Prod, sigT = 0.1, sigF = 0.1) No entanto, a analise nao 'e finalizada pois o teste de shapiro nao 'e aplicado para amostras maiores que 5000. Por favor o que devo fazer para utilizar esse pacote? Att -- Natália da Silva Martins Contato: (19) 8306-4743

4 8