Re: [R-br] Usando o R em grandes bases de dados

2012-11-20 Por tôpico Benilton Carvalho
bigmemory não suporta (ou não suportava, pelo menos) caracteres.
On 20 Nov 2012 05:07, Vinicius Brito Rocha viniciusbri...@gmail.com
wrote:


 fiquei brincando com o bigmemory e carreguei uma tabela mas os campos
 character vieram como NA.

 Quando usava a opção type=char piorava e alterava os valores dos campos
 numéricos. O melhor ( ou menos pior resultado) foi carregar a tabela com a
 opção type=double

 Alguem sabe o que pode ser?


 Abs.

 Vinicius Brito Rocha

 Em 19 de novembro de 2012 15:40, Richard Santos 
 jamesrichardsan...@gmail.com escreveu:

 Valeu!! Toda e qq ajuda nesse sentido é bem vinda!
 Abç.
 Richard

 Em 19/11/12, Daniel C Bezerradanielcbeze...@gmail.com escreveu:
  Eu encontrei estas apresentações e achei legal dividir com vcs:
 
  http://goo.gl/l6fuC
 
  Tem link para os slides, código do R e tudo mais por lá.
 
  Abs,
 
  Daniel
 
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 *Vinicius Brito Rocha.*
 *Estatístico e Atuário (IM / UFRJ)**
 Mestre em Pesquisa Operacional (COPPE / UFRJ)*

 www.aplicademic.blogspot.com
 http://twitter.com/viniciusbritor

 Não se preocupe muito com as suas dificuldades em Matemática, posso
 assegurar-lhe que as minhas são ainda maiores. - Albert Einstein.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Experimento fatorial. Desdobramento

2012-11-20 Por tôpico Walmes Zeviani
Esses dados e tantos outros (de consultorias com permissão dos
proprietários e de livros/materiais didáticos) estão disponíveis no
diretório

http://www.leg.ufpr.br/~walmes/data/

À disposição.
Walmes.

==
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: wal...@ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] excluir mar no ETOPO

2012-11-20 Por tôpico Eder Comunello
Carlos e colegas, bom dia!

Vou  aproveitar um script anterior, fazendo algumas alterações. Como da
outra vez, creio que dá pra adaptar facilmente pro seu caso.

Vou trabalhar extraindo o vetor de interesse, substituindo este vetor no
objeto origem ao final...

# begin #
require(sp); require(rgdal)

### Download do exemplo (399 Kb)
download.file('https://dl.dropbox.com/u/117618178/dourados/DDOS.tif',
  file.path(getwd(),'DDOS.tif'), mode='wb')

alt = readGDAL(DDOS.tif); gridded(alt) ### 'importa' o geotiff
image(alt, col=terrain.colors(21)) ## visualização
slotNames(alt) ### nome dos slots (S4)
proj4string(alt) ### projeçãoo definida no geotiff
head(alt@data) ### observe os valores dos dados
head(alt@data$band1) ### observe os valores dos dados

### Opção para fatiar (0 ou 1)

alt2 - alt ### para manipular uma cópia do objeto
z - alt2@data$band1 ### para manipular apenas o vetor de interesse
n - length(z); n

for(i in 1:n){if (z[i]400){z[i]=0} else{z[i]=1}}; head(z) ### 0 ou 1
alt2@data$band1 - z ### substitui o vetor manipulado no objeto
image(alt2) ## visualização

### Opção para mascarar valores abaixo de critério (com 0 ou NA)

alt3 - alt ### para manipular outra cópia do objeto
z - alt3@data$band1 ### para manipular apenas o vetor de interesse
n - length(z); n

for(i in 1:n){if (z[i]400){z[i]=0}}; head(z) ### 0 ou 1
#for(i in 1:n){if (z[i]400){z[i]=NA}}; head(z) ### 0 ou 1
alt3@data$band1 - z ### substitui o vetor manipulado no objeto
image(alt3) ## visualização
# end #



Éder Comunello

Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)
Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]

Researcher at Embrapa Western Region Agriculture
Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]

UTC-03:00
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] (sem assunto)

2012-11-20 Por tôpico Eder Comunello
Edson e colegas, bom dia!

Ao invés de deixar o campo 'motivo' em branco (NA), experimente atribuir um
valor, p. ex. 'sem resposta',  e rode novamente. Assim estes casos serão
considerados e aparecerão na saída. Vale pros outros campos.


Éder Comunello

Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)
Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]

Researcher at Embrapa Western Region Agriculture
Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]

UTC-03:00
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] excluir mar no ETOPO

2012-11-20 Por tôpico Thiago Veloso
Carlos,

Eu tentaria algo assim:

library (raster) #carrega pacote exigido
r - raster('ETOPO1_Ice_g_geotiff.tif') #carrega arquivo de trabalho
r [ r0 ] - NA #atribui NA a todos os valores negativos no raster.

E caso você queira atribuir NA a um escopo de valores, pode ser mais fácil 
criar uma função:


vals - c(0,24,32,56) # estes são os valores que você quer eliminar da tabela 
de dados do raster
fun - function (x) {
  x - ifelse (x %in% vals, 1, NA) # aqui você substitui os valores acima por 
1 e atribui NA para todo o resto. Você pode substituir o 1 por vals
}

E então aplicá-la sobre o raster:
r.new - calc (r, fun)


Uma vantagem de usar funções é pode realizar uma tarefa para diversos arquivos 
sem precisar repetir comandos para cada um deles.

Espero que ajude,
Thiago.



 From: D'Apolito carlosdapol...@gmail.com
To: r-br@listas.c3sl.ufpr.br 
Sent: Monday, November 19, 2012 8:15 PM
Subject: [R-br] excluir mar no ETOPO
 

Prezados,

Estou graficando mapas com os dados do ETOPO01 de elevação e gostaria de 
excluir o mar, para que os dados iniciem em zero ao invés de números negativos 
(profundidades dos oceanos).
Qualquer ajuda é bem vinda, desde já agradeço

Att
Carlos

#código:

library(raster) 
library(rgdal)

etopo-ETOPO1_Ice_g_geotiff.tif
r.etopo-raster(etopo)

ext - extent(-95,-35,-31,31)

zoom(r.etopo, ext=ext,col=rev(gray.colors(12)))





___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] excluir mar no ETOPO

2012-11-20 Por tôpico D'Apolito
Ola Eder,

O arquivo de origem dos dados (ETOPO1) `e muito grande, de maneira que
tenho obtido erro constantemente, Tentei aumentar a capacidade da memoria
com memory.limit() mas tambem nao funcionou.


att
Carlos


Em 20 de novembro de 2012 11:30, Eder Comunello ecomu...@gmail.comescreveu:

 Carlos e colegas, bom dia!

 Vou  aproveitar um script anterior, fazendo algumas alterações. Como da
 outra vez, creio que dá pra adaptar facilmente pro seu caso.

 Vou trabalhar extraindo o vetor de interesse, substituindo este vetor no
 objeto origem ao final...

 # begin #
 require(sp); require(rgdal)

 ### Download do exemplo (399 Kb)
 download.file('https://dl.dropbox.com/u/117618178/dourados/DDOS.tif',
   file.path(getwd(),'DDOS.tif'), mode='wb')

 alt = readGDAL(DDOS.tif); gridded(alt) ### 'importa' o geotiff
 image(alt, col=terrain.colors(21)) ## visualização
 slotNames(alt) ### nome dos slots (S4)
 proj4string(alt) ### projeçãoo definida no geotiff
 head(alt@data) ### observe os valores dos dados
 head(alt@data$band1) ### observe os valores dos dados

 ### Opção para fatiar (0 ou 1)

 alt2 - alt ### para manipular uma cópia do objeto
 z - alt2@data$band1 ### para manipular apenas o vetor de interesse
 n - length(z); n

 for(i in 1:n){if (z[i]400){z[i]=0} else{z[i]=1}}; head(z) ### 0 ou 1
 alt2@data$band1 - z ### substitui o vetor manipulado no objeto
 image(alt2) ## visualização

 ### Opção para mascarar valores abaixo de critério (com 0 ou NA)

 alt3 - alt ### para manipular outra cópia do objeto
 z - alt3@data$band1 ### para manipular apenas o vetor de interesse
 n - length(z); n

 for(i in 1:n){if (z[i]400){z[i]=0}}; head(z) ### 0 ou 1
 #for(i in 1:n){if (z[i]400){z[i]=NA}}; head(z) ### 0 ou 1
 alt3@data$band1 - z ### substitui o vetor manipulado no objeto
 image(alt3) ## visualização
 # end #


 
 Éder Comunello

 Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)
 Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]

 Researcher at Embrapa Western Region Agriculture
 Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
 
 UTC-03:00


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




-- 
_
*Carlos D'Apolito - PhD candidate - ID:1239480*
Earth Sciences
School of Geography  Environmental Sciences
Aston Webb Building
University of Birmingham
Birmingham
B15 2TT
UK

phone: 44 - 07 919 673 173
cxd...@bham.ac.uk
CV: http://lattes.cnpq.br/4557754746424026
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Graficos de decomposição

2012-11-20 Por tôpico kaue veras
Boa tarde prezados,

Estou realizando um trabalho em que utilizo a decomposição, análises de
previsão usando o modelo ARIMA e Holt-Winters e goataria de mudar os
títulos do gráfico e titulos dos eixos. Com os comandos main, xlab e ylab
não esta dando certo, poderiam me dar uma luz ou me ajudar como faço para
gerar essas modificações?

Desde já agradeço.

Atenciosamente,

Kauê P. Veras da Cunha
Estatística - UERJ
MSN: kaueve...@hotmail.com
E-mail: kaueve...@gmail.com
Cel: (21) 8254-9601
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Gráfico de distribuiçoes estimadas

2012-11-20 Por tôpico Matheus Monteiro
mais alguem pode me ajudar?

Em 10 de novembro de 2012 13:21, Edson Lira edinhoes...@yahoo.com.brescreveu:

 Mateus, no pacote MASS acho que tem alguns gráficos empírica x teórica.

 [. ]'s.
 Edson Lira
 Estatístico
 Ma-Am

 Em 10/11/2012, às 02:20, Matheus Monteiro matheusm...@gmail.com
 escreveu:

 Boa noite, obtive as densidades aproximadas da 't' e da 'normal', gostaria
 de compara-las através de graficos dessas densidades estimadas.
 Alguma sugestão, como faria isso?

  *t *

  fitdistr(x[,1],t)

 msdf

   -0.0018744140.0091227383.522064232

  ( 0.001726835) ( 0.001862692) ( 2.234067399)*

 **normal*

 fitdistr(x[,1],normal)

mean sd

  0.0001040816 0.0128201843

  (0.0018314549)  (0.0012950342)


 x são dados simulados.

 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] (sem assunto)

2012-11-20 Por tôpico Edson Lira
Mas Eder a idéia do NA é não considerar estas respostas na contagem e/ou no 
percentual. 

Mas valeu e muito obrigado. 

Peço desculpas a todos, esqueci de colocar o assunto na mensagem
 
Edson Lira
Estatístico
Manaus-Amazonas



 De: Eder Comunello ecomu...@gmail.com
Para: r-br@listas.c3sl.ufpr.br 
Enviadas: Terça-feira, 20 de Novembro de 2012 7:39
Assunto: Re: [R-br] (sem assunto)
 

Edson e colegas, bom dia!

Ao invés de deixar o campo 'motivo' em branco (NA), experimente atribuir um 
valor, p. ex. 'sem resposta',  e rode novamente. Assim estes casos serão 
considerados e aparecerão na saída. Vale pros outros campos.


Éder Comunello

Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)
Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]

Researcher at Embrapa Western Region Agriculture
Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]

UTC-03:00

___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Usando o R em grandes bases de dados

2012-11-20 Por tôpico Leonard de Assis
Bigmemory eu acredito que lide somente com matrizes, não com data 
frames, resumo: todos os campos devem ter o mesmo tipo.


[]s
Leonard de Assis
http://about.me/ldeassis

Em 20/11/2012 03:06, Vinicius Brito Rocha escreveu:


fiquei brincando com o bigmemory e carreguei uma tabela mas os campos 
character vieram como NA.


Quando usava a opção type=char piorava e alterava os valores dos 
campos numéricos. O melhor ( ou menos pior resultado) foi carregar a 
tabela com a opção type=double


Alguem sabe o que pode ser?


Abs.

Vinicius Brito Rocha

Em 19 de novembro de 2012 15:40, Richard Santos 
jamesrichardsan...@gmail.com mailto:jamesrichardsan...@gmail.com 
escreveu:


Valeu!! Toda e qq ajuda nesse sentido é bem vinda!
Abç.
Richard

Em 19/11/12, Daniel C Bezerradanielcbeze...@gmail.com
mailto:danielcbeze...@gmail.com escreveu:
 Eu encontrei estas apresentações e achei legal dividir com vcs:

 http://goo.gl/l6fuC

 Tem link para os slides, código do R e tudo mais por lá.

 Abs,

 Daniel

___
R-br mailing list
R-br@listas.c3sl.ufpr.br mailto:R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e
forneça código mínimo reproduzível.




--
/Vinicius Brito Rocha./
/Estatístico e Atuário (IM / UFRJ)//
Mestre em Pesquisa Operacional (COPPE / UFRJ)/

www.aplicademic.blogspot.com http://www.aplicademic.blogspot.com
http://twitter.com/viniciusbritor

Não se preocupe muito com as suas dificuldades em Matemática, posso 
assegurar-lhe que as minhas são ainda maiores. - Albert Einstein.





___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.


___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] matrizes de correlação

2012-11-20 Por tôpico alanarocha

Olá, a actualização do R resultou para o meu problema.
agora preciso de ajuda pois tou a tentar fazer uma matriz de  
correlação e outra de covariancias de 4 variaveis de um dataframe e  
dame estes erros:

##matriz de correlação

cor.p-cor.ps(dados[,2:4], Pearson=TRUE,NAs=TRUE)
cor.p
dados=na.omit(dados[,2:4])

cat(valores NA retirados\n)
 x11()
 pairs(dados)
 x11()
 boxplot(dados)
dados=dados[,2:4]
 x11()
 boxplot(dados)
correl=cor(dados[2:4])
result=correl
return(list(sumario=summary(dados),resultado=result))
cor(cbind(X1s09, X1s10, X1s11, X1s12))
cor.test(dados[,2:4], method='pearson')$p.value
cor_with_p_test(cbind(X1s09, X1s10, X1s11, X1s12))
ct - cor.test(X1s10, X1s12)
ct$p
ct$v

#matriz de covariancias
cov(dados[,2:4], na.rm=TRUE)



cor.p-cor.ps(dados[,2:4], Pearson=TRUE,NAs=TRUE)

Error: could not find function cor.ps

cor.p

Error: object 'cor.p' not found

dados=na.omit(dados[,2:4])

cat(valores NA retirados\n)

valores NA retirados





correl=cor(dados[2:4])

Error in `[.data.frame`(dados, 2:4) : undefined columns selected

result=correl

Error: object 'correl' not found

return(list(sumario=summary(dados),resultado=result))

Error: object 'result' not found

cor(cbind(X1s09, X1s10, X1s11, X1s12))

Error in cbind(X1s09, X1s10, X1s11, X1s12) : object 'X1s09' not found

cor.test(dados[,2:4], method='pearson')$p.value

Error in `[.data.frame`(dados, , 2:4) : undefined columns selected

cor_with_p_test(cbind(X1s09, X1s10, X1s11, X1s12))

Error: could not find function cor_with_p_test

ct - cor.test(X1s10, X1s12)

Error in cor.test(X1s10, X1s12) : object 'X1s10' not found

ct$p

Error: object 'ct' not found

ct$v

Error: object 'ct' not found


#matriz de covariancias
cov(dados[,2:4], na.rm=TRUE)

Error in cov(dados[, 2:4], na.rm = TRUE) :
  unused argument(s) (na.rm = TRUE)



cumprimentos e obrigada Ana Rocha



___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.


[R-br] ExpDes - Comparações múltiplas

2012-11-20 Por tôpico Andre Oliveira
Boa tarde, 
estou aqui tentando reproduzir algumas análises do ExpDes. A analise abaixo não 
me deu as comparações múltiplas como faço para tais comparações?

require(ExpDes)
Carregando pacotes exigidos: ExpDes
data(ex1)
attach(ex1)
dic(trat, ig, quali = FALSE, sigF = 0.05)


 
André Oliveira Souza___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dúvida sobre IC simultâneos para modelagem de variáveis binarias no R (multcomp)

2012-11-20 Por tôpico Walmes Zeviani
Seu exemplo não é reproduzível. O guia de postagem não recomenda envio de
dados como anexo das mensagens. Hospede os dados na núvem e acrescente
comandos para leitura do arquivo. Assim rapidamente podemos reproduzir sua
análise e ajudar você.

À disposição.
Walmes.

==
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: wal...@ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux Ubuntu 10.04

2012-11-20 Por tôpico Andre Oliveira
Versão ExpDes_1.1.1.tar.gz agora sim consegui instalar no linux.


 
André Oliveira Souza



 De: Andre Oliveira andreolso...@yahoo.com.br
Para: Eric Ferreira ericbferre...@gmail.com 
Enviadas: Quinta-feira, 5 de Janeiro de 2012 11:05
Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux Ubuntu 
10.04
 

O  Walmes me disse, aqui mesmo na lista,  ha algum tempo que realmente o pacote 
para linux em português ele não consegui colocar para funcionar. E no caso do 
windows eu não sabia que estava em inglês! Obrigado 


 
André Oliveira Souza



 De: Eric Ferreira ericbferre...@gmail.com
Para: Andre Oliveira andreolso...@yahoo.com.br 
Enviadas: Quinta-feira, 5 de Janeiro de 2012 10:55
Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux Ubuntu 
10.04
 

Olá Andre, 

Que bom que não desistiu. Vamos por partes que solucionaremos todos os seus 
problemas.
Quanto ao Windows,
 não se esqueça que a versão do CRAN está em ingles. Portanto, após carregar, 
voce deve pedir o help da função cdr (completely randomized design) e não dic, 
ou seja
?crd
Veja se funciona.


Quanto ao linux, talvez voce deva escrever para o Walmes, que faz parte dessa 
lista. Ele usa o ExpDes no Linux e diz que funciona perfeitamente.

abs

Eric.


Em 5 de janeiro de 2012 10:42, Andre Oliveira andreolso...@yahoo.com.br 
escreveu:

Eric,

eu acompanhei o pacote no seu site e nunca consegui coloca-lo para funcionar 
no Linux a versão em português. Agora com ele no CRAN  voltei a tentar colocar 
o pacote para funcionar no Linux,  mas sem sucesso. E para minha surpresa no 
windows não está funcionado também. Olha a minha tentativa no windows!
require(ExpDes)
Carregando pacotes exigidos: ExpDes
 ?dic
No documentation for ‘dic’ in specified packages and libraries:
you could try ‘??dic’


 Estou ansioso para colocar este pacote para funcionar por aqui. Já li o TCC e 
o manual sobre o mesmo e vi a importância do mesmo. Visto que antes tinha que 
usar muitos pacotes para a mesma análise e com o ExpDes é muito mais prático. 










André Oliveira Souza







 De: Eric Ferreira ericbferre...@gmail.com
Para: r-br@listas.c3sl.ufpr.br; Andre Oliveira andreolso...@yahoo.com.br 
Enviadas: Quarta-feira, 4 de Janeiro de 2012 18:25
Assunto: Re: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu 
10.04
 

Olá André,

Se quiser trocar uma ideia em particular, posso tentar te ajudar. Sou o 
criador do ExpDes, mas não usu Linux. Não sei se serei útil. 

A princípio, posso te dizer que existe o pacote no CRAN e na minha página 
pessoal: https://sites.google.com/site/ericbferreira/unifal/downloads-1

Tente pegar o pacote de lá para ver se o erro se repete.
Peço gentilmente aos usuárioas de (R no) Linux que me ajudem a exclarecer o 
André.

Atenciosamente,

Eric.


Em 4 de janeiro de 2012 16:39, Andre Oliveira andreolso...@yahoo.com.br 
escreveu:

Opa vacilo meu erro no nome do pacote 


 install.packages(ExpDes, dep=T) 






André Oliveira Souza



 De: Andre Oliveira andreolso...@yahoo.com.br
Para: r-br@listas.c3sl.ufpr.br r-br@listas.c3sl.ufpr.br 
Enviadas: Quarta-feira, 4 de Janeiro de 2012 16:36
Assunto: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu 10.04



Pessoal estou com problemas para reinstalar o EspDes !!  Alguém saberia dizer 
o que pode ser?




install.packages(EspDes, dep=T) 
Installing package(s) into ‘/home/andre/R/i486-pc-linux-gnu-library/2.14’
(as ‘lib’ is unspecified)
--- Please select a CRAN mirror for use in this session ---
Carregando a interface Tcl/Tk ... pronto
Mensagens de aviso perdidas:
In getDependencies(pkgs, dependencies, available, lib) :
  package ‘EspDes’ is not available (for R version 2.14.1)






 
André Oliveira Souza
obrigado 


___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.


___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.



-- 
Dr Eric B Ferreira
Exact Sciences Department
Federal University of Alfenas
Brazil





-- 
Dr Eric B Ferreira
Exact Sciences Department
Federal University of Alfenas
Brazil___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux Ubuntu 10.04

2012-11-20 Por tôpico Andre Oliveira
Em português no linux!
 
André Oliveira Souza



 De: Andre Oliveira andreolso...@yahoo.com.br
Para: Eric Ferreira ericbferre...@gmail.com; r-br@listas.c3sl.ufpr.br 
r-br@listas.c3sl.ufpr.br 
Enviadas: Terça-feira, 20 de Novembro de 2012 16:30
Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux Ubuntu 
10.04
 

Versão ExpDes_1.1.1.tar.gz agora sim consegui instalar no linux.


 
André Oliveira Souza



 De: Andre Oliveira andreolso...@yahoo.com.br
Para: Eric Ferreira ericbferre...@gmail.com 
Enviadas: Quinta-feira, 5 de Janeiro de 2012 11:05
Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux Ubuntu 
10.04
 

O  Walmes me disse, aqui mesmo na lista,  ha algum tempo que realmente o pacote 
para linux em português ele não consegui colocar para funcionar. E no caso do 
windows eu não sabia que estava em inglês! Obrigado 


 
André Oliveira Souza



 De: Eric Ferreira ericbferre...@gmail.com
Para: Andre Oliveira andreolso...@yahoo.com.br 
Enviadas: Quinta-feira, 5 de Janeiro de 2012 10:55
Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux Ubuntu 
10.04
 

Olá Andre, 

Que bom que não desistiu. Vamos por partes que solucionaremos todos os seus 
problemas.
Quanto ao Windows,
 não se esqueça que a versão do CRAN está em ingles. Portanto, após carregar, 
voce deve pedir o help da função cdr (completely randomized design) e não dic, 
ou seja
?crd
Veja se funciona.


Quanto ao linux, talvez voce deva escrever para o Walmes, que faz parte dessa 
lista. Ele usa o ExpDes no Linux e diz que funciona perfeitamente.

abs

Eric.


Em 5 de janeiro de 2012 10:42, Andre Oliveira andreolso...@yahoo.com.br 
escreveu:

Eric,

eu acompanhei o pacote no seu site e nunca consegui coloca-lo para funcionar 
no Linux a versão em português. Agora com ele no CRAN  voltei a tentar colocar 
o pacote para funcionar no Linux,  mas sem sucesso. E para minha surpresa no 
windows não está funcionado também. Olha a minha tentativa no windows!
require(ExpDes)
Carregando pacotes exigidos: ExpDes
 ?dic
No documentation for ‘dic’ in specified packages and libraries:
you could try ‘??dic’


 Estou ansioso para colocar este pacote para funcionar por aqui. Já li o TCC e 
o manual sobre o mesmo e vi a importância do mesmo. Visto que antes tinha que 
usar muitos pacotes para a mesma análise e com o ExpDes é muito mais prático. 










André Oliveira Souza







 De: Eric Ferreira ericbferre...@gmail.com
Para: r-br@listas.c3sl.ufpr.br; Andre Oliveira andreolso...@yahoo.com.br 
Enviadas: Quarta-feira, 4 de Janeiro de 2012 18:25
Assunto: Re: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu 
10.04
 

Olá André,

Se quiser trocar uma ideia em particular, posso tentar te ajudar. Sou o 
criador do ExpDes, mas não usu Linux. Não sei se serei útil. 

A princípio, posso te dizer que existe o pacote no CRAN e na minha página 
pessoal: https://sites.google.com/site/ericbferreira/unifal/downloads-1

Tente pegar o pacote de lá para ver se o erro se repete.
Peço gentilmente aos usuárioas de (R no) Linux que me ajudem a exclarecer o 
André.

Atenciosamente,

Eric.


Em 4 de janeiro de 2012 16:39, Andre Oliveira andreolso...@yahoo.com.br 
escreveu:

Opa vacilo meu erro no nome do pacote 


 install.packages(ExpDes, dep=T) 






André Oliveira Souza



 De: Andre Oliveira andreolso...@yahoo.com.br
Para: r-br@listas.c3sl.ufpr.br r-br@listas.c3sl.ufpr.br 
Enviadas: Quarta-feira, 4 de Janeiro de 2012 16:36
Assunto: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu 10.04



Pessoal estou com problemas para reinstalar o EspDes !!  Alguém saberia dizer 
o que pode ser?




install.packages(EspDes, dep=T) 
Installing package(s) into ‘/home/andre/R/i486-pc-linux-gnu-library/2.14’
(as ‘lib’ is unspecified)
--- Please select a CRAN mirror for use in this session ---
Carregando a interface Tcl/Tk ... pronto
Mensagens de aviso perdidas:
In getDependencies(pkgs, dependencies, available, lib) :
  package ‘EspDes’ is not available (for R version 2.14.1)






 
André Oliveira Souza
obrigado 


___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.


___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.



-- 
Dr Eric B Ferreira
Exact Sciences Department
Federal University of Alfenas
Brazil





-- 
Dr Eric B Ferreira
Exact Sciences Department
Federal University of Alfenas
Brazil





___
R-br mailing list

Re: [R-br] Dúvida em summary em regressão logística

2012-11-20 Por tôpico Alexandre Santos
Obrigado Benilton,

      Mas o modelo parece estar correto, estou seguindo a metodologia empregada 
por Farhad et al. 2011 pag 2., Equação 1. (Foraging behavior of Praon Volucre 
 doi: 10.1155/2011/868546).

      Achei um erro no meu CMR que ficaria:

 #
 ###
 #Regressão logística entre a proporção de herbivoros predados e a densidade 
 de herbivoros oferecidos
 #
 pred-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2
+ ,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26
+ ,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros 
predados
 dens-sort(rep(2^(2:6),10)) Densidade de herbivoros oferecidos
 ##
 ##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ exp(P0+P1*N+P2*N^2+P3*N3)
 p.model1-glm(pred/dens~dens+I(dens^2)+I(dens^3),family=binomial)
Mensagens de aviso perdidas:
In eval(expr, envir, enclos) : #sucessos não-inteiro em um glm binomial!
 summary(p.model1)

Call:
glm(formula = pred/dens ~ dens + I(dens^2) + I(dens^3), family = binomial)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
-0.85240  -0.17292  -0.05812   0.19453   1.10008  

Coefficients:
              Estimate Std. Error z value Pr(|z|)
(Intercept) -5.744e-01  1.231e+00  -0.467    0.641
dens         2.238e-01  2.244e-01   0.997    0.319
I(dens^2)   -8.846e-03  9.018e-03  -0.981    0.327
I(dens^3)    8.671e-05  9.154e-05   0.947    0.344

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 8.1347  on 49  degrees of freedom
Residual deviance: 6.8751  on 46  degrees of freedom
AIC: 67.087

Number of Fisher Scoring iterations: 4

 #

Mas é exatamente o seu exemplo que eu procurava, pois eu queria os coeficientes 
de P0, P1, P2 e P3, seguindo sua ajuda:

 naoPred-dens-pred
 p.model2-glm(cbind(pred, naoPred)~poly(dens, 3, raw=TRUE), family='binomial')
 summary(p.model2)

Call:
glm(formula = cbind(pred, naoPred) ~ poly(dens, 3, raw = TRUE), 
    family = binomial)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-3.2393  -0.9997  -0.0732   0.9515   4.2567  

Coefficients:
                             Estimate Std. Error z value Pr(|z|)    
(Intercept)                -8.399e-01  4.949e-01  -1.697 0.089672 .  
poly(dens, 3, raw = TRUE)1  2.688e-01  7.444e-02   3.611 0.000305 ***
poly(dens, 3, raw = TRUE)2 -1.055e-02  2.755e-03  -3.828 0.000129 ***
poly(dens, 3, raw = TRUE)3  1.033e-04  2.681e-05   3.853 0.000117 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 148.29  on 49  degrees of freedom
Residual deviance: 121.77  on 46  degrees of freedom
AIC: 280.61

Number of Fisher Scoring iterations: 4

 #


Obrigado pelas dicas e correções,


Alexandre




 De: Benilton Carvalho beniltoncarva...@gmail.com
Para: r-br r-br@listas.c3sl.ufpr.br; Alexandre Santos 
alexandresanto...@yahoo.com.br 
Enviadas: Terça-feira, 20 de Novembro de 2012 16:30
Assunto: Re: [R-br] Dúvida em summary em regressão logística
 

Seu exemplo nao e' reproduzivel e seu modelo nao esta' correto.

Consulte seu estatistico local para esclarecimentos mais detalhados.

Uma regressao logistica modela a probabilidade de sucesso dado um conjunto de 
covariaveis. No seu caso, sucesso (para o predador) parece ser o herbivoro 
ser predado.

Dito isso, se naoPred fosse o numero de herbivoros que nao foram predados, a 
especificacao do seu modelo seria

glm(cbind(pred, naoPred)~poly(dens, 3, raw=TRUE), family='binomial')

Se vc usar a representacao na escala probalistica (note que faltam uns 
parenteses na sua representacao):

Prob(Sucesso) = exp(P0+P1*N+P2*N^2+P3*N^3)/(1+exp(P0+P1*N+P2*N^2+P3*N^3))

o resultado do summary (os coeficiences) mostrado(s) representa(m) 
respectivamente P0, P1, P2 e P3.

b



2012/11/20 Alexandre Santos alexandresanto...@yahoo.com.br

Boa tarde Pessoal,
       Estou ajustando uma regressão logística e me deparei com a seguinte 
dúvida:
      
###
#Regressão logística entre a proporção de herbívoros predados e a densidade 
de herbívoros oferecidos
#
pred-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2
,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26
,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros 
predados
dens-sort(rep(2^(2:6),10)) Densidade de herbivoros oferecidos
##
##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ exp(P0+P1*N+P2*N^2+P3*N3)
p.model1-glm(pred/par~dens+I(dens^2)+I(dens^3),family=binomial)
summary(p.model1)
#




Minha dúvida é se o coeficiente linear que aparece no summary esta 
transformado em exp(x)/1+exp(x)
ou trata-se do valor sem transformação?


Obrigado,

-- 
==
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e 

Re: [R-br] Usando o R em grandes bases de dados

2012-11-20 Por tôpico Vinicius Brito Rocha
Benilton,

acho que  a nova versão suporta.

Tem opção CHAR, mas tem que ser todo o objeto(tabela) de char.

Att.

Em 20 de novembro de 2012 07:29, Benilton Carvalho 
beniltoncarva...@gmail.com escreveu:

 bigmemory não suporta (ou não suportava, pelo menos) caracteres.
 On 20 Nov 2012 05:07, Vinicius Brito Rocha viniciusbri...@gmail.com
 wrote:


 fiquei brincando com o bigmemory e carreguei uma tabela mas os campos
 character vieram como NA.

 Quando usava a opção type=char piorava e alterava os valores dos campos
 numéricos. O melhor ( ou menos pior resultado) foi carregar a tabela com a
 opção type=double

 Alguem sabe o que pode ser?


 Abs.

 Vinicius Brito Rocha

 Em 19 de novembro de 2012 15:40, Richard Santos 
 jamesrichardsan...@gmail.com escreveu:

 Valeu!! Toda e qq ajuda nesse sentido é bem vinda!
 Abç.
 Richard

 Em 19/11/12, Daniel C Bezerradanielcbeze...@gmail.com escreveu:
  Eu encontrei estas apresentações e achei legal dividir com vcs:
 
  http://goo.gl/l6fuC
 
  Tem link para os slides, código do R e tudo mais por lá.
 
  Abs,
 
  Daniel
 
 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




 --
 *Vinicius Brito Rocha.*
 *Estatístico e Atuário (IM / UFRJ)**
 Mestre em Pesquisa Operacional (COPPE / UFRJ)*

 www.aplicademic.blogspot.com
 http://twitter.com/viniciusbritor

 Não se preocupe muito com as suas dificuldades em Matemática, posso
 assegurar-lhe que as minhas são ainda maiores. - Albert Einstein.



 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.


 ___
 R-br mailing list
 R-br@listas.c3sl.ufpr.br
 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
 Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
 código mínimo reproduzível.




-- 
*Vinicius Brito Rocha.*
*Estatístico e Atuário (IM / UFRJ)**
Mestre em Pesquisa Operacional (COPPE / UFRJ)*

www.aplicademic.blogspot.com
http://twitter.com/viniciusbritor

Não se preocupe muito com as suas dificuldades em Matemática, posso
assegurar-lhe que as minhas são ainda maiores. - Albert Einstein.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dúvida em summary em regressão logística

2012-11-20 Por tôpico Alexandre Santos
Obrigado Vinícius,

       Mas o Benilton já havia chamado a atenção para uma questão crucial: 
regressão logística modela a probabilidade de sucesso dado um conjunto de 
covariaveis, então no modelo original que estava trabalhando existem erros 
conceituais, uma aproximação do que eu queria que era a obtenção dos 
coeficientes do intercepto, linear, quadrático e cúbico, ficou:

 pred-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2
,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26
,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de mumias
dens-sort(rep(2^(2:6),10)) Densidade de pulgoes oferecidos
##
##Prob(Sucesso) = exp(P0+P1*N+P2*N^2+P3*N^3)/(1+exp(P0+P1*N+P2*N^2+P3*N^3))
naoPred-dens-pred
#
p.model2-glm(cbind(pred, naoPred)~poly(dens, 3, raw=TRUE), family=binomial)
summary(p.model2)
#

Abs

Alexandre


 De: Vinicius Brito Rocha viniciusbri...@gmail.com
Para: r-br@listas.c3sl.ufpr.br; Alexandre Santos 
alexandresanto...@yahoo.com.br 
Enviadas: Terça-feira, 20 de Novembro de 2012 21:32
Assunto: Re: [R-br] Dúvida em summary em regressão logística
 

Alexandre,

acho que vc chamou o modelo de forma errada.

deveria ser 
glm(cbind(evento,1-evento)~VA1+VA2+...VAp,data=dados,family=binomial(link = 
logit))

normalmente nos softwares o coeficiente resposta é o do modelo linear, como vc 
esta usando a funçao logit

log(P(X=evento)/P(X=nao evento))= beta0+beta1*X1+beta2*X2+...+betap*Xp


Abs.

Mas o Walmes provavelmente vai poder te responder melhor.

Abs

Vinicius Brito Rocha


Em 20 de novembro de 2012 16:49, Alexandre Santos 
alexandresanto...@yahoo.com.br escreveu:

Boa tarde Pessoal,
       Estou ajustando uma regressão logística e me deparei com a seguinte 
dúvida:
      
###
#Regressão logística entre a proporção de herbívoros predados e a densidade 
de herbívoros oferecidos
#
pred-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2
,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26
,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros 
predados
dens-sort(rep(2^(2:6),10)) Densidade de herbivoros oferecidos
##
##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ exp(P0+P1*N+P2*N^2+P3*N3)
p.model1-glm(pred/par~dens+I(dens^2)+I(dens^3),family=binomial)
summary(p.model1)
#




Minha dúvida é se o coeficiente linear que aparece no summary esta 
transformado em exp(x)/1+exp(x)
ou trata-se do valor sem transformação?


Obrigado,

-- 
==
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso - 
Campus Cáceres
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:alexandresanto...@yahoo.com.br alexandre.san...@cas.ifmt.edu.br 
==


___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.



-- 
Vinicius Brito Rocha.
Estatístico e Atuário (IM / UFRJ)
Mestre em Pesquisa Operacional (COPPE / UFRJ)

www.aplicademic.blogspot.com
http://twitter.com/viniciusbritor

Não se preocupe muito com as suas dificuldades em Matemática, posso 
assegurar-lhe que as minhas são ainda maiores. - Albert Einstein.___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Preencher um serie temporal

2012-11-20 Por tôpico Daniel Marcelino
Caros,
Eu gostaria de expandir os dados abaixo para realizar uma análise
temporal de base diária.
O meu período inicial é '2012-07-01' até '2012-10-06'. É possível eu
fazer isso com meu data frame ou eu terei que criar um novo data frame
com a sequencia de dias e juntar as duas tabelas?
Meu data frame final deveria repetir os valores  apresentados nas
colunas até que uma nova medida surgisse já que os pontos de medida
que tenho são exporádicos.  Alguém já expandiu um data frame assim pra
me dizer como fazer?


structure(list(Date.beg = structure(c(15635, 15617, 15615, 15610,
15609, 15605, 15604, 15601, 15593, 15593, 15586, 15581, 15580,
15577, 15572, 15565, 15552, 15540, 15530, 15516), class = Date),
Date.end = structure(c(15619, 15619, 15616, 15615, 15610,
15607, 15604, 15602, 15595, 15594, 15587, 15582, 15581, 15579,
15572, 15567, 15554, 15541, 15533, 15517), class = Date),
Pollster = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 4L, 1L, 2L,
1L, 1L, 2L, 1L, 4L, 1L, 2L, 2L, 1L, 3L, 1L), .Label = c(Datafolha,
Ibope, Veritá, Vox Populi), class = factor), Serra.PSDB = c(24,
22, 23, 19, 22, 17, 17, 21, 19, 20, 21, 20, 22, 22, 27, 26,
26, 30, 31.4, 31), Russomanno.PRB = c(23, 22, 25, 27, 30,
34, 34, 35, 35, 32, 35, 31, 31, 31, 31, 26, 25, 26, 17.7,
24), Haddad.PT = c(20, 22, 19, 18, 18, 18, 17, 15, 15, 17,
16, 16, 14, 14, 8, 9, 6, 7, 9.5, 6), Chalita.PMDB = c(11,
11, 11, 10, 9, 7, 5, 8, 6, 8, 7, 5, 7, 5, 6, 5, 5, 6, 4.3,
6), Others = c(8, 7, 7, 4, 7, 6, 4, 5, 6, 7, 9, 5, 8, 6,
10, 12, 15, 15, 8.6, 16), Null = c(8, 8, 8, NA, 8, 10, 10,
NA, 13, 9, 8, 12, 10, NA, 10, 12, 14, 11, 13.1, 11), Undecided = c(6,
8, 6, NA, 6, 8, 13, NA, 6, 7, 4, 9, 7, 13, 6, 10, 9, 6, 13.3,
5), Round = structure(c(3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L), .Label = c(year before,
off campaign, first, second), class = factor), Stage =
structure(c(2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L,
1L, 1L, 1L, 1L), .Label = c(0, 1, 2), class = factor),
Serra.Haddad = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c(,
1, 0), class = factor), N = c(3959, 1204, 2099, 1204,
1799, 1204, 2000, 1802, 1001, 1221, 1078, 1001, 1069, 1200,
1077, 805, 805, 1075, 1331, 1081), Err = c(2, 3, 2, 3, 2,
3, 2.2, 2, 3, 3, 3, 3, 3, 2.8, 3, 3, 3, 3, 2.7, 3)), .Names =
c(Date.beg,
Date.end, Pollster, Serra.PSDB, Russomanno.PRB, Haddad.PT,
Chalita.PMDB, Others, Null, Undecided, Round, Stage,
Serra.Haddad, N, Err), row.names = c(6, 7, 8, 9,
10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20,
21, 22, 23, 24, 25), class = data.frame)



--
Small steps toward a much better world


-- 
--
View my research on these pages:
http://www.crcee.umontreal.ca
http://danielmarcelino.com
http://ssrn.com/author=1640522
--
Small steps toward a much better world

\begin{signature}
Daniel Marcelino
Land Phone 1+514 343 6111 #3799
3200 Jean Brillant,  Office C5071
Montreal, QC; H3T 1N8
Canada
\end{signature}
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.