- R-br - DINF UFPR

Sorting
by David Feitosa 06 Abr '15

06 Abr '15

Caros, boa noite. Uma dúvida de ordenação. Suponha que eu tenha 3 list compostas: um = list(id=1, m=matrix(1:9, nrow=3)) dois = list(id=2, m=matrix(1:4, nrow=2)) tres = list(id=3, m=matrix(1, nrow=1)) e crie um vetor com as 3 list: lista = c(um, dois, tres) Haveria alguma forma de ordená-las pelo tamanho do elemento m, ou seja, de forma que a lista final fosse equivalente a: c(tres, dois, um) ? Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

1 0

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Oi David, cheguei a tentar instalar por fora mas não consegui... como te disse pensei que seria talvez incompatibilidade com a versão mai recente do R... Vc lembra em qual versão do R usou o pacote? Tem como disponibilizar o link de onde baixou o arquivo zip com o pacote?? Obrigado!

2 1

Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Caros amigos, estou com uma dúvida simples entretanto não consigo solucionar o problema... Já procurei solucionar usando os pacotes R2Cuba e o cubature, entretanto sem sucesso... Quero resolver uma integral dupla entretanto o limite superior de uma das variáveis depende da outra variável, e nos pacotes listados acima consigo apenas colocar limites numéricos... Existe algo que eu possa fazer para resolver esse tipo de problema no R?? A função que desejo integrar é a seguinte: f=(2*exp(-2*x))*(3*exp(-3*y)) Resumindo quero resolver o seguinte: Inf x ∫ ∫ (2*exp(-2*x))*(3*exp(-3*y)) dy dx 0 0 Agradeço desde Já!!! Romero.

2 1

Como fazer esse grafico no ggplot2
by Pedro Emmanuel Alvarenga Americano do Brasil 06 Abr '15

06 Abr '15

Amigos de R, Essa semana eu decidi que iria aprender a usar o ggplot2. Logo na primeira tarefa empaquei uns 40 min. Como fazer esse grafico que faço no graphics com o ggplot2? > rec.cum b2 r3 r24 r25 1 1 2014-04-09 Sim 2014-04-09 10 2 2014-04-15 Não 2014-04-15 12 3 2014-04-16 Sim 2014-04-16 21 4 2014-05-16 Sim 2014-05-16 28 5 2014-05-28 Sim 2014-05-28 35 6 2014-07-29 Sim 2014-07-29 42 7 2014-08-27 Sim 2014-08-27 49 8 2014-10-07 Sim 2014-10-08 55 9 2014-10-22 Sim 2014-10-22 61 10 2014-10-27 Não 2014-10-27 78 11 2014-10-29 Sim 2014-10-29 84 12 2014-11-10 Sim 2014-11-10 89 13 2014-11-10 Sim 2014-11-10 93 14 2014-11-19 Não 2014-11-19 94 15 2014-12-02 Sim 2014-12-02 95 16 2014-12-08 Sim 2014-12-08 99 17 2015-02-25 Sim 2015-02-25 100 18 2015-03-03 Sim 2015-03-03 101 19 2015-03-10 Sim 2015-03-10 102 20 2015-03-17 Sim 2015-03-17 plot.recrut <- function(){ plot(rec.cum$r3,1:nrow(rec.cum),type='b',xlab='',ylab='Número de pacientes rastreados/recrutados',lty=2,col=2,pch=6,cex=1,lwd=3,xlim=c(min(rec.cum$r3,na.rm=T),max(rec.cum$r25,na.rm=T)),axes=F);grid();box();axis(2) tmp <- rec.cum[order(rec.cum$r25),] tmp <- tmp[which(tmp$r24=='Sim'),] lines(tmp$r25,1:nrow(tmp),type='b',lty=3,col=4,pch=2,cex=1,lwd=3) rm(tmp) axis.Date(1,rec.cum$r3,format='%d') # axis.Date(1,rec.cum$r3,format='%b',line=1.5,tick=F) axis.Date(1,rec.cum$r3,format='%Y',line=2.5,tick=F) legend('bottomrigh',c(paste('Primeiro rastreado:',min(rec.cum$r3,na.rm=T)),paste('Último rastreado:',max(rec.cum$r3,na.rm=T))),bty='n') legend('topleft',c(paste('Primeiro recrutado:',min(rec.cum$r25,na.rm=T)),paste('Último recrutado:',max(rec.cum$r25,na.rm=T))),bty='n') legend('top',c('Rastreamento','Recrutamento'),lty=c(2,3),col=c(2,4),pch=c(6,2),pt.cex=1.5,lwd=3,bty='n',xpd=T,inset=-.15,horiz=T) } plot.recrut() Abraço forte, Pedro Brasil

3 5

successive projection algorithm (SPA)
by Ari Clecius 01 Abr '15

01 Abr '15

Boa noite estou realizando uma seleção de variáveis para análise discriminante linear e gostaria de usar o algoritmo SPA, alguém saberia me dizer se existe implementação em R? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)84174333 (085)33669042

3 4

Mudar formatação de uma timestamp separando em colunas.
by Michelle Bau Graczyk 31 Mar '15

31 Mar '15

Boa tarde, Eu tenho um arquivo nesse formato e queria que a coluna de timestamp fosse separada em uma de data e outra de hora pois preciso calcular a media do volume por horário. Para isso eu tentei usar o programa abaixo mas a saida foi NA: > INTC<-read.table("/Users/bau/Documents/Mestrado/Dadosbaixadosdothebonnotgang/INTC_1m.txt", header=TRUE, sep=";", dec=",") > INTC[1:25,] symbol nr timestamp open high low close volume 1 INTC 1 2011-09-09 09:30:59 19.7500 19.86 19.75 19.8400 615100 2 INTC 2 2011-09-09 09:31:59 19.8350 19.86 19.79 19.8600 242900 3 INTC 3 2011-09-09 09:32:58 19.8501 19.86 19.80 19.8001 314400 4 INTC 4 2011-09-09 09:33:56 19.8068 19.81 19.77 19.7800 183500 5 INTC 5 2011-09-09 09:34:59 19.7900 19.81 19.78 19.7901 150900 6 INTC 6 2011-09-09 09:35:59 19.8000 19.81 19.74 19.7700 264600 7 INTC 7 2011-09-09 09:36:57 19.7700 19.81 19.77 19.8001 174300 8 INTC 8 2011-09-09 09:37:57 19.8100 19.81 19.76 19.7700 146200 9 INTC 9 2011-09-09 09:38:59 19.7600 19.83 19.76 19.8300 192100 10 INTC 10 2011-09-09 09:39:59 19.8400 19.87 19.83 19.8600 172500 11 INTC 11 2011-09-09 09:40:59 19.8700 19.87 19.76 19.7700 275000 12 INTC 12 2011-09-09 09:41:58 19.7700 19.77 19.72 19.7500 269200 13 INTC 13 2011-09-09 09:42:59 19.7500 19.82 19.75 19.8100 209600 14 INTC 14 2011-09-09 09:43:57 19.8197 19.82 19.79 19.7900 303200 15 INTC 15 2011-09-09 09:44:57 19.7900 19.82 19.78 19.7900 148000 16 INTC 16 2011-09-09 09:45:58 19.7900 19.79 19.75 19.7550 168700 17 INTC 17 2011-09-09 09:46:59 19.7600 19.85 19.75 19.8500 224200 18 INTC 18 2011-09-09 09:47:57 19.8594 19.88 19.85 19.8700 158600 19 INTC 19 2011-09-09 09:48:01 19.8800 19.91 19.87 19.9000 219500 20 INTC 20 2011-09-09 09:49:59 19.8900 19.92 19.89 19.9200 121000 21 INTC 21 2011-09-09 09:50:59 19.9200 19.93 19.90 19.9200 241400 22 INTC 22 2011-09-09 09:51:59 19.9100 19.95 19.90 19.9300 214000 23 INTC 23 2011-09-09 09:52:59 19.9300 19.94 19.92 19.9300 74200 24 INTC 24 2011-09-09 09:53:59 19.9292 19.93 19.88 19.9000 274400 25 INTC 25 2011-09-09 09:54:59 19.9000 19.92 19.87 19.9200 75800 > strptime(INTC$timestamp," %H:%M:%S") [1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [31] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [61] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [91] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [121] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [151] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [181] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [211] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [241] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [271] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [ reached getOption("max.print") -- omitted 335637 entries ] Muito obrigada!!!

2 3

Manipulação matriz com laço
by Michelle Bau Graczyk 31 Mar '15

31 Mar '15

Boa noite, Eu estou normalizando autovetores de uma matriz e colocando eles em outra. Para isso eu montei o programa a seguir. A primeira forma, ele só calcula o ultimo elemento da coluna 4 da matriz, e da segunda forma os valores do resultado estão errados. Primeira forma: > a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE) > b<-as.matrix(a) > matriznorm<-matrix(0,30,4) > matriznorm[,1]<-b[,1]# só vou usar este valores do arquivo acima > show(matriznorm[,1]) [1] 10.373173554 6.125265914 3.148949951 2.722401273 2.045041946 1.498078581 1.350122780 [8] 1.335457110 1.190051106 0.999610112 0.850212256 0.716720583 0.600704889 0.474433269 [15] 0.431826856 0.388576986 0.266709936 0.176960029 0.045935755 0.020126910 -0.003897135 [22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973 [29] -1.222537293 -1.273419528 > for( i in 1:30){ + matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1] + } > > somaQuadrados<-sum(matriznorm[,2]) > raiz<-sqrt(somaQuadrados) > matriznorm[,3]<-raiz > > for(j in i:30){ + matriznorm[j,4]<-matriznorm[j,1]/raiz + } > > write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt") > show(matriznorm) [,1] [,2] [,3] [,4] [1,] 10.373173554 1.076027e+02 13.46876 0.00000000 [2,] 6.125265914 3.751888e+01 13.46876 0.00000000 [3,] 3.148949951 9.915886e+00 13.46876 0.00000000 [4,] 2.722401273 7.411469e+00 13.46876 0.00000000 [5,] 2.045041946 4.182197e+00 13.46876 0.00000000 [6,] 1.498078581 2.244239e+00 13.46876 0.00000000 [7,] 1.350122780 1.822832e+00 13.46876 0.00000000 [8,] 1.335457110 1.783446e+00 13.46876 0.00000000 [9,] 1.190051106 1.416222e+00 13.46876 0.00000000 [10,] 0.999610112 9.992204e-01 13.46876 0.00000000 [11,] 0.850212256 7.228609e-01 13.46876 0.00000000 [12,] 0.716720583 5.136884e-01 13.46876 0.00000000 [13,] 0.600704889 3.608464e-01 13.46876 0.00000000 [14,] 0.474433269 2.250869e-01 13.46876 0.00000000 [15,] 0.431826856 1.864744e-01 13.46876 0.00000000 [16,] 0.388576986 1.509921e-01 13.46876 0.00000000 [17,] 0.266709936 7.113419e-02 13.46876 0.00000000 [18,] 0.176960029 3.131485e-02 13.46876 0.00000000 [19,] 0.045935755 2.110094e-03 13.46876 0.00000000 [20,] 0.020126910 4.050925e-04 13.46876 0.00000000 [21,] -0.003897135 1.518766e-05 13.46876 0.00000000 [22,] -0.030136936 9.082349e-04 13.46876 0.00000000 [23,] -0.074815367 5.597339e-03 13.46876 0.00000000 [24,] -0.145410630 2.114425e-02 13.46876 0.00000000 [25,] -0.294553778 8.676193e-02 13.46876 0.00000000 [26,] -0.453450244 2.056171e-01 13.46876 0.00000000 [27,] -0.551343912 3.039801e-01 13.46876 0.00000000 [28,] -0.710794973 5.052295e-01 13.46876 0.00000000 [29,] -1.222537293 1.494597e+00 13.46876 0.00000000 [30,] -1.273419528 1.621597e+00 13.46876 -0.09454616 > Da segunda forma: > a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE) > b<-as.matrix(a) > matriznorm<-matrix(0,30,4) > matriznorm[,1]<-b[,1] > show(matriznorm[,1]) [1] 10.373173554 6.125265914 3.148949951 2.722401273 2.045041946 1.498078581 1.350122780 [8] 1.335457110 1.190051106 0.999610112 0.850212256 0.716720583 0.600704889 0.474433269 [15] 0.431826856 0.388576986 0.266709936 0.176960029 0.045935755 0.020126910 -0.003897135 [22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973 [29] -1.222537293 -1.273419528 > for( i in 1:30){ + matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1] + + somaQuadrados<-sum(matriznorm[,2]) + raiz<-sqrt(somaQuadrados) + matriznorm[,3]<-raiz + + for(j in i:30){ + matriznorm[j,4]<-matriznorm[j,1]/raiz + } + } > write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt") > show(matriznorm) [,1] [,2] [,3] [,4] [1,] 10.373173554 1.076027e+02 13.46876 1.0000000000 [2,] 6.125265914 3.751888e+01 13.46876 0.5084624670 [3,] 3.148949951 9.915886e+00 13.46876 0.2528991576 [4,] 2.722401273 7.411469e+00 13.46876 0.2135962692 [5,] 2.045041946 4.182197e+00 13.46876 0.1584251428 [6,] 1.498078581 2.244239e+00 13.46876 0.1152793173 [7,] 1.350122780 1.822832e+00 13.46876 0.1033376904 [8,] 1.335457110 1.783446e+00 13.46876 0.1016853664 [9,] 1.190051106 1.416222e+00 13.46876 0.0902440218 [10,] 0.999610112 9.992204e-01 13.46876 0.0755856426 [11,] 0.850212256 7.228609e-01 13.46876 0.0641564605 [12,] 0.716720583 5.136884e-01 13.46876 0.0540043425 [13,] 0.600704889 3.608464e-01 13.46876 0.0452163558 [14,] 0.474433269 2.250869e-01 13.46876 0.0356888678 [15,] 0.431826856 1.864744e-01 13.46876 0.0324667092 [16,] 0.388576986 1.509921e-01 13.46876 0.0292025272 [17,] 0.266709936 7.113419e-02 13.46876 0.0200398900 [18,] 0.176960029 3.131485e-02 13.46876 0.0132951406 [19,] 0.045935755 2.110094e-03 13.46876 0.0034511674 [20,] 0.020126910 4.050925e-04 13.46876 0.0015121392 [21,] -0.003897135 1.518766e-05 13.46876 -0.0002927926 [22,] -0.030136936 9.082349e-04 13.46876 -0.0022641887 [23,] -0.074815367 5.597339e-03 13.46876 -0.0056207917 [24,] -0.145410630 2.114425e-02 13.46876 -0.0109238800 [25,] -0.294553778 8.676193e-02 13.46876 -0.0221227473 [26,] -0.453450244 2.056171e-01 13.46876 -0.0340370871 [27,] -0.551343912 3.039801e-01 13.46876 -0.0413498303 [28,] -0.710794973 5.052295e-01 13.46876 -0.0532327969 [29,] -1.222537293 1.494597e+00 13.46876 -0.0911767942 [30,] -1.273419528 1.621597e+00 13.46876 -0.0945461625 Muito obrigada!

2 4

Sub-matrizes em R
by Tiago Souza Marçal 31 Mar '15

31 Mar '15

Boa noite, alguém sabe como selecionar 2 ou mais sub-matrizes de uma matriz maior de forma direta sem usar um loop???? CMR: mat <- matrix(1:25,nc=5) z = c(1,4)lz = c(2,5)k <- numeric(0) for(j in 1:2) print( k <- mat[z[j]:lz[j],z[j]:lz[j]] ) # Gostaria de chegar a este resultado sem o loop Estou tentando evitar este loop, pois, ele esta dentro de while. Att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

2 2

Analise fatorial
by Bruce Kelly 30 Mar '15

30 Mar '15

Boa tarde colegas, Estou tentando usar a função 'factanal' para fazer uma analise fatorial ja utilizei mas a base de dados era um pouco, a mensagem de erro é essa Erro em cov.wt(z) : 'x' deve conter apenas valores finitos Alguém já se deparou com esse problema, ou sabe como resolver? Obrigado Bruce Kelly

2 2

Látice quadrado parcialmente balanceado
by Rodrigo Lima 30 Mar '15

30 Mar '15

Prezados listeiros Bom dia. Eu tenho um experimento em látice 8x8 com duas repetições. Fiz a análise de variância utilizando a função "PBIB" - Analysis of the Partially Balanced Incomplete Block Design. Script: exp01.gy<-PBIB.test(bloco,gen,rep,GY,k=8,group=TRUE, console=TRUE, method="REML") Minha dúvida é como fazer para estimar os resíduos desse modelo, uma vez que não estou conseguindo com a função "residuals". O que você me sugerem fazer? Desde já agradeço. Rodrigo. -- Prof. Rodrigo Oliveira de Lima Federal University of Viçosa Department of Plant Sciences Viçosa, MG, Brazil, 36570-900, Phone: +55(31)3899-2781 Fax: +55(31)3899-2614 Email: rodrigoodelima(a)ufv.br CV: http://lattes.cnpq.br/2761299616351806

2 4