Novembro 2015 - R-br

Re: [R-br] RES: Laço para gerar Subconjuntos
by sznelwar 12 Nov '15

12 Nov '15

1 0

[OFF-TOPIC] Teste sobre funcionamento da lista
by Walmes Zeviani 11 Nov '15

11 Nov '15

Mensagem de teste para avaliar se voltou a funcionar o serviço de email que hospeda a lista de discussão. W. ======================================================================== Walmes Marques Zeviani LEG - Laboratório de Estatística e Geoinformação [25.450418S 49.231759W] Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 linux user number: 531218 pages & profiles: http://www.leg.ufpr.br/~walmes | http://lattes.cnpq.br/4410617539281650 http://git.leg.ufpr.br/walmes | https://github.com/walmes https://www.mendeley.com/profiles/walmes-zeviani1/ ========================================================================

10 9

Re: [R-br] Laço para gerar Subconjuntos
by Robert Iquiapaza 11 Nov '15

11 Nov '15

Pode fazer com laço ou com funções da família apply: matriz<-matrix(NA,3,4)## numero de linhas é o numero de dias. colnames(matriz)<-c("open","close","low","higth") i<-1 for( x in unique(AA$X37988)){# x=37988 AA1<-subset(AA,AA$X37988==x) matriz[i,1]<-AA1[1,7] matriz[i,2]<-AA1[nrow(AA1),7] matriz[i,3]<-min(AA1[,7]) matriz[i,4]<-max(AA1[,7]) i=i+1 } matriz #ou t(sapply(split(AA,AA$X37988),function(AA1)c(AA1[1,7], AA1[nrow(AA1),7], min(AA1[,7]), max(AA1[,7])))) [,1] [,2] [,3] [,4] 37988 38.0 38.05 38.0 38.06 37991 37.9 37.90 37.9 37.90 37992 38.9 38.90 38.9 38.90 Att Em 11 de novembro de 2015 13:00, Michelle Bau Graczyk <mbgraczyk(a)gmail.com> escreveu: > Caros, bom dia… > > Estou com uma dúvida bem boba mas não consigo resolver! > > Eu tenho um data frame que a primeira coluna é a data, a segunda é o > horário e as outras seis são dados como preço e volume. > Os dias estão em forma numérica e eu os deixei assim porque achei que > seria mais fácil fazer o laço. O objetivo é que para cada daia eu ache o > valor inicial e final do preço assim como o seu máximo e seu mínimo. Estou > tendo problemas somente em trocar de dia. > Abaixo segue um pedaço reduzido do data frame e o código. > Alguém poderia me dar uma dica? > Muito obrigada! > AA<-read.table("/Users/bau/Dropbox/AA.N_ExemploReduzido.csv",header=TRUE > + ,dec=".",as.is=TRUE, sep=",") > > #AA > > dput(AA) > structure(list(X37988 = c(37988L, 37988L, 37988L, 37988L, 37988L, > 37988L, 37991L, 37991L, 37991L, 37991L, 37992L, 37992L, 37992L, > 37992L), X34244.059 = c(34251.038, 34263.36, 34296.73, 34301.369, > 34318.91, 34324.707, 34398.295, 34400.585, 34416.511, 34425.892, > 34250.443, 34252.127, 34266.757, 34273.928), X.1 = c(-1, -1, > 38.05, 38.06, 38.06, 38.04, 37.8, 37.8, 37.85, 37.85, -1, -1, > 38.85, 38.85), X.1.1 = c(-1L, -1L, 32L, 2L, 4L, 1L, 50L, 50L, > 1L, 1L, -1L, -1L, 1L, 1L), X.1.2 = c(-1, -1, 38.08, 38.08, 38.07, > 38.07, 38.3, 38.3, 37.95, 37.95, -1, -1, 38.95, 38.95), X.1.3 = c(-1L, > -1L, 1L, 1L, 1L, 1L, 203L, 203L, 1L, 1L, -1L, -1L, 1L, 1L), X38 = c(38, > 38, 38.05, 38.06, 38.06, 38.05, 37.9, 37.9, 37.9, 37.9, 38.9, > 38.9, 38.9, 38.9), X55900 = c(600L, 5000L, 6000L, 1500L, 500L, > 4500L, 161800L, 500L, 3700L, 1100L, 125000L, 1100L, 400L, 1000L > )), .Names = c("X37988", "X34244.059", "X.1", "X.1.1", "X.1.2", > "X.1.3", "X38", "X55900"), class = "data.frame", row.names = c(NA, > -14L)) > > > > matriz<-matrix(NA,230,4)## numero de linhas é o numero de dias. > > colnames(matriz)<-c("open","close","low","higth") > > i<-1 > > for( x in 37988:38000){ > + > + > + #x<-37988 > + > + if(AA$data==x){ > + AA1<-subset(AA,AA$data==x) > + matriz[i,1]<-AA1[1,7] > + matriz[i,2]<-AA1[nrow(AA1),7] > + matriz[i,3]<-min(AA1$last) > + matriz[i,4]<-max(AA1$last) > + > + i=i+1 > + } > + AA<-read.table("/Users/bau/Dropbox/AA.N_ExemploReduzido.csv",header=TRUE > + ,dec=".",as.is=TRUE, sep=",") > + #AA > + dput(AA) > + > + matriz<-matrix(NA,3,4)## numero de linhas é o numero de dias. > + colnames(matriz)<-c("open","close","low","higth") > + i<-1 > + for( x in 37988:37992){ > + > + > + > + > + if(AA$data==x){ > + AA1<-subset(AA,AA$data==x) > + matriz[i,1]<-AA1[1,7] > + matriz[i,2]<-AA1[nrow(AA1),7] > + matriz[i,3]<-min(AA1$last) > + matriz[i,4]<-max(AA1$last) > + > + i=i+1 > + } > + } > + matriz > + write.table(matriz,"Matriz4ParametrosPreço.txt") > + > > _______________________________________________ > 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

Intervalo de confiança das estimativas dos parâmetros das EDO geradas pela função modFit, do pacote FME
by Adriele Giaretta Biase 11 Nov '15

11 Nov '15

Como pode ser calculado o intervalo de confiança das estimativas dos parâmetros das EDO (ou modelos não lineares) geradas pela função modFit, do pacote FME? Abaixo segue o summay obtido Estimate Std. Error t value Pr(>|t|) Beta1 4.285 0.0041 0.411 0.00687 Beta2 0.082 0.0769 0.530 0.00605 Beta3 2.041 0.0900 0.411 0.00068 Eu devo usar a distribuição t-student retornado pela função, fazendo: Estimate ± Std. Error * t value Em que CI_95%(Beta1)= (4.285- 0.0041* 0.411 ; 4.285+ 0.0041* 0.411 ) CI_95%(Beta2)= (0.082- 0.0769* 0.530 ; 0.082+ 0.0769* 0.530) CI_95%(Beta3)= (2.041- 0.0900* 0.411; 2.041+ 0.0900* 0.411) Por que os valores da distribuição t estão variando, os graus de liberdade não deveriam ser os mesmos (n-p)? Eem que o n é o tamanho da amostra e p o número de parâmetros? Ou devo usar a distribuição normal Z(α/2): CI_95%(Beta1)= (4.285- 0.0041* 1.96 ; 4.285+ 0.0041* 1.96 ) Onde eu posso conseguir um material falando sobre o cálculo dos Intervalos de Confiança produzidos pela função modFit do pacote FME. O manual do pacote não deixa claro esse cálculo. O cálculo dos intervalos sempre serão baseados em distribuições simétricas? Adriele -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutoranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8861-0619.

1 0

Meetup R São Paulo Users Group (IME-USP) - 27 de novembro do 2015
by Luis Benites Sánchez 11 Nov '15

11 Nov '15

Caros, A sexta 27 novembro teremos nossa reunião número 7 do Meetup de usuários de R de São Paulo, no auditório Jacy Monteiro no Bloco B. Segue a programação desse dia: Programação da sexta *27 de novembro* (14h-16:00h) - [14:00 - 14:55] Uma introdução aos modelos mistos em R - Ana Paula Zerbeto (IME-USP) - [Presencial] - [14:55 - 15:05] Coffee Break - [15:05 - 16:00] Web scraping usando o R - Julio Trecenti (IME-USP) - [Presencial] As palestras serão transmitidas ao vivo pela conexão IPTV da USP e YouTube (*https://plus.google.com/u/0/events/cpjg070m652k6su1gkeuivrr6bc <https://plus.google.com/u/0/events/cpjg070m652k6su1gkeuivrr6bc>*) . Para poder participar, pode confirmar sua participação gratuitamente no site do Meetup :) Para mais informações: Site: www.meetup.com/pt/useR-SP Canal de YouTube: www.youtube.com/channel/UCb9PdJYYaElE6Uqb1KxhE8w e-mail: meetuprsp(a)gmail.com Um abraço, Luis Benites Sánchez Rocío Paola Maehara Organizadores do São Paulo R Users Group --

1 0

Trabalhar com subsets de data frame
by Michelle Bau Graczyk 11 Nov '15

11 Nov '15

Caros, bom dia, Eu estou com um data frame da forma : "data" "hora" "bid" "bidQ" "ask" "askQ" "last" "volume" "1" 2004-01-04 "09:30" -1 -1 -1 -1 38 55900 "2" 2004-01-04 "09:30" -1 -1 -1 -1 38 600 "3" 2004-01-04 "09:31" -1 -1 -1 -1 38 5000 "4" 2004-01-04 "09:31" 38.05 32 38.08 1 38.05 6000 "5" 2004-01-04 "09:31" 38.06 2 38.08 1 38.06 1500 "6" 2004-01-04 "09:31" 38.06 4 38.07 1 38.06 500 "7" 2004-01-04 "09:32" 38.04 1 38.07 1 38.05 4500 "8" 2004-01-04 "09:32" 38 1 38.2 1 38.03 1600 "9" 2004-01-04 "09:32" 38 1 38.2 1 38.03 3800 "10" 2004-01-04 "09:32" 38 1 38.04 1 38.04 500 "11" 2004-01-04 "09:32" 38 1 38.07 1 38.03 200 "12" 2004-01-04 "09:32" 38 1 38.08 2 38.07 100 "13" 2004-01-04 "09:32" 38 1 38.08 2 38.01 5000 "14" 2004-01-04 "09:32" 38 21 38.08 2 38.05 200 "15" 2004-01-04 "09:33" 38.03 9 38.13 8 38.08 700 "16" 2004-01-04 "09:33" 38.08 15 38.13 8 38.08 100 "17" 2004-01-04 "09:33" 38.11 16 38.13 10 38.12 900 onde as datas vão de 2004-01-04 até 2013-06-30 e os horários de 09:30 até 16:10. Eu tenho que achar o volume médio para cada minuto e depois normalizar esses volumes para cada minuto pelo volume médio. Eu tento montar um programa cheio de for() para data e for() para horário mas não está dando certo… Acho que deve ter uma lógica mais simples. Alguém saberia me dar uma dica? Muito obrigada! Michelle

2 2

equacao de S.E. para regressao
by Cleber Borges 10 Nov '15

10 Nov '15

Prezados UseRs, Usando o livro do Faraway (practical regression and anova using r) adaptei um trecho de código (abaixo) e consigo traçar um gráfico de modelo de regressão (reta) e com os intervalos de confiança (IC) da regressão. Com uso da função *predicit* eu obtenho os IC... Minha perguntinha é: Tem alguma formulazinha bonitinha e elegante para se obter esses *s.e.* da regressão tal como são bonitinhas e elegantes as formulazinhas de estimativas de coeficientes dos parâmetros e erros segundo OLS?? tipo: 1) b = ( X^t * X )^1 * X^t * y 2) Var( b ) = ( X^t * X )^1 * Var_residual 3) Algo do gênero de 1) e 2) ??? Agradeço pelo help com antecedência! Cleber ps.: vi uma formula mas é muita feia! ############ x <- rnorm( 100 ) y <- x + rnorm( 100, s=5 ) reg <- lm( y ~ x ) grid <- seq( -3, 3, 0.1 ) p <- predict( reg, data.frame( x=grid ), se=T ) IC <- cbind( p$fit - 1.96*p$se, p$fit, p$fit + 1.96*p$se ) matplot( grid, IC, lty=1, col=c(2,1,2), type="l", lwd=3 ) --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

2 2

Índice de concordância Kappa
by André Lucas de Oliveira Moreira 10 Nov '15

10 Nov '15

Olá pessoal, tudo bem? Aqui no trabalho, pela primeira vez me deparei com uma situação em que é necessário comparar a eficiência de técnicas que detectam que uma mesma espécie de parasita. A principio não sabia qual seria o melhor teste para verificar isso, mas me falaram pela primeira vez do kappa. Então tenho buscado sobre o assunto. Porém, nos livros que eu tive acesso fala disso muito rapidamente e não me deram segurança suficiente para executar no R e interpretar os resultados. Pensando que tenho resultados categóricos (presença e ausência dos parasitos para ambas técnicas), qual seria a melhor forma de analisar isso, melhor pacote e função? Também preciso do valor de p. Seria muito útil se alguém indicasse um tutorial e também algum material teórico que me dê mais segurança no uso desta técnica. Agradecido, André -- *MSc. André Lucas de O. Moreira* http://lattes.cnpq.br/7258065668864153 <http://www.wikiaves.com.br/perfil_andrelukinhas> http://www.wikiaves.com.br/perfil_andrelukinhas 79 8837-3562 79 9132-9093

5 10

Como aumentar o espaço de strip, para caber um nome de 2 linhas. bwplot (lattice)
by Augusto Ribas 09 Nov '15

09 Nov '15

Ola, estou com dificuldades de escrever o nome o titulo dos strips numa figura usando o lattice. Meu problema é o seguinte, eu tenho os seguintes dados, por exemplo resposta<-rnorm(90) preditor1<-rep(rep(c("a","b"),each=15),3) preditor2<-rep(c("sp1","sp2","sp3"),each=30) Dai preciso fazer uma figura, mas os nomes nos strips são espécies biologicas, então precisam estar em italico, resolvi da seguinte forma. library(lattice) bwplot(resposta~preditor1|preditor2, strip = strip.custom(factor.levels = c( expression(italic("P. paradoxa")), expression(italic("H. raniceps")), expression(italic("P. azurea"))) ) ) Mas além disso, gostaria de pular uma linha, e escrever uma segunda informação, mas essa segunda informação não é em italico, mas a primeira precisa. Pesquisando uma pouco, vi a possibilidade de usar o atop do expression, mas ai o texto não cabe no strip, como eu posso aumentar o tamanho do strip, para caber tudo? bwplot(resposta~preditor1|preditor2, strip = strip.custom(factor.levels = c( expression(atop(italic("P. paradoxa"),"outra info")), expression(atop(italic("H. raniceps"),"outra info")), expression(atop(italic("P. azurea"),"outra info"))) ) ) Bem é isso, agradeço desde ja a atenção. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Auxilio com Funções de Estatistica
by Yury Duarte 09 Nov '15

09 Nov '15

Bom dia! Meu nome é Yury Duarte e sou pós graduando em Engenharia de Biossistemas. Estou com uma necessidade de aplicar alguns índices estatísticos a uma serie de dados, entretanto, não consegui encontrar uma função no R que as faça. Os índices são: - R² - r - c - D - RMSE - Erro médio - Erro médio absoluto Dessa forma, gostaria de saber se já existem comandos para calcular esses índices ou se terei que desenvolver o código para cada um deles. Desde já, agradeço pela colaboração! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 3