- R-br - DINF UFPR

sobre abóboras
by Antonio Silva 09 Nov '16

09 Nov '16

Olá Hoje atializei o R. A nova versão é a Sincere Pumpkin Patch Vem de um diálogo entre o Linus e a Sally Brown. Outras versões anteriores também faziam referência às tirinhas do Snoopy e do Charlie Brown (não sei quando começou) Linus: He'll come here because I have the most *sincere pumpkin patch* and he respects sincerity. Sally Brown: Do you really think he will come? Linus: Tonight the Great Pumpkin will rise out of the pumpkin patch. He flies through the air and brings toys to all the children of the world. Sally Brown: That's a good story. Linus: You don't believe the story of the Great Pumpkin? I thought little girls always believed everything that was told to them. I thought little girls were innocent and trusting. Sally Brown: Welcome to the 20th century! Procurei na rede algo sobre o motivo desta homenagem, mas não achei. Alguém sabe algo sobre? Obrigado -- Antônio Olinto

2 1

Re: [R-br] Heterocedasticidade e teste de média
by Luiz Henrique Leal 09 Nov '16

09 Nov '16

Prezado César, toda ajuda é bem vinda pois estou com muita dificuldade em modelar esses dados por causa das violações. Endereçei aos outros colegas porque eles deram sugestões. O tratamento 0 é o controle e quero compará-lo com os outros tratamentos l. E ainda tem o problema da heteroscedasticidade. Todos os dados que trabalho tem essa característica. Abs Luiz Conectado pela Motorola Cesar Rabak <cesar.rabak(a)gmail.com> escreveu: >Luiz, > > >Malgrado você tenha endereçado a pergunta, dou meu "pitaco" enquanto os outros colegas não se manifestam. > > >Todas as regressões necessitam que em última análise os erros aleatórios possam ser modelados por um processo cuja melhor representação é a distribuição gaussiana com média zero e alguma dispersão em torno desse zero (que gera o "ruído" nas respostas versus a modelagem matemática da regressão). > > >A grande contribuição dos modelos generalizados é identificação de diversas funções de ligação que permitem modelar processos com respostas diferentes da média pura (ligação == identidade) de tal forma que quando você analisa os resíduos eles "ficam" normais (gaussianos). > > >No script que você enviou em anexo ao e-mail, você está aplicando testes com a suposição ("modelagem") de um processo que deveria ter uma relação "linear" entre x e y. > > >Ora, você tem três amostras para todos os grupos "x" com exceção do grupo "0" que tem 14 casos, como você espera obter alguma coisa razoável em relação aos resíduos? > > >Uma análise exploratória dos seus dados mostra algo interessante, veja no seu ambiente: > > >> library(gplots) > >> plotmeans(y ~ x, data=dados) > > >Uma outra instrutiva visualização dos seus dados é a seguinte: > > >> dotchart(dados$y, groups=dados$x) > > >HTH > >-- > >Cesar Rabak > > > >2016-11-08 11:11 GMT-02:00 Luiz Leal via R-br <r-br(a)listas.c3sl.ufpr.br>: > >Prezados Walmes e Wagner > >Costumo trabalhar com conjuntos de dados em que todos os pressupostos são violados (os dados e os resíduos desviam muito da normalidade além de heterocedasticidade que conforme o Walmes mencionou a solução transformação Box-Cox não resolve). Além disso não posso excluir ou agrupar os tratamentos porque os mesmos representam resultados (medições) de laboratórios. Preciso compará-los com um tratamento (laboratório) considerado de "referência". > >O modelo gls() é capaz de abarcar esses problemas? > >Desde já agradeço > >Luiz > > > >On Monday, November 7, 2016 10:34 PM, Mauro Sznelwar via R-br <r-br(a)listas.c3sl.ufpr.br> wrote: > > > >Ele não reconhece estes comandos 'mcglm','gof', tem alguma bibliotéca para isto? > > > >Caros, > > >Alguém postou esse conjunto de dados com problema de pressupostos, principalmente heterocedasticidade. Agora a pouco veio outro e-mail com um problema similar. Fiz um exemplo um pouco mais detalhado de como isso pode > >ser facilmente resolvido e mostrando o efeito disso no modelo. > ># Example 2 ------------------------------------------------------------ >Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, > 328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, > 342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, > 88.571, 88.571, 90.252, 41.513, 52.437, 49.076, 88.571, > 88.571, 90.252, 64.202, 60.000, 61.681) >Cor = factor(c(rep("ambar",6), rep("ambar_claro",3), rep("ambar",6), > rep("ambar_claro",6),rep("branco",6), > rep("extra_ambar_claro",3),rep("branco",3))) > ># Exploratory analysis >boxplot(Fenois ~ Cor) >tapply(Fenois, Cor, sd) >dados <- data.frame(Fenois, Cor) >dados$id <- 1 > ># Fitting >fit1 <- mcglm(c(Fenois ~ Cor), list(mc_id(dados)), data = dados) >fit2 <- mcglm(c(Fenois ~ Cor), list(mc_dglm(~ Cor, id = "id", data = dados)), > covariance = "expm", data = dados) ># Goodness-of-fit >gof(fit1) >gof(fit2) > ># Comparing estimates and standard errors >coef(fit1, type = "beta", std.error = TRUE) >coef(fit2, type = "beta", std.error = TRUE) > >O interessante é que a estimativa pontual é exatamente a mesma, porém >olha a enorme diferença nos erros padrões dos betas. > > >-- > >Wagner Hugo Bonat >---------------------------------------------------------------------------------------------- >Department of Mathematics and Computer Science (IMADA) >University of Southern Denmark (SDU) and >Laboratório de Estatística e Geoinformação (LEG) >Universidade Federal do Paraná (UFPR) > >_______________________________________________ > > >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne� c�igo m�imo reproduz�el. > > > > >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. > >

1 0

Re: [R-br] Heterocedasticidade e teste de média
by sznelwar＠uol.com.br 08 Nov '16

08 Nov '16

3 2

Re: [R-br] Regressão linear - Excluindo outliers
by sznelwar＠uol.com.br 07 Nov '16

07 Nov '16

2 1

Re: [R-br] heterocedasticidade
by sznelwar＠uol.com.br 07 Nov '16

07 Nov '16

1 0

Re: [R-br] Heterocedasticidade e teste de média
by sznelwar＠uol.com.br 05 Nov '16

05 Nov '16

2 1

Heterocedasticidade e teste de média
by Wagner Bonat 04 Nov '16

04 Nov '16

Caros, Alguém postou esse conjunto de dados com problema de pressupostos, principalmente heterocedasticidade. Agora a pouco veio outro e-mail com um problema similar. Fiz um exemplo um pouco mais detalhado de como isso pode ser facilmente resolvido e mostrando o efeito disso no modelo. # Example 2 ------------------------------------------------------------ Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, 328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, 342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, 88.571, 88.571, 90.252, 41.513, 52.437, 49.076, 88.571, 88.571, 90.252, 64.202, 60.000, 61.681) Cor = factor(c(rep("ambar",6), rep("ambar_claro",3), rep("ambar",6), rep("ambar_claro",6),rep("branco",6), rep("extra_ambar_claro",3),rep("branco",3))) # Exploratory analysis boxplot(Fenois ~ Cor) tapply(Fenois, Cor, sd) dados <- data.frame(Fenois, Cor) dados$id <- 1 # Fitting fit1 <- mcglm(c(Fenois ~ Cor), list(mc_id(dados)), data = dados) fit2 <- mcglm(c(Fenois ~ Cor), list(mc_dglm(~ Cor, id = "id", data = dados)), covariance = "expm", data = dados) # Goodness-of-fit gof(fit1) gof(fit2) # Comparing estimates and standard errors coef(fit1, type = "beta", std.error = TRUE) coef(fit2, type = "beta", std.error = TRUE) O interessante é que a estimativa pontual é exatamente a mesma, porém olha a enorme diferença nos erros padrões dos betas. -- Wagner Hugo Bonat ---------------------------------------------------------------------------------------------- Department of Mathematics and Computer Science (IMADA) University of Southern Denmark (SDU) and Laboratório de Estatística e Geoinformação (LEG) Universidade Federal do Paraná (UFPR)

3 2

heterocedasticidade
by Luiz Leal 04 Nov '16

04 Nov '16

PrezadosNos casos em que a transformação box-cox não é capaz remover a heterocedasticidade há alguma outra solução? Eu utilizei a sugestão do Leg/ UFPR nos meus dados, mas não funcionou.(exemplo: <http://www.leg.ufpr.br/Rpira/Rpira/node17.html>) Segue em anexo o script e os dados. Desde já agradeço Luiz

2 2

Sorteios aleatórios de um cadastro (AAS ou amostra estratificada)
by Andre Oliveira 02 Nov '16

02 Nov '16

Bom dia, estou buscando informações de como realizar sorteios aleatórios de um cadastro (AAS ou amostra estratificada) com objetivo de inferir para média ou para a proporção. Alguma dica de pacote para tal trabalho? obg 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

mutate vs within
by Leonard Mendonça de Assis 02 Nov '16

02 Nov '16

Boa tarde a todos Surgiu uma dúvida conceitual aqui: Eu tenho um data.frame que precisa unir com vários outros, além de fazer várias transformações. Estou utilizando a seguinte sequencia para fazer tudo bhz.data <- read_csv2( "Belo Horizonte 20161027.csv", col_types=colunas) %>% filter(Filtro == 1 & Filtro.Votar==1) %>% select(-Filtro, -Filtro.Votar) %>% left_join(label.Sexo,by=c("Sexo"="value")) %>% select(-Sexo) %>% left_join(label.Faixa.Idade, by=c("Faixa.Idade"="value")) %>% select(-Faixa.Idade) %>% left_join(label.Renda, by=c("Renda"="value")) %>% select(-Renda) %>% left_join(label.classe, by=c("Renda.Classe"="value")) %>% select(-Renda.Classe) %>% left_join(label.Escolaridade, by=c("Escolaridade"="value")) %>% select(-Escolaridade) %>% left_join(label.Etinia, by=c('Etinia'="value")) %>% select(-Etinia) %>% left_join(label.Religião, by=c("Religião"="value")) %>% select(-Religião) %>% left_join(label.Trabalha, by=c("Trabalha"="value")) %>% select(-Trabalha) %>% mutate(Area.1 = ifelse(is.na(Area.1),0,Area.1), Area.2 = ifelse(is.na(Area.2),0,Area.2), Area.3 = ifelse(is.na(Area.3),0,Area.3)) %>% mutate(area.Educação = ifelse(Area.1 == 1,1,0) + ifelse(Area.2 == 1,1,0) + ifelse(Area.3 == 1,1,0), area.Habitação = ifelse(Area.1 == 2,1,0) + ifelse(Area.2 == 2,1,0) + ifelse(Area.3 == 2,1,0), area.Limpeza = ifelse(Area.1 == 3,1,0) + ifelse(Area.2 == 3,1,0) + ifelse(Area.3 == 3,1,0), area.Obras.Mobilidade = ifelse(Area.1 == 4,1,0) + ifelse(Area.2 == 4,1,0) + ifelse(Area.3 == 4,1,0), area.Obras.Pavimentação = ifelse(Area.1 == 5,1,0) + ifelse(Area.2 == 5,1,0) + ifelse(Area.3 == 5,1,0), area.Saneamento = ifelse(Area.1 == 6,1,0) + ifelse(Area.2 == 6,1,0) + ifelse(Area.3 == 6,1,0), area.Saúde = ifelse(Area.1 == 7,1,0) + ifelse(Area.2 == 7,1,0) + ifelse(Area.3 == 7,1,0), area.Transpoorte = ifelse(Area.1 == 8,1,0) + ifelse(Area.2 == 8,1,0) + ifelse(Area.3 == 8,1,0), area.Conservação = ifelse(Area.1 == 9,1,0) + ifelse(Area.2 == 9,1,0) + ifelse(Area.3 == 9,1,0), area.Acessibilidade = ifelse(Area.1 == 10,1,0) + ifelse(Area.2 == 10,1,0) + ifelse(Area.3 == 10,1,0), area.Combate.Mosquito = ifelse(Area.1 == 11,1,0) + ifelse(Area.2 == 11,1,0) + ifelse(Area.3 == 11,1,0), area.MeioAmbiente = ifelse(Area.1 == 12,1,0) + ifelse(Area.2 == 12,1,0) + ifelse(Area.3 == 12,1,0), area.Social = ifelse(Area.1 == 13,1,0) + ifelse(Area.2 == 13,1,0) + ifelse(Area.3 == 13,1,0), area.Esporte.Lazer = ifelse(Area.1 == 14,1,0) + ifelse(Area.2 == 14,1,0) + ifelse(Area.3 == 14,1,0), area.Geração.Empregos = ifelse(Area.1 == 15,1,0) + ifelse(Area.2 == 15,1,0) + ifelse(Area.3 == 15,1,0)) %>% select(-Area.1, -Area.2, -Area.3) A questão é a seguinte: No caso da variável área, que é recodificada de 3 para 15 variáveis, quando eu substituo esta sequencia enorme de mutate por um within, as variáveis ficam em uma ordem diferente da que eu preciso. A questão é: como deixar o código mais limpo nesta situação, visto que esta alternativa funciona, mas parece ser meio gambiarra? Leonard

3 3

Teste de independência dos resíduos
by Maurício Sangiogo 02 Nov '16

02 Nov '16

Bom dia listeiros! Estou iniciando no R, e gostaria de pedir ajuda a vocês sobre uma dúvida... Como eu faço no R para testar o pressuposto da anova de "independência dos resíduos"? Eu li algumas citações ao uso do Durbin-watson test para análise de resíduos em regressões, porém não sei se o mesmo se encaixaria (ou se ao menos está correto) para análise de experimentos com fatores qualitativos (em DBC ou DIC). Também vi algo em plotar um gráfico dos resíduos, mas não entendi o como. .. Se alguma boa alma estatística puder me ajudar, muito obrigado, se não puderem agradeço igualmente. Att .................. Eng. Agr. Maurício Sangiogo TEL:(55) 96822030

4 3

comparação de médias - problema com os pressupostos
by Maurício Lordêlo 02 Nov '16

02 Nov '16

Caros, Preciso fazer um teste de comparação de médias mas estou tendo problemas com os pressupostos. Tentei uma transformação Box Cox mas não consegui êxito. Alguma sugestão? Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, 328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, 342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, 88.571, 88.571, 90.252, 41.513, 52.437, 49.076, 88.571, 88.571, 90.252, 64.202, 60.000, 61.681) Cor = factor(c(rep("ambar",6),rep("ambar_claro",3),rep("ambar",6),rep("ambar_claro",6),rep("branco",6), rep("extra_ambar_claro",3),rep("branco",3))) Fenois; Cor boxplot(Fenois~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) tapply(Fenois, Cor, mean) tapply(Fenois, Cor, sd) tapply(Fenois, Cor, var) #Teste de homogeneidade de variâncias bartlett.test(Fenois ~ Cor) #Gráficos para avaliar a variável resposta qqnorm(Fenois) qqline(Fenois) hist(Fenois) boxplot(Fenois) #Transformação Box-Cox require(MASS) boxcox(Fenois~Cor, lambda = seq(-2, 2, 1/10)) boxcox(Fenois~Cor, lambda = seq(0.5, 1.5, 1/100)) #Optando por lambda = 0.6 Fenois2 = ((Fenois^0.6) - 1)/0.6 par(mfrow=c(1,2)) boxplot(Fenois~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) boxplot(Fenois2~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) bartlett.test(Fenois2 ~ Cor)

3 2

Sorteios aleatórios de um cadastro (AAS ou amostra estratificada)
by Andre Oliveira 02 Nov '16

02 Nov '16

Bom dia, estou buscando informações de como realizar sorteios aleatórios de um cadastro (AAS ou amostra estratificada) com objetivo de inferir para média ou para a proporção. Alguma dica de pacote para tal trabalho? obg 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

1 0

Barplot
by Emerson Cotta Bodevan 01 Nov '16

01 Nov '16

Prezados, bom dia. Estou tentando construir um barplot horizontal tal que as barras representem a frequencia de uso de dez medicamentos. Gostaria de colocar os nomes de todos os medicamentos a esquerda das barras. Não estou conseguindo colocar todos os nomes. Somente alguns deles aparecem. Alguma sugestão? Agradeço qualquer ajuda. Att., *Emerson*

4 4

Dúvida para gerar gráfico
by Paulo Henrique Pimenta 31 Out '16

31 Out '16

Olá Listeiros, Tenho uma dúvida bem simples da qual gostaria de pedir ajuda a vocês. Criei um código que gera as evapotranspirações pelo método de Penman, Thronthwait e Blaney-Criddle. A partir desses resultados gostaria de gerar gráficos em dois modos: O primeiro seria pegando o resultado de Penman mostrando sua evolução com os meses do ano O segundo é pegando os métodos de Thornthwait e Blaney-Criddle juntos mostrando, também suas evoluções no ano. Sempre que eu tento plotar algo aparece no meu IDE (RStudio) uma série de caixas. Com este conjunto de dados, os .csv que possuo não conisgo criar os dados desejados. Os dados estão em: http://s000.tinyupload.com/?file_id=49893669914373744636 O código está em: http://s000.tinyupload.com/?file_id=96211791048541413747 Abraços a todos e obrigado desde já. Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435.

3 3

Re: [R-br] Teste de independência dos resíduos
by Walmes Zeviani 31 Out '16

31 Out '16

Maurício, O primeiro gráfico até pode mostrar não independência quando ela é induzida por falta de ajuste. O código abaixo ilustra minha afirmação. library(lattice) library(latticeExtra) data(wtloss, package = "MASS") str(wtloss) xyplot(Weight ~ Days, data = wtloss) # Ajuste de um modelo linear. m0 <- lm(Weight ~ Days, data = wtloss) # Indica falta de ajuste e não independência. plot(m0, which = 1) # Mas eles exibem uma não independência devido à falta de ajuste. plot(acf(residuals(m0))) # Ajuste do modelo quadrático. m1 <- lm(Weight ~ poly(Days, degree = 2), data = wtloss) # Ausência de falta de ajuste. plot(m1, which = 1) # E sem correlação seria. plot(acf(residuals(m1))) start <- list(th0 = 185, th1 = -100, th2 = 140) with(start, { plot(Weight ~ Days, data = wtloss) curve(th0 + th1 * (1 - 2^(-Days/th2)), xname = "Days", add = TRUE, col = 2) }) # Ajuste do modelo não linear. m2 <- nls(Weight ~ th0 + th1 * (1 - 2^(-Days/th2)), data = wtloss, start = start) with(as.list(coef(m2)), { plot(Weight ~ Days, data = wtloss) curve(th0 + th1 * (1 - 2^(-Days/th2)), xname = "Days", add = TRUE, col = 2) }) # Valores preditos de cada modelo. pred <- with(wtloss, data.frame(Days = seq(min(Days), max(Days), length.out = 30))) pred$y0 <- predict(m0, newdata = pred) pred$y1 <- predict(m1, newdata = pred) pred$y2 <- predict(m2, newdata = pred) # Resíduos e valores ajustados de cada modelo. L <- lapply(list(m0, m1, m2), function(m) { data.frame(r = residuals(m), f = fitted(m)) }) names(L) <- 0:2 rf <- plyr::ldply(L, .id = "model") # Gráfico da curva predita com os valores observados. p1 <- xyplot(y0 + y1 + y2 ~ Days, outer = TRUE, type = "l", layout = c(NA, 1), data = pred) + layer(panel.points(x = wtloss$Days, y = wtloss$Weight)) # Gráfico dos resíduos vs valores ajustados. p2 <- xyplot(r ~ f | model, data = rf, layout = c(NA, 1), type = c("p", "smooth")) # Todos numa mesma "página". c(p1, p2, layout = c(NA, 2)) À disposição. Walmes.

1 0

Passar linhas de comandos do DOS para o script do R
by Adriele Giaretta Biase 28 Out '16

28 Out '16

Como executar um script do “R” com uma linha de comando do DOS passando parâmetros para o script. Tenho duas perguntas para essa dúvida: 1) Como executar uma linha de comando do DOS, passando parâmetros para o script do R? 2) Como reconhecer as linhas de comando com os parâmetros dentro do script do R? -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

4 3

boxplot comparativos
by Elisa Henning 23 Out '16

23 Out '16

Boa noite Gostaria de construir boxplots múltiplos comparativos. *Exemplificando:* tenho dois arquivos com dados climáticos referentes a duas cidades. Gostaria de comparar, mês a mês, por meio de boxplots, as temperaturas destas duas cidades. Alguém poderia indicar algum material que auxilie a fazer estes gráficos? Desde já agradeço Elisa

4 4

não foi possível encontrar a função "%!in%"
by Elias Carvalho 23 Out '16

23 Out '16

Pessoal Executando uma função do pacote bn ocorreu o seguinte erro: check.bn.strength(strength) Error in check.bn.strength(strength) : não foi possível encontrar a função "%!in%" Alguém sabe porque dá esse erro no R ? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1

servidor cloud para rodar o R + shiny
by Marcio B 20 Out '16

20 Out '16

Olá pessoal boa tarde. Estou pesquisando um serviço de cloud para utilizar o R + Shiny. Minha duvida esta no custo e por isso gostaria de saber se alguém pode me ajudar para entender melhor: Quem é a empresa (cloud) que esta usando? AWS, Azure, Google ou IBM Qual o tipo de serviço para configurar o R + Shiny? Nestes servidores oferecem vários tipos, o que mais encontro na web é o EC2 - AWS, como vários tutoriais de instalação e configuração Se puder nos informar o uso e o custo deste serviço? Principalmente no uso de Memoria e Storage de dados, algum parâmetro para temos uma noção. Somos uma startup e queremos oferecer nos serviços utilizando o R + Shiny. Porem gostaríamos de um parâmetro para tomar a decisão de configurar este serviço e começar a utilizar. Espero que fui claro e toda a ajuda será bem vinda. Obrigado desde já. Marcio. -- Marcio B. Cel.: 19 99648.3949 tim 99798.0104 vivo skype: maborbaa http://www.linkedin.com/in/marciobar https://www.facebook.com/marciobar Empresa http://www.techagr.com/ ==== PAZ E BEM ====

6 8

geração de variáveis aleatórias normais com estrutura de correlação
by Adriele Giaretta Biase 19 Out '16

19 Out '16

Olá pessoal, tenho uma dúvida com relação à geração de variáveis aleatórias normais multivariadas no R. Eu gostaria de gerar duas variáveis aleatórias (x1, x2) usando distribuição normal multivariada com estrutura de correlação entre elas. A Estrutura da matriz de correlação foi estimada antes com base num banco de dados reais (com correlação fraca de - 0.2). Porém, essas variáveis não podem ser negativas, Ex. x1 tem média e variância, respectivamente, iguais a 0.0067 e 0.0017; x2 tem média e variância, respectivamente, iguais a 0.1374 e 0.0024. Eu gero as variáveis da seguinte forma: *Programa usado no R:* # Simulando os parâmetros com estrutura de correlação entre eles n <- 1000 # tamanho da amostra gerada p <- 2 # numero de variáveis a serem geradas library(MASS) # construíndo a matriz de correlação para usar nas simulações, baseadas na característica da amostra coletada mu<-rep(0, times = p) rho <- - 0.2 # correlacao negativa - 0.2 sigma2 <- 1 Sigma <- sigma2 * ((1-rho)*diag(p)+rho*matrix(1, p, p)) X1 <- 0.0096 # Media da variavel x1 X2 <- 0.1203 # Media da variavel x2 media <- c(X1,X2) y <- mvrnorm(n, media, Sigma) cor(y) apply(y, 2, mean) y [1,] 0.309910452 1.0642521083 [2,] -0.251583312 1.8909032058 [3,] 1.330362012 -1.1239501814 [4,] -0.793399464 -1.5433056284 [5,] 2.165144843 -0.2645184534 [6,] 0.532777085 1.0910864562 [7,] -1.612135390 2.0489354648 [8,] -0.430529913 1.0312062602 ... Quando gero as simulações para x1 e x2 usando a função *mvrnorm*, o resultado me retorna alguns valores negativos para as variáveis, isso não poderia ocorrer. Teria alguma outra função em que eu possa fornecer a variância de cada uma dessas variáveis, além da estrutura de correlação? Como poderia contornar essa situação, usando o R? *Obs:* No SAS, a seguinte programação funciona perfeitamente: *proc* *iml*; wrksize=*100000*; K=*2*; N=*1000*; /* tamanho da amostra */ M={*0* *0*}; S={ *1* -*0.5283*, -*0.5283* *1* }; X=shape(*0*,K,N); ME=*0*; SI=*1*; DO I=*1* TO K; DO J=*1* TO N; if I>*1* then do; ME=M[I]+(S[*1*:I-*1*,I])`*(inv(S[*1*:I-*1*,*1*:I-*1*])*(X[*1*:I-*1* ,J]-M[*1*:I-*1*])); SI=S[I,I]-(S[*1*:I-*1*,I])`*(inv(S[*1*:I-*1*,*1*:I-*1*])*(S[*1*:I- *1*,I])); end; X[I,J]=ME+NORMAL(*0*)*SQRT(SI); END; END; Z=t(X); varnames='X1':'X2'; create NOVO from Z[colname=varnames]; append from Z; *quit*; *data* MEXT; set NOVO; options ps=*66* ls=*75*; X1=*0.0067*+*0.0017**X1; /* normal */ X2=*0.1374*+*0.0024**X2; /* normal */ *proc* *corr* data=MEXT; var X1 X2; *run*; -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

5 5

Re: [R-br] Digest R-br, volume 70, assunto 18
by Adriele Giaretta Biase 18 Out '16

18 Out '16

Amigos, gostaria de agradecer muito a ajuda de vocês (Walmes, Eder e Luiz Roberto). Compreendi que vocês construíram as variáveis normais multivariadas com a estrutura de correlação e variâncias com as características biológicas que eu precisava, usando o R. Muito obrigada mesmo pela contribuição! Adriele. Em 18 de outubro de 2016 17:09, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: geração de variáveis aleatórias normais com estrutura > de correlação (Éder Comunello) > 2. Re: servidor cloud para rodar o R + shiny (Marcio B) > 3. Re: servidor cloud para rodar o R + shiny (Fernando Gama) > 4. Re: geração de variáveis aleatórias normais com estrutura > de correlação (Luiz Roberto Martins Pinto) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 18 Oct 2016 14:52:22 -0300 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: Adriele Giaretta Biase <adrielegbiase(a)gmail.com>, a lista > Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] geração de variáveis aleatórias normais com > estrutura de correlação > Message-ID: > <CABmC8gn0eGcp2ROBMvO+H184T2n1aLNgLsAN92PhXj39RRt+ > UA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Adriele, boa tarde! > > Tentei adaptar seu código R com base no que entendi do script em SAS. Veja > se pode ser útil... > > ### <code r> > library(MASS) > { > set.seed(1357) > n <- 1000 # tamanho da amostra gerada (N) > p <- 2 # numero de variáveis a serem geradas (K) > ME <- rep(1, p) > rho <- -0.2 # correlacao negativa - 0.2 > sigma2 <- 1 > sigma <- sigma2 * ((1-rho)*diag(p)+rho*matrix(1, p, p)) > SI <- mvrnorm(n, ME, sigma) > } > apply(SI, 2, mean) # ~ 1 > # [1] 1.037774 1.017846 > > cor(SI) # ~ sigma > # [,1] [,2] > # [1,] 1.0000000 -0.1973338 > # [2,] -0.1973338 1.0000000 > > mu <- c(0.0067, 0.1374) > mean(0.0017*SI[,1]) # ~ 0.0017 > mean(0.0024*SI[,2]) # ~ 0.0024 > > RES <- data.frame(X1=0.0067+0.0017*SI[,1], X2=0.1374+0.0024*SI[,2]) > > round(cbind(RES=colMeans(RES), INI=mu),4) > # RES INI > # X1 0.0085 0.0067 > # X2 0.1398 0.1374 > > cor(RES) > # X1 X2 > # X1 1.0000000 -0.1973338 > # X2 -0.1973338 1.0000000 > ### </code> > > > ================================================ > Éder Comunello > Researcher at Brazilian Agricultural Research Corporation (Embrapa) > DSc in Agricultural Systems Engineering (USP/Esalq) > MSc in Environ. Sciences (UEM), Agronomist (UEM) > --- > Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| > ================================================ > GEO, -22.2752, -54.8182, 408m > UTC-04:00 / DST: UTC-03:00 > > > > > Em 17 de outubro de 2016 16:05, Adriele Giaretta Biase via R-br < > r-br(a)listas.c3sl.ufpr.br> escreveu: > > > Olá pessoal, > > > > > > tenho uma dúvida com relação à geração de variáveis aleatórias normais > > multivariadas no R. > > > > Eu gostaria de gerar duas variáveis aleatórias (x1, x2) usando > > distribuição normal multivariada com estrutura de correlação entre elas. > A > > Estrutura da matriz de correlação foi estimada antes com base num banco > de > > dados reais (com correlação fraca de - 0.2). Porém, essas variáveis não > > podem ser negativas, Ex. x1 tem média e variância, respectivamente, > iguais > > a 0.0067 e 0.0017; x2 tem média e variância, respectivamente, iguais a > > 0.1374 e 0.0024. Eu gero as variáveis da seguinte forma: > > > > > > > > *Programa usado no R:* > > > > # Simulando os parâmetros com estrutura de correlação entre eles > > > > n <- 1000 # tamanho da amostra gerada > > > > p <- 2 # numero de variáveis a serem geradas > > > > > > > > library(MASS) > > > > # construíndo a matriz de correlação para usar nas simulações, baseadas > na > > característica da amostra coletada > > > > mu<-rep(0, times = p) > > > > rho <- - 0.2 # correlacao negativa - 0.2 > > > > sigma2 <- 1 > > > > Sigma <- sigma2 * ((1-rho)*diag(p)+rho*matrix(1, p, p)) > > > > X1 <- 0.0096 # Media da variavel x1 > > > > X2 <- 0.1203 # Media da variavel x2 > > > > media <- c(X1,X2) > > > > y <- mvrnorm(n, media, Sigma) > > > > cor(y) > > > > apply(y, 2, mean) > > > > y > > > > > > [1,] 0.309910452 1.0642521083 > > > > [2,] -0.251583312 1.8909032058 > > > > [3,] 1.330362012 -1.1239501814 > > > > [4,] -0.793399464 -1.5433056284 > > > > [5,] 2.165144843 -0.2645184534 > > > > [6,] 0.532777085 1.0910864562 > > > > [7,] -1.612135390 2.0489354648 > > > > [8,] -0.430529913 1.0312062602 > > > > ... > > > > > > > > Quando gero as simulações para x1 e x2 usando a função *mvrnorm*, o > > resultado me retorna alguns valores negativos para as variáveis, isso não > > poderia ocorrer. Teria alguma outra função em que eu possa fornecer a > > variância de cada uma dessas variáveis, além da estrutura de correlação? > > Como poderia contornar essa situação, usando o R? > > > > > > > > *Obs:* No SAS, a seguinte programação funciona perfeitamente: > > > > > > *proc* *iml*; > > > > wrksize=*100000*; > > > > K=*2*; > > > > N=*1000*; /* tamanho da amostra */ > > > > M={*0* *0*}; > > > > S={ *1* -*0.5283*, > > > > -*0.5283* *1* }; > > > > X=shape(*0*,K,N); > > > > ME=*0*; SI=*1*; > > > > DO I=*1* TO K; > > > > DO J=*1* TO N; > > > > if I>*1* > > > > then > > > > do; > > > > ME=M[I]+(S[*1*:I-*1*,I])`*(inv(S[*1*:I-*1*,*1*:I-*1*])*( > X[*1*:I- > > *1*,J]-M[*1*:I-*1*])); > > > > SI=S[I,I]-(S[*1*:I-*1*,I])`*(inv(S[*1*:I-*1*,*1*:I-*1*])*(S[*1* > > :I-*1*,I])); > > > > end; > > > > X[I,J]=ME+NORMAL(*0*)*SQRT(SI); > > > > > > > > END; > > > > > > > > END; > > > > Z=t(X); > > > > > > > > varnames='X1':'X2'; > > > > create NOVO from Z[colname=varnames]; > > > > append from Z; > > > > *quit*; > > > > > > > > *data* MEXT; set NOVO; > > > > options ps=*66* ls=*75*; > > > > X1=*0.0067*+*0.0017**X1; /* normal */ > > > > X2=*0.1374*+*0.0024**X2; /* normal */ > > > > > > > > *proc* *corr* data=MEXT; > > > > var X1 X2; > > > > *run*; > > > > -- > > Adriele Giaretta Biase. > > Mestre em Estatística e Experimentação Agropecuária - UFLA. > > Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP > > Contato: (19) 98861-0619. > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > > código mínimo reproduzível. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20161018/8311c3ef/attachment-0001.html> > > ------------------------------ > > Message: 2 > Date: Tue, 18 Oct 2016 16:23:33 -0200 > From: Marcio B <marciobar(a)gmail.com> > To: Raphael Saldanha <rfsaldanha(a)gmail.com> > Cc: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] servidor cloud para rodar o R + shiny > Message-ID: > <CAHL8Jho88ReiZabFeB3PuMwftD-9e8ts1_UWJ3p68s2NLPrXpQ@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > Raphael, obrigado > > Observei a Digital Ocean, mas ela não fornece todas as vantagens que > buscamos. Somos uma empresa de desenvolvimento de software. > > No momento montei um servidor R na Amazon e vamos testar outros ambientes, > como a IBM. > > Mesmo assim agradeço. > > Abs., > > Em 29 de setembro de 2016 14:18, Raphael Saldanha <rfsaldanha(a)gmail.com> > escreveu: > > > Olá Márcio, > > > > Tenho uma experiência muito boa com a Digital Ocean: https://www. > > digitalocean.com/ > > > > 2016-09-29 12:53 GMT-03:00 Marcio B via R-br <r-br(a)listas.c3sl.ufpr.br>: > > > >> Olá pessoal boa tarde. > >> > >> Estou pesquisando um serviço de cloud para utilizar o R + Shiny. > >> > >> Minha duvida esta no custo e por isso gostaria de saber se alguém pode > me > >> ajudar para entender melhor: > >> > >> Quem é a empresa (cloud) que esta usando? AWS, Azure, Google ou IBM > >> Qual o tipo de serviço para configurar o R + Shiny? Nestes servidores > >> oferecem vários tipos, o que mais encontro na web é o EC2 - AWS, como > >> vários tutoriais de instalação e configuração > >> Se puder nos informar o uso e o custo deste serviço? Principalmente no > >> uso de Memoria e Storage de dados, algum parâmetro para temos uma noção. > >> > >> Somos uma startup e queremos oferecer nos serviços utilizando o R + > >> Shiny. Porem gostaríamos de um parâmetro para tomar a decisão de > configurar > >> este serviço e começar a utilizar. > >> > >> Espero que fui claro e toda a ajuda será bem vinda. > >> > >> Obrigado desde já. > >> > >> Marcio. > >> -- > >> Marcio B. > >> > >> Cel.: 19 99648.3949 tim 99798.0104 vivo > >> skype: maborbaa > >> http://www.linkedin.com/in/marciobar > >> https://www.facebook.com/marciobar > >> > >> Empresa > >> http://www.techagr.com/ > >> > >> ==== PAZ E BEM ==== > >> > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br > >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > >> código mínimo reproduzível. > >> > > > > > > > -- > Marcio B. > > Cel.: 19 99648.3949 tim 99798.0104 vivo > skype: maborbaa > http://www.linkedin.com/in/marciobar > https://www.facebook.com/marciobar > > Empresa > http://www.techagr.com/ > > ==== PAZ E BEM ==== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20161018/a5513393/attachment-0001.html> > > ------------------------------ > > Message: 3 > Date: Tue, 18 Oct 2016 15:49:22 -0300 > From: Fernando Gama <f.fabiogama88(a)gmail.com> > To: Marcio B <marciobar(a)gmail.com>, a lista Brasileira oficial de > discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] servidor cloud para rodar o R + shiny > Message-ID: > <CAOwqRysOgYTBitciwDr3VhEpmALmNBEwV8c8jr=LDrqEYskrNw@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > Marcio por questão de curiosidade, > > tive um problema na empresa que trabalhava quanto a uso do shiny haja vista > que os servidores disponíveis não forneciam (na sua versão gratuita) a > disponibilização de um serviço multithread para que os usuários das > aplicações pudessem acessar de forma simultânea os serviços disponíveis... > > Como pretendes lidar com essa limitação já que o valor para habilitar esse > serviço é extremamente elevado? > > Em 18 de outubro de 2016 15:23, Marcio B via R-br < > r-br(a)listas.c3sl.ufpr.br> > escreveu: > > > Raphael, obrigado > > > > Observei a Digital Ocean, mas ela não fornece todas as vantagens que > > buscamos. Somos uma empresa de desenvolvimento de software. > > > > No momento montei um servidor R na Amazon e vamos testar outros > ambientes, > > como a IBM. > > > > Mesmo assim agradeço. > > > > Abs., > > > > Em 29 de setembro de 2016 14:18, Raphael Saldanha <rfsaldanha(a)gmail.com> > > escreveu: > > > >> Olá Márcio, > >> > >> Tenho uma experiência muito boa com a Digital Ocean: > >> https://www.digitalocean.com/ > >> > >> 2016-09-29 12:53 GMT-03:00 Marcio B via R-br <r-br(a)listas.c3sl.ufpr.br > >: > >> > >>> Olá pessoal boa tarde. > >>> > >>> Estou pesquisando um serviço de cloud para utilizar o R + Shiny. > >>> > >>> Minha duvida esta no custo e por isso gostaria de saber se alguém pode > >>> me ajudar para entender melhor: > >>> > >>> Quem é a empresa (cloud) que esta usando? AWS, Azure, Google ou IBM > >>> Qual o tipo de serviço para configurar o R + Shiny? Nestes servidores > >>> oferecem vários tipos, o que mais encontro na web é o EC2 - AWS, como > >>> vários tutoriais de instalação e configuração > >>> Se puder nos informar o uso e o custo deste serviço? Principalmente no > >>> uso de Memoria e Storage de dados, algum parâmetro para temos uma > noção. > >>> > >>> Somos uma startup e queremos oferecer nos serviços utilizando o R + > >>> Shiny. Porem gostaríamos de um parâmetro para tomar a decisão de > configurar > >>> este serviço e começar a utilizar. > >>> > >>> Espero que fui claro e toda a ajuda será bem vinda. > >>> > >>> Obrigado desde já. > >>> > >>> Marcio. > >>> -- > >>> Marcio B. > >>> > >>> Cel.: 19 99648.3949 tim 99798.0104 vivo > >>> skype: maborbaa > >>> http://www.linkedin.com/in/marciobar > >>> https://www.facebook.com/marciobar > >>> > >>> Empresa > >>> http://www.techagr.com/ > >>> > >>> ==== PAZ E BEM ==== > >>> > >>> _______________________________________________ > >>> R-br mailing list > >>> R-br(a)listas.c3sl.ufpr.br > >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > >>> código mínimo reproduzível. > >>> > >> > >> > > > > > > -- > > Marcio B. > > > > Cel.: 19 99648.3949 tim 99798.0104 vivo > > skype: maborbaa > > http://www.linkedin.com/in/marciobar > > https://www.facebook.com/marciobar > > > > Empresa > > http://www.techagr.com/ > > > > ==== PAZ E BEM ==== > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > > código mínimo reproduzível. > > > > > > -- > Att, > > | Fernando Gama da Mata | > | Database Specialist | Master's Degree UFPA | > > | Contacts: +55 91 99150 0365 | f.fabiogama88(a)gmail.com | Social > Networks: [ > <https://www.facebook.com/fernando.gama.13>][ > <https://plus.google.com/+FernandoGama13>][ > <https://www.linkedin.com/in/fernandogama>] | > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20161018/7b9f2b7d/attachment-0001.html> > > ------------------------------ > > Message: 4 > Date: Tue, 18 Oct 2016 16:09:12 -0300 > From: Luiz Roberto Martins Pinto <luizroberto.uesc(a)gmail.com> > To: comunello.eder(a)gmail.com, a lista Brasileira oficial de > discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] geração de variáveis aleatórias normais com > estrutura de correlação > Message-ID: > <CAK4=-FBYwLbRMFfh3KYBOAJ+ktXi0MQm+xWDRKVL6SkXjE-4Ow@ > mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Se você desejar entender o significado de correlação experimente: > > n=10000000 > cor_x1x2=-0.2 > medx1=0.0067 > varx1=0.0017 > sdx1=varx1^0.5 > medx2=0.1374 > varx2=0.0024 > sdx2=varx2^0.5 > cov_x1x2=-cor_x1x2*(sdx1*sdx2) > > sda=(varx1-cov_x1x2)^0.5 > sdb=(varx2-cov_x1x2)^0.5 > a=rnorm(n,medx1,sda);var(a) > b=rnorm(n,medx2,sdb);var(b) > c=rnorm(n,0,cov_x1x2^0.5);var(c) > > x1=a+c > x2=b-c > var(x1) > var(x2) > cor(x1,x2) > > > Abraços, > > 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 > > "*The s**cience exists because there are patterns. The patterns exist > because God created them*. > * The statistic exists to research the patterns that God created.*" > > > > Em 18 de outubro de 2016 14:52, Éder Comunello via R-br < > r-br(a)listas.c3sl.ufpr.br> escreveu: > > > Adriele, boa tarde! > > > > Tentei adaptar seu código R com base no que entendi do script em SAS. > Veja > > se pode ser útil... > > > > ### <code r> > > library(MASS) > > { > > set.seed(1357) > > n <- 1000 # tamanho da amostra gerada (N) > > p <- 2 # numero de variáveis a serem geradas (K) > > ME <- rep(1, p) > > rho <- -0.2 # correlacao negativa - 0.2 > > sigma2 <- 1 > > sigma <- sigma2 * ((1-rho)*diag(p)+rho*matrix(1, p, p)) > > SI <- mvrnorm(n, ME, sigma) > > } > > apply(SI, 2, mean) # ~ 1 > > # [1] 1.037774 1.017846 > > > > cor(SI) # ~ sigma > > # [,1] [,2] > > # [1,] 1.0000000 -0.1973338 > > # [2,] -0.1973338 1.0000000 > > > > mu <- c(0.0067, 0.1374) > > mean(0.0017*SI[,1]) # ~ 0.0017 > > mean(0.0024*SI[,2]) # ~ 0.0024 > > > > RES <- data.frame(X1=0.0067+0.0017*SI[,1], X2=0.1374+0.0024*SI[,2]) > > > > round(cbind(RES=colMeans(RES), INI=mu),4) > > # RES INI > > # X1 0.0085 0.0067 > > # X2 0.1398 0.1374 > > > > cor(RES) > > # X1 X2 > > # X1 1.0000000 -0.1973338 > > # X2 -0.1973338 1.0000000 > > ### </code> > > > > > > ================================================ > > Éder Comunello > > Researcher at Brazilian Agricultural Research Corporation (Embrapa) > > DSc in Agricultural Systems Engineering (USP/Esalq) > > MSc in Environ. Sciences (UEM), Agronomist (UEM) > > --- > > Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| > > ================================================ > > GEO, -22.2752, -54.8182, 408m > > UTC-04:00 / DST: UTC-03:00 > > > > > > > > > > Em 17 de outubro de 2016 16:05, Adriele Giaretta Biase via R-br < > > r-br(a)listas.c3sl.ufpr.br> escreveu: > > > >> Olá pessoal, > >> > >> > >> tenho uma dúvida com relação à geração de variáveis aleatórias normais > >> multivariadas no R. > >> > >> Eu gostaria de gerar duas variáveis aleatórias (x1, x2) usando > >> distribuição normal multivariada com estrutura de correlação entre > elas. A > >> Estrutura da matriz de correlação foi estimada antes com base num banco > de > >> dados reais (com correlação fraca de - 0.2). Porém, essas variáveis não > >> podem ser negativas, Ex. x1 tem média e variância, respectivamente, > iguais > >> a 0.0067 e 0.0017; x2 tem média e variância, respectivamente, iguais a > >> 0.1374 e 0.0024. Eu gero as variáveis da seguinte forma: > >> > >> > >> > >> *Programa usado no R:* > >> > >> # Simulando os parâmetros com estrutura de correlação entre eles > >> > >> n <- 1000 # tamanho da amostra gerada > >> > >> p <- 2 # numero de variáveis a serem geradas > >> > >> > >> > >> library(MASS) > >> > >> # construíndo a matriz de correlação para usar nas simulações, baseadas > >> na característica da amostra coletada > >> > >> mu<-rep(0, times = p) > >> > >> rho <- - 0.2 # correlacao negativa - 0.2 > >> > >> sigma2 <- 1 > >> > >> Sigma <- sigma2 * ((1-rho)*diag(p)+rho*matrix(1, p, p)) > >> > >> X1 <- 0.0096 # Media da variavel x1 > >> > >> X2 <- 0.1203 # Media da variavel x2 > >> > >> media <- c(X1,X2) > >> > >> y <- mvrnorm(n, media, Sigma) > >> > >> cor(y) > >> > >> apply(y, 2, mean) > >> > >> y > >> > >> > >> [1,] 0.309910452 1.0642521083 > >> > >> [2,] -0.251583312 1.8909032058 > >> > >> [3,] 1.330362012 -1.1239501814 > >> > >> [4,] -0.793399464 -1.5433056284 > >> > >> [5,] 2.165144843 -0.2645184534 > >> > >> [6,] 0.532777085 1.0910864562 > >> > >> [7,] -1.612135390 2.0489354648 > >> > >> [8,] -0.430529913 1.0312062602 > >> > >> ... > >> > >> > >> > >> Quando gero as simulações para x1 e x2 usando a função *mvrnorm*, o > >> resultado me retorna alguns valores negativos para as variáveis, isso > não > >> poderia ocorrer. Teria alguma outra função em que eu possa fornecer a > >> variância de cada uma dessas variáveis, além da estrutura de correlação? > >> Como poderia contornar essa situação, usando o R? > >> > >> > >> > >> *Obs:* No SAS, a seguinte programação funciona perfeitamente: > >> > >> > >> *proc* *iml*; > >> > >> wrksize=*100000*; > >> > >> K=*2*; > >> > >> N=*1000*; /* tamanho da amostra */ > >> > >> M={*0* *0*}; > >> > >> S={ *1* -*0.5283*, > >> > >> -*0.5283* *1* }; > >> > >> X=shape(*0*,K,N); > >> > >> ME=*0*; SI=*1*; > >> > >> DO I=*1* TO K; > >> > >> DO J=*1* TO N; > >> > >> if I>*1* > >> > >> then > >> > >> do; > >> > >> ME=M[I]+(S[*1*:I-*1*,I])`*(inv(S[*1*:I-*1*,*1*:I-*1*])*( > X[*1*:I- > >> *1*,J]-M[*1*:I-*1*])); > >> > >> SI=S[I,I]-(S[*1*:I-*1*,I])`*(inv(S[*1*:I-*1*,*1*:I-*1*])*( > S[*1* > >> :I-*1*,I])); > >> > >> end; > >> > >> X[I,J]=ME+NORMAL(*0*)*SQRT(SI); > >> > >> > >> > >> END; > >> > >> > >> > >> END; > >> > >> Z=t(X); > >> > >> > >> > >> varnames='X1':'X2'; > >> > >> create NOVO from Z[colname=varnames]; > >> > >> append from Z; > >> > >> *quit*; > >> > >> > >> > >> *data* MEXT; set NOVO; > >> > >> options ps=*66* ls=*75*; > >> > >> X1=*0.0067*+*0.0017**X1; /* normal */ > >> > >> X2=*0.1374*+*0.0024**X2; /* normal */ > >> > >> > >> > >> *proc* *corr* data=MEXT; > >> > >> var X1 X2; > >> > >> *run*; > >> > >> -- > >> Adriele Giaretta Biase. > >> Mestre em Estatística e Experimentação Agropecuária - UFLA. > >> Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP > >> Contato: (19) 98861-0619. > >> > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br > >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > >> código mínimo reproduzível. > >> > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > > código mínimo reproduzível. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20161018/d54ee426/attachment.html> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 70, assunto 18 > ***************************************** > -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

1 0

Ajuste de modelo não linear
by Andre Oliveira 17 Out '16

17 Out '16

Boa noite pessoal,o que eu poderia modificar neste ajuste para obter exito? Já variei os chutes iniciais e não consegui. Brasil<- c(13.63 ,13.93,14.09,15.12,15.98,16.08,16.45,17.96,18.83,19.28,19.18,20.89,21.87,21.85,23.43,25.38,25.64,26.99,28.68,31.92,32.92,33.32,34.37,34.99,36.08,37.72,40.41,42.81,44.18,44.91,48.97,51.59,53.44,56.49,62.2,66.49,67.01,68.41,74.14,77.58,80.39,85.71,89.58,96.3,100.9,102.8,107.69,111.51,111.42,116.59,118.85,121.42,127.53); anos <- c(1961,1962,1963,1964,1965,1966,1967,1968,1969,1970,1971,1972,1973,1974,1975,1976,1977,1978,1979,1980,1981,1982,1983,1984,1985,1986,1987,1988,1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001,2002,2003,2004,2005,2006,2007,2008,2009,2010,2011,2012,2013); fit<-nls(Brasil~a*(b^(anos)), start = list(a =18, b=1.4)) Error in nls(Brasil ~ a * (b^(anos)), start = list(a = 18, b = 1.4)) : número de iterações excedeu o máximo de 50 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

5 6

(sem assunto)
by Wilmar Freitas 17 Out '16

17 Out '16

1 0

(sem assunto)
by Wilmar Freitas 17 Out '16

17 Out '16

1 0

Regressão linear - Excluindo outliers
by Mac David S. Pinto 14 Out '16

14 Out '16

Pessoal, bom dia Estou fazendo algumas regressões lineares entre Comprimento x Peso de peixes com um conjunto de dados bastante extenso (+- 1000 linhas). Acontece que alguns destes valores são outliers (+- 50) e na analise, eles não são importantes. Tem alguma função no R que eu consiga "selecionar" os outliers e pedir que não seja considerado na analise? Sem precisa mudar minha planilha original (já que vou utilizar para outras coisas futuramente) Tem alguns valores que estão no meio dos planilhas de dados. Eu consegui resolver os extremos selecionado um conjunto de linhas Exemplo: Local.F.1<-lm(pt[4:970,7]~ls[4:970,6]) ## Com este comando eu retire os outliers dos extremos (linhas 1 a 3 e linhas 971 a 1000), mais ainda tem alguns valores no meio do analise (+- entre as linhas 251 a 260) De já, eu agradeço a ajuda David

6 6

Re: [R-br] Contents of Digest R-br, volume 70, assunto 13 Regressão linear - Excluindo outliers
by sirio＠ig.com.br 14 Out '16

14 Out '16

David, a função boxplot oferece a possibilidade de identificar outliers e logo após excluir Como exemplo: boxplot.stats(meus_dados)$out outliers <- which(meus_dados%in%boxplot.stats(meus_dados)$out) Você pode usar o outliers para filtrar os dados e trabalhar com eles limpo. Espero ter ajudado Humberto =========================================== Walter Humberto Subiza Piña - Eng. Cart., Dr. Gerencia de Imageamento- Coordenação de Cartografia Diretoria de Geociências IBGE - Instituto Brasileiro de Geografia e Estatística "In other words, maps hold a clue to what makes us human" Simon Garfield - On the map =========================================== Em 14/10/2016 12:00, r-br-request(a)listas.c3sl.ufpr.br escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br [1] > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > Tópicos de Hoje: > > 1. Pacotes do R (Edmar Caldas) > 2. Re: Pacotes do R (Tito Conte) > 3. Re: Pacotes do R (Daniela Petruzalek) > 4. Regressão linear - Excluindo outliers (Mac David S. Pinto) > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 14 Oct 2016 12:57:57 +0000 (UTC) > From: Edmar Caldas <edmar_caldas(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Pacotes do R > Message-ID: <1259226042.208586.1476449877464(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Bom Dia! > Preciso baixar vários pacotes do R. São 120 pacotes onde posso encontrar porque ficar pesquisando um por um não dá. Edmar > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/d53aaf45/atta… [3]> > > ------------------------------ > > Message: 2 > Date: Fri, 14 Oct 2016 10:33:36 -0300 > From: Tito Conte <tito.conte(a)gmail.com> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Pacotes do R > Message-ID: > <CACqq46zTwEcC=2zzWGjOf7vCEV1+_EPc_90WF+hzbEGkez59rA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Edmar, o que você precisa fazer exatamente? É mais fácil buscar o que você > precisa fazer na internet do que ficar consultando pacote por pacote. > > Tito Conte > > Em 14 de outubro de 2016 09:57, Edmar Caldas via R-br < > r-br(a)listas.c3sl.ufpr.br> escreveu: > >> Bom Dia! Preciso baixar vários pacotes do R. São 120 pacotes onde posso encontrar porque ficar pesquisando um por um não dá. Edmar _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br [1] Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia [2]) e forneça código mínimo reproduzível. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/47d8fb25/atta… [4]> > > ------------------------------ > > Message: 3 > Date: Fri, 14 Oct 2016 13:51:43 +0000 > From: Daniela Petruzalek <daniela.petruzalek(a)gmail.com> > To: Edmar Caldas <edmar_caldas(a)yahoo.com.br>, a lista Brasileira > oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Pacotes do R > Message-ID: > <CA+=iygtmtjXTpkvPEvt3i8nYnZtccVwuj0wLdvuyEAi1DHLVuw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > A melhor resposta para esta pergunta é "depende" :) > > Se estes pacotes estiverem na CRAN você não precisa baixá-los de um site, > basta rodar a função install.packages passando o nome do pacote para a > função. > > Se não estiverem na CRAN você vai ter que descobrir onde os autores > disponibilizam eles e seguir as instruções para cada pacote. Aí a internet > é o caminho. > > Em sex, 14 de out de 2016 às 09:58, Edmar Caldas via R-br < > r-br(a)listas.c3sl.ufpr.br> escreveu: > >> Bom Dia! Preciso baixar vários pacotes do R. São 120 pacotes onde posso encontrar porque ficar pesquisando um por um não dá. Edmar _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br [1] Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia [2]) e forneça código mínimo reproduzível. > > -- > -- > Daniela Petruzalek > Consultora de TI, OCP DBA, OCP PL/SQL > daniela.petruzalek(a)gmail.com > +55 42 9116-1104 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/af340a63/atta… [5]> > > ------------------------------ > > Message: 4 > Date: Fri, 14 Oct 2016 14:26:08 +0000 (UTC) > From: "Mac David S. Pinto" <macdavidsp(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Regressão linear - Excluindo outliers > Message-ID: <186111342.346980.1476455168933(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Pessoal, bom dia > Estou fazendo algumas regressões lineares entre Comprimento x Peso de peixes com um conjunto de dados bastante extenso (+- 1000 linhas). > Acontece que alguns destes valores são outliers (+- 50) e na analise, eles não são importantes. > Tem alguma função no R que eu consiga "selecionar" os outliers e pedir que não seja considerado na analise? Sem precisa mudar minha planilha original (já que vou utilizar para outras coisas futuramente) > Tem alguns valores que estão no meio dos planilhas de dados. > Eu consegui resolver os extremos selecionado um conjunto de linhas > Exemplo: > Local.F.1<-lm(pt[4:970,7]~ls[4:970,6]) > ## Com este comando eu retire os outliers dos extremos (linhas 1 a 3 e linhas 971 a 1000), mais ainda tem alguns valores no meio do analise (+- entre as linhas 251 a 260) > De já, eu agradeço a ajuda > David > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/723ff501/atta… [6]> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br [1] > > ------------------------------ > > Fim da Digest R-br, volume 70, assunto 13 > ***************************************** Links: ------ [1] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br [2] http://www.leg.ufpr.br/r-br-guia [3] http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/d53aaf45/atta… [4] http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/47d8fb25/atta… [5] http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/af340a63/atta… [6] http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161014/723ff501/atta…

1 0

Pacotes do R
by Edmar Caldas 14 Out '16

14 Out '16

Bom Dia! Preciso baixar vários pacotes do R. São 120 pacotes onde posso encontrar porque ficar pesquisando um por um não dá. Edmar

3 2

Aprendendo a manipular dados netcdf no R
by Paulo Henrique Pimenta 13 Out '16

13 Out '16

Olá Listeiros, Estou começando a manipular dados netcdf no R, antes usava o software grads para manipular dados netcdf de tempo e clima, por isso gostaria de opinião e ajuda de vcs sobre tutoriais ou materiais para quem está iniciando. Gostaria de aprender a manipular e criar os mapas. Caso alguém tenha alguma dica ou até mesmo códigos por favor, peço que me informe. Desde já agradeço a todos! Paulo Henrique Pimenta Graduando em Meteorologia - IAG USP Celular/WhatsApp/Telegram: +5511981418435 Este email foi enviado de um telemóvel Samsung Galaxy S5 mini Obter o Outlook para Android

1 0

Problemas para instalar RNetCDF no Rstudio (linux)
by Paulo Henrique Pimenta 13 Out '16

13 Out '16

Ol? listeiros! Usando as dicas do link do stack over flow tive ?xito na instala??o dos pacotes RNetCDF e ncdf4. Muito obrigado pela ajuda de todos vcs! Abra?os a todos. Paulo Henrique Pimenta Graduando em Meteorologia - IAG USP Celular/WhatsApp/Telegram: +5511981418435 Este email foi enviado de um telem?vel Samsung Galaxy S5 mini Obter o Outlook para Android

1 0

Problemas para instalar RNetCDF no Rstudio (linux)
by Paulo Henrique Pimenta 12 Out '16

12 Out '16

Olá a todos, Estou tentando instalar o RNetCDF no R do qual uso, num ubunutu 12.04 no RStudio, porém não consigo instalar o pacote aparecendo o seguinte erro: ####### > install.packages("RNetCDF") Installing package into '/home/notemaster/R/i686-pc-linux-gnu-library/3.0' (as 'lib' is unspecified) --2016-10-11 17:49:52-- https://cran.rstudio.com/src/contrib/RNetCDF_1.8-2.tar.gz Resolvendo cran.rstudio.com (cran.rstudio.com) 54.192.56.107 Conectando-se a cran.rstudio.com (cran.rstudio.com)|54.192.56.107|:443... conectado. A requisição HTTP foi enviada, aguardando resposta... 200 OK Tamanho: 95693 (93K) [application/x-gzip] Salvando em: "/tmp/Rtmpbh1Jqr/downloaded_packages/RNetCDF_1.8-2.tar.gz" 0K .......... .......... .......... .......... .......... 53% 68,7K 1s 50K .......... .......... .......... .......... ... 100% 50,2K=1,6s 2016-10-11 17:49:54 (58,7 KB/s) - "/tmp/Rtmpbh1Jqr/downloaded_packages/RNetCDF_1.8-2.tar.gz" salvo [95693/95693] * installing *source* package 'RNetCDF' ... ** package 'RNetCDF' successfully unpacked and MD5 sums checked checking for nc-config... yes checking for gcc... gcc -std=gnu99 checking whether the C compiler works... yes checking for C compiler default output file name... a.out checking for suffix of executables... checking whether we are cross compiling... no checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc -std=gnu99 accepts -g... yes checking for gcc -std=gnu99 option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -std=gnu99 -E checking for grep that handles long lines and -e... /bin/grep checking for egrep... /bin/grep -E checking for ANSI C header files... yes checking for sys/types.h... yes checking for sys/stat.h... yes checking for stdlib.h... yes checking for string.h... yes checking for memory.h... yes checking for strings.h... yes checking for inttypes.h... yes checking for stdint.h... yes checking for unistd.h... yes checking netcdf linker flags... -L/usr/lib -lnetcdf checking netcdf pre-processor and compiler flags... -I/usr/include -DgFortran checking for library containing XML_StopParser... no checking for library containing utFree... no checking for library containing utScan... no configure: error: "unable to use udunits2 or udunits" ERROR: configuration failed for package 'RNetCDF' * removing '/home/notemaster/R/i686-pc-linux-gnu-library/3.0/RNetCDF' Warning in install.packages : installation of package 'RNetCDF' had non-zero exit status The downloaded source packages are in '/tmp/Rtmpbh1Jqr/downloaded_packages' > library(RNetCDF) Erro em library(RNetCDF) : there is no package called 'RNetCDF' #################### Tentei usar o pacote netcdf, mas não consegui devido ao pacote ter sido descontinuado. Espero que alguém possa me ajudar, pois sou novo no uso de netcdf no R. Abraços a todos. Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435.

3 2

Re: [R-br] Extração de dados climáticos
by sznelwar＠uol.com.br 10 Out '16

10 Out '16

2 1

pacote Hyfo
by Elisa Henning 10 Out '16

10 Out '16

Bom dia a todos Estóu com um problema ao usar o pacote Hyfo do R. Precisamos que os dados estejam na forma de lista (list); Gostaria de saber como converter em lista um data.frame contendo na primeira coluna a data e nas demais, as observações (precipitação) de cada estação. Obrigado Elisa

2 1

Re: [R-br] Regexp no grep para linhas que não iniciam com sequência
by Rafael Tieppo 05 Out '16

05 Out '16

Eder, tudo bem?Tentando ajudar, caso tenha alguma máquina com linux, tenho a seguinte sugestão. obs.: parece que o update do windows 10 tem ou terá uma layer do shell do linux test <- c("SPA100", "MSA200", "MSB300", "MSC400", "MSC500", "PRA100", "PRC200", "MGV100", "MTJ400", "MTK500") write.csv(test, "eder.csv", row.names = FALSE) ## indo para o Shell grep -v "MS" < eder.csv ## para salvar em outro arquivo grep -v "MS" < eder.csv > eder2.csv cat eder2.csv espero ter ajudadoSaudações Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://docente.unemat.br/rafaeltieppo/ blog: https://fuidebicicleta.wordpress.com"Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente". On Tuesday, October 4, 2016 11:00 AM, "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> wrote: Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. "alternativa não-paramétrica" para o teste de Dunnett (Luiz Leal) 2. Re: Processamento paralelo para um modelo de treinamento (Tito Conte) 3. Re: "alternativa não-paramétrica" para o teste de Dunnett (Tito Conte) 4. Identificar variáveis dicotomicas (Elias Carvalho) 5. Re: Identificar variáveis dicotomicas (Eduardo Junior) 6. Re: Identificar variáveis dicotomicas (Felipe) 7. Re: "alternativa não-paramétrica" para o teste de Dunnett (Luiz Henrique Leal) 8. Regexp no grep para linhas que não iniciam com sequência (Éder Comunello) 9. Re: Regexp no grep para linhas que não iniciam com sequência (Marcus Nunes) ---------------------------------------------------------------------- Message: 1 Date: Mon, 3 Oct 2016 16:24:09 +0000 (UTC) From: Luiz Leal <richfield1974(a)yahoo.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] "alternativa não-paramétrica" para o teste de Dunnett Message-ID: <318158849.4552633.1475511849793(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Boa tarde a todos.Alguém sabe se existe uma "alternativa não-paramétrica" para o teste de Dunnett?Desde já agradeçoLuiz ...Dunnett's test is a multiple comparison procedure to compare each of a number of treatments with a single control. ... -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161003/2222ae8c/atta…> ------------------------------ Message: 2 Date: Mon, 3 Oct 2016 17:37:36 -0300 From: Tito Conte <tito.conte(a)gmail.com> To: Fernando Gama <f.fabiogama88(a)gmail.com>, a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Processamento paralelo para um modelo de treinamento Message-ID: <CACqq46womHtn9eivqKSQG2K_SdSD9zuSF=9JQG82wWoJuJQAYQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Ele está rodando paralelo? Você verificou? Se unix use o comando top e aperte 1 Se windows abra o gerenciador de tarefas e verifique os núcleos de trabalho. Aproveite para ver se é a CPU do seu PC o problema. Pode ser a memória ou a escrita do disco. Uma solução de baixissimo nível, mas que funciona. é você quebrar o código. Em partes e rodar vários Rs, cada um com uma etapa, gerar uns arquivos intermediários e depois concatena tudo. Tito Conte Em 2 de outubro de 2016 08:41, Fernando Gama via R-br < r-br(a)listas.c3sl.ufpr.br> escreveu: > Pessoal, > > Estou tentando treinar um modelo com 192 atributos e meu objetivo é, no > pós-treinamento, identificar os atributos mais importantes. (seleção de > features). > > A questão é que estou tendo problemas para treinar o modelo porque o > processamento é extremamente lento. Pesquisei sobre o parallel e o > doParallel e coloquei no meu código mas aparentemente não obtive resultados > segue um trecho do código: > > library(caret) > library(doParallel) > > myControl <- trainControl(method = "repeatedcv", number = 10, repeats = 3, > allowParallel = TRUE) > > t<-proc.time() > > cl <- makeCluster(detectCores()) > > registerDoParallel(cl) > > model <- train(GENRE~., data=dtset_genres, method="lvq", preProcess = > "scale", trControl = myControl) > > stopCluster(cl) > > proc.time()-t > > > Alguma sugestão? > > -- > Att, > > | Fernando Gama da Mata | > | Database Specialist | Master's Degree UFPA | > > | Contacts: +55 91 99150 0365 | f.fabiogama88(a)gmail.com | Social > Networks: [ <https://www.facebook.com/fernando.gama.13>][ > <https://plus.google.com/+FernandoGama13>][ > <https://www.linkedin.com/in/fernandogama>] | > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161003/6d0dada5/atta…> ------------------------------ Message: 3 Date: Mon, 3 Oct 2016 17:38:38 -0300 From: Tito Conte <tito.conte(a)gmail.com> To: Luiz Leal <richfield1974(a)yahoo.com>, a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] "alternativa não-paramétrica" para o teste de Dunnett Message-ID: <CACqq46zura6UmNb=YD0RG84RF4hucLL-MT2aRKQQC8aXVbqnsQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Achei isso: "I'm pretty sure Dunnett's Test is for parametric data only. You can use Mann Whitney follow up tests and just divide alpha by how many follow ups you need to run to take a conservative approach." Tito Conte 2016-10-03 13:24 GMT-03:00 Luiz Leal via R-br <r-br(a)listas.c3sl.ufpr.br>: > Boa tarde a todos. > Alguém sabe se existe uma "alternativa não-paramétrica" para o teste de > Dunnett? > Desde já agradeço > Luiz > > ... > *Dunnett's test* is a multiple comparison > <https://en.wikipedia.org/wiki/Multiple_comparisons> procedure to compare > each of a number of treatments with a single control. > ... > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161003/1a2c2fb0/atta…> ------------------------------ Message: 4 Date: Mon, 3 Oct 2016 20:18:27 -0300 From: Elias Carvalho <ecacarva(a)gmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Identificar variáveis dicotomicas Message-ID: <CAH8E6UvJE4FJdg3sqPSyEJuUy1QR45KEW0B3kTS4i6QUCyypEw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Boa noite pessoal Existe alguma função do R que permite identificar variáveis dicotomicas em um dataframe ? Como por exemplo: 0 e 1 ou 1 e 2 no mesmo dataframe Obrigado -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161003/4c256e4c/atta…> ------------------------------ Message: 5 Date: Mon, 3 Oct 2016 20:48:23 -0300 From: Eduardo Junior <edujrrib(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Identificar variáveis dicotomicas Message-ID: <CAO7eVPz+Nbop7a6e6=yvmMcT4ab2jJOp7KWQoHYvcbT_omCeDA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" As funções nativas. ## Mostra quantos valores unicos tem em cada coluna sapply(dados, function(x) length(unique(x))) ## Os iguais a dois ... sapply(dados, function(x) length(unique(x))) == 2L Att, Eduardo E. R, Junior <http://jreduardo.github.io/> -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161003/bf2f4bad/atta…> ------------------------------ Message: 6 Date: Mon, 3 Oct 2016 20:57:38 -0300 From: Felipe <felipe.e.barletta(a)gmail.com> To: Elias Carvalho <ecacarva(a)gmail.com>, a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Identificar variáveis dicotomicas Message-ID: <efd9e8d7-87a5-26cf-bca9-3bf1576561b0(a)gmail.com> Content-Type: text/plain; charset=utf-8; format=flowed Se as variáveis dicotômicas tiverem valores 0 ou 1, pode tentar da seguinte maneira: dados <- data.frame(x1=rbinom(100,1,.5),x2=rnorm(100)) sapply(dados, function(x) all(x%in%c(0,1))) x1 x2 TRUE FALSE ------------------------------ Message: 7 Date: Mon, 03 Oct 2016 22:50:14 -0300 From: Luiz Henrique Leal <richfield1974(a)yahoo.com> To: Tito Conte <tito.conte(a)gmail.com> Cc: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] "alternativa não-paramétrica" para o teste de Dunnett Message-ID: <owgspse5irq36a1succbe6m8.1475545814367(a)email.android.com> Content-Type: text/plain; charset="utf-8" Muito obrigado Conectado pela Motorola Tito Conte <tito.conte(a)gmail.com> escreveu: >Achei isso: > >"I'm pretty sure Dunnett's Test is for parametric data only. You can use Mann Whitney follow up tests and just divide alpha by how many follow ups you need to run to take a conservative approach." > > >Tito Conte > > >2016-10-03 13:24 GMT-03:00 Luiz Leal via R-br <r-br(a)listas.c3sl.ufpr.br>: > >Boa tarde a todos. > >Alguém sabe se existe uma "alternativa não-paramétrica" para o teste de Dunnett? > >Desde já agradeço > >Luiz > > >... > >Dunnett's test is a multiple comparison procedure to compare each of a number of treatments with a single control. > >... > > >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161003/8ca1874d/atta…> ------------------------------ Message: 8 Date: Tue, 4 Oct 2016 08:49:15 -0400 From: Éder Comunello <comunello.eder(a)gmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Regexp no grep para linhas que não iniciam com sequência Message-ID: <CABmC8gmXyf0Cv1movTQqUktHucnSUSd-rxLWv0CLDZLcG6aEcw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Senhores, bom dia! Considerando o vetor: test <- c("SPA100", "MSA200", "MSB300", "MSC400", "MSC500", "PRA100", "PRC200", "MGV100", "MTJ400", "MTK500") Posso obter facilmente os valores iniciados com a sequência "MS": grep("^MS", test, val=T) # [1] "MSA200" "MSB300" "MSC400" "MSC500" Uma primeira tentativa, sem sucesso, para obter linhas que não iniciam com a sequência foi: grep("^!(MS)", test, val=T) # character(0) Consegui uma solução estranha com: grep("^[^M]|^M[^S]", test, val=T) # [1] "SPA100" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" Mas essa solução é limitada quando usando strings maiores, tal como "MSC". Embora existam alternativas como: test[-grep("^MSC", test)] # [1] "SPA100" "MSA200" "MSB300" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" Gostaria de saber se há uma forma fácil de fazer isso usando diretamente uma regexp. Grato, ================================================ Éder Comunello Researcher at Brazilian Agricultural Research Corporation (Embrapa) DSc in Agricultural Systems Engineering (USP/Esalq) MSc in Environ. Sciences (UEM), Agronomist (UEM) --- Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| ================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161004/ed7972df/atta…> ------------------------------ Message: 9 Date: Tue, 4 Oct 2016 10:14:36 -0300 From: Marcus Nunes <marcus.nunes(a)gmail.com> To: comunello.eder(a)gmail.com, a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Regexp no grep para linhas que não iniciam com sequência Message-ID: <CA+QGQvvMACE6i_RA46nvpY-T34QS9hxZadhbNsaezwhXTg4SbQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Se entendi corretamente o teu problema, não é necessário utilizar uma regex: test <- c("SPA100", "MSA200", "MSB300", "MSC400", "MSC500", "PRA100", "PRC200", "MGV100", "MTJ400", "MTK500") grep("MS", test, val=T, invert=T) 2016-10-04 9:49 GMT-03:00 Éder Comunello via R-br <r-br(a)listas.c3sl.ufpr.br> : > Senhores, bom dia! > > Considerando o vetor: > > test <- c("SPA100", "MSA200", "MSB300", "MSC400", "MSC500", > "PRA100", "PRC200", "MGV100", "MTJ400", "MTK500") > > Posso obter facilmente os valores iniciados com a sequência "MS": > grep("^MS", test, val=T) > > # [1] "MSA200" "MSB300" "MSC400" "MSC500" > > > Uma primeira tentativa, sem sucesso, para obter linhas que não iniciam com > a sequência foi: > grep("^!(MS)", test, val=T) > # character(0) > > Consegui uma solução estranha com: > grep("^[^M]|^M[^S]", test, val=T) > > # [1] "SPA100" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" > > > Mas essa solução é limitada quando usando strings maiores, tal como "MSC". > > Embora existam alternativas como: > test[-grep("^MSC", test)] > > # [1] "SPA100" "MSA200" "MSB300" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" > > > Gostaria de saber se há uma forma fácil de fazer isso usando diretamente > uma regexp. > > Grato, > > ================================================ > Éder Comunello > Researcher at Brazilian Agricultural Research Corporation (Embrapa) > DSc in Agricultural Systems Engineering (USP/Esalq) > MSc in Environ. Sciences (UEM), Agronomist (UEM) > --- > Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| > ================================================ > GEO, -22.2752, -54.8182, 408m > UTC-04:00 / DST: UTC-03:00 > > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -- Marcus Nunes http://marcusnunes.me/ -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161004/585080f4/atta…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 70, assunto 6 ****************************************

4 3

Re: [R-br] Regexp no grep para linhas que não iniciam com sequência
by Éder Comunello 04 Out '16

04 Out '16

Pessoal, obrigado pelas dicas! Tive sucesso com a sugestão do Marcus e acabei encontrando mais uma forma, transcrita abaixo. test <- c("SPA100", "MSA200", "MSB300", "MSC400", "MSC500", "PRA100", "PRC200", "MGV100", "MTJ400", "MTK500") grep("MSC", test, val=T, invert=T) # [1] "SPA100" "MSA200" "MSB300" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" grep("^(?!MSC)", test, val=T, perl=T) # [1] "SPA100" "MSA200" "MSB300" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" ================================================ Éder Comunello Researcher at Brazilian Agricultural Research Corporation (Embrapa) DSc in Agricultural Systems Engineering (USP/Esalq) MSc in Environ. Sciences (UEM), Agronomist (UEM) --- Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| ================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00 Em 4 de outubro de 2016 10:25, Tiago Fragoso <mingote(a)gmail.com> escreveu: > Oi > > Estou longe do R no momento, mas na sua segunda solução, você tentou usar > a regex "^[^MSC]{0,3}"? > > Em princípio, ela força o início da string e identifica todos os > caracteres fora do conjunto MSC, para valores entre 0 e 3 ocorrências. > > Note que nesse caso a regex vai identificar strings iniciadas em M,S ou C. > Me parece que você tem os padrões M, MS ou MSC, então talvez a regex > > "^[^M|^MS|^MSC]" > > resolva. > > On Oct 4, 2016 9:49 AM, "Éder Comunello via R-br" < > r-br(a)listas.c3sl.ufpr.br> wrote: > >> Senhores, bom dia! >> >> Considerando o vetor: >> >> test <- c("SPA100", "MSA200", "MSB300", "MSC400", "MSC500", >> "PRA100", "PRC200", "MGV100", "MTJ400", "MTK500") >> >> Posso obter facilmente os valores iniciados com a sequência "MS": >> grep("^MS", test, val=T) >> >> # [1] "MSA200" "MSB300" "MSC400" "MSC500" >> >> >> Uma primeira tentativa, sem sucesso, para obter linhas que não iniciam >> com a sequência foi: >> grep("^!(MS)", test, val=T) >> # character(0) >> >> Consegui uma solução estranha com: >> grep("^[^M]|^M[^S]", test, val=T) >> >> # [1] "SPA100" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" >> >> >> Mas essa solução é limitada quando usando strings maiores, tal como >> "MSC". >> >> Embora existam alternativas como: >> test[-grep("^MSC", test)] >> >> # [1] "SPA100" "MSA200" "MSB300" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" >> >> >> Gostaria de saber se há uma forma fácil de fazer isso usando diretamente >> uma regexp. >> >> Grato, >> >> ================================================ >> Éder Comunello >> Researcher at Brazilian Agricultural Research Corporation (Embrapa) >> DSc in Agricultural Systems Engineering (USP/Esalq) >> MSc in Environ. Sciences (UEM), Agronomist (UEM) >> --- >> Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| >> ================================================ >> GEO, -22.2752, -54.8182, 408m >> UTC-04:00 / DST: UTC-03:00 >> >> >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> >

1 0

Regexp no grep para linhas que não iniciam com sequência
by Éder Comunello 04 Out '16

04 Out '16

Senhores, bom dia! Considerando o vetor: test <- c("SPA100", "MSA200", "MSB300", "MSC400", "MSC500", "PRA100", "PRC200", "MGV100", "MTJ400", "MTK500") Posso obter facilmente os valores iniciados com a sequência "MS": grep("^MS", test, val=T) # [1] "MSA200" "MSB300" "MSC400" "MSC500" Uma primeira tentativa, sem sucesso, para obter linhas que não iniciam com a sequência foi: grep("^!(MS)", test, val=T) # character(0) Consegui uma solução estranha com: grep("^[^M]|^M[^S]", test, val=T) # [1] "SPA100" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" Mas essa solução é limitada quando usando strings maiores, tal como "MSC". Embora existam alternativas como: test[-grep("^MSC", test)] # [1] "SPA100" "MSA200" "MSB300" "PRA100" "PRC200" "MGV100" "MTJ400" "MTK500" Gostaria de saber se há uma forma fácil de fazer isso usando diretamente uma regexp. Grato, ================================================ Éder Comunello Researcher at Brazilian Agricultural Research Corporation (Embrapa) DSc in Agricultural Systems Engineering (USP/Esalq) MSc in Environ. Sciences (UEM), Agronomist (UEM) --- Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| ================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00

2 1

Re: [R-br] "alternativa não-paramétrica" para o teste de Dunnett
by Luiz Henrique Leal 03 Out '16

03 Out '16

Muito obrigado Conectado pela Motorola Tito Conte <tito.conte(a)gmail.com> escreveu: >Achei isso: > >"I'm pretty sure Dunnett's Test is for parametric data only. You can use Mann Whitney follow up tests and just divide alpha by how many follow ups you need to run to take a conservative approach." > > >Tito Conte > > >2016-10-03 13:24 GMT-03:00 Luiz Leal via R-br <r-br(a)listas.c3sl.ufpr.br>: > >Boa tarde a todos. > >Alguém sabe se existe uma "alternativa não-paramétrica" para o teste de Dunnett? > >Desde já agradeço > >Luiz > > >... > >Dunnett's test is a multiple comparison procedure to compare each of a number of treatments with a single control. > >... > > >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. > >

1 0

Identificar variáveis dicotomicas
by Elias Carvalho 03 Out '16

03 Out '16

Boa noite pessoal Existe alguma função do R que permite identificar variáveis dicotomicas em um dataframe ? Como por exemplo: 0 e 1 ou 1 e 2 no mesmo dataframe Obrigado -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho*

3 2

"alternativa não-paramétrica" para o teste de Dunnett
by Luiz Leal 03 Out '16

03 Out '16

Boa tarde a todos.Alguém sabe se existe uma "alternativa não-paramétrica" para o teste de Dunnett?Desde já agradeçoLuiz ...Dunnett's test is a multiple comparison procedure to compare each of a number of treatments with a single control. ...

2 1

Processamento paralelo para um modelo de treinamento
by Fernando Gama 03 Out '16

03 Out '16

Pessoal, Estou tentando treinar um modelo com 192 atributos e meu objetivo é, no pós-treinamento, identificar os atributos mais importantes. (seleção de features). A questão é que estou tendo problemas para treinar o modelo porque o processamento é extremamente lento. Pesquisei sobre o parallel e o doParallel e coloquei no meu código mas aparentemente não obtive resultados segue um trecho do código: library(caret) library(doParallel) myControl <- trainControl(method = "repeatedcv", number = 10, repeats = 3, allowParallel = TRUE) t<-proc.time() cl <- makeCluster(detectCores()) registerDoParallel(cl) model <- train(GENRE~., data=dtset_genres, method="lvq", preProcess = "scale", trControl = myControl) stopCluster(cl) proc.time()-t Alguma sugestão? -- Att, | Fernando Gama da Mata | | Database Specialist | Master's Degree UFPA | | Contacts: +55 91 99150 0365 | f.fabiogama88(a)gmail.com | Social Networks: [ <https://www.facebook.com/fernando.gama.13>][ <https://plus.google.com/+FernandoGama13>][ <https://www.linkedin.com/in/fernandogama>] |

2 1

Bando de dado .txt
by Andre Oliveira 02 Out '16

02 Out '16

Boa noite,não estou conseguindo identificar o motivo de não abrir este arquivo .txt corretamente, ele tem em torno de 6.3 GBAlguém teria uma dica para compartilhar? Acredito que não estou trabalhando corretamente! library(data.table) readLines("DADOS_ENEM_2011.TXT", n = 1) [1] "3000000000012011 2803548708SAO BERNARDO DO CAMPO SP110000000000000000000011 . . ...3548708SAO BERNARDO DO CAMPO SP1111 577.10 576.00 575.50 737.10CEBDCACBEAABBAABAEBDCABDCDEAAADCBBEBADDBEBECDDEECABABEBDDEBAEADCACABCBBECADDAADBADDAABABCAEEAEAEDDDAAEEBDECADCEABDDBAECAAABDAEBCCACBAECBAEBECCAAEBABECECCDDAACDADADBACCBAADCAECCCEDD1241201251290CECDDACBCAEBEBBACEEDCAEBDAEDADBEBBDDADCBCBACDDEECEBABAADDEBAECDCACEDCBBEBCDCAADBADDDBAEACAEEDBDBDDCAEDDDACDEBBECAACDABADBAECCEADDCEBCBACBAEEBAEEECCEEEBCBECEBADCBDBEADCDBDCCCBADCAECCCCDDP 140.00 140.00 140.00 140.00 60.00 620.000 . . . require(bit64) enem2011<- fread("DADOS_ENEM_2011.TXT", nrow = 1, header = TRUE, sep=";") head(enem2011) 3000000000012011 2803548708SAO BERNARDO DO CAMPO SP110000000000000000000011 . . ...3548708SAO BERNARDO DO CAMPO SP1111 577.10 576.00 575.50 737.10CEBDCACBEAABBAABAEBDCABDCDEAAADCBBEBADDBEBECDDEECABABEBDDEBAEADCACABCBBECADDAADBADDAABABCAEEAEAEDDDAAEEBDECADCEABDDBAECAAABDAEBCCACBAECBAEBECCAAEBABECECCDDAACDADADBACCBAADCAECCCEDD1241201251290CECDDACBCAEBEBBACEEDCAEBDAEDADBEBBDDADCBCBACDDEECEBABAADDEBAECDCACEDCBBEBCDCAADBADDDBAEACAEEDBDBDDCAEDDDACDEBBECAACDABADBAECCEADDCEBCBACBAEEBAEEECCEEEBCBECEBADCBDBEADCDBDCCCBADCAECCCCDDP 140.00 140.00 140.00 140.00 60.00 620.000 . . . ... 1: 3000000000022011 2214125506SAO JOSE DOS PINHAIS PR110000000000000000000011 . . ...4125506SAO JOSE DOS PINHAIS PR0000 . . . . 1211171251290DCEACDBECBCAEBACBEAECEDEDAABDBEDBBDDDCABCBDCAEDEBCEAEDBDAACACECDBBCACEDCDCBAEADADBBDDEAABAEEDBDBDDCAEDDDACDEBBECAACDABADBAECCEADDCEBCBACBAEEBAEEECCEEEBCBECEBADCBDBEADCDBDCCCBADCAECCCCDDF 0.00 0.00 0.00 0.00 0.00 0.000 . . . obrigado. 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

3 8

Re: [R-br] Digest R-br, volume 70, assunto 3rrrregghtgrgfffdf
by Josmar Mazucheli 02 Out '16

02 Out '16

8ttgdftfese Em 02/10/2016 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Bando de dado .txt (Andre Oliveira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 2 Oct 2016 12:03:27 +0000 (UTC) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: Marcos Silva <marcosfs2006(a)gmail.com> > Cc: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Bando de dado .txt > Message-ID: <1862678262.297135.1475409807537(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Oi Marcos,obrigado pelas sugestões. Eu estava indo muito bem com a > data.table, tentei library(ff) e library(ffbase), mas sem sucesso. > > 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 > > Em Domingo, 2 de Outubro de 2016 9:00, Marcos Silva < > marcosfs2006(a)gmail.com> escreveu: > > > Dá uma olhada nesse post: > http://rcoster.blogspot.com.br/2014/02/lendo-grandes-bancos-de-dados.html > > Em 2 de outubro de 2016 08:36, Marcos Silva <marcosfs2006(a)gmail.com> > escreveu: > > André, > Este tipo de arquivo não tem um caractere delimitador, tal como virgula, > ponto e virgula, pipe ("|"), Você não vai conseguir importar com > read.table() e derivações (read.csv(), read.csv2()) ou mesmo com fread(). > Mas para vc "ler" esse arquivo, vc precisa saber qual o comprimento de > cada variável (quantidade de caracteres que ocupa no arquivo), informação > essa que deveria estar no dicionário de variáveis. > > Este tipo de arquivo vc importa com read.fwf() ou read_fwf() do pacote > readr. Mas como essa base é muito grande, recomendo utilizar algum pacote, > a exemplo do pacote LaF. Existem outras possibilidades talvez melhores. > > > > > > Em 2 de outubro de 2016 03:13, Andre Oliveira <andreolsouza(a)yahoo.com.br> > escreveu: > > Oi Marcos,neste dicionário fala que é formato em formato ASCII. Não tem > mais instruções pelo que percebi. Descreve os nomes da varáveis também! > 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 > > Em Domingo, 2 de Outubro de 2016 0:13, Marcos Silva < > marcosfs2006(a)gmail.com> escreveu: > > > Acho que este e um arquivo de formato fixo. Vc vai precisar do dicionario > de dados para importar. Normalmente este dicionario e disponibilizado com > os dados pelo inep. > Em 01/10/2016 22:31, "Andre Oliveira via R-br" <r-br(a)listas.c3sl.ufpr.br> > escreveu: > > Boa noite,não estou conseguindo identificar o motivo de não abrir este > arquivo .txt corretamente, ele tem em torno de 6.3 GBAlguém teria uma dica > para compartilhar? Acredito que não estou trabalhando corretamente! > > library(data.table) > readLines("DADOS_ENEM_2011. TXT", n = 1) > > [1] "3000000000012011 2803548708SAO BERNARDO DO > CAMPO > > SP110000000000000000000011 . . > > > ...3548708SAO BERNARDO DO CAMPO > > SP1111 577.10 576.00 > 575.50 737. 10CEBDCACBEAABBAABAEBDCABDCDEA AADCBBEBADDBEBECDDEECABABEBDDE > BAEADCACABCBBECADDAADBADDAABAB CAEEAEAEDDDAAEEBDECADCEABDDBAE > CAAABDAEBCCACBAECBAEBECCAAEBAB ECECCDDAACDADADBACCBAADCAECCCE > DD1241201251290CECDDACBCAEBEBB ACEEDCAEBDAEDADBEBBDDADCBCBACD > DEECEBABAADDEBAECDCACEDCBBEBCD CAADBADDDBAEACAEEDBDBDDCAEDDDA > CDEBBECAACDABADBAECCEADDCEBCBA CBAEEBAEEECCEEEBCBECEBADCBDBEA > DCDBDCCCBADCAECCCCDDP 140.00 140.00 140.00 140.00 60.00 > 620.000 . . . > > > require(bit64) > enem2011<- fread("DADOS_ENEM_2011.TXT", nrow = 1, header = TRUE, sep=";") > head(enem2011) > > 3000000000012011 2803548708SAO BERNARDO DO CAMPO > > SP110000000000000000000011 > . . > > ...3548708SAO BERNARDO DO > CAMPO > > SP1111 577.10 576.00 575.50 737. 10CEBDCACBEAABBAABAEBDCABDCDEA > AADCBBEBADDBEBECDDEECABABEBDDE BAEADCACABCBBECADDAADBADDAABAB > CAEEAEAEDDDAAEEBDECADCEABDDBAE CAAABDAEBCCACBAECBAEBECCAAEBAB > ECECCDDAACDADADBACCBAADCAECCCE DD1241201251290CECDDACBCAEBEBB > ACEEDCAEBDAEDADBEBBDDADCBCBACD DEECEBABAADDEBAECDCACEDCBBEBCD > CAADBADDDBAEACAEEDBDBDDCAEDDDA CDEBBECAACDABADBAECCEADDCEBCBA > CBAEEBAEEECCEEEBCBECEBADCBDBEA DCDBDCCCBADCAECCCCDDP 140.00 140.00 > 140.00 140.00 60.00 620.000 . . > . > > ... > 1: 3000000000022011 2214125506SAO JOSE DOS PINHAIS > > > PR110000000000000000000011 . . > > > ...4125506SAO JOSE DOS PINHAIS > > PR0000 . > . . . > > > 1211171251290DCEACDBECBCAEBACB > EAECEDEDAABDBEDBBDDDCABCBDCAED EBCEAEDBDAACACECDBBCACEDCDCBAE > ADADBBDDEAABAEEDBDBDDCAEDDDACD EBBECAACDABADBAECCEADDCEBCBACB > AEEBAEEECCEEEBCBECEBADCBDBEADC DBDCCCBADCAECCCCDDF 0.00 0.00 > 0.00 0.00 0.00 0.000 . . > . obrigado. > > > 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 > ______________________________ _________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/ cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br- guia) e forneça > código mínimo reproduzível. > > > > > > > > -- > Marcos F. Silva > http://sites.google.com/site/ marcosfs2006 > > > > -- > Marcos F. Silva > http://sites.google.com/site/marcosfs2006 > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20161002/24cbb345/attachment-0001.html> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 70, assunto 3 > **************************************** >

1 0

Título de gráficos em mais de uma linha
by Guilherme Biz 30 Set '16

30 Set '16

Boa tarde. Estou tentando fazer um gráfico com título extenso e gostaria que fosse dividido em duas ou mais linhas de forma automática, sem ser pelo comando \n. Ex: a=rnorm(100) b=rnorm(100) plot(a,b,main="Este título é muito grande para ficar só em uma linha e gostaria de dividi-lo em duas ou mais linhas de forma automática") Sei que usando "\n" eu passo para a linha de baixo e no Rstudio se eu escrever o título em duas linhas também ocorre isso, mas preciso de forma automática, pois são muitos gráficos e estou fazendo de forma iterativa. Alguém tem uma solução? Agradeço desde já. Att, Guilherme Biz

3 3

Re: [R-br] Parcela Subdividida
by Angélica Ricarte 27 Set '16

27 Set '16

Cesar Rabak e demais colegas, boa tarde. Agradeço o retorno e peço desculpas pela demora em responder esta postagem. Respondendo os questionamentos: - Por que então que você constrói (especifica) o seu termo Error como PARC:REP? R: o modelo da anova que eu utilizei [mod<-aov(var~PARC*SubP+Error(PARC:REP))] retirei de outras postagens feitas no fórum e também o vi em outras consultas, portanto, não sei o porquê do "Error como PARC:REP" (e gostaria de entender), mas entendo que seria o erro referente às parcelas ('Erro a' do resultado da anova). Em todas as minhas consultas o modelo para parcela subdividida é sempre este, mas comigo não funcionou (ou melhor, funcionou, mas apareceu uma mensagem de erro, como mostrado anteriormente). Entendo que se eu não inserir o termo do 'Error', a análise será gerada considerando um esquema fatorial. - Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer (uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar em conta na sua ANOVA? R: como dito anteriormente, o erro referente à Parcela (Erro a), visto que no resultado da anova devo ter dois Erros (a e b). - Outra perguntinha: você entende que ao especificar o termo erro você reduziu o número de pontos para ajuste para dez e esse é o número de fatores PARC vezes o número de parcelas SubP? R: não consigo compreender isto. Se você, ou alguém aqui do fórum, puder me ajudar, agradeço muito, estou ansiosa para resolver este problema. Se eu estiver errada em algum ponto, por favor, corrijam-me. Se o problema está no modelo, qual seria o que eu poderia utilizar para parcela subdividida (PS) (tenho um experimento em DIC, com parcelas subdividas (2 x 5) e 4 repetições)? Em suma, meu objetivo é gerar o resíduo (que não sei se utilizei o comando correto, visto que é uma PS) para poder realizar o teste de Shapiro-Wilk (SW), pois o quadro da anova que utilizo é aquele do pacote 'ExpDes.pt' (e, infelizmente, para o caso de PS este pacote não apresenta o resultado do teste de SW). Assim, se houver uma solução ainda mais simples para que possa atingir este objetivo, peço, por gentileza, que me apresentem. Muito obrigada. Aguardo retorno. Att., Angélica. ________________________________ De: Cesar Rabak <cesar.rabak(a)gmail.com> Enviado: sexta-feira, 9 de setembro de 2016 16:14 Para: Angélica Ricarte Assunto: Re: [R-br] Parcela Subdividida Olá Angélica, CMR (que está descrito no FAQ deste grupo) é "código mínimo reproduzível". A ideia é permitir que a gente refaça seus passos e acelere a ajuda. Seu código acima não é exatamente reproduzível porque esta linha: > PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T) diz respeito a um arquivo que está na sua máquina. Contudo, como você imprimiu o total do dataframe importado co o comando seguinte, está quase! OK, então você tem uma resposta "var" para o fator PARC e aninhado em PARC você tem quatro subparcelas, com quatro repetições por medida. Por que então que você constrói (especifica) o seu termo Error como PARC:REP? Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer (uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar em conta na sua ANOVA? Outra perguntinha: você entende que ao especificar o termo erro você reduziu o número de pontos para ajuste para dez e esse é o número de fatores PARC vezes o número de parcelas SubP? HTH -- Cesar Rabak 2016-09-08 14:27 GMT-03:00 Angélica Ricarte <angelica.ricarte(a)hotmail.com<mailto:angelica.ricarte@hotmail.com>>: Boa tarde, Cezar Rabak. Entendi melhor, porém ainda não por completo. Mas, o meu problema agora é que não consigo mais obter o mesmo resultado. Estou importando exatamente os mesmos dados e utilizando exatamente os mesmos comandos, porém não consigo mais o mesmo resultado que enviei no fórum dia 04/09/16. Quando realizo a ANOVA o resultado gerado está sendo diferente (e está errado) daquele que apresentei anteriormente (este está correto). Gostaria de saber como posso obter o "resíduo" para inserir no comando para o teste de Shapiro Willk: "shapiro.test(resíduo)". Seguem abaixo os comandos utilizados (acho que é o 'CMR' a que você se referiu) e o dados, alguém poderia me ajudar a encontrar o meu erro e mostrar uma solução? Já procurei muito, mas não encontro. _____________________________________________________________ > PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T);PS PARC SubP REP var 1 C 30 1 12 2 C 30 2 13 3 C 30 3 14 4 C 30 4 12 5 C 45 1 18 6 C 45 2 17 7 C 45 3 16 8 C 45 4 17 9 C 60 1 22 10 C 60 2 21 11 C 60 3 24 12 C 60 4 23 13 C 75 1 24 14 C 75 2 23 15 C 75 3 22 16 C 75 4 24 17 C 90 1 21 18 C 90 2 20 19 C 90 3 21 20 C 90 4 19 21 S 30 1 22 22 S 30 2 21 23 S 30 3 24 24 S 30 4 23 25 S 45 1 21 26 S 45 2 22 27 S 45 3 23 28 S 45 4 21 29 S 60 1 20 30 S 60 2 18 31 S 60 3 19 32 S 60 4 20 33 S 75 1 17 34 S 75 2 15 35 S 75 3 16 36 S 75 4 15 37 S 90 1 12 38 S 90 2 13 39 S 90 3 12 40 S 90 4 11 > attach(PS) > SubP<-as.factor(SubP);is.factor(SubP) [1] TRUE > mod<-aov(var~PARC*SubP+Error(PARC:REP));summary(mod) Error: PARC:REP Df Sum Sq Mean Sq F value Pr(>F) PARC 1 6.453 6.453 80.67 0.0706 . Residuals 1 0.080 0.080 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) PARC 1 1.7 1.67 1.535 0.226 SubP 4 108.9 27.23 25.076 6.64e-09 *** PARC:SubP 4 496.9 124.22 114.418 < 2e-16 *** Residuals 28 30.4 1.09 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > resíduo<-residuals(mod);resíduo NULL ___________________________________________________________________ Aguardo retorno. Obrigada. Att., Angélica. ________________________________ From: cesar.rabak(a)gmail.com<mailto:cesar.rabak@gmail.com> Date: Tue, 6 Sep 2016 21:49:25 -0300 Subject: Re: [R-br] Parcela Subdividida To: angelica.ricarte(a)hotmail.com<mailto:angelica.ricarte@hotmail.com>; r-br(a)listas.c3sl.ufpr.br<mailto:r-br@listas.c3sl.ufpr.br> Angélica, Como você deu nomes "crípticos" para seu dataframe, e não postou um CMR, dá só para fazer uma 'interpretação' das mensagens. Aparentemente você quer a resposta PS[,4] para a interação entre as variáveis (grupos[?]) PS[,1] e PS[,2] e considera a combinação dos fatores PS[,1] e PS[,3] como medidas repetidas e tem 28 casos ao todo. Além do mais você considera o fator PS[,1] como resposta e Error() ou seja parte da especificação das medidas repetidas. A mensagem sobre o modelo Error() ser singular ocorre quando você tenta ajustar mais parâmetros que seus dados permitem... Se você continuar com dificuldades, você precisa explicar melhor o que deseja fazer ou consultar uma obra sobre ANOVA com medidas repetidas. HTH -- Cesar Rabak 2016-09-05 21:05 GMT-03:00 Angélica Ricarte via R-br <r-br(a)listas.c3sl.ufpr.br<mailto:r-br@listas.c3sl.ufpr.br>>: Prezado Cesar, Não consegui compreender a mensagem de erro e nem o porquê que minha matriz é nula. Você pode ajudar? Agradeço e aguardo retorno. Obrigada. Att., Angélica ________________________________ Date: Mon, 5 Sep 2016 12:52:41 -0700 From: ml-node+s2285057n4666617h13(a)n4.nabble.com<mailto:ml-node%2Bs2285057n4666617h13@n4.nabble.com> To: angelica.ricarte(a)hotmail.com<mailto:angelica.ricarte@hotmail.com> Subject: Re: [R-br] Parcela Subdividida Angélica, Você entendeu a mensagem de erro? Não é possível analisar uma matriz nula.... 2016-09-04 1:22 GMT-03:00 Angélica Ricarte via R-br <[hidden email]>: Prezados, Estou tentando realizar uma análise de variância para um experimentos em parcelas subdivididas em DIC, mas, apesar de ser gerado um resultado, aparece uma mensagem de erro. Quando tento realizar o cálculo dos resíduos também dá erro e, consequentemente, não consigo fazer o teste de normalidade. Alguém pode me ajudar? Vejam como realizei as análises e os seus resultados: #ANOVA: > mod<-aov(PS[,4]~PS[,1]*PS[,2]+Error(PS[,1]:PS[,3]));summary(mod) Warning message: In aov(PS[, 4] ~ PS[, 1] * PS[, 2] + Error(PS[, 1]:PS[, 3])) : modelo Error() é singular Error: PS[, 1]:PS[, 3] Df Sum Sq Mean Sq F value Pr(>F) PS[, 1] 1 8.1 8.100 11.3 0.0152 * Residuals 6 4.3 0.717 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) PS[, 2] 4 108.9 27.23 24.94 3.02e-08 *** PS[, 1]:PS[, 2] 4 496.9 124.23 113.79 3.09e-15 *** Residuals 24 26.2 1.09 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Resíduos: > residuals(mod) NULL #Teste de Shapiro-Wilk: > res<-residuals(mod) > shapiro.test(res) Error: is.numeric(x) is not TRUE Obrigada. Att., Angélica. _______________________________________________ R-br mailing list [hidden email] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. _______________________________________________ R-br mailing list [hidden email] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne?a c?digo m?nimo reproduz?vel. ________________________________ If you reply to this email, your message will be added to the discussion below: http://r-br.2285057.n4.nabble.com/R-br-Parcela-Subdividida-tp4666613p466661… To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br<mailto:R-br@listas.c3sl.ufpr.br> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Salvar gráfico no Linux
by Alexandre Loures 26 Set '16

26 Set '16

Boa tarde! Para conseguir visualizar os gráficos preciso maximizar a janela. Porém, não estou sabendo salvar a janela maximizada no Linux. Segue pdf (salvo no windows) em anexo para exemplificar. Alguém poderia me ajudar? Desde já muito obrigado! Alexandre Rodrigues Loures Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE Site: www.ccsa.ufpb.br/ppge www.orcid.org/0000-0002-1288-0135

3 2

Precipitação diária baseada em uma rede neural
by Alisson Lucrécio 25 Set '16

25 Set '16

Caro Colegas, Boa tarde. Estou tentando calcular a ocorrência de dias secos e chuvosos e a quantidade de precipitação diária através de uma rede neural. Algum colega teria alguns exemplo que poderia me ajudar? Ou algum livro que poderia recorrer? Obrigado. -- Alisson Lucrecio da Costa

2 1

Extração de dados climáticos
by ana paula coelho madeira 23 Set '16

23 Set '16

________________________________ Prezados colegas, bom dia! Então, é o seguinte, no site http://berkeleyearth.org/ eu tenho dados de temperatura média de toda a terra. Eu preciso extrair de acordo com os meus pontos geográficos, a temperatura média anual para cada ponto. Só que eu tenho os específicos que eu preciso. Comecei o script de forma que eu pudesse extrair para um ponto específico e para todos os anos, para depois arrumar o script para fazer para todos os meus pontos e anos de interesse. Segue o script utilizado! library(raster) b <- brick("D:/Documentos/Desktop/analises_cap2/dados_climaticos/precipt.nc", varname <- "precip") print(b) b idx <- getZ(b) coords <- matrix(c(-55, -10), ncol = 2) vals <- extract(b, coords, df = T) df <- data.frame(idx, t(vals)[-1,]) rownames(df) <- NULL names(df) <- c('date','value') head(df) Seguindo o script no local onde deveria apresentar as variáveis extraídas , aparece NA. Alguém já fez este tipo de extração? Peço ajuda. Desde já, muito obrigada. Ana Paula

2 1

"Append" arquivos csv de diversos subdiretórios
by Alexandre Loures 22 Set '16

22 Set '16

Bom dia! Tenho um diretório com 50 subdiretórios sendo que cada um possui 11 arquivos csv. O que preciso é "empilhar" esses arquivos (em um único arquivo) e para isso estou fazendo o seguinte: filenames <- list.files () mydata <- do.call ("rbind", lapply (filenames, read.csv, header = TRUE)) Entretanto, gostaria de "automatizar" esse script, ou seja, o próprio comando fosse em cada subdiretório e fizesse o "append" (dos 550 arquivos csv) dando o resultado final. Alguém poderia me ajudar? Desde já muito obrigado!

4 4

Salvar uma tabela/matriz como imagem
by Elias Carvalho 20 Set '16

20 Set '16

Boa tarde a todos. Eu tenho uma matriz, chamada "matrix.results", como resultado de algumas análises: Consigo imprimir uma tabela no console com o comando: kable(head(matrix.results ), format = "rst", align = "r") Agora eu gostaria de salvar uma imagem desta tabela para inserir em um documento. No entanto eu gostaria de salvar como imagem uma tabela mais bem formatada, como esta: Alguém poderia ajudar ? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho*

3 2

Re: [R-br] Distribuição exponencial mista
by Alisson Lucrécio 20 Set '16

20 Set '16

Jaildo Rocha, Muito obrigado por sua ajuda. Alisson 2016-09-20 20:53 GMT-03:00 Alisson Lucrécio < alisson.lucrecio(a)ifgoiano.edu.br>: > Jaildo Rocha, > > Muito obrigado por sua ajuda. > > Alisson > > 2016-09-20 19:27 GMT-03:00 Jaildo Rocha <j.vieirarocha(a)gmail.com>: > >> Prezado Alisson, >> >> Para realizar as estimativas dos parâmetros via o package “Renext” >> utilize a função “EM.mixexp”. Apesar disso, é importante notar dois >> detalhes desta distribuição; >> >> >> - O primeiro é o caso do coeficiente de variação maior ou igual a um >> que tende a gerar “erros” (convergir para uma exponencial simples) quando a >> amostra apresenta tal coeficiente menor que um. >> >> >> >> - Um segundo detalhe é que as estimativas dos parâmetros,no que diz >> respeito a convergência em relação a amostra, é bem mais “lenta” com o >> crescimento da amostra do que distribuições usuals como a gama, weibull, >> gumbel e outras (o segundo detalhe é apenas um feeling particular). >> >> >> Caso não queira usar o “Renext”, outra possibilidade é utilizar o pacote >> “fitdistplus” para a realização das estimativas dos parâmetros, o único >> inconveniente é a necessidade de inserir um valor inicial.Por fim, segue os >> exemplos de estimativas utilizando os dois pacotes diferentes. >> >> >> ##### >> >> # Via Renext >> >> library("Renext") >> mix <- rmixexp2(n = 1000, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0) >> fit_mixexp <- EM.mixexp(x = mix) >> >> # Parâmetros >> fit_mixexp[[1]] >> >> >> # Via fitdistrplus >> >> library("Renext") >> library("fitdistrplus") >> >> fit_mixexp2 <- fitdist(mix, "mixexp2", start = list(prob1 = 0.8,rate1 = >> 1, rate2 = 2)) >> fit_mixexp2[[1]] >> >> >> #### >> >> Obs: No segundo exemplo, apesar de implícito, utilizei as funções >> pmixexp2, dmixexp2, qmixexp2 do pacote “Renext” mas caso seja de sua >> preferência é possível você mesmo defini-las. >> >> Att, >> ------- >> Jaildo Rocha >> MSc. pelo Programa Pós Graduação em Saneamento Meio Ambiente e Recursos >> Hídricos – UFMG >> Área de concentração: Recursos Hídricos >> >> >> Em 20 de set de 2016, à(s) 00:21, Alisson Lucrécio via R-br < >> r-br(a)listas.c3sl.ufpr.br> escreveu: >> >> Cesar, >> >> Preciso solucionar essa equação f(x) = (alpha/beta1*exp(-x/beta1)) + >> ((1-alpha)/beta2*exp(-x/beta2)) para obter os paramentos da >> distribuição. Assim como solucionado pelo fitdistr do pacote MASS para >> as ditribuções gamma, webull e exponencial. Observe que fiz uma tentativa >> com o pacote Renext, mas sem sucesso. >> >> Muito obrigado. >> >> 2016-09-19 16:12 GMT-03:00 Alisson Lucrécio <alisson.lucrecio@ifg >> oiano.edu.br>: >> >>> Como solucionar essa equação f(x) = (alpha/beta1*exp(-x/beta1)) + >>> ((1-alpha)/beta2*exp(-x/beta2)) para obter os paramentos da >>> distribuição. Assim como é feito pelo fitdistr do pacote MASS para as >>> ditribuções gamma, webull e exponencial. Observe que fiz uma tentativa com >>> o pacote Renext, mas sem sucesso. >>> >>> Muito obrigado. >>> >>> >>> On 19 September 2016 at 16:05, Cesar Rabak <cesar.rabak(a)gmail.com> >>> wrote: >>> >>>> . . . e sua dúvida seria? >>>> >>>> On Mon, Sep 19, 2016 at 10:54 AM, Alisson Lucrécio via R-br < >>>> r-br(a)listas.c3sl.ufpr.br> wrote: >>>> >>>>> Caro Colegas, >>>>> >>>>> Bom dia. >>>>> >>>>> Algum colega da poderia me ajudar na solução de uma distribução >>>>> exponencial mista. >>>>> >>>>> # Script >>>>> >>>>> library(MASS) >>>>> >>>>> # Distribuição Gamma >>>>> >>>>> rg <- rgamma(100, shape = 1.5, rate = 1) >>>>> >>>>> fit_gamma <- fitdistr(rg,"gamma") >>>>> >>>>> par(bg="white", las=1, cex=1.1) >>>>> plot(density(rg, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> >>>>> # Distribuição Weibull >>>>> >>>>> rw <- rweibull(100, scale=1, shape=1.5) >>>>> >>>>> par(bg="white", las=1, cex=1.1) >>>>> plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> >>>>> fit_weibull <- fitdistr(rw,"weibull") >>>>> >>>>> # Distribuição Exponential >>>>> >>>>> ex <- rexp(100, rate = 1) >>>>> >>>>> par(bg="white", las=1, cex=1.1) >>>>> plot(density(ex, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> >>>>> fit_exp <- fitdistr(ex, "exponential") >>>>> >>>>> plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> lines(density(rg, bw=0.5, cut=0),col="green") >>>>> lines(density(ex, bw=0.5, cut=0),col="red") >>>>> >>>>> # Distribuição Exponencial mista >>>>> >>>>> library(Renext) >>>>> >>>>> # f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)) >>>>> >>>>> mixexp <- rmixexp2(n = 100, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0) >>>>> >>>>> fit_mixexp <- mom.mixexp2(mixexp) >>>>> >>>>> plot(density(mixexp, bw=0.5, cut=1), lwd=2, >>>>> xlim=c(0, 5),col="steelblue") >>>>> lines(density(ex, bw=0.5, cut=0),col="red") >>>>> >>>>> mod <- nls(mixexp ~ prob1 * dexp(x, rate1) + (1 - prob1) * dexp(x, >>>>> rate2), >>>>> start = list(prob1 = 1, rate1 = 1, rate2 = 3.0), control = >>>>> nls.control(maxiter = 100)) >>>>> >>>>> Obrigado. >>>>> -- >>>>> Alisson Lucrecio da Costa >>>>> >>>>> _______________________________________________ >>>>> R-br mailing list >>>>> R-br(a)listas.c3sl.ufpr.br >>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>> código mínimo reproduzível. >>>>> >>>> >>>> >>> >>> >>> -- >>> Alisson Lucrecio da Costa >>> >> >> >> >> -- >> Alisson Lucrecio da Costa >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> >> >> > > > -- > Alisson Lucrecio da Costa > -- Alisson Lucrecio da Costa

1 0