- R-br - DINF UFPR

Re: [R-br] SOMAR LINHAS
by sznelwar＠uol.com.br 08 Mai '17

08 Mai '17

2 1

Re: [R-br] Reproduzir gráfico 3D no R
by sznelwar＠uol.com.br 07 Mai '17

07 Mai '17

1 0

Reproduzir gráfico 3D no R
by Fernando de Souza Bastos 07 Mai '17

07 Mai '17

Prezados, alguém tem alguma sugestão de como resolver essa questão: https://pt.stackoverflow.com/questions/202827/reproduzir-gr%C3%A1fico-3d-co… Desde já agradeço! Att, Fernando de Souza Bastos Professor Assistente Universidade Federal de Viçosa (UFV) Campus UFV - Florestal Doutorando em Estatística Universidade Federal de Minas Gerais (UFMG) Cel: (31) 99751-6586

2 2

Re: [R-br] SOMAR LINHAS
by Fernando Antonio de souza 07 Mai '17

07 Mai '17

Mauro Realmente havia um erro no comando que enviei. Peço desculpas por isso Faltou fechar o parenteses do for (Em vermelho). Agora deve funcionar Obrigado pelo retorno Abrcs dados <- structure(list(id = c(49L, 49L, 49L, 64L, 64L, 64L, 23L, 23L, 23L, 23L), evento1 = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L), .Label = c("A", "D", "E"), class = "factor"), nomes = structure(c(3L, 3L, 3L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("Anicuns", "Goiânia", "Santo Antônio do Descoberto"), class = "factor"), cons_jan_15 = c(9033.36, NA, NA, NA, NA, 4091.54, 6833.28, NA, NA, NA), valor_jan_15 = c(4848.67, NA, NA, NA, NA, 3076.74, 3772.8, NA, NA, NA), cons_fev_15 = c(13131.13, NA, NA, NA, NA, 5482.15, 8706.81, NA, NA, NA), valor_fev_15 = c(4320.57, NA, NA, NA, NA, 2243.94, 2961.72, NA, NA, NA), cons_mar_15 = c(8121.07, NA, NA, NA, NA, 3380.04, 5524.23, NA, NA, NA), valor_mar_15 = c(4410.8, NA, NA, NA, NA, 2033.13, 3285.14, NA, NA, NA), cons_abr_15 = c(8562.85, NA, NA, NA, NA, 3003.66, 6157.68, NA, NA, NA), valor_abr_15 = c(6241.59, NA, NA, NA, NA, 2250.64, 4381.01, NA, NA, NA), cons_mai_15 = c(8261.5, NA, NA, NA, NA, 3601.44, 5997.78, NA, NA, NA), valor_mai_15 = c(5586.92, NA, NA, NA, NA, 2550.23, 4178.95, NA, NA, NA)), .Names = c("id", "evento1", "nomes", "cons_jan_15", "valor_jan_15", "cons_fev_15", "valor_fev_15", "cons_mar_15", "valor_mar_15", "cons_abr_15", "valor_abr_15", "cons_mai_15", "valor_mai_15"), class = "data.frame", row.names = c(NA, -10L)) install.package("plyr") library(plyr) dados <- transform(dados,id=factor(id),evento1=factor(evento1),nomes= factor(nomes)) adicLinha <- function(x) { soma<-numeric() for(i in 1:ncol(x)*)* { if(is.numeric(x[,i])){ soma[i] <-sum(x[,i],na.rm=TRUE) } } return(rbind(x,soma)) } ddply(dados,.(id,evento1,nomes),adicLinha) Em 7 de maio de 2017 14:01, Fernando Antonio de souza < nandodesouza(a)gmail.com> escreveu: > Estranho Mauro, porque aqui funcionou completamente > > Enviou um vídeo (anexo) para comprovação > > O erro está indicando um colchete ({), de uma olhada se ao copiar e colar > o último colchete, foi copiado. > > Se puder gravar sua tela , como fiz, e enviar podemos tentar ver o que > está acontecendo. > > Mas aqui, funcionou perfeitamente > > A disposição > > > > Em 7 de maio de 2017 13:15, <sznelwar(a)uol.com.br> escreveu: > >> Tentei rodar e não consegui: >> >> > dados <- transform(dados,id=factor(id),evento1=factor(evento1),nomes= >> factor(nomes)) >> > >> > adicLinha <- function(x) { >> + soma<-numeric() >> + for(i in 1:ncol(x){ >> Erro: '{' inesperado in: >> " soma<-numeric() >> for(i in 1:ncol(x){" >> > >> >> >> >> Utilizando os dados organizados pelo Rodrigo. Segue a minha solução para >> seu problema de forma reproduzível. >> >> Abraço a todos >> >> dados <- structure(list(id = c(49L, 49L, 49L, 64L, 64L, 64L, 23L, 23L, >> 23L, 23L), evento1 = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L), >> .Label = c("A", "D", "E"), class = "factor"), nomes = structure(c(3L, 3L, >> 3L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("Anicuns", "Goiânia", "Santo >> Antônio do Descoberto"), class = "factor"), cons_jan_15 = c(9033.36, NA, >> NA, NA, NA, 4091.54, 6833.28, NA, NA, NA), valor_jan_15 = c(4848.67, NA, >> NA, NA, NA, 3076.74, 3772.8, NA, NA, NA), cons_fev_15 = c(13131.13, NA, NA, >> NA, NA, 5482.15, 8706.81, NA, NA, NA), valor_fev_15 = c(4320.57, NA, NA, >> NA, NA, 2243.94, 2961.72, NA, NA, NA), cons_mar_15 = c(8121.07, NA, NA, NA, >> NA, 3380.04, 5524.23, NA, NA, NA), valor_mar_15 = c(4410.8, NA, NA, NA, NA, >> 2033.13, 3285.14, NA, NA, NA), cons_abr_15 = c(8562.85, NA, NA, NA, NA, >> 3003.66, 6157.68, NA, NA, NA), valor_abr_15 = c(6241.59, NA, NA, NA, NA, >> 2250.64, 4381.01, NA, NA, NA), cons_mai_15 = c(8261.5, NA, NA, NA, NA, >> 3601.44, 5997.78, NA, NA, NA), valor_mai_15 = c(5586.92, NA, NA, NA, NA, >> 2550.23, 4178.95, NA, NA, NA)), .Names = c("id", "evento1", "nomes", >> "cons_jan_15", "valor_jan_15", "cons_fev_15", "valor_fev_15", >> "cons_mar_15", "valor_mar_15", "cons_abr_15", "valor_abr_15", >> "cons_mai_15", "valor_mai_15"), class = "data.frame", row.names = c(NA, >> -10L)) >> >> install.package("plyr") >> >> >> library(plyr) >> >> dados <- transform(dados,id=factor(id),evento1=factor(evento1),nomes= >> factor(nomes)) >> >> adicLinha <- function(x) { >> soma<-numeric() >> for(i in 1:ncol(x){ >> >> if(is.numeric(x[,i])){ >> >> soma[i] <-sum(x[,i],na.rm=TRUE) >> >> >> } >> >> >> } >> return(rbind(x,soma)) >> } >> >> ddply(dados,.(id,evento1,nomes),adicLinha) >> >> >> >> Em 3 de maio de 2017 14:03, jose gomes de souza filho Gomes via R-br < >> r-br(a)listas.c3sl.ufpr.br >> <http://../../../undefined//compose?to=r-br@listas.c3sl.ufpr.br>> >> escreveu: >> >>> Boa Tarde, >>> >>> Tenho a seguinte base ( uma parte, são várias linhas): >>> >>> bd<-"id;evento1;nomes;cons_jan_15 ;valor_jan_15;cons_fev_15;valo >>> r_fev_15;cons_mar_15;valor_mar_15 ;cons_abr_15;valor_abr_15;cons >>> _mai_15;valor_mai_15 >>> 49;A;Santo Antônio do Descoberto;9033.36;4848.67;131 >>> 31.13;4320.57;8121.07;4410.8;8562.85;6241.59;8261.5;5586.92 >>> 49;A;Santo Antônio do Descoberto;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA >>> 49;A;Santo Antônio do Descoberto;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA >>> 64;D;Goiânia;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA >>> 64;D;Goiânia;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA >>> 64;D;Goiânia;4091.54;3076.74;5482.15;2243.94;3380.04;2033.13 >>> ;3003.66;2250.64;3601.44;2550.23 >>> 23;E;Anicuns;6833.28;3772.8;8706.81;2961.72;5524.23;3285.14; >>> 6157.68;4381.01;5997.78;4178.95 >>> 23;E;Anicuns;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA >>> 23;E;Anicuns;NA;NA ;NA;NA;NA;NA;NA;NA;NA;NA >>> 23;E;Anicuns;NA;NA ;NA;NA;NA;NA;NA;NA;NA;NA" >>> dads <- read.table(header = TRUE, text=bd, sep=";") >>> dads >>> >>> Quero somar as linhas, de modo que: >>> >>> id >>> >>> evento1 >>> >>> nomes >>> >>> cons_jan_15 >>> >>> valor_jan_15 >>> >>> cons_fev_15 >>> >>> valor_fev_15 >>> >>> cons_mar_15 >>> >>> valor_mar_15 >>> >>> cons_abr_15 >>> >>> valor_abr_15 >>> >>> cons_mai_15 >>> >>> valor_mai_15 >>> >>> 49 >>> >>> A >>> >>> Santo Antônio do Descoberto >>> >>> 9033,36 >>> >>> 4848,67 >>> >>> 13131,13 >>> >>> 4320,57 >>> >>> 8121,07 >>> >>> 4410,8 >>> >>> 8562,85 >>> >>> 6241,59 >>> >>> 8261,5 >>> >>> 5586,92 >>> >>> 49 >>> >>> A >>> >>> Santo Antônio do Descoberto >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 49 >>> >>> A >>> >>> Santo Antônio do Descoberto >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 49 >>> >>> A >>> >>> Santo Antônio do Descoberto >>> >>> 9033,36 >>> >>> 4848,67 >>> >>> 13131,13 >>> >>> 4320,57 >>> >>> 8121,07 >>> >>> 4410,8 >>> >>> 8562,85 >>> >>> 6241,59 >>> >>> 8261,5 >>> >>> 5586,92 >>> >>> 64 >>> >>> D >>> >>> Goiânia >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 64 >>> >>> D >>> >>> Goiânia >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 64 >>> >>> D >>> >>> Goiânia >>> >>> 4091,54 >>> >>> 3076,74 >>> >>> 5482,15 >>> >>> 2243,94 >>> >>> 3380,04 >>> >>> 2033,13 >>> >>> 3003,66 >>> >>> 2250,64 >>> >>> 3601,44 >>> >>> 2550,23 >>> >>> 64 >>> >>> D >>> >>> Goiânia >>> >>> 4091,54 >>> >>> 3076,74 >>> >>> 5482,15 >>> >>> 2243,94 >>> >>> 3380,04 >>> >>> 2033,13 >>> >>> 3003,66 >>> >>> 2250,64 >>> >>> 3601,44 >>> >>> 2550,23 >>> >>> 23 >>> >>> E >>> >>> Anicuns >>> >>> 6833,28 >>> >>> 3772,8 >>> >>> 8706,81 >>> >>> 2961,72 >>> >>> 5524,23 >>> >>> 3285,14 >>> >>> 6157,68 >>> >>> 4381,01 >>> >>> 5997,78 >>> >>> 4178,95 >>> >>> 23 >>> >>> E >>> >>> Anicuns >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 23 >>> >>> E >>> >>> Anicuns >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 23 >>> >>> E >>> >>> Anicuns >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 0 >>> >>> 23 >>> >>> E >>> >>> Anicuns >>> >>> 6833,28 >>> >>> 3772,8 >>> >>> 8706,81 >>> >>> 2961,72 >>> >>> 5524,23 >>> >>> 3285,14 >>> >>> 6157,68 >>> >>> 4381,01 >>> >>> 5997,78 >>> >>> 4178,95 >>> >>> >>> Obrigado pela ajuda, >>> >>> >>> >>> >>> >>> >>> _______________________________________________ >> >> > > > -- > ========================================= > Fernando Souza > Zootecnista, DSc. Produção e Alimentação Animal > Celular: (31)99796-8781 (Vivo) > E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> > Lattes: http://lattes.cnpq.br/6519538815038307 > Blog: https://producaoanimalcomr.wordpress.com/ > ========================================== > -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

1 0

SOMAR LINHAS
by jose gomes de souza filho Gomes 06 Mai '17

06 Mai '17

Boa Tarde, Tenho a seguinte base ( uma parte, são várias linhas): bd<-"id;evento1;nomes;cons_jan_15 ;valor_jan_15;cons_fev_15;valor_fev_15;cons_mar_15;valor_mar_15 ;cons_abr_15;valor_abr_15;cons_mai_15;valor_mai_15 49;A;Santo Antônio do Descoberto;9033.36;4848.67;13131.13;4320.57;8121.07;4410.8;8562.85;6241.59;8261.5;5586.92 49;A;Santo Antônio do Descoberto;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA 49;A;Santo Antônio do Descoberto;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA 64;D;Goiânia;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA 64;D;Goiânia;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA 64;D;Goiânia;4091.54;3076.74;5482.15;2243.94;3380.04;2033.13;3003.66;2250.64;3601.44;2550.23 23;E;Anicuns;6833.28;3772.8;8706.81;2961.72;5524.23;3285.14;6157.68;4381.01;5997.78;4178.95 23;E;Anicuns;NA;NA;NA;NA;NA;NA;NA;NA;NA;NA 23;E;Anicuns;NA;NA ;NA;NA;NA;NA;NA;NA;NA;NA 23;E;Anicuns;NA;NA ;NA;NA;NA;NA;NA;NA;NA;NA" dads <- read.table(header = TRUE, text=bd, sep=";") dads Quero somar as linhas, de modo que: id evento1 nomes cons_jan_15 valor_jan_15 cons_fev_15 valor_fev_15 cons_mar_15 valor_mar_15 cons_abr_15 valor_abr_15 cons_mai_15 valor_mai_15 49 A Santo Antônio do Descoberto 9033,36 4848,67 13131,13 4320,57 8121,07 4410,8 8562,85 6241,59 8261,5 5586,92 49 A Santo Antônio do Descoberto 0 0 0 0 0 0 0 0 0 0 49 A Santo Antônio do Descoberto 0 0 0 0 0 0 0 0 0 0 49 A Santo Antônio do Descoberto 9033,36 4848,67 13131,13 4320,57 8121,07 4410,8 8562,85 6241,59 8261,5 5586,92 64 D Goiânia 0 0 0 0 0 0 0 0 0 0 64 D Goiânia 0 0 0 0 0 0 0 0 0 0 64 D Goiânia 4091,54 3076,74 5482,15 2243,94 3380,04 2033,13 3003,66 2250,64 3601,44 2550,23 64 D Goiânia 4091,54 3076,74 5482,15 2243,94 3380,04 2033,13 3003,66 2250,64 3601,44 2550,23 23 E Anicuns 6833,28 3772,8 8706,81 2961,72 5524,23 3285,14 6157,68 4381,01 5997,78 4178,95 23 E Anicuns 0 0 0 0 0 0 0 0 0 0 23 E Anicuns 0 0 0 0 0 0 0 0 0 0 23 E Anicuns 0 0 0 0 0 0 0 0 0 0 23 E Anicuns 6833,28 3772,8 8706,81 2961,72 5524,23 3285,14 6157,68 4381,01 5997,78 4178,95 Obrigado pela ajuda,

4 6

R-Ladies Rio
by Celina Rebello 03 Mai '17

03 Mai '17

Pessoal, Aproveitando para divulgar: Qui, 18 de mai, 19:00 Faculdades Senac Rio Rua Santa Luzia, 735 - 4° andar - sala 402 - Próximo ao Metrô da Cinelândia., Rio de Janeiro Evento é gratuito, temos poucas vagas. Inscrições em: https://goo.gl/AOsCQp Quem estiver pelo Rio, vale muito ir.

1 0

Re: [R-br] Ajuda para elaboração de boxplot em painel
by Gilenio Borges Fernandes 03 Mai '17

03 Mai '17

Prezad@s Estou com dificuldade de elaborar um boxplot por grupo (ou painel), no formato do que está no arquivo pdf em anexo a esse E-mail. Segue o CMR do banco de dados (Documentação do SAS) O arquivo pdf pode ser acessado em https://www.dropbox.com/s/5q65geuq7e9m5ey/Dados%26grafico-DocSAS.pdf?dl=0 Desde já agradeço Parts=structure(list(Machine = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("A386", "A455", "C334"), class = "factor"), Sample = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L, 20L, 20L, 20L, 20L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 23L, 23L, 23L, 23L, 24L, 24L, 24L, 24L), Diam = c(4.32, 4.55, 4.16, 4.44, 4.49, 4.3, 4.52, 4.61, 4.44, 4.32, 4.25, 4.5, 4.55, 4.15, 4.42, 4.49, 4.21, 4.3, 4.29, 4.63, 4.56, 4.61, 4.29, 4.56, 4.63, 4.3, 4.41, 4.58, 4.38, 4.65, 4.43, 4.44, 4.12, 4.49, 4.3, 4.36, 4.45, 4.56, 4.38, 4.51, 4.62, 4.67, 4.7, 4.58, 4.33, 4.23, 4.34, 4.58, 4.29, 4.38, 4.28, 4.41, 4.15, 4.35, 4.28, 4.23, 4.21, 4.3, 4.32, 4.38, 4.16, 4.28, 4.31, 4.59, 4.14, 4.18, 4.08, 4.21, 4.51, 4.2, 4.28, 4.19, 4.1, 4.33, 4.37, 4.47, 3.99, 4.09, 4.47, 4.25, 4.24, 4.54, 4.43, 4.38, 4.23, 4.48, 4.31, 4.57, 4.27, 4.4, 4.32, 4.56, 4.7, 4.65, 4.49, 4.38)), .Names = c("Machine", "Sample", "Diam"), class = "data.frame", row.names = c(NA, -96L)) Gilênio Borges Fernandes Professor Associado IV (Aposentado) Professor Adjunto A (Substituto) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 Skype: gilenio.fernandes Lattes: http://lattes.cnpq.br/6764860618464860 Em 3 de maio de 2017 15:58, Gilenio Borges Fernandes <gilenio(a)ufba.br> escreveu: > Prezad@s > > Estou com dificuldade de elaborar um boxplot por grupo (ou painel), no > formato do que está no arquivo pdf em anexo a esse E-mail. Segue o CMR do > banco de dados (Documentação do SAS) > > Desde já agradeço > > > Parts=structure(list(Machine = structure(c(1L, 1L, 1L, 1L, 1L, 1L, > > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, > > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > > 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, > > 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, > > 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("A386", "A455", > > "C334"), class = "factor"), Sample = c(1L, 1L, 1L, 1L, 2L, 2L, > > 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, > > 6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 10L, > > 10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L, > > 13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, > > 16L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L, > > 20L, 20L, 20L, 20L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 23L, > > 23L, 23L, 23L, 24L, 24L, 24L, 24L), Diam = c(4.32, 4.55, 4.16, > > 4.44, 4.49, 4.3, 4.52, 4.61, 4.44, 4.32, 4.25, 4.5, 4.55, 4.15, > > 4.42, 4.49, 4.21, 4.3, 4.29, 4.63, 4.56, 4.61, 4.29, 4.56, 4.63, > > 4.3, 4.41, 4.58, 4.38, 4.65, 4.43, 4.44, 4.12, 4.49, 4.3, 4.36, > > 4.45, 4.56, 4.38, 4.51, 4.62, 4.67, 4.7, 4.58, 4.33, 4.23, 4.34, > > 4.58, 4.29, 4.38, 4.28, 4.41, 4.15, 4.35, 4.28, 4.23, 4.21, 4.3, > > 4.32, 4.38, 4.16, 4.28, 4.31, 4.59, 4.14, 4.18, 4.08, 4.21, 4.51, > > 4.2, 4.28, 4.19, 4.1, 4.33, 4.37, 4.47, 3.99, 4.09, 4.47, 4.25, > > 4.24, 4.54, 4.43, 4.38, 4.23, 4.48, 4.31, 4.57, 4.27, 4.4, 4.32, > > 4.56, 4.7, 4.65, 4.49, 4.38)), .Names = c("Machine", "Sample", > > "Diam"), class = "data.frame", row.names = c(NA, -96L)) > > > > Gilênio Borges Fernandes > Professor Associado IV (Aposentado) > Professor Adjunto A (Substituto) > Universidade Federal da Bahia > Instituto de Matemática e Estatística > Departamento de Estatística > Av. Adhemar de Barros, s/n – Ondina. > 40.170-110 - Salvador - BA, Brasil > Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 > Skype: gilenio.fernandes > Lattes: http://lattes.cnpq.br/6764860618464860 >

2 1

geoR > Krige.bayes com trend e como obter DIC
by Wagner Wolff 03 Mai '17

03 Mai '17

Olá pessoal da lista! Estou com dificuldade de implementar na função krige.bayes() o atributo trend. Primeiramente, por uma abordagem frequentista escolhi um modelo de trend, quando tento usar o mesmo na função krige.bayes() ocorre alguns erros com questões matriciais. Alguem tem alguma dica para corrigir isso? Outra questão cé omo obtendo o DIC para os modelos que testarei na função krige.bayes(). Por exemplo, modelo com trend e componentes espaciais, modelo só com trend e modelo só com os componentes espacias. A figura abaixo mostra um exemplo de um artigo que fez essas comparações utilizando binom.krige.bayes(). Fonte: Mapping, Bayesian Geostatistical Analysis and Spatial Prediction of Lymphatic Filariasis Prevalence in Africa -- *Wagner Wolff, **PhD* "*Luiz de Queiroz**" College of Agriculture,* University of São Paulo Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil Phone: +55 19 982385582 <+55%2019%2098238-5582> http://orcid.org/0000-0003-3426-308X https://github.com/wwolff7 http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1

1 0

Comparações Múltiplas
by Isabela M. 03 Mai '17

03 Mai '17

Boa tarde, Meu nome é Isabela, faço iniciação científica e trabalho com análise multivariada, mais especificamente com a MANOVA. Alguém poderia me ajudar com script do R, para realizar a comparação múltipla multivariada, no caso para 5 variáveis? Desde já agradeço!

5 7

Banco de Dados
by jose gomes de souza filho Gomes 02 Mai '17

02 Mai '17

Boa tarde, Tem o seguinte banco de dados *cidade* *parcela2* *parcela3* *ano* *item* *quantidade* *GOIÂNIA* *Ana* *serv* *2011* *baixados* *82331* *GOIÂNIA* *Paulo* *terc* *2011* *baixados* *15733* *ANÁPOLIS* *Roberto* *serv* *2011* *baixados* *9989* *GOIÂNIA* *Ana* *serv* *2011* *Pendentes* *7282* *GOIÂNIA* *Paulo* *terc* *2011* *Pendentes* *6819* *ANÁPOLIS* *Roberto* *serv* *2011* *Pendentes* *5420* *GOIÂNIA* *Ana* *serv* *2011* *novos* *4635* *GOIÂNIA* *Paulo* *terc* *2011* *novos* *4270* *ANÁPOLIS* *Roberto* *serv* *2011* *novos* *4125* *GOIÂNIA* *Ana* *serv* *2012* *baixados* *82332* *GOIÂNIA* *Paulo* *terc* *2012* *baixados* *15734* *ANÁPOLIS* *Roberto* *serv* *2012* *baixados* *9990* *GOIÂNIA* *Ana* *serv* *2012* *Pendentes* *7283* *GOIÂNIA* *Paulo* *terc* *2012* *Pendentes* *6820* *ANÁPOLIS* *Roberto* *serv* *2012* *Pendentes* *5421* *GOIÂNIA* *Ana* *serv* *2012* *novos* *4636* *GOIÂNIA* *Paulo* *terc* *2012* *novos* *4271* *ANÁPOLIS* *Roberto* *serv* *2012* *novos* *4126* Gostaria que ficasse assim: *cidade* *parcela2* *parcela3* *baixados_2011* *baixados_2012* *Pendentes_2011* *Pendentes_2012* *novos_2011* *novos_2012* *GOIÂNIA* *Ana* *serv* *82331* *82332* *7282* *7283* *4635* *4636* *GOIÂNIA* *Paulo* *terc* *15733* *15734* *6819* *6820* *4636* *4271* *ANÁPOLIS* *Roberto* *serv* *9989* *9990* *5420* *5421* *4125* *4126* Obrigado.

3 2