- R-br - DINF UFPR

Distribuir o resultado de duas condições de ifelse em duas colunas
by ASANTOS 04 Jun '15

04 Jun '15

Bom dia Pessoal, Tenho um vetor de dados artificiais chamado dados e gostaria de testar se os valores de dados são <=0.05767 e jogar o resultado (sendo que se true =1, senão 0) em uma coluna e ao mesmo tempo testar se x<=0.89000 e colocar o resultados (sendo que se true =1, senão 0) em outra coluna, para isso fiz: dados=rnorm(10) a1=0.05767 a2=0.89000 RES<-NULL for (i in 1:length(dados)){ I1 <- x[i] x[i] <- as.vector(x[i]) ifelse(x[i] <= a1, 1, 0) I2 <- x[i] x[i] <- as.vector(x[i]) ifelse (x[i] <= a2, 1, 0) RES<-rbind(RES,c(I1,I2)) } RES Mas estou fazendo algo errado, pois são estou obtendo os resultados esperados nas colunas I1 e I2, alguém poderia me ajudar? Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

5 7

Shapefile sobre image.kriging
by Felinto COSTA 03 Jun '15

03 Jun '15

Prezados. Sabem alguma maneira de plotar um arquivo em formato shape sobre a imagem gerada pela função image.kriging (geoR)? A imagem que disponho já foi delimitada anteriormente na krigagem (em krige.cov, borders) mas o arquivo que desejo sobrepor possui delimitações internas (bairros de uma cidade) e o anterior apenas o perímetro externo da idade. Grato Felinto

2 1

equação no scatter3d
by Natalia Martins 03 Jun '15

03 Jun '15

Prezados boa tarde. Estou tentando colocar um modelo de regressao juntamente com o R^2 em um scatter3d, mas a equação sai como texto. Usei dois comandos "mtext3d" e "text3d" , mas não consegui inserir o modelo com nenhum. Se alguem puder me ajudar, desde já meu obrigada! scatter3d(y ~ x1 + x2) mtext3d(expression(paste(hat(y),'=7.5 + 1.59x[1] + 8.68x[2], ', R^2,'=0,9747')) text3d(expression(paste(hat(y),'=7.5 + 1.59x[1] + 8.68x[2], ', R^2,'=0,9747')) Natalia Martins

1 0

Otimizando desempenho do Servidor - dados grandes
by Pedro Emmanuel Alvarenga Americano do Brasil 02 Jun '15

02 Jun '15

Amigos de R, Eu venho com um dúvida de iniciante, porque nunca trabalhei assim. Eu não entendo muito de informática por isso estou sem luz no caminho. Bom, a minha questão é que me envolvi num projeto em que os dado são na caso de centenas de milhares de linhas, com umas 150 colunas. Estou trabalhando numa área de trabalho remota num servidor windows 2008. Eu não sei quais as especificaçõs do servidor, mas essa área esta servidor está servindo somente para esse fim. Então gostaria de uma luz ou do que ler para otimizar o tempo de processamento nas ediçoes e nas análises. Algum tutorial ou alguma dica que aproveite melhor o meu tempo de espera. O que ja fucei que é pouco... > memory.size(max = FALSE) [1] 1044.11 > memory.size(max = T) [1] 1814.94 > memory.limit(size = NA) [1] 4095 > Abraço, Pedro Brasil

3 4

Método monte carlo
by Fernando Antonio de souza 02 Jun '15

02 Jun '15

Caros amigos, Ajustei um modelo de regressão através de uma meta analise e gostaria de avaliar a sua capacidade preditiva. Infelizmente não tenho um banco de dados independente(não utilizado para ajustar modelo) para poder compara valores preditos e observados. Penso que o método de monte Carlo pode me ajudar a resolver este problema (conforme algumas leituras aleatorias pela internet e) . Gostaria de uma opinião sobre o uso desse método para esse fim. E se possível dicas ou compartilhamento de material didático( for dummies) para que eu possa introduzir ao estudo de monte Carlo. Att

2 6

Impressão do objeto Date
by David Feitosa 01 Jun '15

01 Jun '15

Boa noite! Estou deparando-me com um pequeno problema: Tenho um data.frame que desejo cortar de acordo com o valor de uma coluna. Esta coluna está em formato de data com hora: YYYY-MM-DD HH:MM:SS Minha intenção é de criar um diretório com a data e um arquivo dentro de cada diretório. Ex: <dir>/data/2015-02-02/saida.csv , <dir>/data/2015-02-03/saida.csv , <dir>/data/2015-02-04/saida.csv etc Entretanto, o que está sendo impresso, e, criado como diretório, é a representação interna do objeto Date Ex: <dir>/data/16467/saida.csv , <dir>/data/16468/saida.csv , <dir>/data/16469/saida.csv etc Já tentei usar as.character(), format(), mas estou ficando sem ideias de como "corrigir". Alguém já passou por isso ou teria alguma sugestão por favor? Segue trecho de código pertinente com um pedaço dos dados: library(R.utils) lines <- " 5928671,14,-3.831290,-38.50146,2015-02-01 00:00:01 5926098,44,-3.709390,-38.58986,2015-02-02 00:01:13 5915217,145,-3.743671,-38.48788,2015-02-03 00:01:16 5444641,174,-3.766271,-38.57452,2015-02-01 00:01:38 2250561,977,-3.738654,-38.48400,2015-02-02 00:01:47 5671754,14,-3.831290,-38.50146,2015-02-03 00:02:43 " con <- textConnection(lines) df <- read.csv(con, header = F) close(con) for (dia in unique(as.Date(df[,5]))){ dfi <- df[as.Date(df[,5]) == dia,] diretorio <- paste(getwd(), "rsc/data", format(dia, format = "%Y-%m-%d"), sep = "/", collapse = "") mkdirs(diretorio) write.csv(dfi, paste(diretorio, "saida.csv", sep = "/", collapse = "")) } Dados do ambiente: R version 3.2.0 (2015-04-16) -- "Full of Ingredients" Platform: i386-w64-mingw32/i386 (32-bit) Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

5 6

Re: [R-br] Área foliar
by Andre Oliveira 30 Mai '15

30 Mai '15

Continua reclamando da tiff, alguma sugestão? .......................................................... tentando a URL 'http://bioconductor.org/packages/3.1/bioc/src/contrib/EBImage_4.10.0.tar.gz'Content type 'application/x-gzip' length 4114806 bytes (3.9 MB)==================================================downloaded 3.9 MB * installing *source* package ‘tiff’ ...** package ‘tiff’ successfully unpacked and MD5 sums checked** libsgcc -std=gnu99 -I/usr/share/R/include -DNDEBUG -fpic -g -O2 -fstack-protector --param=ssp-buffer-size=4 -Wformat -Werror=format-security -D_FORTIFY_SOURCE=2 -g -c common.c -o common.oIn file included from common.c:1:0:common.h:5:18: fatal error: tiff.h: Arquivo ou diretório não encontrado #include <tiff.h> ^compilation terminated.make: ** [common.o] Erro 1ERROR: compilation failed for package ‘tiff’* removing ‘/home/andre/R/x86_64-pc-linux-gnu-library/3.2/tiff’ERROR: dependency ‘tiff’ is not available for package ‘EBImage’* removing ‘/home/andre/R/x86_64-pc-linux-gnu-library/3.2/EBImage’ 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 Sábado, 30 de Maio de 2015 9:18, Augusto Ribas <ribas.aca(a)gmail.com> escreveu: Opa, ele está reclamando que não tem todas as livrarias que precisa para compitar. Ta faltando do ffw3h.h Ta usando linux certo? No ubuntu, abra um terminal e instale os pacotes necessários com apt-get sudo apt-get install fftw3 fftw3-dev pkg-config depois tenta instalar denovo que vai funcionar. Em 30 de maio de 2015 07:48, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: Pois é o erro que ocorre ..mas consegui identificar com clareza os motivos! > source("http://bioconductor.org/biocLite.R")Bioconductor version 3.1 (BiocInstaller 1.18.2), ?biocLite for help> biocLite("EBImage")BioC_mirror: http://bioconductor.orgUsing Bioconductor version 3.1 (BiocInstaller 1.18.2), R version 3.2.0.Installing package(s) ‘EBImage’also installing the dependencies ‘tiff’, ‘fftwtools’ tentando a URL 'http://cran.rstudio.com/src/contrib/tiff_0.1-5.tar.gz'Content type 'application/x-gzip' length 28925 bytes (28 KB)==================================================downloaded 28 KB tentando a URL 'http://cran.rstudio.com/src/contrib/fftwtools_0.9-7.tar.gz'Content type 'application/x-gzip' length 145998 bytes (142 KB)==================================================downloaded 142 KB tentando a URL 'http://bioconductor.org/packages/3.1/bioc/src/contrib/EBImage_4.10.0.tar.gz'Content type 'application/x-gzip' length 4114806 bytes (3.9 MB)==================================================downloaded 3.9 MB * installing *source* package ‘tiff’ ...** package ‘tiff’ successfully unpacked and MD5 sums checked** libsgcc -std=gnu99 -I/usr/share/R/include -DNDEBUG -fpic -g -O2 -fstack-protector --param=ssp-buffer-size=4 -Wformat -Werror=format-security -D_FORTIFY_SOURCE=2 -g -c common.c -o common.oIn file included from common.c:1:0:common.h:5:18: fatal error: tiff.h: Arquivo ou diretório não encontrado #include <tiff.h> ^compilation terminated.make: ** [common.o] Erro 1ERROR: compilation failed for package ‘tiff’* removing ‘/home/andre/R/x86_64-pc-linux-gnu-library/3.2/tiff’* installing *source* package ‘fftwtools’ ...** package ‘fftwtools’ successfully unpacked and MD5 sums checked** libsgcc -std=gnu99 -I/usr/share/R/include -DNDEBUG -fpic -g -O2 -fstack-protector --param=ssp-buffer-size=4 -Wformat -Werror=format-security -D_FORTIFY_SOURCE=2 -g -c fftwtools.c -o fftwtools.offtwtools.c:28:18: fatal error: fftw3.h: Arquivo ou diretório não encontrado #include<fftw3.h> ^compilation terminated.make: ** [fftwtools.o] Erro 1ERROR: compilation failed for package ‘fftwtools’* removing ‘/home/andre/R/x86_64-pc-linux-gnu-library/3.2/fftwtools’ERROR: dependencies ‘tiff’, ‘fftwtools’ are not available for package ‘EBImage’* removing ‘/home/andre/R/x86_64-pc-linux-gnu-library/3.2/EBImage’ The downloaded source packages are in ‘/tmp/RtmpDc7Umf/downloaded_packages’Old packages: 'car', 'tcltk2', 'MASS', 'spatial'Update all/some/none? [a/s/n]: 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 Sexta-feira, 29 de Maio de 2015 21:03, Augusto Ribas <ribas.aca(a)gmail.com> escreveu: O pacote é do biocondutor, você tentou instalar pela função biocLite? #################################################### source("http://bioconductor.org/biocLite.R") biocLite("EBImage") Em 29 de maio de 2015 18:54, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: Boa noite,existe alguma forma de fazer o cálculo da área da folha de uma árvore no R? Tentei esta dica abaixo, mas não achei o pacote referido! https://ridiculas.wordpress.com/page/2/ 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. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/Github: https://github.com/SquiercgLattes: http://lattes.cnpq.br/7355685961127056 -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/Github: https://github.com/SquiercgLattes: http://lattes.cnpq.br/7355685961127056

2 3

Área foliar
by Andre Oliveira 29 Mai '15

29 Mai '15

Boa noite,existe alguma forma de fazer o cálculo da área da folha de uma árvore no R? Tentei esta dica abaixo, mas não achei o pacote referido! https://ridiculas.wordpress.com/page/2/ 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

2 1

Gerar observações
by Alexandre Loures 29 Mai '15

29 Mai '15

Boa noite! Tenho uma base, 1989-2010, sem observações para alguns anos. Por exemplo, ano ido idd code valor 2008 bra arg 0101 10 2007 bra pry 8508 20 2010 bra ury 1011 5 1992 bra swe 7989 30 1995 bra pry 8508 15 Sendo assim, gostaria de gerar tais observações. Ou seja, ano ido idd code valor 1989 bra arg 0101 0 1989 bra pry 8508 0 1989 bra ury 1011 0 1989 bra swe 7989 0 . . . . . . . . . . . . . . . 1995 bra arg 0101 0 1995 bra ury 1011 0 1995 bra swe 7989 0 Isto é, para o ano de 1995 há informações para: ido=bra e idd=pry. Então, não precisa ser gerada. 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

3 4

Multivariate Statistics: Exercises and Solutions
by Luiz Roberto Martins Pinto 28 Mai '15

28 Mai '15

Ola. Alguém tem o pdf deste livro? Multivariate Statistics: Exercises and Solutions Wolfgang K. Härdle <http://quantlet.de/index.php?p=searchResults&w=author&id=125>, Zdenek Hlávka <http://quantlet.de/index.php?p=searchResults&w=author&id=310> 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

2 2

Re: [R-br] Fatorial tiplo com dados adicionais
by Alisson Lucrécio 27 Mai '15

27 Mai '15

Tiago, Obrigado por sua dica. Entretanto, o pacote ExpDes não analisa apenas um único tratamento adicional, mas quando tenho dentro de um fator. No meu caso tenho um tratamento adicional dentro de solo. Att. 2015-05-10 0:28 GMT-03:00 Alisson Lucrécio <alisson.lucrecio(a)ifgoiano.edu.br >: > Tiago, > > Obrigado por sua dica. Entretanto, o pacote ExpDes não análisa apenas um > único tratamento adicional, mas quando tenho dentro de um fator. No meu > caso tenho um tratamento adicional dentro de solo. > > Att. > > 2015-05-09 23:42 GMT-03:00 Tiago de Souza [via R-br] < > ml-node+s2285057n4664460h96(a)n4.nabble.com>: > >> Te indico o pacote ExpDes se você quer uma rotina simples. >> >> Att. >> >> Tiago. >> >> Alisson Lucrécio <[hidden email] >> <http:///user/SendEmail.jtp?type=node&node=4664460&i=0>> escreveu: >> >> Caro Colegas da lista r-br, >> >> Preciso de uma ajudar para montar uma análise de variância para um >> fatorial tripo (Solo x Solução de Extração x Concentração) e um fator >> adicional água. >> >> Os dados fatoriais são compostos por 5 níveis de solos, 2 níveis de >> solução e 3 níveis de concentração. Os dados adicionais são compostos por 5 >> níveis de solo e 1 nível de solução (Água). >> >> Segue um exemplo abaixo. >> >> fat_data <- data.frame(Solo = gl(5, 18, la = LETTERS[1:5]), >> Sol = gl(2, 9, length = 90, la = c("PO4", "SO4")), >> Conc = gl(3, 3, length = 90, la = c(0.10, 0.25, >> 0.50))) >> >> Rep <- apply(data.frame(Solo = rep(1:5, each = 18), >> Sol = rep(rep(c(0.5, 1), each = 9), 5), >> Conc = rep(rep(c(0.4, 0.9, 1.5), each = 3), 10)), >> 1, sum) >> >> fat_data$Rep <- Rep + rnorm(length(Rep),0,0.5) >> >> adi_data <- data.frame(Solo = gl(5, 3, la = LETTERS[1:5]), >> Sol = rep("Água", 15)) >> >> Rep <- rep(1:5, each = 3) >> adi_data$Rep <- Rep + rnorm(length(Rep),0,0.5) >> >> Muito obrigado. >> >> -- >> Alisson Lucrecio da Costa >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664460&i=1> >> 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] <http:///user/SendEmail.jtp?type=node&node=4664460&i=2> >> 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-Fatorial-tiplo-com-dados-adicionais-… >> 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…> >> > > > > -- > Alisson Lucrecio da Costa > -- Alisson Lucrecio da Costa

2 3

Modelos Mistos
by David Matta 27 Mai '15

27 Mai '15

Prezados, Estou utilizando da função lme (pacote nlme) para ajustar um modelo misto. Gostaria de saber porém se é possível flexibilizar a estrutura de correlação dos efeitos aleatórios para modelos fora da função "pdMat" e suas respectivas classes ( "pdClass"). Posso (vide código a seguir) propor uma estrutura bloco diagonal, com blocos na classe "corClasses" ou outras formas e/ou funções que eu mesmo descreveria? lme(Y~X1+X2-1, data=dados.frame, random=list(subject=pdBlocked(list(corARMA (~Z1+Z2-1), corARMA(~ Z3+Z4-1))))) Grato David

2 2

Alinhamento de eixos 'y'
by Leonardo Monteiro 26 Mai '15

26 Mai '15

Prezados, bom dia. tenho o seguinte: ### barras <- 1:10 linhas <- 1/barras grafico<- barplot(barras, beside = T, xaxt = 'n', ylim = c(0,20), las = 2, ylab = expression (Harvested~area~(x~10^{6}~ha)), xlab = "years") axis(1, at = grafico, labels = 1990:1999, las = 3) par(new = T) with(dados, plot(linhas, ty = 'l', lwd = 2, axes = F, ann = F, ylim = c(0,1))) mtext('test number', side = 4, line = 3) x = seq(0,1, 0.1) axis(4,x, las = 2) ### entretanto, nao consigo alinhar os eixos 'y'.. alguem teria uma sugestão? obrigado, Leonardo.

2 1

Re: [R-br] Fwd: Pacote "survival"
by sznelwar 25 Mai '15

25 Mai '15

1 0

Cálculo da correlação entre vetores que apresentam Nas
by Michelle Bau Graczyk 25 Mai '15

25 Mai '15

Boa tarde, Eu estou calculando a correlação entre vetores que apresentam NAs. Eu usei a função use="pairwise.complete.obs" mas da o erro: Mensagens de aviso perdidas: In cor(AA[, i], get(Symbols[j])[, i], use = "pairwise.complete.obs") : the standard deviation is zero Eu preciso que ele calcule a correlação entre os pares de vectores que contém valores sem propagar os NAs. Alguém saberia me dizer se preciso usar outra função ou se estou fazendo algo errado? Abaixo segue o código e um exemplo de um dos banco de dados usados. Muito obrigada! AA<-read.table("/Users/bau/MatrizCorrelação25%/AA_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) AXP<-read.table("/Users/bau/MatrizCorrelação25%/AIG_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) AIG<-read.table("/Users/bau/MatrizCorrelação25%/AXP_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) BA<-read.table("/Users/bau/MatrizCorrelação25%/BA_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) CAT<-read.table("/Users/bau/MatrizCorrelação25%/C_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) C<-read.table("/Users/bau/MatrizCorrelação25%/CAT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) DD<-read.table("/Users/bau/MatrizCorrelação25%/DD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) DIS<-read.table("/Users/bau/MatrizCorrelação25%/DIS_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) GE<-read.table("/Users/bau/MatrizCorrelação25%/GE_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) GM<-read.table("/Users/bau/MatrizCorrelação25%/GM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) HD<-read.table("/Users/bau/MatrizCorrelação25%/HD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) HON<-read.table("/Users/bau/MatrizCorrelação25%/HON_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) HPQ<-read.table("/Users/bau/MatrizCorrelação25%/HPQ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) IBM<-read.table("/Users/bau/MatrizCorrelação25%/IBM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) INTCc<-read.table("/Users/bau/MatrizCorrelação25%/INTC_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) INTC<-INTCc[,91:485] JNJ<-read.table("/Users/bau/MatrizCorrelação25%/JNJ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) JPM<-read.table("/Users/bau/MatrizCorrelação25%/JPM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) KO<-read.table("/Users/bau/MatrizCorrelação25%/KO_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MCD<-read.table("/Users/bau/MatrizCorrelação25%/MCD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MMM<-read.table("/Users/bau/MatrizCorrelação25%/MMM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MO<-read.table("/Users/bau/MatrizCorrelação25%/MO_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) MRK<-read.table("/Users/bau/MatrizCorrelação25%/MRK_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MSFTt<-read.table("/Users/bau/MatrizCorrelação25%/MSFT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) MSFT<-MSFTt[,91:485] PFE<-read.table("/Users/bau/MatrizCorrelação25%/PFE_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt", header=TRUE) PG<-read.table("/Users/bau/MatrizCorrelação25%/PG_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) SBC<-read.table("/Users/bau/MatrizCorrelação25%/SBC_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) UTX<-read.table("/Users/bau/MatrizCorrelação25%/UTX_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) VZ<-read.table("/Users/bau/MatrizCorrelação25%/VZ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) WMT<-read.table("/Users/bau/MatrizCorrelação25%/WMT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) XOM<-read.table("/Users/bau/MatrizCorrelação25%/XOM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE) matrizAA<-matrix(nrow=395, ncol=31) Symbols<-c("AA","AIG","AXP","BA","C","CAT","DD","DIS","GE","GM","HD","HON","HPQ","IBM","INTC","JNJ","JPM","KO","MCD","MMM","MO","MRK","MSFT","PFE","PG","SBC","UTX","VZ","WMT","XOM") for(j in 1:30) { empresa<-Symbols[j] for (i in 2:395){ matrizAA[i-1,j+1]<-cor(AA[,i],get(Symbols[j])[,i], use="pairwise.complete.obs") ##show(get(Symbols[j])[,i]) } } matrizAA[,1]<-c(1:395) write.table(matrizAA,"/Users/bau/MatrizCorrelação25%/AAcorrelaçãoParaMatriz25.txt") > show(INTC) X09.31 X09.32 X09.33 X09.34 X09.35 X09.36 X09.37 X09.38 X09.39 X09.40 1 95991 NA NA NA NA NA NA NA NA NA 2 NA NA NA NA NA NA NA NA NA NA 3 165142 288968 NA NA 108902 NA NA NA NA NA 4 118675 NA NA 106849 NA NA NA NA NA NA 5 143731 NA NA NA NA NA NA NA NA NA 6 182431 NA 320616 121284 102200 110047 177755 NA 89746 NA 7 385596 233253 102160 98434 NA 107502 125117 NA NA NA 8 164819 NA NA NA 82725 NA NA NA NA NA 9 1037347 1391429 1168113 2144977 996927 555185 399196 665563 911069 594260 10 375979 229744 123207 142073 389064 164749 NA NA NA 208184 11 734943 258833 133831 228149 173972 148795 108966 182753 109337 211194 12 302175 95693 NA 131569 NA NA 148767 185680 103596 NA 13 301830 NA NA NA NA NA NA NA 96945 111397 14 380442 175182 NA NA 154135 107285 NA 94698 82688 NA 15 257121 122878 NA NA 113273 NA 170090 NA NA NA 16 135805 NA NA NA NA 96105 230759 178512 98562 152584 17 143954 105240 NA 135897 NA 162295 NA 89958 175599 171491 18 110088 NA NA 232589 91542 NA 90570 118529 NA 176981 19 133450 NA NA NA 103523 NA NA NA NA NA 20 151424 535320 124938 202910 150950 228580 NA NA NA NA 21 174497 146076 379007 230754 NA 259161 262723 141953 NA 209095 22 202101 125722 178210 241326 222106 NA NA NA NA 234913 23 193474 81358 126753 157524 NA 140525 NA 132483 122650 94888 24 226467 NA 106886 NA NA NA NA NA NA NA 25 94991 NA NA NA NA 100226 NA 138565 NA NA X09.41 X09.42 X09.43 X09.44 X09.45 X09.46 X09.47 X09.48 X09.49 X09.50 X09.51 1 NA NA NA NA NA NA NA 125230 NA NA NA 2 NA NA NA NA NA NA NA NA NA NA NA 3 95210 NA 83615 NA NA NA NA NA 81950 NA NA 4 NA NA NA NA 525228 NA NA NA NA 122300 NA 5 NA NA NA NA NA NA NA NA NA NA NA 6 NA NA NA 118085 NA NA NA NA NA NA NA 7 101348 NA NA 85359 NA NA NA 84274 95450 90829 NA 8 NA NA NA NA NA NA NA NA NA 92065 NA 9 475210 1588546 828340 591364 849922 604366 592396 500498 231407 622671 633772 10 NA NA 82904 115703 133982 600720 NA NA 139664 142184 92034 11 91429 84681 288650 NA 182264 134730 NA 234712 82140 NA 138098 12 NA NA NA 103916 130810 164900 94492 282866 117227 94789 NA 13 NA NA NA NA NA NA NA NA NA NA NA 14 NA 290261 82617 NA NA 91485 NA NA NA NA NA 15 NA NA NA NA NA NA 114630 120390 NA 150594 NA 16 191759 NA NA NA NA 111460 212266 135440 NA NA 143878 17 163178 NA NA 181344 94675 NA 124844 NA NA NA 89650 18 NA NA NA NA NA NA NA NA NA NA NA 19 168241 NA 204800 NA NA NA NA NA NA NA NA 20 382726 197013 261650 NA 193270 191328 162645 NA NA 187996 134218 21 257324 248162 391465 212684 125445 101610 228834 NA NA 226410 128451 22 NA NA 209705 NA NA NA NA NA NA NA 88528 23 NA 145010 NA NA NA 86450 106323 NA 127235 131685 NA 24 NA NA NA NA NA NA NA NA NA NA NA 25 NA NA NA NA NA NA NA NA NA NA NA X09.52 X09.53 X09.54 X09.55 X09.56 X09.57 X09.58 X09.59 X10.00 X10.01 X10.02 1 NA NA NA NA NA NA 189095 NA NA NA NA 2 NA NA 212714 124657 NA 125427 NA NA NA NA NA 3 NA NA NA NA NA NA 108377 NA NA NA NA 4 NA NA NA NA NA NA NA NA NA NA NA 5 NA NA NA NA NA NA NA NA 111310 NA NA 6 135987 NA NA NA NA NA NA NA NA NA NA 7 NA 143816 86290 91866 NA NA 186791 NA NA NA NA 8 NA NA NA NA NA NA NA NA NA NA NA 9 398724 348269 209850 192490 427251 419138 879114 231852 289164 164084 560886 10 152704 NA 115500 NA 91505 NA 211354 107150 135024 NA 110234 11 88621 212824 194126 147968 NA NA 113810 NA NA NA NA 12 NA 85738 NA 84984 NA NA NA NA NA NA NA 13 94627 NA NA NA NA NA NA NA NA NA NA 14 107388 NA NA NA 160841 97684 291855 NA NA NA NA 15 NA 148828 NA 370458 NA NA NA NA 135246 NA NA 16 NA 187592 NA NA 87650 NA NA NA NA NA NA 17 NA 138967 117446 NA NA NA NA NA NA NA NA 18 230329 261154 NA NA 153873 NA NA NA NA NA NA 19 NA NA NA NA NA NA NA NA NA NA NA 20 295116 163870 145028 231132 154587 166770 286393 143550 190173 NA 146736 21 NA NA NA NA NA NA 91722 325090 NA NA NA 22 NA NA NA NA NA NA 112804 NA NA NA NA 23 NA NA NA NA NA NA NA NA NA NA 81763 24 NA 125616 NA NA NA NA NA NA NA NA NA 25 NA NA NA NA NA NA 84216 NA NA NA NA

1 0

Pacote "survival"
by Stephan M. CARVALHO 25 Mai '15

25 Mai '15

Boa tarde a todos, Gostaria de saber se alguém tem experiência com o pacote "survival". Estou com dificuldades para extrair os valores de "forma e escala" quando uso a função /survreg /assumindo a distribuição de Weibull. Obrigado! Stephan -- *************************************************** STEPHAN M. CARVALHO "Chemical control and ecotoxicology" Universidade Federal de Uberlândia Instituto de Ciências Agrárias Avenida Amazonas s/n – Bloco 2E Caixa Postal 593 - Umuarama 38.400-902 Uberlândia/MG Brazil stephan.carvalho(a)iciag.ufu.br +55 34 3225-8444 (university) +55 34 3218-2225 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! ***************************************************

2 3

Re: [R-br] Fwd: Pacote "survival"
by Robert Iquiapaza 24 Mai '15

24 Mai '15

http://www.leg.ufpr.br/doku.php/software:rbr?&#guia_de_postagem_da_r-br_em_… http://r-br.2285057.n4.nabble.com/ Em 24 de maio de 2015 09:26, Stephan M. CARVALHO < stephan.carvalho(a)iciag.ufu.br> escreveu: > Bom dia a todos! > > Semana passada fiz o post abaixo e, como não obtive retorno e mesmo > por não saber se foi da maneira correta, estou reencaminhando. > > Obrigado, > > Stephan > > > > -------- Forwarded Message -------- Assunto: Pacote "survival" Data: Sun, > 17 May 2015 14:12:11 -0300 De: Stephan M. CARVALHO > <stephan.carvalho(a)iciag.ufu.br> <stephan.carvalho(a)iciag.ufu.br> Para: Lista > R <r-br(a)listas.c3sl.ufpr.br> <r-br(a)listas.c3sl.ufpr.br> > > Boa tarde a todos, > > Gostaria de saber se alguém tem experiência com o pacote "survival". > > Estou com dificuldades para extrair os valores de "forma e escala" > quando uso a função *survreg *assumindo a distribuição de Weibull. > > Obrigado! > > Stephan > > > -- > > > *************************************************** > > STEPHAN M. CARVALHO > > > "Chemical control and ecotoxicology" > > > Universidade Federal de Uberlândia > Instituto de Ciências Agrárias > > Avenida Amazonas s/n – Bloco 2E > Caixa Postal 593 - Umuarama > > 38.400-902 Uberlândia/MG > Brazil > > > stephan.carvalho(a)iciag.ufu.br > > +55 34 3225-8444 (university)+55 34 3218-2225 (fax)+55 35 9103-9822 (mobile)+55 35 9949-0707 (mobile) > > > Antes de imprimir, pense no ambiente! > Before printing, think in the environment! > Avant d'imprimer, pensez à l'environnement! > > *************************************************** > > > > > > > > > > _______________________________________________ > 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. >

2 1

Cotações Históricas BMF e R
by Rodolfo Timoteo da Silva 23 Mai '15

23 Mai '15

Senhores, alguém já fez importação destes dados da BMF para usar com o R? Ou conhece algum "passo a passo"? Antecipadamente. -- --- Zhu

1 0

OFF TOPIC
by Andre Oliveira 23 Mai '15

23 Mai '15

Boa noite,alguém da lista usa e-mail yahoo? Meu e-mail yahoo não esta fazendo busca desde obtém, tenho este e-mail para a lista da Rbr e não consigo fazer buscas, alguém com mesmo problema? Sei que posso buscar no histórico da lista, mas entre os e-mails antigos desde a R_STAT não consigo buscar desde ontém em meu webmail!

3 4

Uso da função aggregate com NAS
by Michelle Bau Graczyk 21 Mai '15

21 Mai '15

Boa noite, Eu tenho um arquivo com NAs na coluna de volume, e quando tento usar a função aggregate ele retorna tudo NAs. Já tente usar o na.action como na.pass e na.omit e não obtive êxito. Eu queria que as linhas que tivessem o volume NA não fossem contatas na hora de agregar e calcular a média. Abaixo segue o programa: > file<-read.table("/Users/bau/MatrizCorrelação25%/AA_teste25.txt",header=TRUE) > file[1:300,] data hora preço volume 1 01.07.2004 09:31:00 32.87 NA 2 01.07.2004 09:32:00 32.89 NA 3 01.07.2004 09:33:00 32.82 NA 4 01.07.2004 09:34:00 32.82 NA 5 01.07.2004 09:35:00 32.87 NA 6 01.07.2004 09:36:00 32.94 NA 7 01.07.2004 09:37:00 32.95 NA 8 01.07.2004 09:38:00 32.88 NA 9 01.07.2004 09:39:00 32.90 NA 10 01.07.2004 09:40:00 32.90 NA 11 01.07.2004 09:41:00 32.89 NA 12 01.07.2004 09:42:00 32.88 NA 13 01.07.2004 09:43:00 32.89 NA 14 01.07.2004 09:44:00 32.86 NA 15 01.07.2004 09:45:00 32.88 NA 16 01.07.2004 09:46:00 32.83 NA 17 01.07.2004 09:47:00 32.81 NA 18 01.07.2004 09:48:00 32.83 NA 19 01.07.2004 09:49:00 32.83 NA 20 01.07.2004 09:50:00 32.91 NA 21 01.07.2004 09:51:00 32.99 NA 22 01.07.2004 09:52:00 32.98 NA 23 01.07.2004 09:53:00 32.91 NA 24 01.07.2004 09:54:00 32.92 NA 25 01.07.2004 09:55:00 32.94 NA 26 01.07.2004 09:56:00 32.93 NA 27 01.07.2004 09:57:00 32.91 NA 28 01.07.2004 09:58:00 32.93 NA 29 01.07.2004 09:59:00 32.94 NA 30 01.07.2004 10:00:00 32.98 NA 31 01.07.2004 10:01:00 32.95 NA 32 01.07.2004 10:02:00 32.95 NA 33 01.07.2004 10:03:00 32.97 NA 34 01.07.2004 10:05:00 33.00 NA 35 01.07.2004 10:06:00 32.97 NA 36 01.07.2004 10:07:00 32.99 NA 37 01.07.2004 10:08:00 32.94 NA 38 01.07.2004 10:09:00 32.90 NA 39 01.07.2004 10:10:00 32.89 NA 40 01.07.2004 10:11:00 32.86 NA 41 01.07.2004 10:12:00 32.85 NA 42 01.07.2004 10:13:00 32.86 NA 43 01.07.2004 10:14:00 32.84 NA 44 01.07.2004 10:15:00 32.81 NA 45 01.07.2004 10:16:00 32.79 NA 46 01.07.2004 10:17:00 32.75 NA 47 01.07.2004 10:18:00 32.74 NA 48 01.07.2004 10:19:00 32.75 NA 49 01.07.2004 10:20:00 32.82 NA 50 01.07.2004 10:21:00 32.84 NA 51 01.07.2004 10:22:00 32.80 NA 52 01.07.2004 10:24:00 32.80 NA 53 01.07.2004 10:26:00 32.82 NA 54 01.07.2004 10:27:00 32.84 NA 55 01.07.2004 10:28:00 32.88 NA 56 01.07.2004 10:29:00 32.88 NA 57 01.07.2004 10:30:00 32.88 NA 58 01.07.2004 10:31:00 32.88 NA 59 01.07.2004 10:32:00 32.87 NA 60 01.07.2004 10:33:00 32.88 NA 61 01.07.2004 10:34:00 32.95 NA 62 01.07.2004 10:35:00 32.94 NA 63 01.07.2004 10:36:00 32.93 NA 64 01.07.2004 10:37:00 32.91 NA 65 01.07.2004 10:38:00 32.93 NA 66 01.07.2004 10:40:00 32.93 NA 67 01.07.2004 10:41:00 33.00 NA 68 01.07.2004 10:42:00 32.99 NA 69 01.07.2004 10:43:00 33.00 NA 70 01.07.2004 10:44:00 33.05 NA 71 01.07.2004 10:45:00 33.03 NA 72 01.07.2004 10:46:00 33.02 NA 73 01.07.2004 10:47:00 33.02 NA 74 01.07.2004 10:48:00 33.02 NA 75 01.07.2004 10:49:00 33.02 NA 76 01.07.2004 10:50:00 33.04 NA 77 01.07.2004 10:51:00 33.03 NA 78 01.07.2004 10:52:00 33.04 NA 79 01.07.2004 10:53:00 33.03 NA 80 01.07.2004 10:54:00 33.03 NA 81 01.07.2004 10:55:00 33.00 NA 82 01.07.2004 10:56:00 32.96 NA 83 01.07.2004 10:57:00 32.96 NA 84 01.07.2004 10:58:00 32.95 NA 85 01.07.2004 10:59:00 32.90 NA 86 01.07.2004 11:00:00 32.90 NA 87 01.07.2004 11:02:00 32.81 32000 88 01.07.2004 11:03:00 32.75 15300 89 01.07.2004 11:04:00 32.72 NA 90 01.07.2004 11:05:00 32.73 NA 91 01.07.2004 11:06:00 32.74 NA 92 01.07.2004 11:07:00 32.70 11100 93 01.07.2004 11:08:00 32.69 NA 94 01.07.2004 11:09:00 32.63 NA 95 01.07.2004 11:10:00 32.64 NA 96 01.07.2004 11:11:00 32.63 NA 97 01.07.2004 11:12:00 32.55 NA 98 01.07.2004 11:13:00 32.50 NA 99 01.07.2004 11:14:00 32.56 NA 100 01.07.2004 11:15:00 32.57 NA 101 01.07.2004 11:16:00 32.55 NA 102 01.07.2004 11:17:00 32.56 NA 103 01.07.2004 11:18:00 32.62 NA 104 01.07.2004 11:19:00 32.64 NA 105 01.07.2004 11:20:00 32.60 NA 106 01.07.2004 11:21:00 32.59 NA 107 01.07.2004 11:22:00 32.56 13800 108 01.07.2004 11:23:00 32.62 NA 109 01.07.2004 11:24:00 32.65 NA 110 01.07.2004 11:25:00 32.61 NA 111 01.07.2004 11:26:00 32.67 NA 112 01.07.2004 11:27:00 32.71 NA 113 01.07.2004 11:28:00 32.71 NA 114 01.07.2004 11:29:00 32.74 NA 115 01.07.2004 11:30:00 32.75 NA 116 01.07.2004 11:31:00 32.75 NA 117 01.07.2004 11:32:00 32.75 NA 118 01.07.2004 11:33:00 32.75 NA 119 01.07.2004 11:34:00 32.75 NA 120 01.07.2004 11:35:00 32.76 NA 121 01.07.2004 11:36:00 32.76 NA 122 01.07.2004 11:37:00 32.76 NA 123 01.07.2004 11:38:00 32.75 NA 124 01.07.2004 11:39:00 32.76 NA 125 01.07.2004 11:40:00 32.77 NA 126 01.07.2004 11:41:00 32.76 NA 127 01.07.2004 11:42:00 32.75 NA 128 01.07.2004 11:43:00 32.76 NA 129 01.07.2004 11:44:00 32.75 NA 130 01.07.2004 11:45:00 32.71 45700 131 01.07.2004 11:46:00 32.68 NA 132 01.07.2004 11:47:00 32.65 NA 133 01.07.2004 11:48:00 32.67 NA 134 01.07.2004 11:49:00 32.67 NA 135 01.07.2004 11:50:00 32.64 NA 136 01.07.2004 11:51:00 32.64 NA 137 01.07.2004 11:52:00 32.67 NA 138 01.07.2004 11:53:00 32.68 NA 139 01.07.2004 11:54:00 32.68 NA 140 01.07.2004 11:55:00 32.70 NA 141 01.07.2004 11:56:00 32.69 NA 142 01.07.2004 11:57:00 32.69 NA 143 01.07.2004 11:58:00 32.69 NA 144 01.07.2004 11:59:00 32.72 NA 145 01.07.2004 12:00:00 32.71 NA 146 01.07.2004 12:01:00 32.66 NA 147 01.07.2004 12:02:00 32.68 15300 148 01.07.2004 12:03:00 32.67 NA 149 01.07.2004 12:04:00 32.65 NA 150 01.07.2004 12:05:00 32.61 14300 151 01.07.2004 12:06:00 32.61 NA 152 01.07.2004 12:07:00 32.61 NA 153 01.07.2004 12:08:00 32.61 NA 154 01.07.2004 12:09:00 32.63 NA 155 01.07.2004 12:10:00 32.62 11000 156 01.07.2004 12:11:00 32.62 NA 157 01.07.2004 12:12:00 32.62 NA 158 01.07.2004 12:13:00 32.62 NA 159 01.07.2004 12:14:00 32.62 NA 160 01.07.2004 12:15:00 32.63 NA 161 01.07.2004 12:16:00 32.62 NA 162 01.07.2004 12:17:00 32.62 NA 163 01.07.2004 12:18:00 32.63 NA 164 01.07.2004 12:19:00 32.65 NA 165 01.07.2004 12:20:00 32.61 NA 166 01.07.2004 12:21:00 32.61 NA 167 01.07.2004 12:22:00 32.61 NA 168 01.07.2004 12:23:00 32.62 NA 169 01.07.2004 12:24:00 32.61 NA 170 01.07.2004 12:25:00 32.62 NA 171 01.07.2004 12:26:00 32.62 NA 172 01.07.2004 12:27:00 32.61 NA 173 01.07.2004 12:28:00 32.62 NA 174 01.07.2004 12:29:00 32.62 12500 175 01.07.2004 12:30:00 32.64 NA 176 01.07.2004 12:31:00 32.63 NA 177 01.07.2004 12:32:00 32.63 NA 178 01.07.2004 12:33:00 32.64 10800 179 01.07.2004 12:34:00 32.66 25100 180 01.07.2004 12:35:00 32.69 NA 181 01.07.2004 12:36:00 32.70 NA 182 01.07.2004 12:37:00 32.70 NA 183 01.07.2004 12:38:00 32.70 NA 184 01.07.2004 12:39:00 32.68 NA 185 01.07.2004 12:40:00 32.68 NA 186 01.07.2004 12:41:00 32.67 NA 187 01.07.2004 12:42:00 32.65 NA 188 01.07.2004 12:43:00 32.63 NA 189 01.07.2004 12:44:00 32.62 NA 190 01.07.2004 12:45:00 32.62 NA 191 01.07.2004 12:46:00 32.62 NA 192 01.07.2004 12:47:00 32.62 NA 193 01.07.2004 12:48:00 32.62 21600 194 01.07.2004 12:49:00 32.60 NA 195 01.07.2004 12:50:00 32.58 NA 196 01.07.2004 12:51:00 32.56 NA 197 01.07.2004 12:52:00 32.55 NA 198 01.07.2004 12:53:00 32.55 NA 199 01.07.2004 12:54:00 32.56 30100 200 01.07.2004 12:55:00 32.58 NA 201 01.07.2004 12:56:00 32.57 NA 202 01.07.2004 12:57:00 32.57 NA 203 01.07.2004 12:58:00 32.53 NA 204 01.07.2004 12:59:00 32.51 NA 205 01.07.2004 13:00:00 32.52 NA 206 01.07.2004 13:01:00 32.52 NA 207 01.07.2004 13:02:00 32.52 NA 208 01.07.2004 13:03:00 32.53 NA 209 01.07.2004 13:04:00 32.51 NA 210 01.07.2004 13:05:00 32.47 NA 211 01.07.2004 13:06:00 32.43 NA 212 01.07.2004 13:07:00 32.40 NA 213 01.07.2004 13:08:00 32.37 NA 214 01.07.2004 13:09:00 32.39 11900 215 01.07.2004 13:10:00 32.39 NA 216 01.07.2004 13:11:00 32.40 NA 217 01.07.2004 13:12:00 32.41 NA 218 01.07.2004 13:13:00 32.46 NA 219 01.07.2004 13:14:00 32.44 NA 220 01.07.2004 13:15:00 32.41 NA 221 01.07.2004 13:16:00 32.42 NA 222 01.07.2004 13:17:00 32.42 NA 223 01.07.2004 13:18:00 32.42 NA 224 01.07.2004 13:19:00 32.41 NA 225 01.07.2004 13:20:00 32.41 NA 226 01.07.2004 13:21:00 32.41 NA 227 01.07.2004 13:22:00 32.41 NA 228 01.07.2004 13:23:00 32.36 NA 229 01.07.2004 13:24:00 32.34 22800 230 01.07.2004 13:25:00 32.34 NA 231 01.07.2004 13:26:00 32.33 NA 232 01.07.2004 13:27:00 32.33 NA 233 01.07.2004 13:28:00 32.34 NA 234 01.07.2004 13:29:00 32.33 NA 235 01.07.2004 13:30:00 32.33 NA 236 01.07.2004 13:31:00 32.33 NA 237 01.07.2004 13:32:00 32.33 NA 238 01.07.2004 13:33:00 32.30 NA 239 01.07.2004 13:34:00 32.33 16000 240 01.07.2004 13:35:00 32.32 12500 241 01.07.2004 13:36:00 32.33 NA 242 01.07.2004 13:37:00 32.31 19100 243 01.07.2004 13:38:00 32.34 22600 244 01.07.2004 13:39:00 32.34 NA 245 01.07.2004 13:40:00 32.34 NA 246 01.07.2004 13:41:00 32.33 NA 247 01.07.2004 13:42:00 32.31 NA 248 01.07.2004 13:43:00 32.31 NA 249 01.07.2004 13:44:00 32.31 NA 250 01.07.2004 13:45:00 32.30 NA 251 01.07.2004 13:46:00 32.30 NA 252 01.07.2004 13:47:00 32.28 NA 253 01.07.2004 13:48:00 32.25 19200 254 01.07.2004 13:49:00 32.25 21200 255 01.07.2004 13:50:00 32.23 25300 256 01.07.2004 13:51:00 32.20 NA 257 01.07.2004 13:52:00 32.20 15400 258 01.07.2004 13:53:00 32.20 17500 259 01.07.2004 13:54:00 32.17 11700 260 01.07.2004 13:55:00 32.16 NA 261 01.07.2004 13:56:00 32.14 17900 262 01.07.2004 13:57:00 32.14 16300 263 01.07.2004 13:58:00 32.17 19000 264 01.07.2004 13:59:00 32.19 15600 265 01.07.2004 14:00:00 32.23 25100 266 01.07.2004 14:01:00 32.20 NA 267 01.07.2004 14:02:00 32.24 NA 268 01.07.2004 14:03:00 32.26 12500 269 01.07.2004 14:04:00 32.31 19800 270 01.07.2004 14:05:00 32.31 NA 271 01.07.2004 14:06:00 32.30 NA 272 01.07.2004 14:07:00 32.33 NA 273 01.07.2004 14:08:00 32.33 14000 274 01.07.2004 14:09:00 32.33 NA 275 01.07.2004 14:10:00 32.33 NA 276 01.07.2004 14:11:00 32.34 13300 277 01.07.2004 14:12:00 32.34 11200 278 01.07.2004 14:13:00 32.34 NA 279 01.07.2004 14:14:00 32.33 NA 280 01.07.2004 14:15:00 32.33 NA 281 01.07.2004 14:16:00 32.33 NA 282 01.07.2004 14:17:00 32.32 NA 283 01.07.2004 14:18:00 32.28 NA 284 01.07.2004 14:19:00 32.26 NA 285 01.07.2004 14:20:00 32.26 NA 286 01.07.2004 14:21:00 32.27 NA 287 01.07.2004 14:22:00 32.28 NA 288 01.07.2004 14:23:00 32.29 NA 289 01.07.2004 14:24:00 32.30 NA 290 01.07.2004 14:25:00 32.31 NA 291 01.07.2004 14:26:00 32.32 NA 292 01.07.2004 14:27:00 32.34 NA 293 01.07.2004 14:28:00 32.35 10100 294 01.07.2004 14:29:00 32.37 NA 295 01.07.2004 14:30:00 32.36 NA 296 01.07.2004 14:31:00 32.37 NA 297 01.07.2004 14:32:00 32.35 NA 298 01.07.2004 14:33:00 32.37 14100 299 01.07.2004 14:34:00 32.31 NA 300 01.07.2004 14:35:00 32.30 NA > agg <- aggregate.data.frame(x=file, by=list(file$hora), FUN=mean) Muito obrigada! Michelle

2 2

Definição da linha com os nomes das variaveis
by SILVIA PATRICIA BARRETO SANTANA 21 Mai '15

21 Mai '15

Boa noite a todos! Tenho uma base de dados com extensão .txt obtida a partir de uma planilha do Excel, separadas as colunas por espaço. Segue uma pequena amostra abaixo Data Nivel1 Nivel2 01/01/2014 00:00 29 26 01/01/2014 00:01 29 26 01/01/2014 00:23 29 26 01/01/2014 00:49 29 26 01/01/2014 00:50 29 26 01/01/2014 00:57 29 26 01/01/2014 00:58 30 26 01/01/2014 00:59 30 26 01/01/2014 01:00 30 26 01/01/2014 01:01 30 26 05/01/2014 14:42 29 46 05/01/2014 14:43 I/O_Timeout I/O_Timeout 05/01/2014 14:44 I/O_Timeout I/O_Timeout 05/01/2014 14:45 I/O_Timeout I/O_Timeout 05/01/2014 14:46 I/O_Timeout I/O_Timeout 05/01/2014 14:47 I/O_Timeout I/O_Timeout 05/01/2014 14:48 29 45 05/01/2014 14:49 29 45 Para carregar o banco no R, estou usando o seguinte comando: dados<-read.delim("Dados.txt", header=TRUE, sep ="", na.strings=c("Bad","I/O_Timeout","Error","Failed","Intf_Shut"), dec=".", row.names = NULL) Problema: a variavel Data retorna apenas com as horas e minutos e a informação do dia aparece como row.names. Gostaria que a Data mantivesse o formato '%d/%m/%Y %H:%M'. Se eu remover 'row.names = NULL' aparece a seguinte informação de erro Error in read.table(file = file, header = header, sep = sep, quote = quote, : duplicate 'row.names' are not allowed Alguém pode me ajudar? Grata. Silvia

1 1

Ajuste Regressão Robusta
by Manoel Vítor de Souza Veloso 20 Mai '15

20 Mai '15

Boa tarde. Preciso fazer diagnósticos de resíduos para uma regressão robusta. Como por exemplo, plotar a distância de cook para cada observação amostral. No ajuste com "lm", basta rodar: plot(cooks.distance(lm(y~x))) Já na Regressão Robusta, usando o comando "rlm" da biblioteca MASS, o comando gráfico não roda: plot(cooks.distance(rlm(y~x))) Tem algum comando no R? -- ______________________________________________ *MANOEL VITOR DE SOUZA VELOSO----------------------------------------* *ProfessorInstituto de Ciências Sociais AplicadasUniversidade Federal de Alfenas - MGCampus Avançado de Varginha-MGhttp://www.unifal-mg.edu.br/icsa/ <http://www.unifal-mg.edu.br/icsa/>-------Doutor* * em Estatística e ExperimentaçãoUniversidade Federal de Lavras - MG----------------------------------------Tel: 35 - 8865.6116 35 - 3219.8705* P Antes de imprimir pense em sua responsabilidade e compromisso com o* MEIO AMBIENTE* !!!

2 1

Re: [R-br] Definição da linha com os nomes das variaveis
by Robert Iquiapaza 20 Mai '15

20 Mai '15

Seria isso dados<-read.delim("Dados.txt", header=TRUE, sep ="", na.strings=c("Bad","I/O_Timeout","Error","Failed","Intf_Shut"), dec=".", row.names = NULL,as.is=c(2)) str(dados) colnames(dados)[1:2]=c("Dia","Hora") dados$Data=as.POSIXct(paste(dados$Dia, dados$Hora), format="%d/%m/%Y %H:%M") dados Boa sorte Em 20 de maio de 2015 18:09, SILVIA PATRICIA BARRETO SANTANA < silpat(a)ufrb.edu.br> escreveu: > Boa noite a todos! > > Tenho uma base de dados com extensão .txt obtida a partir de uma planilha > do Excel, separadas as colunas por espaço. > Segue uma pequena amostra abaixo > > Data Nivel1 Nivel2 > 01/01/2014 00:00 29 26 > 01/01/2014 00:01 29 26 > 01/01/2014 00:23 29 26 > 01/01/2014 00:49 29 26 > 01/01/2014 00:50 29 26 > 01/01/2014 00:57 29 26 > 01/01/2014 00:58 30 26 > 01/01/2014 00:59 30 26 > 01/01/2014 01:00 30 26 > 01/01/2014 01:01 30 26 > 05/01/2014 14:42 29 46 > 05/01/2014 14:43 I/O_Timeout I/O_Timeout > 05/01/2014 14:44 I/O_Timeout I/O_Timeout > 05/01/2014 14:45 I/O_Timeout I/O_Timeout > 05/01/2014 14:46 I/O_Timeout I/O_Timeout > 05/01/2014 14:47 I/O_Timeout I/O_Timeout > 05/01/2014 14:48 29 45 > 05/01/2014 14:49 29 45 > > Para carregar o banco no R, estou usando o seguinte comando: > > dados<-read.delim("Dados.txt", header=TRUE, sep ="", > na.strings=c("Bad","I/O_Timeout","Error","Failed","Intf_Shut"), dec=".", > row.names = NULL) > > Problema: a variavel Data retorna apenas com as horas e minutos e a > informação do dia aparece como row.names. Gostaria que a Data mantivesse o > formato > > '%d/%m/%Y %H:%M'. > > Se eu remover 'row.names = NULL' aparece a seguinte informação de erro > > Error in read.table(file = file, header = header, sep = sep, quote = > quote, : > duplicate 'row.names' are not allowed > > > Alguém pode me ajudar? > > > Grata. > > Silvia > > _______________________________________________ > 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

outro erro likfit
by Hélio Gallo Rocha 18 Mai '15

18 Mai '15

Caros da lista, Com certeza tem algo nos dados, pois além do erro na krigagem apareceu este ao estimar por verossimillhança, com comando normal, sempre usei: LIKFIT.reml <- likfit(gd, ini = ini.grid ,lik.method ="REML",nug =grid.pepita ) onde: ini.grid é uma matriz com combinações do paramar e alcance grid.pepita são várias "chutes" de nugge gera este erro: Erro em optim(par = c(1000, 0.04), fn = function (pars, fp, ip, temp.list) : valor não finito fornecido por optim SE alguém puder ter uma ideia... -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 3

Dúvida ao juntar duas bases de dados
by Victor Eduardo 18 Mai '15

18 Mai '15

Bom dia, pessoal! Estou com um problema na hora de juntar a base de dados. Estou com duas bases com as seguintes dimensões: > dim(dados) [1] 2056 11 > dim(dados1) [1] 1782 2 Ambas as bases possuem uma variável chamada "id". O que queria é formar um subconjunto a partir de dados1 baseado nesses ids, ou seja, aproveitaria as variáveis da base "dados" com as variáveis em "dados1". Eu tentei com o merge e não está dando certo. Existe alguma função que faça isso? Abraços e ótima semana!

8 9

Instalar R no slackware 64 bits
by Carlos Andrade 18 Mai '15

18 Mai '15

Prezados amigos da lista. Fiz um upgrade no notebook para 4 GB de RAM. Instalei em dual boot Windows 7 64bits e Slackware 64 bits. Minha dúvida é sobre a versão do R para instalar no Slack. No repositório do SlackBuilds tem a versão 3.1.2 do R , mas não informa se é 64bits. No Windows 7 o instalador do R deu a opção de 32 bits e 64 bits. Depois da instalação apareceram dois ícones do R na área de trabalho: um 32 bits e outro 64bits. Agradeço qualquer dica pois estou iniciando no mundo 64 bits -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

3 3

Maximizar função
by Andre Oliveira 17 Mai '15

17 Mai '15

Pessoal,como devo escrever a função fxy para que o optim pegue fxy e maximize mesma? Estou pesquisando aqui, mas não consegui. Este é um exemplo apenas para fins didáticos! x = seq(-2,2,length=20) y = seq(-2,2,length=20) f = function(x,y){ fxy = -x^2 - y^2 return(fxy) } z = outer(x,y,f) persp(x,y,z,theta=30,phi=30, col="green",xlab="x",ylab="Y",zlab="Produção") # Gostaria de acrescentar a escala no eixo! fxy=function(x,y) -x^2 - y^2 optim(c(2,2),fxy, method="BFGS", control=list(fnscale=-1)) obrigado

2 1

Mudança da fonte de grafico
by Leonardo Monteiro 17 Mai '15

17 Mai '15

Prezados, eu havia mandado um email anteriormente sobre as fontes e tal.. entretanto acho que nao fui bem compreendido. Eu estou usando o expression() para adicionar caracteres especiais no xlab do grafico.. entretanto, eu nao estou conseguindo modificar a fonte. Modificar a fonte seria o tipo da fonte - da fonte padrão dos graficos (que eu acho que é curie news' para "times new roman', por exemplo. alguem teria alguma dica? obrigado, Leonardo

3 3

Modificar fonte do expressio()
by Leonardo Monteiro 17 Mai '15

17 Mai '15

Prezados, estou fazendo uns graficos e usando o expression() para xlab e ylab. Das outras partes do gráfico (legendas, etc..) eu já consegui. Entretanto nao consigo modificar a fonte.. (fica a fonte padrão) Alguem poderia me ajudar? obrigado.

2 1

Re: [R-br] Dúvida sobre escrito em xlab
by Tiago Souza Marçal 16 Mai '15

16 Mai '15

expression ("t ha"^-1) Att. Tiago. Cesar Rabak <cesar.rabak(a)gmail.com> escreveu: Veja no help (e no demo do R) "expression" 2015-05-15 17:08 GMT-03:00 Leonardo Monteiro <monteiroleonardo6(a)gmail.com>: > Senhores, boa noite. > > Como modifico o formato da legenda: t/ha para t ha-1? > Sendo que o -1 deve ficar sobreescrito. > > obrigado, > > Leonardo > > _______________________________________________ > 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.

1 0

erro na krigagem
by Hélio Gallo Rocha 16 Mai '15

16 Mai '15

Caros listeiros. Rodando krigagem de clima, com retirada de tendência polinomial de segunda ordem: KC.2nd <- krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0)) deu esse erro: Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 Além disso: Mensagens de aviso perdidas: 'overlayPointsWithPolygons' is deprecated. Use 'over' instead. See help("Deprecated") and help("sp-deprecated"). Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

4 6

Dúvida sobre escrito em xlab
by Leonardo Monteiro 15 Mai '15

15 Mai '15

Senhores, boa noite. Como modifico o formato da legenda: t/ha para t ha-1? Sendo que o -1 deve ficar sobreescrito. obrigado, Leonardo

3 2

eficiência de mem?==?UTF-8?Q?ória no R
by Fatima do Nascimento Silva 15 Mai '15

15 Mai '15

prezados, boa noite. Venho por meio deste solicitar uma dica de como fazer a leitura(importação) de uma base de dados em formato txt, com caracteres separados por ; e com 29 colunas e 1.520.171 linhas no R, de uma maneira que não perca eficiência de memória. estou fazendo a leitura usando o comando: (dados <- read.table("SAJ_WIN1252.txt",header=T,dec=",",sep=";")) Como o arquivo tem 3Gb, o R demora uns 10 min para abri-lo, mas não consigo fazer mais nada porque o programa trava. Alguém teria dicas de como me ajudar? agradeço a atenção. Atenciosamente, Fátima Nascimento Estatística - UFRN Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

3 6

Como fazer para dois graficos do ggplot ficarem lado a lado
by Pedro Emmanuel Alvarenga Americano do Brasil 12 Mai '15

12 Mai '15

Amigos de R, Estou mais uma vez batendo cabeça com o ggplot. Eu gostaria que doi graficos ficasse lado a lado na mesma janela, mas infelizmente ainda estou pensando com o plot.default() e não estou conseguindo. Alguma luz? d <- data.frame(ID=rep(1:10,100),outcome=c(rep(0,500),rep(1,500)),Sex=c(sample(c("F","M"),500,T,c(.8,.2)),sample(c("M","F"),500,T,c(.8,.2))),score=c(rnorm(500,5,10),rnorm(500,20,10)),grave=c("Sim","Não")) library("lme4") fit <- glmer(outcome ~ Sex + score + (1 | ID), data = d, family = binomial,subset=grave=="Sim") summary(fit) fit2 <- glmer(outcome ~ Sex + score + (1 | ID), data = d, family = binomial,subset=grave=="Não") summary(fit2) library("sjPlot") sjp.glmer(fit,type='fe') sjp.glmer(fit2,type='fe') Que a força esteja sempre com voces! Pedro Brasil

2 2

Primeiro contato
by Kamioka 12 Mai '15

12 Mai '15

Olá pessoal, boa noite. Já estudo data science/analysis há alguns meses usando R e python. Fiquei feliz em saber que temos uma lista brasileira. Eduardo

2 1

[OFF-TOPIC] Avaliação geral de modelos em SEM
by Vinícius Lionel Mateus 12 Mai '15

12 Mai '15

Olá a todos, Alguém tem uma boa sugestão de bibliografia para avaliação de "model fit" em SEM? Tenho uma amostra relativamente grande (454 x 25) e que não segue a multinormal. Minhas referências básicas têm sido Kline (2014) e Ralph (1996). Gostaria de investigar a magnitude e fontes de /misfit/ entre minhas matrizes de covariância, oriundas da amostra e impostas pelo modelo. Desde já, muito obrigado. Abs, Vinícius -- Best regards, Vinícius Lionel Mateus, MSc (http://lattes.cnpq.br/6501001637020665) Chemistry PhD Student Atmospheric Chemistry Laboratory - Dep. Chemistry PUC-Rio - Pontifical Catholic University of Rio de Janeiro Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brazil CEP: 22451-900 Telefone: (+45) 26 28 28 51 (+55) (21) 3527-1327 (+55) (21) 993 - 588 - 051 Skype: vinicius.lionel http://www.qui.puc-rio.br/index.html

1 0

Knitr iniciante
by Carlos Andrade 12 Mai '15

12 Mai '15

Prezados amigos da lista. Baixei e instalei o pacote knitr. Estou usando o RStudio versão 99.9.9 no Slackware 14.1, o R é da versão 3.1.2. Em uma pasta que fiz alguns exercicios e salvei um gráfico no formato png iniciei o RStudio e abri um arquivo testeKnitr.Rmd. Executei knitHTML e funcionou. Substitui o plot(cars) do arquivo demo que funcionou por plot(luglikp1) com um arquivo que criei. Executei knitHTML e recebi a mensagem de erro: ## Error in plot(luglikp1): objeto 'luglikp1' não encontrado Também deu erro no Windows 7, R version 3.2.0 Que ajustes estão faltando para funcionar corretamente? -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

3 2

javareconf + windows + R 3.2
by Leonard de Assis 11 Mai '15

11 Mai '15

Olá, O manual do R (http://cran.r-project.org/doc/manuals/r-release/R-admin.html#Java-support) diz o seguinte: <PRINT> Note that this unfortunately depends on the exact version of the JRE/JDK installed, and so may need updating if the Java installation is updated. This can be done by running R CMD javareconf which updates settings in both etc/Makeconf and R_HOME/etc/ldpaths. See R CMD javareconf --help for details. </PRINT> Contudo, ao executar tal comando no Windows, recebo a seguinte msg: <PRINT> C:\Program Files\R\R-3.2.0\bin>R CMD javareconf 'javareconf' não é reconhecido como um comando interno ou externo, um programa operável ou um arquivo em lotes. </PRINT> Ao executar o Help do comando, deparo com o seguinte: <PRINT> C:\Program Files\R\R-3.2.0\bin>R --arch x64 --help Usage: Rterm [options] [< infile] [> outfile] [EnvVars] Start R, a system for statistical computation and graphics, with the specified options EnvVars: Environmental variables can be set by NAME=value strings Options: -h, --help Print usage message and exit --version Print version info and exit --encoding=enc Specify encoding to be used for stdin --encoding enc ditto --save Do save workspace at the end of the session --no-save Don't save it --no-environ Don't read the site and user environment files --no-site-file Don't read the site-wide Rprofile --no-init-file Don't read the .Rprofile or ~/.Rprofile files --restore Do restore previously saved objects at startup --no-restore-data Don't restore previously saved objects --no-restore-history Don't restore the R history file --no-restore Don't restore anything --vanilla Combine --no-save, --no-restore, --no-site-file, --no-init-file and --no-environ --max-mem-size=N Set limit for memory to be used by R --max-ppsize=N Set max size of protect stack to N -q, --quiet Don't print startup message --silent Same as --quiet --slave Make R run as quietly as possible --verbose Print more information about progress --internet2 Use Internet Explorer settings for proxies etc. --args Skip the rest of the command line --ess Don't use getline for command-line editing and assert interactive use -f file Take input from 'file' --file=file ditto -e expression Use 'expression' as input One or more -e options can be used, but not together with -f or --file An argument ending in .RData (in any case) is taken as the path to the workspace to be restored (and implies --restore) Or: R CMD command args where 'command' is one of: INSTALL Install add-on packages REMOVE Remove add-on packages SHLIB Make a DLL for use with dynload BATCH Run R in batch mode build Build add-on packages check Check add-on packages Rprof Post process R profiling files Rdconv Convert Rd format to various other formats Rdiff difference R output files Rd2pdf Convert Rd format to PDF Rd2txt Convert Rd format to pretty text Stangle Extract S/R code from vignette Sweave Process vignette documentation config Obtain configuration information about R open Open a file via Windows file associations texify Process a latex file Use R CMD command --help for usage information for each command. </PRINT> A questão é a seguinte: Foi abandonado o suporte R ao java na versão 3.2, visto que tal comando agora não existe? Leonard

3 6

Java 8 e bibliotecas
by Andre Oliveira 11 Mai '15

11 Mai '15

Pessoal,configurei tudo aqui a xlsx, gdata, XLConnect entre outras que dependem do java, deu tudo certo, mas tive que colocar java 8 por causa do site da caixa e ai voltou tudo a zero e não consigo instalar de novo. Alguma dica? 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

2 2

Re: [R-br] Estimativas de parâmetros de uma Weibull utilizando duas metodologias: Newton-Raphson "vs" Nelder Mead (estimativas muito diferentes) o que pode estar acontecendo??
by Romero Luiz M. Sales Filho 11 Mai '15

11 Mai '15

Olá Wagner, segui o seu conselho e coloquei para o caso Nelder-Mead as estimativas iniciais como sendo as estimativas encontradas no caso de Newton-Raphson, o valor foi realmente quase em cima! Entaõ o problema realmente era esse, a seleção dos pontos iniciais. Fica então uma dúvida: Quando tenho dist. de probabilidades que não conseguimos estimar via Newton Raphson, por exemplo, como podemos ter uma idéias dos melhores pontos iniciais a serem utilizados?? Tem como encontrar alguma forma de se encontrar algo interessante ou é somente por tentativa e erro? Elias, muito obrigado também por compartilhar informação! Acho que vc deu uma dica de como encontrar os melhores pontos iniciais, né? Infelizmente não entendi... vc poderia me explicar melhor... Ah, outra coisa, com relação ao ks.test tbm tinha dúvidas se ele realmente era adequado, uma vez que não tenho os parâmetros e sim estimativas desses, aí gero uma amostra utilizando essas estimativas e isso talvez não fosse o mais indicado pq se o algoritomo tivesse convergido para valores errados de parâmetros e mesmo assim a distribuição osse Weibull o ajuste seria excelente quando na verdade os parâmetros estariam errados, né isso (por favor, me corrija se eu estiver errado)??... Então vc colocou uma abordagem utilizando bootstrap, essa seria uma forma alternativa de utilização do ks.test sendo mais apropriada para esse tipo de problema em que não conheço os parâmetros?? Muito obrigado à todos!! Abraço, Romero. -- *Profº Romero Luiz Mendonça Sales Filho* *Estatístico / Mestre em Engª de ProduçãoUFRPE - UAGCel.: (81) 9667-4477* *romero.sfilho(a)gmail.com* <romero.sfilho(a)gmail.com> *romero_sfilho(a)yahoo.com.br* <romero_sfilho(a)yahoo.com.br>

3 3

Instalação - Ubuntu 64 bits
by Andre Oliveira 10 Mai '15

10 Mai '15

Pessoal apesar de não usar sempre instalo o Rcmdr para apresentar a quem tenha resistência ao uso do R. O fato de ser 64 meu linux e R pode haver problemas na instalação do Rcmdr? 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

Re: [R-br] Fatorial tiplo com dados adicionais
by Tiago Souza Marçal 09 Mai '15

09 Mai '15

Te indico o pacote ExpDes se você quer uma rotina simples. Att. Tiago. Alisson Lucrécio <alisson.lucrecio(a)ifgoiano.edu.br> escreveu: Caro Colegas da lista r-br, Preciso de uma ajudar para montar uma análise de variância para um fatorial tripo (Solo x Solução de Extração x Concentração) e um fator adicional água. Os dados fatoriais são compostos por 5 níveis de solos, 2 níveis de solução e 3 níveis de concentração. Os dados adicionais são compostos por 5 níveis de solo e 1 nível de solução (Água). Segue um exemplo abaixo. fat_data <- data.frame(Solo = gl(5, 18, la = LETTERS[1:5]), Sol = gl(2, 9, length = 90, la = c("PO4", "SO4")), Conc = gl(3, 3, length = 90, la = c(0.10, 0.25, 0.50))) Rep <- apply(data.frame(Solo = rep(1:5, each = 18), Sol = rep(rep(c(0.5, 1), each = 9), 5), Conc = rep(rep(c(0.4, 0.9, 1.5), each = 3), 10)), 1, sum) fat_data$Rep <- Rep + rnorm(length(Rep),0,0.5) adi_data <- data.frame(Solo = gl(5, 3, la = LETTERS[1:5]), Sol = rep("Água", 15)) Rep <- rep(1:5, each = 3) adi_data$Rep <- Rep + rnorm(length(Rep),0,0.5) Muito 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.

1 0

Fatorial tiplo com dados adicionais
by Alisson Lucrécio 09 Mai '15

09 Mai '15

Caro Colegas da lista r-br, Preciso de uma ajudar para montar uma análise de variância para um fatorial tripo (Solo x Solução de Extração x Concentração) e um fator adicional água. Os dados fatoriais são compostos por 5 níveis de solos, 2 níveis de solução e 3 níveis de concentração. Os dados adicionais são compostos por 5 níveis de solo e 1 nível de solução (Água). Segue um exemplo abaixo. fat_data <- data.frame(Solo = gl(5, 18, la = LETTERS[1:5]), Sol = gl(2, 9, length = 90, la = c("PO4", "SO4")), Conc = gl(3, 3, length = 90, la = c(0.10, 0.25, 0.50))) Rep <- apply(data.frame(Solo = rep(1:5, each = 18), Sol = rep(rep(c(0.5, 1), each = 9), 5), Conc = rep(rep(c(0.4, 0.9, 1.5), each = 3), 10)), 1, sum) fat_data$Rep <- Rep + rnorm(length(Rep),0,0.5) adi_data <- data.frame(Solo = gl(5, 3, la = LETTERS[1:5]), Sol = rep("Água", 15)) Rep <- rep(1:5, each = 3) adi_data$Rep <- Rep + rnorm(length(Rep),0,0.5) Muito obrigado. -- Alisson Lucrecio da Costa

1 0

Re: [R-br] Emparelhamento de arquivos
by Éder Comunello 09 Mai '15

09 Mai '15

Desculpe, Mauro! Troque xts2 por xts3: merge.xts(xts1, xts2, join="inner") # Data: # numeric(0) # # Index: # integer(0) merge.xts(xts1, xts3, join="inner") # Open High Low Close Open.1 High.1 Low.1 Close.1 # 2007-01-04 50.42096 50.42096 50.26414 50.33236 50.42096 50.42096 50.26414 50.33236 # 2007-01-05 50.37347 50.37347 50.22103 50.33459 50.37347 50.37347 50.22103 50.33459 # 2007-01-06 50.24433 50.24433 50.11121 50.18112 50.24433 50.24433 50.11121 50.18112 Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 8 de maio de 2015 23:42, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > Por que eu não estou conseguindo rodar este seu exemplo? > > > > >

1 0

Emparelhamento de arquivos
by Edilson Santos 08 Mai '15

08 Mai '15

Caros usuários, Preciso emparelhar dados de dois arquivos "xts", por data, sendo que em cada um dos arquivos somente algumas datas são coincidentes. São essas que eu gostaria de emparelhadas. Alguém idéia como poderia fazer. Utilizando o código abaixo não tive sucesso: OS dois arquivos: Shanghai.xts e vale.xts ShVale=merge.xts(Shanghai.xts,vale.xts,all = FALSE,suffixes=NULL, join="outer",retside=TRUE,retclass="xts") ShVale = na.exclude(object=ShVale) o arquivo resultante (ShVale) tem não tem dados, embora existam datas coincidentes entre os dois arquivos originais. Edilson Floriano dos Santos

2 6

Re: [R-br] Obtendo efeitos INDIVIDUAIS de múltiplas variáveis qualitativas.
by Rafael Costa 08 Mai '15

08 Mai '15

Prezados colegas, Esclareço que antes de clicar em DOWNLOAD, por favor desmarquem a opção"Use our download manager and get recommended downloads" para baixar o arquivo diretamente sem problemas. Qualquer contribuição será muito bem vinda! Em 4 de maio de 2015 22:51, Rafael Costa <rafael_costa(a)caen.ufc.br> escreveu: > Mauro, obrigado pela atenção. Segue em anexo o arquivo. > > Em 4 de maio de 2015 21:18, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > >> Poderia me enviar este arquivo anexado, pois meu antivírus não deixa >> baixar! >> >> >> >> Prezados colegas, >> >> >> >> Possuo dados de um questionário e quero estimar o efeito individual de >> cada variável explicativa (todas são qualitativas) sobre a variável >> dependente (contínua). Contudo o padrão é considerar os coeficientes >> estimados como a diferença entre o o grupo de referência (valor estimado do >> intercepto) e o grupo do coeficiente. Cada variável qualitativa diz >> respeito a uma característica de uma determinada atividade e a variável >> contínua é o tempo tomado para executar tal atividade. Destaco que o nível >> de referência de cada fator diz respeito ao caso em que nenhuma das opções >> para aquele fator foi marcado. Os dados estão em " >> http://www.datafilehost.com/d/c7f0d342". Não os coloquei no CMR, pois >> ainda não consegui (peço minhas sinceras desculpas), mas espero que isto >> não seja empecilho.Também não coloquei apenas uma amostra dos dados, pois >> estava havendo problemas de matriz singular. >> >> >> >> 1º (e principal) Problema - Para obter o efeito individual dos níveis de >> cada fator, considerei que o grupo de referência possui efeito zero e fiz >> os seguintes passos: >> >> >> >> #Como o arquivo não está sendo carregado no CMR, supõe-se aqui que ele >> foi baixado da internet e já se encontra carregado no R. >> >> #Utilizarei uma regressão quantílica, por isso segue o pacote abaixo. >> >> install.packages(quantreg) >> >> library(quantreg) >> >> #Transformando as variáveis categóricas em objetos únicos >> >> p_1 = table(1:length(tabela1.1$p1),as.ordered(tabela1.1$p1)) >> >> p_21 = table(1:length(tabela1.1$p21),as.factor(tabela1.1$p21)) >> >> p_22 = table(1:length(tabela1.1$p22),as.factor(tabela1.1$p22)) >> >> p_23 = table(1:length(tabela1.1$p23),as.factor(tabela1.1$p23)) >> >> p_24 = table(1:length(tabela1.1$p24),as.factor(tabela1.1$p24)) >> >> p_25 = table(1:length(tabela1.1$p25),as.factor(tabela1.1$p25)) >> >> p_34 = table(1:length(tabela1.1$p34),as.ordered(tabela1.1$p34)) >> >> p_5 = table(1:length(tabela1.1$p5),as.ordered(tabela1.1$p5)) >> >> p_6 = table(1:length(tabela1.1$p6),as.ordered(tabela1.1$p6)) >> >> p_7 = table(1:length(tabela1.1$p7),as.ordered(tabela1.1$p7)) >> >> p_8 = table(1:length(tabela1.1$p8),as.ordered(tabela1.1$p8)) >> >> p_9 = table(1:length(tabela1.1$p9),as.ordered(tabela1.1$p9)) >> >> #Regressão do modelo sem intercepto, mas também considerando que o grupo >> de referência tem valor 0, haja vista que o grupo de referência significa >> que nenhum dos fatores foi marcado (se nenhum item foi marcado, considero >> que o tempo tomado para executar a atividade é praticamente zero). >> >> qrModel=rq(data=tabela1.1, pontoefetivo ~ 0 + p_1[,-1] + p_21[,-1] + >> p_22[,-1] + p_23[,-1] + p_24[,-1] + p_25[,-1] + p_34[,-1] + p_5[,-1] + >> p_6[,-1] + p_7[,-1] + p_8[,-1] + p_9[,-1], tau=0.5) >> >> summary(qrModel) >> >> >> >> A minha ideia foi que, dado que o grupo de referência é zero, o >> coeficiente estimado é justamente o efeito individual do nível escolhido da >> variável. A minha ideia está certa? Se não, o que faço para obter esses >> efeitos individuais? >> >> >> >> 2º Problema: Supondo que fiz tudo certinho acima, os fatores ordenados >> não possuem efeitos crescentes [Ver summary(qrModel)]. Mas não deveriam >> possuir? Se deveriam, o que faço para garantir tal efeito? >> >> >> >> 3º Problema: Mas uma vez supondo que tudo está correto, esperava que >> todos os efeitos sobre o tempo de execução das atividades fossem não >> negativos. Ou estou errado a respeito disso? Se estou certo, o que faço >> para garantir que todos os valores sejam não negativos. >> >> >> >> Desde já agradeço a todos e aguardo qualquer ajuda. >> <http://www.datafilehost.com/get.php?file=c7f0d342>-- >> >> >> >> >> *Rafael Carneiro da Costa Doutorando em Economia Centro de Pós-Graduação >> em Economia - CAEN Universidade Federal do Ceará - UFC* >> >> >> ------------------------------ >> [image: Avast logo] <http://www.avast.com/> >> >> Este email foi escaneado pelo Avast antivírus. >> www.avast.com >> >> > > > -- > > > > *Rafael Carneiro da CostaDoutorando em EconomiaCentro de Pós-Graduação em > Economia - CAENUniversidade Federal do Ceará - UFC* > -- *Rafael Carneiro da CostaDoutorando em EconomiaCentro de Pós-Graduação em Economia - CAENUniversidade Federal do Ceará - UFC*

1 0

Criar mapa com Image
by Wagner Wolff 08 Mai '15

08 Mai '15

Olá pessoal, quero gerar um mapa em que a variável em questão é obtida por meio de uma função, na qual as variáveis de entrada são dados espacializados em um grid. Entretanto estou me deparando com este erro: Error in image.default(x, y, z) : increasing 'x' and 'y' values expected Vocês teriam uma solução? Segue o código exemplo Qp <- function(mu,sigma,area,p){ Qp<-(exp(mu+(sigma*qnorm(1-p))))*area } grid <- pred_grid(range(Longitude),range(Latitude), by=0.5) gr <- locations.inside(grid, limite) x <- gr$Var1 y <- gr$Var2 z <- Qp(kc.mu.ANO$predict,kc.sigma.ANO$predict,1,0.9) image(x,y,z) contour(z, add=T, nlev=5)

2 3

R erro Rstudio
by Andre Oliveira 08 Mai '15

08 Mai '15

Pessoal,não uso Windows, mas hoje precisei de usar e encontrei este problema ao iniciar o Rstudio. Alguém poderia ajudar resolver este erro? É sobre permissão mas não sei resolver! R version 3.2.0 (2015-04-16) -- "Full of Ingredients" Copyright (C) 2015 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. Error installing package: Erro: ERROR: no permission to install to directory 'C:/Program Files/R/R-3.2.0/library' Error installing package: Erro: ERROR: no permission to install to directory 'C:/Program Files/R/R-3.2.0/library' 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 7

Banco de dados
by Andre Oliveira 07 Mai '15

07 Mai '15

Pessoal,gostaria de deixa aqui meus agradecimentos nas sugestões do grupo sobre banco de dados públicos e no problema do Rstudio não citarei ninguém para não cometer injustiça de esquecimento. No Rstudio de fato era permissão que o aluno da escola adicionou permissão daquela pasta do Rstudio. Só para constar voltarei para linux semana que vem hj foi algo pontual. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espírito Santo. IFES

1 0

Sequência e repetição de série temporal (minutos)
by Aline Lipsky 07 Mai '15

07 Mai '15

Boa tarde prezados colegas, Eu preciso criar um fator dena série temporal entre 0 a 10 minutos, com intervalo de 15 segundos e repetir essa sequência 21 vezes em ordem crescente. Para isso usei o seguinte script: library ( chron ) sort (rep (seq ( chron ( times = c("00:00:00") ), chron ( times = c("00:10:00") ), chron ( times = c("00:00:15"))), 21 )) Mas quando faço isso os valores que retornam não são o intervalo de tempo descrito.Consigo apenas criar a sequência sem repetição com o seguinte script: seq ( chron ( times = c ( "00:00:00" ) ), chron ( times = c ( "00:10:00" ) ), chron ( times = c ( "00:00:15" ))) Alguém saberia o que posso fazer? Devo usar outro pacote? Grata pela atenção. Atenciosamente, Aline Lipsky

2 3