- R-br - DINF UFPR

[Dúvida] Maximização e Minimização de funções.
by Pedro Rafael 14 Nov '12

14 Nov '12

Pessoal qual a diferença entre a função nlminb para a optim? O método de Nelder-Mead pelo que percebo é um pouco mais lento que o método BFGS contudo tem menos problemas com os chutes iniciais. Em geral, quando quero maximizar uma função eu sempre faço uma inspeção gráfica para escolher meus chutes iniciais. Existe alguma rotina de métodos fazem essas escolhas dos chutes de forma automática como nos algoritmos evolutivos? Obrigado à todos. [ ], Pedro Rafael Diniz Marinho.

1 0

(sem assunto)
by Carlos pombo sonderblohm 13 Nov '12

13 Nov '12

http://www.topcouponsdeals.com/wp-content/themes/twentyten/ugoogle.html

1 0

Re: [R-br] Fwd: como extrair dados de arquivos GEOtiff
by Thiago Veloso 13 Nov '12

13 Nov '12

Olá Carlos, Estou respondendo também para a lista para com vistas a guardar a discussão que no futuro pode ajudar alguém. Fica difícil analisar o problema sem ter uma amostra real dos dados, mas pelo que você relata eu tentaria mudar a projeção da imagem. A mais comum, e que eu tentaria logo de cara, é a lat/lon WGS84. Você pode usar a função "projection" do próprio pacote raster, em um comando parecido com: projection(r) <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 " Saudações, Thiago. ________________________________ From: D'Apolito <carlosdapolito(a)gmail.com> To: Thiago Veloso <thi_veloso(a)yahoo.com.br> Sent: Tuesday, November 13, 2012 5:17 PM Subject: Re: [R-br] Fwd: como extrair dados de arquivos GEOtiff Olá Thiago, Realmente, dividindo-se tudo por 10 os valores condizem perfeitamente com o esperado. Eles devem ter subido os dados assim para pesar menos, ou sei lá por que. Agora meu problema é o seguinte, as coordenadas geográficas aparecem de 0-3600 pontos correspondendo a 30 graus = 120/grau. Mas quando uso plot(r) a escala de Lat e Lon aparecem normais, haverá uma maneira fácil de transformar esses dados, do contrário as tabelas em excel ficam demasiado grandes. Muito obrigado Carlos Em 13 de novembro de 2012 14:05, Thiago Veloso <thi_veloso(a)yahoo.com.br> escreveu: Carlos, > > > O procedimento parece estar correto. Se você está recebendo valores que não esperava, pode ser algum problema com o arquivo tif de entrada. Você tem certeza de que ele já expressa temperatura? Dependendo da fonte desse arquivo, você pode ter que fazer uma calibração dos números digitais para obter valores de temperatura fisicamente consistentes. > > > Como se trata de uma imagem, a primeira coisa que eu faria seria plotar o objeto raster (imagem) para fazer uma análise exploratória: > > > plot(r) > > > Thiago. > > > >________________________________ > From: D'Apolito <carlosdapolito(a)gmail.com> >To: r-br(a)listas.c3sl.ufpr.br >Sent: Tuesday, November 13, 2012 12:21 PM >Subject: [R-br] Fwd: como extrair dados de arquivos GEOtiff > > > >Prezados, > > > >Tenho arquivos no formato GEOtiff (ex. anexo) e quero extrair dados de temperatura média para cada coordenada geográfica. > > >Estou usando isso até o momento: > > >library(raster) >f <- "tmean1_23.tif" >r <- raster(f) >v = getValues(r, format='matrix') > > > >Mas o que obtenho com getValues() é apenas uma matriz enorme com números que não tem nada a ver com temperatura. >Meu nível de R não é muito avançado, por isso não estou conseguindo sair desse passo, agradeço qualquer ajuda. > > >Att >Carlos > > > > > > > > > >_______________________________________________ >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. > > -- _____________________________________ 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 CXD280(a)bham.ac.uk CV: http://lattes.cnpq.br/4557754746424026

1 0

codigo para graficos de dispersão
by alanarocha＠sapo.pt 13 Nov '12

13 Nov '12

caro luciano e todos eu tentei fazer o seguinte codigo: e mostro o erro que me dá: setwd("G:/SIN/Users/Ana Rua/Ana_trabalho") #objectivo: #avaliar a correlação de pthi com DOSE_Cinacalcet, controlando os efeitos da variável GRUPO, que é do tipo fator com dois estados, H e C. dados<-read.table("correlacaocinacalcet.csv",header=TRUE,sep=";") dados Doente Grupo Pthi Dose_Cinacalcet 1 54020637 H 221.60 1080 2 54020698 H 317.00 2790 3 54020681 H 412.60 930 4 54131227 H 462.00 900 5 54030338 H 474.90 930 6 54220370 H 547.00 930 7 54110496 H 547.50 930 8 54110160 H 589.00 930 9 54140258 H 709.00 930 10 54110468 H 804.90 930 11 54130380 H 842.00 930 12 54300217 H 855.40 930 13 54110542 H 911.50 930 14 54300021 H 979.30 930 15 54020732 H 1000.00 930 16 54300160 H 1054.00 840 17 54300040 H 1178.00 900 18 54131239 H 1409.00 930 19 54140376 H 1579.00 420 20 54020303 H 2794.00 1860 21 54030135 H 204.00 420 22 54140295 H 236.00 2790 23 54140011 H 253.00 930 24 54210252 H 334.00 930 25 54130307 H 358.00 390 26 54300160 H 360.80 930 27 54210085 H 426.00 930 28 54030451 H 494.50 1860 29 54020681 H 535.90 900 30 54030344 H 594.60 1860 31 54110542 H 615.70 900 32 54110255 H 659.00 1860 33 54110468 H 666.10 900 34 54020637 H 705.70 480 35 54300040 H 765.20 930 36 54300217 H 875.80 930 37 54110496 H 953.20 900 38 54131239 H 982.00 930 39 54300021 H 1275.00 930 40 54110160 H 1297.00 900 41 54110070 H 1312.00 1440 42 54130380 H 1330.00 930 43 54020732 H 1388.00 900 44 54020303 H 2060.00 2790 45 54110386 H 2432.00 900 46 54210217 H 63.50 1860 47 54110255 H 86.19 1710 48 54210104 H 236.00 2040 49 54210333 H 249.00 1860 50 54300040 H 256.00 390 51 54030135 H 265.90 60 52 54140295 H 322.70 2790 53 54030451 H 354.60 1860 54 54110542 H 436.00 930 55 54060515 H 458.80 510 56 54110468 H 460.40 930 57 54210252 H 480.10 930 58 54300160 H 480.90 900 59 54210085 H 627.80 930 60 54110496 H 700.50 930 61 54110160 H 700.90 930 62 54020732 H 739.90 870 63 54030344 H 760.70 1860 64 54020681 H 808.60 870 65 54300021 H 830.10 900 66 54210233 H 883.00 930 67 54140326 H 1033.00 390 68 54300217 H 1074.00 900 69 54210181 H 1138.00 2820 70 54130380 H 1145.00 900 71 54110070 H 1174.00 1350 72 54020637 H 1308.00 870 73 54220660 H 1401.00 1860 74 54300248 H 1462.00 1860 75 54020303 H 1536.00 2700 76 54131031 C 76.50 780 77 54210217 C 159.30 390 78 54210252 C 161.10 780 79 54210181 C 268.00 2070 80 54300021 C 280.90 390 81 54030451 C 364.10 720 82 54300160 C 574.30 390 83 54110468 C 603.90 390 84 54030135 C 651.80 390 85 54300255 C 727.10 360 86 54110542 C 738.60 390 87 54110255 C 745.00 390 88 54210233 C 764.50 1170 89 54020681 C 771.20 420 90 54300217 C 772.30 840 91 54140011 C 788.70 540 92 54030344 C 807.50 780 93 54210085 C 813.60 900 94 54110160 C 856.30 780 95 54300268 C 865.10 840 96 54110070 C 882.70 810 97 54140326 C 883.90 540 98 54140295 C 916.30 600 99 54110496 C 923.70 780 100 54210333 C 984.70 780 101 54300068 C 1044.00 420 102 54210104 C 1084.00 810 103 54020637 C 1287.00 2970 104 54020732 C 1332.00 390 105 54130380 C 1334.00 600 106 54210180 C 1401.00 420 107 54300248 C 1452.00 390 108 54110386 C 2664.00 780 109 54020303 C 2980.00 1620 110 54300068 C 176.70 390 111 54210217 C 217.00 390 112 54210104 C 221.60 1080 113 54110160 C 369.70 720 114 54300160 C 373.00 420 115 54030451 C 484.40 840 116 54300156 C 512.10 390 117 54030135 C 574.50 390 118 54210252 C 613.00 780 119 54210181 C 613.30 1080 120 54130380 C 626.30 840 121 54060515 C 706.10 840 122 54140326 C 777.30 720 123 54030344 C 793.20 780 124 54110468 C 819.00 360 125 54300255 C 822.60 390 126 54110255 C 855.20 390 127 54140011 C 861.50 960 128 54300021 C 948.70 960 129 54110496 C 984.50 660 130 54300268 C 1069.00 780 131 54140295 C 1130.00 720 132 54110070 C 1177.00 780 133 54210180 C 1250.00 1170 134 54300217 C 1318.00 780 135 54210333 C 1382.00 870 136 54300248 C 1404.00 570 137 54110542 C 1446.00 360 138 54210085 C 2230.00 1170 139 54110386 C 2577.00 720 140 54020303 C 2924.00 1440 141 54210180 C 116.80 1080 142 54210181 C 229.70 960 143 54210217 C 267.90 270 144 54210252 C 285.80 690 145 54300068 C 311.80 390 146 54300156 C 470.20 420 147 54110255 C 485.90 390 148 54030451 C 526.40 780 149 54110374 C 628.10 720 150 54110070 C 687.10 780 151 54210333 C 702.70 1170 152 54030344 C 722.20 840 153 54140326 C 765.30 840 154 54130380 C 811.80 780 155 54300217 C 845.20 720 156 54030135 C 853.80 420 157 54300255 C 949.30 570 158 54210104 C 978.10 780 159 54140295 C 1076.00 780 160 54300268 C 1240.00 780 161 54210085 C 1391.00 1260 162 54300248 C 1738.00 780 grupo <- dados$Grupo pthi <- dados$Pthi dose <- dados$Dose_Cinacalcet grupo [1] H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H [38] H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H [75] H C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C [112] C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C [149] C C C C C C C C C C C C C C Levels: C H pthi [1] 221.60 317.00 412.60 462.00 474.90 547.00 547.50 589.00 709.00 [10] 804.90 842.00 855.40 911.50 979.30 1000.00 1054.00 1178.00 1409.00 [19] 1579.00 2794.00 204.00 236.00 253.00 334.00 358.00 360.80 426.00 [28] 494.50 535.90 594.60 615.70 659.00 666.10 705.70 765.20 875.80 [37] 953.20 982.00 1275.00 1297.00 1312.00 1330.00 1388.00 2060.00 2432.00 [46] 63.50 86.19 236.00 249.00 256.00 265.90 322.70 354.60 436.00 [55] 458.80 460.40 480.10 480.90 627.80 700.50 700.90 739.90 760.70 [64] 808.60 830.10 883.00 1033.00 1074.00 1138.00 1145.00 1174.00 1308.00 [73] 1401.00 1462.00 1536.00 76.50 159.30 161.10 268.00 280.90 364.10 [82] 574.30 603.90 651.80 727.10 738.60 745.00 764.50 771.20 772.30 [91] 788.70 807.50 813.60 856.30 865.10 882.70 883.90 916.30 923.70 [100] 984.70 1044.00 1084.00 1287.00 1332.00 1334.00 1401.00 1452.00 2664.00 [109] 2980.00 176.70 217.00 221.60 369.70 373.00 484.40 512.10 574.50 [118] 613.00 613.30 626.30 706.10 777.30 793.20 819.00 822.60 855.20 [127] 861.50 948.70 984.50 1069.00 1130.00 1177.00 1250.00 1318.00 1382.00 [136] 1404.00 1446.00 2230.00 2577.00 2924.00 116.80 229.70 267.90 285.80 [145] 311.80 470.20 485.90 526.40 628.10 687.10 702.70 722.20 765.30 [154] 811.80 845.20 853.80 949.30 978.10 1076.00 1240.00 1391.00 1738.00 dose [1] 1080 2790 930 900 930 930 930 930 930 930 930 930 930 930 930 [16] 840 900 930 420 1860 420 2790 930 930 390 930 930 1860 900 1860 [31] 900 1860 900 480 930 930 900 930 930 900 1440 930 900 2790 900 [46] 1860 1710 2040 1860 390 60 2790 1860 930 510 930 930 900 930 930 [61] 930 870 1860 870 900 930 390 900 2820 900 1350 870 1860 1860 2700 [76] 780 390 780 2070 390 720 390 390 390 360 390 390 1170 420 840 [91] 540 780 900 780 840 810 540 600 780 780 420 810 2970 390 600 [106] 420 390 780 1620 390 390 1080 720 420 840 390 390 780 1080 840 [121] 840 720 780 360 390 390 960 960 660 780 720 780 1170 780 870 [136] 570 360 1170 720 1440 1080 960 270 690 390 420 390 780 720 780 [151] 1170 840 840 780 720 420 570 780 780 780 1260 780 modelo1 <- lm(pthi~dose) modelo1 summary(modelo1) anova(modelo1) #criar e abrir arquivo jpeg("modelo1.jpg") plot(pthi~dose,xlab="dose",ylab="pthi","red") #Adicionar a reta da regressão no gráfico: abline(modelo1) #Adicionar a equação da reta no gráfico: text("y=738.3511+0.1177dose, p=0.1279") erro: plot(pthi~grupo,xlab="grupo",ylab="pthi","blue") Error in FUN(X[[1L]], ...) : invalid 'envir' argument > #Adicionar a reta da regressão no gráfico: > abline(modelo2) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet > #Adicionar a equação da reta no gráfico: > text("y=88.49-84.84grupo, p=0.1279") Error in text.default("y=88.49-84.84grupo, p=0.1279") : plot.new has not been called yet In addition: Warning message: In xy.coords(x, y, recycle = TRUE) : NAs introduced by coercion modelo2 <- lm(pthi~grupo) modelo2 summary(modelo2) plot(pthi~grupo,xlab="grupo",ylab="pthi","blue") Error in FUN(X[[1L]], ...) : invalid 'envir' argument > #Adicionar a reta da regressão no gráfico: > abline(modelo2) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet > #Adicionar a equação da reta no gráfico: > text("y=88.49-84.84grupo, p=0.1279") Error in text.default("y=88.49-84.84grupo, p=0.1279") : plot.new has not been called yet In addition: Warning message: In xy.coords(x, y, recycle = TRUE) : NAs introduced by coercion obrigada desde ja. cumprimentos Ana

2 1

Fwd: como extrair dados de arquivos GEOtiff
by D'Apolito 13 Nov '12

13 Nov '12

Prezados, Tenho arquivos no formato GEOtiff (ex. anexo) e quero extrair dados de temperatura média para cada coordenada geográfica. Estou usando isso até o momento: library(raster) f <- "tmean1_23.tif" r <- raster(f) v = getValues(r, format='matrix') Mas o que obtenho com getValues() é apenas uma matriz enorme com números que não tem nada a ver com temperatura. Meu nível de R não é muito avançado, por isso não estou conseguindo sair desse passo, agradeço qualquer ajuda. Att Carlos

2 1

dúvida gráficos de dispersão
by alanarocha＠sapo.pt 13 Nov '12

13 Nov '12

nBom dia, eu tenho um script onde teho tres modelos lineares (lm) e queria saber como faço um grafico que mostre as rectas de regressão de cada um. obrigada desde já. Cumprimentos Ana Rocha

2 1

Record Linkage - compare.dedup
by Fátima Lima Paula 12 Nov '12

12 Nov '12

Prezados, estou começando a estudar o artigo enviado como sugestão de leitura sobre o RecordLinkage. Logo no início diz que a função compare.dedup serve para deduplicação de um banco simples. Ou seja, pelo que entendi por essa função conseguimos os pares do banco. Bem, tentei e retornou um erro. Aí vai o CMR (é isso mesmo?). nome1=c("ANITA","ULISSES","CREUSA","HILDA","ANITA") nome3=c("COSTA","SILVA","RAMALHO","ALVES","COSTA") datanas=c("1947-06-01","1925-08-09","1946-02-03","1931-09-23","1947-06-01") teste=data.frame(nome1,nome3,datanas) teste nome1 nome3 datanas 1 ANITA COSTA 1947-06-01 2 ULISSES SILVA 1925-08-09 3 CREUSA RAMALHO 1946-02-03 4 HILDA ALVES 1931-09-23 5 ANITA COSTA 1947-06-01 rpairs=compare.dedup(teste, identity=identity.teste) Erro em identical(identity, NA) : objeto 'identity.teste' não encontrado Alguém que já tenha feito um linkage usando esse pacote pode me ajudar? Obrigada Fátima

2 2

Identificar cadastros com dados diferentes
by Rodrigo Coster 12 Nov '12

12 Nov '12

Caros, venho mais uma vez solicitar a ajuda de vocês. O sistema de onde trabalho me fornece um relatório que tem 13 informações, mas só 4 informações primordiais: ID, Contrato, Info1 e Info2. Os relatórios são mensais e cada ID pode ter mais de um contrato (mas somente um por mês). Estou montando um relatório geral, juntando informações dos últimos 5 anos. Era esperado que, para cada combinação ID-Contrato, Info1 e Info2 fossem iguais em todos os meses, mas não é o que está acontecendo, e preciso identificar os pares ID-Contrato que isso não ocorre, só que não me veio a cabeça nenhuma maneira prática de fazer isso. Por exemplo: (dados <- structure(list(ID = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), CONTRATO = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L), INFO1 = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L), INFO2 = c(1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L)), .Names = c("ID", "CONTRATO", "INFO1", "INFO2" ), class = "data.frame", row.names = c(NA, -29L))) Nesse exemplo ele deveria identificar: ID: 1 Contrato: 2 ID: 2 Contrato: 2 ID: 3 Contrato: 1 O banco total (a junção de todos relatórios) possui 200.000 linhas, portanto estou buscando alternativas antes de vasculhar linha por linha num for. Valeu!

2 2

sample condicional
by Daniel Marcelino 12 Nov '12

12 Nov '12

Bom dia, Alguém conhece uma maneira de viesar a função sample? Talvez a palavra correta seria 'condicionar' as amostras seguintes com base nas amostras anteriores. Eu gostaria de gerar varias distribuições categoricas mas não quero que o mesmo nível do fator se repita na linha. set.seed(51) data <- data.frame( id=as.factor(1:100), a=as.factor(sample(1:10, size=100, replace=TRUE)), b=as.factor(sample(1:10, size=100, replace=TRUE)), c=as.factor(sample(1:10, size=100, replace=TRUE)), d=as.factor(sample(1:10, size=100, replace=TRUE)), e=as.factor(sample(1:10, size=100, replace=TRUE)) ) -- "Small steps toward a much better world" Daniel Marcelino Land Phone 1+514 343 6111 #3799 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada

2 2

Ler .xls no Rcmdr no Ubuntu
by Andre Oliveira 11 Nov '12

11 Nov '12

Boa tarde, alguns dias atrás tive aqui com mesmo assunto na lista e não consegui resolver a questão. Pesquisei muito e não encontrei como carregar um .xls pelo menu (importar arquivos) do Rcmdr o que é possivel pelo Win$$. Atualmente estou trabalhando no curso de administração e gostaria de trabalhar com Ubuntu+R+Rcmdr. Na escola tem ubuntu em todas máquinas este o motivo de voltar ao assunto. Agora estou com dúvidas será que esta é uma tarefa impossível para Ubuntu+R+Rcmdr?

4 8

altera ano na data
by Fátima Lima Paula 11 Nov '12

11 Nov '12

Prezados, tenho uma banco com uma variável datao, da seguinte forma: datao=c("2008-03-25", "2008-03-24", "2008-03-25", "2008-03-25", "0008-12-20", "2008-03-19", "0011-04-23", "0009-04-01", "2009-12-26") Como se vê, existem alguns anos que foram digitados com erro: "0008-12-20", "0011-04-23" e "0009-04-01". Eu quero corrigir o erro. Consegui passar "0008" para "2008" usando os comandos: abc=as.POSIXlt(datao)$year+1900 e depois fazendo abc[abc==8]=2008, abc[abc==9]=2009, por exemplo, corrijo o erro de digitação, para o ano certo. Mas, não é isso que eu quero. O que eu quero é alterar o ano dentro da variável datao. Tem como fazer isso? Alguém pode me ajudar? Obrigada Fátima

2 4

Erro no rbind
by Fátima Lima Paula 11 Nov '12

11 Nov '12

Prezados, estou tentando juntar 4 bancos de dados com rbind e está retornando os seguintes erros. Os bancos têm as mesmas colunas, mas linhas diferentes. Tive um problema parecido, mas foi porque não estava sendo reconhecido o banco como data.frame. Dessa vez já usei as.data.frame(banco) para todos. Alguém pode me ajudar? Obrigada. sim=rbind(sim2008,sim2009,sim2010,sim2011) Mensagens de aviso perdidas: 1: In `[<-.factor`(`*tmp*`, ri, value = c(700L, 2030L, 602L, 240L, : nível de fator inválido, NAs gerados 2: In `[<-.factor`(`*tmp*`, ri, value = c(21320120L, 22745270L, 22775040L, : nível de fator inválido, NAs gerados 3: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 4: In `[<-.factor`(`*tmp*`, ri, value = c(1445L, 1240L, 640L, 1300L, : nível de fator inválido, NAs gerados 5: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, 22060020L, NA, 22775113L, : nível de fator inválido, NAs gerados 6: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 7: In `[<-.factor`(`*tmp*`, ri, value = c(230L, 515L, 1505L, 910L, : nível de fator inválido, NAs gerados 8: In `[<-.factor`(`*tmp*`, ri, value = c(23825400L, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 9: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados

8 14

Edição de arquivos .RD
by tiago souza marçal 10 Nov '12

10 Nov '12

Boa tarde prezados, Eu estou editando um arquivo .RD do arquivo man de um pacote e quando eu escrevi 5% não consegui gerar o pacote, portanto, tive que mudar para cinco por cento. Alguem sabe se existe algum comando para se adicionar o termo 5%?

3 2

Re: [R-br] Regressão logistica (conceito)
by viniciusbritor 09 Nov '12

09 Nov '12

amostrar com reposição toda a informação da classe alvo , Y=1 para que tenha o mesmo tamanho da classe Y=0. apenas na amostra de treinamento Enviado por Samsung MobileLeonard Mendonça de Assis <assis.leonard(a)gmail.com> escreveu:Vinícius o que você está chamando de reamostragem? []s Leonard de Assis http://about.me/ldeassis Em 06/11/2012 19:14, Vinicius Brito Rocha escreveu: Pessoal, preciso tirar uma dúvida a respeito de regressão logistica. Tenho uma conjunto de dados, onde existe um grande desbalanceamento nas classes da variável resposta (meu evento Y=1 a ser modelado, possui poucas observações) O que fiz foi: separei 70% dos dados para amostra treino e 30% para teste. na amostra treino (70%): utilizei uma re-amostragem, apenas na classe do evento (Y=1) e fiz com que as linhas desta classe fossem re-amostradas até que a classe com resposta Y=1 possui-se a mesma quantidade de linhas da classe (Y=0) Com minha base de amostra treino balanceada ajustei um modelo de regressão logistica na amostra teste(30%): calculei minhas probabilidades de respostas a partir das variáveis independentes da amostra teste. arbitrei um ponto de classificação de P_CHAP>= 0,7 para classificar meu evento como classe (Y_CHAP=1) e P_CHAP<0,7 classificar meu evento como Y_CHAP=0 Construi uma tabela de confusão comparando os resultados Y_CHAP e Y para comparar minha Sensitividade e 1-Especificidade. Dúvidas: A questão é que meus resultados estão muito ruins. Acredito que é a estrutura de dados. Estou sendo questionado a respeito da técnica de re-amostragem na amostra treino para equilibrar as classes. Pois o demandante acredita que é necessário algum tipo de correção do modelo(feito a re-amostragem no ajuste da amostra treino) ao aplica-lo no conjunto teste, que não sofreu nenhuma alteração. Alguém tem algum material que justifique o uso de re-amostragem nos dados da amostra treino? Abs. -- 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(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.

5 7

Re: [R-br] Regressão logistica (conceito)
by Roney Fraga 08 Nov '12

08 Nov '12

Vinicius, Se o problema for a quantidade de eventos desproporcionais acredito que os capítulos 4, 5 e 6 dessa dissertação pode ajudar. http://www.dca.fee.unicamp.br/~vonzuben/research/semolini_mest.html Abraço -- Roney

1 0

erro em OpenStreetMap
by D'Apolito 08 Nov '12

08 Nov '12

Olá, alguém poderia me ajudar com esse erro? Windows, R versão 2.15.2, 64-bit > library(OpenStreetMap) Carregando pacotes exigidos: rJava *Error *in get(Info[i, 1], envir = env) : internal error -3 in R_decompress1 *Erro*: package ‘rJava’ could not be loaded obrigado Carlos

3 3

Re: [R-br] Erro no rbind
by Humberto Hazin 08 Nov '12

08 Nov '12

Fátima vc pode mandar esses dois arquivos (ou parte deles) por email para hghazin(a)hotmail.com ou disponibilizar um link de download, de forma que possa olhar melhor ele! Humberto Enviado do Email do Windows De: Fátima Lima Paula Enviado: ‎8‎ de ‎novembro‎ de ‎2012 ‎20‎:‎00 Para: r-br(a)listas.c3sl.ufpr.br, Alisson Lucrecio Assunto: Re: [R-br] Erro no rbind Os nomes são os mesmos e o número de linhas não. Em 8 de novembro de 2012 19:55, Alisson Lucrecio <alissonluc(a)yahoo.com.br> escreveu: Fátima, confere suas data.frame se ela tem o mesmo tamanho de linhas. Att Alisson Lucrécio da Costa From: Fátima Lima Paula <fatima.lima.paula(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Sent: Thursday, November 8, 2012 7:37 PM Subject: Re: [R-br] Erro no rbind Desculpe, sou uma anta. Li o help, mas nem reparei esse detalhe. De qualquer forma tentei, mas reportou a mesma advertência: sim=merge(sim2008,sim2009) Mensagens de aviso perdidas: 1: In `[<-.factor`(`*tmp*`, ri, value = c(700L, 2030L, 602L, 240L, : nível de fator inválido, NAs gerados 2: In `[<-.factor`(`*tmp*`, ri, value = c(21320120L, 22745270L, 22775040L, : nível de fator inválido, NAs gerados 3: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados Na verdade, não sei o que querem dizer essas mensagens, porque se faço outros comandos, depois de rbind ele retorna os valores como se o data.frame "sim" tivesse sido construído. Por exemplo, faço: sim=rbind(sim2008,sim2009,sim2010,sim2011) Mensagens de aviso perdidas: 1: In `[<-.factor`(`*tmp*`, ri, value = c(700L, 2030L, 602L, 240L, : nível de fator inválido, NAs gerados 2: In `[<-.factor`(`*tmp*`, ri, value = c(21320120L, 22745270L, 22775040L, : nível de fator inválido, NAs gerados 3: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 4: In `[<-.factor`(`*tmp*`, ri, value = c(1445L, 1240L, 640L, 1300L, : nível de fator inválido, NAs gerados 5: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, 22060020L, NA, 22775113L, : nível de fator inválido, NAs gerados 6: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 7: In `[<-.factor`(`*tmp*`, ri, value = c(230L, 515L, 1505L, 910L, : nível de fator inválido, NAs gerados 8: In `[<-.factor`(`*tmp*`, ri, value = c(23825400L, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 9: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados > dim(sim) [1] 163842 79 Em 8 de novembro de 2012 19:24, Rodrigo Coster <rcoster(a)gmail.com> escreveu: Tu ao menos leu o help do merge? Pra utilizar ele, tem que juntar 2 a 2... 2012/11/8 Fátima Lima Paula <fatima.lima.paula(a)gmail.com> Tentei, mas não sei se usei corretamente. Coloquei: sim=merge(sim2008,sim2009,sim2010,sim2011) Erro em fix.by(by.x, x) : 'by' deve especificar coluna(s) como números, nomes ou lógico Não entendi como especificar colunas Em 8 de novembro de 2012 17:53, FHRB Toledo <fernandohtoledo(a)gmail.com> escreveu: Fatima, Não seria o caso de usar o merge? att, FH 2012/11/8 Fátima Lima Paula <fatima.lima.paula(a)gmail.com> Prezados, estou tentando juntar 4 bancos de dados com rbind e está retornando os seguintes erros. Os bancos têm as mesmas colunas, mas linhas diferentes. Tive um problema parecido, mas foi porque não estava sendo reconhecido o banco como data.frame. Dessa vez já usei as.data.frame(banco) para todos. Alguém pode me ajudar? Obrigada. sim=rbind(sim2008,sim2009,sim2010,sim2011) Mensagens de aviso perdidas: 1: In `[<-.factor`(`*tmp*`, ri, value = c(700L, 2030L, 602L, 240L, : nível de fator inválido, NAs gerados 2: In `[<-.factor`(`*tmp*`, ri, value = c(21320120L, 22745270L, 22775040L, : nível de fator inválido, NAs gerados 3: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 4: In `[<-.factor`(`*tmp*`, ri, value = c(1445L, 1240L, 640L, 1300L, : nível de fator inválido, NAs gerados 5: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, 22060020L, NA, 22775113L, : nível de fator inválido, NAs gerados 6: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 7: In `[<-.factor`(`*tmp*`, ri, value = c(230L, 515L, 1505L, 910L, : nível de fator inválido, NAs gerados 8: In `[<-.factor`(`*tmp*`, ri, value = c(23825400L, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 9: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados _______________________________________________ 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. _______________________________________________ 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. _______________________________________________ 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

Re: [R-br] Erro no rbind
by Humberto Hazin 08 Nov '12

08 Nov '12

Fátima verifique se as duas planilhas os nomes das colunas são idênticos nos dois arquivos! Humberto Enviado do Email do Windows De: Fátima Lima Paula Enviado: ‎8‎ de ‎novembro‎ de ‎2012 ‎19‎:‎37 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Erro no rbind Desculpe, sou uma anta. Li o help, mas nem reparei esse detalhe. De qualquer forma tentei, mas reportou a mesma advertência: sim=merge(sim2008,sim2009) Mensagens de aviso perdidas: 1: In `[<-.factor`(`*tmp*`, ri, value = c(700L, 2030L, 602L, 240L, : nível de fator inválido, NAs gerados 2: In `[<-.factor`(`*tmp*`, ri, value = c(21320120L, 22745270L, 22775040L, : nível de fator inválido, NAs gerados 3: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados Na verdade, não sei o que querem dizer essas mensagens, porque se faço outros comandos, depois de rbind ele retorna os valores como se o data.frame "sim" tivesse sido construído. Por exemplo, faço: sim=rbind(sim2008,sim2009,sim2010,sim2011) Mensagens de aviso perdidas: 1: In `[<-.factor`(`*tmp*`, ri, value = c(700L, 2030L, 602L, 240L, : nível de fator inválido, NAs gerados 2: In `[<-.factor`(`*tmp*`, ri, value = c(21320120L, 22745270L, 22775040L, : nível de fator inválido, NAs gerados 3: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 4: In `[<-.factor`(`*tmp*`, ri, value = c(1445L, 1240L, 640L, 1300L, : nível de fator inválido, NAs gerados 5: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, 22060020L, NA, 22775113L, : nível de fator inválido, NAs gerados 6: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 7: In `[<-.factor`(`*tmp*`, ri, value = c(230L, 515L, 1505L, 910L, : nível de fator inválido, NAs gerados 8: In `[<-.factor`(`*tmp*`, ri, value = c(23825400L, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 9: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados > dim(sim) [1] 163842 79 Em 8 de novembro de 2012 19:24, Rodrigo Coster <rcoster(a)gmail.com> escreveu: Tu ao menos leu o help do merge? Pra utilizar ele, tem que juntar 2 a 2... 2012/11/8 Fátima Lima Paula <fatima.lima.paula(a)gmail.com> Tentei, mas não sei se usei corretamente. Coloquei: sim=merge(sim2008,sim2009,sim2010,sim2011) Erro em fix.by(by.x, x) : 'by' deve especificar coluna(s) como números, nomes ou lógico Não entendi como especificar colunas Em 8 de novembro de 2012 17:53, FHRB Toledo <fernandohtoledo(a)gmail.com> escreveu: Fatima, Não seria o caso de usar o merge? att, FH 2012/11/8 Fátima Lima Paula <fatima.lima.paula(a)gmail.com> Prezados, estou tentando juntar 4 bancos de dados com rbind e está retornando os seguintes erros. Os bancos têm as mesmas colunas, mas linhas diferentes. Tive um problema parecido, mas foi porque não estava sendo reconhecido o banco como data.frame. Dessa vez já usei as.data.frame(banco) para todos. Alguém pode me ajudar? Obrigada. sim=rbind(sim2008,sim2009,sim2010,sim2011) Mensagens de aviso perdidas: 1: In `[<-.factor`(`*tmp*`, ri, value = c(700L, 2030L, 602L, 240L, : nível de fator inválido, NAs gerados 2: In `[<-.factor`(`*tmp*`, ri, value = c(21320120L, 22745270L, 22775040L, : nível de fator inválido, NAs gerados 3: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 4: In `[<-.factor`(`*tmp*`, ri, value = c(1445L, 1240L, 640L, 1300L, : nível de fator inválido, NAs gerados 5: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, 22060020L, NA, 22775113L, : nível de fator inválido, NAs gerados 6: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 7: In `[<-.factor`(`*tmp*`, ri, value = c(230L, 515L, 1505L, 910L, : nível de fator inválido, NAs gerados 8: In `[<-.factor`(`*tmp*`, ri, value = c(23825400L, NA, NA, NA, NA, : nível de fator inválido, NAs gerados 9: In `[<-.factor`(`*tmp*`, ri, value = c(NA, NA, NA, NA, NA, NA, NA, : nível de fator inválido, NAs gerados _______________________________________________ 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. _______________________________________________ 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

Teste de nº obs
by Carlos Mendonça 08 Nov '12

08 Nov '12

Prezados, como posso fazer para testar se um data tem observações ou não. Por exemplo, se o data não tiver obs eu quero terminar o programa, caso contrário, continue executando. Carlos Mendonça.

2 2

Dúvida gráfico
by Sérgio Henrique almeida da silva ju 08 Nov '12

08 Nov '12

Prezados Estou precisando gerar um gráfico tipo esse: http://www.datafilehost.com/download-a49920d3.html Observem que existem diferenças nas alturas de cada grupo, isso ocorre porque cada grupo, possui quantidade de doenças diferentes Alguém tem alguma ideia como posso gerar isso? Obrigado -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 2

Encaminahndo: ... Comparação entre estimação ...
by FHRB Toledo 08 Nov '12

08 Nov '12

Matheus, Mandando suas dúvidas para a lista do R, suas chances de alguém ajudá-lo são bem maiores. Leia o Guia de Postagens e produza um CMR, consulte: http://www.leg.ufpr.br/doku.php/software:rbr. att, FH 2012/11/8 Matheus Monteiro <matheusmont(a)gmail.com> > boa tarde fernando, sera que pode me ajudar? > > fitdistr(y[,1],"t") > m s df > 26.847104 22.705444 2.502808 > ( 7.253774) ( 9.024268) ( 2.189144) > > minha intenção é plotar essa informação e depois compara-la com a normal > como devo proceder? > > abração >

1 0

o arqui
by alanarocha＠sapo.pt 08 Nov '12

08 Nov '12

arquivo

2 1

dúvida para trabalhar com arquivo importado
by Matheus Monteiro 08 Nov '12

08 Nov '12

é o seguinte: estou trabalhando com a distribuição t quero estimar um conjunto de dados eu chamando esse conjunto de dados de x fazendo y <- c(...........................) tá ok por ex: y=c(49,-67,8,16,6,23,28,41,14,29,56,24,75,60,-48) fitdistr (y, "t") mas se meus dados importados de um arquivo txt (vamos chamar de X) eu leio os dados, aparece tudo certinho porém a estimação não acontece

5 8

Material introdutório à geoestatística
by Ligia Souza 08 Nov '12

08 Nov '12

Boa tarde. Estou começando a estudar geoestatística. Alguém poderia me indicar leituras introdutórias à este assunto? Att., Ligia Pereira de Souza Graduanda em Ecologia UNESP - Rio Claro/SP

5 6

duvida em diferenciação de séries temporais
by Natalia Martins 08 Nov '12

08 Nov '12

Bom dia prezados, Estou finalizando minha dissertação e estou com um problema se alguem puder ajudar fico muito grata. Tenho uma serie temporal de temperatura para essa ajustei um modelo que tem uma diferença sazonal fiz as previsoes 6 passos a frente, no entanto minhas previsoes saem na "escala" da serie diferenciada(-0.45, 1.21, -0.89) e eu gostaria de volta-las a escala original (17,19,16,21....) A classe de modelos que estou estudando não tem implentado em nenhum software. Obrigada. -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

1 0

interpretação dos resultados
by alanarocha＠sapo.pt 07 Nov '12

07 Nov '12

Boa tarde, eu tenho um script cujo o bjectivo é relacionar as duas variaveis com a .variavel de control.. Queria fazer um texto que esplicasse os resultados mas não estou a conseguir interpretalos. Gostava de uma ajuda. Cumprimentos Ana

2 1

norma de matriz complexa
by Cleber N.Borges 07 Nov '12

07 Nov '12

olá a todos, alguém conhece alguma função do R para calcular a norma de matrizes complexas? (ou: existe um jeito mais direto de calcular isso?) Agradeço antecipadamente por qualquer ajuda. Cleber

1 0

Consistência de um estimador...
by andrebvs 07 Nov '12

07 Nov '12

3 3

Tinn-R 2.4.0.1 liberado
by Jose Claudio Faria 07 Nov '12

07 Nov '12

Caros usuários do projeto Tinn-R. Uma nova versão com vários melhoramentos foi liberada hoje no servidor SourceForge. Para conhecer as novidades: http://sourceforge.net/p/tinn-r/news/2012/11/tinn-r-2401-released/ Att, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

teste de hipotese
by Edson Lira 06 Nov '12

06 Nov '12

Quero pedir desculpas aos colegas, pois em um email anterior coloquei um anexo acidentalmente, errei o endereço. Edson Lira Estatístico Manaus-Amazonas

1 0

Regressão logistica (conceito)
by Vinicius Brito Rocha 06 Nov '12

06 Nov '12

Pessoal, preciso tirar uma dúvida a respeito de regressão logistica. Tenho uma conjunto de dados, onde existe um grande desbalanceamento nas classes da variável resposta (meu evento Y=1 a ser modelado, possui poucas observações) O que fiz foi: separei 70% dos dados para amostra treino e 30% para teste. na amostra treino (70%): - utilizei uma re-amostragem, apenas na classe do evento (Y=1) e fiz com que as linhas desta classe fossem re-amostradas até que a classe com resposta Y=1 possui-se a mesma quantidade de linhas da classe (Y=0) - Com minha base de amostra treino balanceada ajustei um modelo de regressão logistica na amostra teste(30%): - calculei minhas probabilidades de respostas a partir das variáveis independentes da amostra teste. - arbitrei um ponto de classificação de P_CHAP>= 0,7 para classificar meu evento como classe (Y_CHAP=1) e P_CHAP<0,7 classificar meu evento como Y_CHAP=0 - Construi uma tabela de confusão comparando os resultados Y_CHAP e Y para comparar minha Sensitividade e 1-Especificidade. Dúvidas: A questão é que meus resultados estão muito ruins. Acredito que é a estrutura de dados. - Estou sendo questionado a respeito da técnica de re-amostragem na amostra treino para equilibrar as classes. Pois o demandante acredita que é necessário algum tipo de correção do modelo(feito a re-amostragem no ajuste da amostra treino) ao aplica-lo no conjunto teste, que não sofreu nenhuma alteração. Alguém tem algum material que justifique o uso de re-amostragem nos dados da amostra treino? Abs. -- *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.

2 1

[Dúvida] Listar todos os comandos.
by Pedro Rafael 06 Nov '12

06 Nov '12

Pessoal uma dúvida besta. Como faço para listar todos os comandos de um determinado pacote no R? Obrigado deste já. [ ], Pedro Rafael Diniz Marinho.

3 4

Gráfico de previsão ARIMA
by kaue veras 06 Nov '12

06 Nov '12

Boa noite Pessoal, Será que vocês poderiam me ajudar, pois estou fazendo meu trabalho de monografia baseado no modelo ARIMA. Alguém saberia como eu poderia gerar um gráfico de previsão do modelo ARIMA? Desde já agradeço. Atenciosamente, Kauê P. Veras da Cunha Estatística - UERJ MSN: kaueveras(a)hotmail.com E-mail: kaueveras(a)gmail.com Cel: (21) 8254-9601

2 1

Dúvida sobre Regressão Logística Condicional
by Victor Eduardo 05 Nov '12

05 Nov '12

Pessoal, estou tentando rodar o script abaixo, mas está apresentando o seguinte erro: clogit(con1$status~con1$tnf+strata(con1$cond),data=con1)Erro em Surv(rep(1, 143L), con1$status) : Invalid status value, must be logical or numeric Tentei rodar sem o argumento data, mas apresentou o mesmo erro. OBS: con1 é o meu banco de dados A classe de cada uma das variáveis são: > class(con$status)[1] "factor"> class(con$tnf)[1] "genotype" "factor" > class(con$cond)[1] "integer" Abraços!

2 3

Duvida, como exibir help-pages do R no navegador usando Ubunto-linux e emacs?
by Augusto Ribas 05 Nov '12

05 Nov '12

Ola colegas, recentemente começei a usar o ubunto e estou tendo algumas dificuldades básicas. Estou tentando usar o emacs como editor para o R e minha duvida é de como fazer para as paginas de documentação, quando se usa o comando ?função aparecam no navegador. Eu não entendi se isso tem que ser configurado no arquivo .emacs ou no arquivo .Rprofile Alias o arquivo .Rprofile não é criado automaticamente quando instalamos o R, ele tem que ser criado manualmente no diretorio home? no arquivo .emacs eu coloquei a seguinte linha: (setq ess-help-own-frame t) Mas somente qd digito por exemplo ?rnorm no buffer que esta o R que abre um novo buffer com o help. Se digito ?rnorm no buffer que estou fazendo um script o help sai na pagina do R. Eu gostaria que de ambas as formas a documentação aparecesse no navegador com hyperlinks. Aproveitando, eu gostaria de saber também como o emacs pode colorir todos os comando do R. Atualmente se eu digito algumas coisas como require no buffer de script e o emacs pinta de azul, mas comando por exemplo rnorm( ele somente mostra os argumentos la em baixo, mas não pinta com uma cor diferente. Vou colar o que ja fiz de modificação no arquivo.emacs: ########################################################################################################################### (require 'ess-eldoc) ;; numeração das linhas na margem esquerda (global-linum-mode 1) ;; inicia Emacs com ctrl-{zxcv} abilitado para desf/recor/cop/colar (cua-mode t) ;; realçador de pareamento de parenteses, chaves, colchetes, aspas... (show-paren-mode 1) ;; tamanho da fonte (120~12pt) ;; comandos C-x C-+ e C-x C-- para aumentar e diminiur a fonte (set-face-attribute 'default nil :height 120) ;;-------------------------------------------------------------------- ;; aux. na identação do código, colunas coloridas após certa distância (require 'whitespace) (setq whitespace-line-column 77) (setq whitespace-style '(face tabs lines-tail trailing)) (global-whitespace-mode t) ;;-------------------------------------------------------------------- ;;;;; create a new frame for each help instance (setq ess-help-own-frame t) ############################################################################################################################ O arquivo .Rprofile eu não achei ainda, ou não criei se eu tenho que criar ele. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.blogspot.com.br/<http://augustoribas.heliohost.org> Lattes: http://lattes.cnpq.br/7355685961127056

5 6

Re: [R-br] Fw: extrair bordas de arquivos .shp
by Eder Comunello 05 Nov '12

05 Nov '12

Mauro e colegas, boa tarde! As operações com polígonos do {maptools} dependem do pacote {gpclib} que tem licença restrita. Ele é desativado por default, por isso usa-se o gpclibPermit(). Verifique a presença do pacote {gpclib}. Pra funcionar gpclibPermit() deve retornar TRUE. PS: peço enviar as dúvidas diretamente para a lista, assim outras pessoas poderão ajudar ou ter ideias melhores. Saudações, -- ================================================ É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 Em 5 de novembro de 2012 00:57, Mauro Sznelwar <sznelwar(a)uol.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=sznelwar@uol.com.br> > escreveu: > ** > Tentei rodar o seu script, mas a partir deste comando dá pau! > > > gpclibPermit() ## para permitir unionSpatialPolygons > [1] FALSE > > borda.spol <- (unionSpatialPolygons(ctba, ID=rep(1, length(ctba)), > avoidGEOS=T)) > Erro: isTRUE(gpclibPermitStatus()) is not TRUE > > borda <- data.frame(borda.spol@polygons[[1]]@Polygons[[1]]@coords<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=borda.spol@polygons[[1]]…> > ) > Erro em data.frame(borda.spol@polygons[[1]]@Polygons[[1]]@coords<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=borda.spol@polygons[[1]]…>) > : > objeto 'borda.spol' não encontrado > > >

1 0

Re: [R-br] Eclipse Juno
by Leandro Marino 05 Nov '12

05 Nov '12

*Izidro,* * * *fiz exatamente o que você falou e é o mesmo procedimento que faço no linux. Mas ao fazer a busca no diretorio ele busca este xml que parece estar faltando. * Não sei o que pode estar acontecendo. O interessante é que algumas horas antes havia feito o mesmo procedimento no linux e tudo funcionou perfeitamente. Obrigado pela ajuda! Um abraço, Leandro 2012/11/3 <izidro3.1416(a)gmail.com> > Olá Leandro > > Verifique o endereço! > O endereço correto é http://download.walware.de/eclipse-3.8 > > e NÃO http://download.walware.de/eclipse-3.8/content.xml > > Veja no site da WalWare > http://www.walware.de/goto/statet > > []' > > Izidro >

1 1

Re: [R-br] Material introdutório à geoestatística
by Eder Comunello 04 Nov '12

04 Nov '12

Olá, Ligia! Tem um material que encontrei hoje que está bem didático. É do Prof. Landim, também de Rio Claro. Você consegue o texto completo no link que consta na referência que segue. LANDIM, Paulo M. Barbosa. Sobre Geoestatística e mapas.* Terrae didat.*, Campinas, v. 2, n. 1, 2006 . Disponível em < http://ppegeo.igc.usp.br/scielo.php?script=sci_arttext&pid=S1980-4407200600…>. acessos em 04 nov. 2012. Boa leitura, -- ================================================ É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

1 0

Re: [R-br] extrair bordas de arquivos .shp
by Eder Comunello 04 Nov '12

04 Nov '12

Olá, pessoal! Tentei dar um 'up' no CMR, seguindo a dica do Walmes... ############################### #install.packages(c("geoR", maptools", "sp"), dep=T) require(geoR); require(maptools); require(sp) ### download direto! ;D ### joga na pasta/dir de trabalho ativa download.file('http://dl.dropbox.com/u/117618178/ctba/bairros.shp', paste(getwd(),'bairros.shp',sep='/'), mode='wb') download.file('http://dl.dropbox.com/u/117618178/ctba/bairros.dbf', paste(getwd(),'bairros.dbf',sep='/'), mode='wb') download.file('http://dl.dropbox.com/u/117618178/ctba/bairros.shx', paste(getwd(),'bairros.shx',sep='/'), mode='wb') ctba <- readShapePoly("bairros.shp", IDvar="CODE") slotNames(ctba) npol <- as.integer(length(ctba@polygons)) ### numero de poligonos for (i in 1:npol){ nsubpol <- length(ctba@polygons[[i]]@Polygons) ### numero de poligonos internos for (ii in 1:nsubpol){ pol <- ctba@polygons[[i]]@Polygons[[ii]]@coords assign(paste("P",formatC(i, flag=0, width=3),ii,sep="-"),pol) } } ### unir todos os poligonos e então exporta uma borda unica gpclibPermit() ## para permitir unionSpatialPolygons borda.spol <- (unionSpatialPolygons(ctba, ID=rep(1, length(ctba)), avoidGEOS=T)) borda <- data.frame(borda.spol@polygons[[1]]@Polygons[[1]]@coords) ############################### Observações sobre o download.file() e Dropbox 1) arquivos locados em um subdiretório da pasta 'Public'; 2) funcionou tanto com http quanto https; 3) precisei definir mode='wb' (modo binário) para não corromper os arquivos na gravação. R version 2.15.2 (2012-10-26) Platform: x86_64-w64-mingw32/x64 (64-bit) -- ================================================ É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

1 0

Re: [R-br] [Dúvida] Função Gama Incompleta.
by Pedro Rafael 04 Nov '12

04 Nov '12

Mesmo assim deu valores diferentes. hy Mathematica: Gamma[1.5, 0, 0.10] = 0.019861 OBS: O argumento zero é para dizer que estou calculando a gamma inferior, integral de 0 à x; R: pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214 Uma outra dúvida que tenho é sobre a função hipergeométrica 2F1. Estou utilizando o pacote hypergeo. A minha função que envolva a função hypergeo é: fda_bllg <- function(par,x){ a <- par[1] b <- par[2] alpha <- par[3] Beta <- par[4] (((x^(-Beta))^(-a))*((1/alpha)^(a*Beta))*Re(hypergeo(a, a+b, 1+a, -(x/alpha)^Beta, tol = 0, maxiter=20000)))/(a*beta(a,b)) } fda_bllg trata-se da função de distribuição acumulada da distribuição Beta Log-Logistica KLL, uma distribuição com 4 parâmetros utilizada em análise de sobrevivência e suporte x>=0. Quando substituo os estimadores de máxima verossimilhança ocorre erro. As EMV da distribuição KLL são: est = c(6.474587e+02,2.043228e+02, 5.867625e-03,1.210073e-01) Seja x = 10. Temos que: fda_bllg(est,10) A mensagem de erro é: [1] Inf Mensagens de aviso perdidas: 1: In i15.3.8(A, B, C) : valor fora de limites em 'gammafn' [ ], Pedro Rafael Diniz Marinho. Em 1 de novembro de 2012 14:41, Pedro Rafael <pedro.rafael.marinho(a)gmail.com > escreveu: > Mesmo assim deu valores diferentes. > hy > > > Mathematica: > Gamma[1.5, 0, 0.10] = 0.019861 > > OBS: O argumento zero é para dizer que estou calculando a gamma inferior, > integral de 0 à x; > > R: > > pgamma(1.5, .10, lower=FALSE)*gamma(.10) = 0.1079214 > > > Uma outra dúvida que tenho é sobre a função hipergeométrica 2F1. Estou > utilizando o pacote hypergeo. A minha função que envolva a função hypergeo > é: > > fda_bllg <- function(par,x){ > a <- par[1] > b <- par[2] > alpha <- par[3] > Beta <- par[4] > (((x^(-Beta))^(-a))*((1/alpha)^(a*Beta))*Re(hypergeo(a, a+b, 1+a, > -(x/alpha)^Beta, tol = 0, maxiter=20000)))/(a*beta(a,b)) > } > > fda_bllg trata-se da função de distribuição acumulada da distribuição Beta > Log-Logistica KLL, uma distribuição com 4 parâmetros utilizada em análise > de sobrevivência e suporte x>=0. Quando substituo os estimadores de máxima > verossimilhança ocorre erro. > > As EMV da distribuição KLL são: > > est = c(6.474587e+02,2.043228e+02, 5.867625e-03,1.210073e-01) > > Seja x = 10. Temos que: > > fda_bllg(est,10) > > A mensagem de erro é: > [1] Inf > Mensagens de aviso perdidas: > 1: In i15.3.8(A, B, C) : valor fora de limites em 'gammafn' > > [ ], > Pedro Rafael Diniz Marinho. > > > Em 1 de novembro de 2012 06:44, beniltoncarvalho [via R-br] < > ml-node+s2285057n4656804h1(a)n4.nabble.com> escreveu: > > o problema e' q vc vai precisar revisitar as definicoes da funcao gama >> incompleta... ler a documentacao detalhada de pgamma() e comparar com o que >> Gamma[] no Mathematica faz... >> >> o que vc quer e': >> >> pgamma(1.5, .10, lower=FALSE)*gamma(.10) >> >> >> 2012/11/1 Pedro Rafael <[hidden email]<http://user/SendEmail.jtp?type=node&node=4656804&i=0> >> > >> >>> Pessoal qual o comando certo ou pacote que trabalha com a função gamma >>> incompleta? Estou obtendo valores muito diferentes da gamma incompleta no R >>> e no mathematica. >>> >>> >>> Por exemplo: >>> >>> No R: pgamma(.10,1.5) = 0.0224107 >>> No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921 >>> >>> Qual seria o problema que não estou enxergando? >>> >>> [ ], >>> Pedro Rafael Diniz Marinho. >>> >>> _______________________________________________ >>> R-br mailing list >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4656804&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=4656804&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-Duvida-Funcao-Gama-Incompleta-tp4656… >> 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…> >> > >

2 6

Re: [R-br] extrair bordas de arquivos .shp
by Eder Comunello 04 Nov '12

04 Nov '12

Boa tarde, Samuel e colegas! Segue uma sugestão de código... Testei com a delimitação dos bairros de Curitiba, do material do Prof. Paulo Justiniano Para testar será necessário baixar os 3 arquivos disponibilizados a seguir: https://dl.dropbox.com/u/117618178/ctba/bairros.shp' https://dl.dropbox.com/u/117618178/ctba/bairros.dbf' https://dl.dropbox.com/u/117618178/ctba/bairros.shx' Ou então baixar e descompactar: https://dl.dropbox.com/u/117618178/ctba/bairros.zip' ######################## #install.packages(c("geoR", maptools", "sp"), dep=T) require(geoR); require(maptools); require(sp) ctba <- readShapePoly("bairros.shp", IDvar="CODE") slotNames(ctba) npol <- as.integer(length(ctba@polygons)) ### numero de poligonos for (i in 1:npol){ nsubpol <- length(ctba@polygons[[i]]@Polygons) ### numero de poligonos internos for (ii in 1:nsubpol){ pol <- ctba@polygons[[i]]@Polygons[[ii]]@coords #assign(paste("P",i,ii,sep="-"),pol) assign(paste("P",formatC(i, flag=0, width=3),ii,sep="-"),pol) } } ######################## Neste caso aí serão originados 75 objetos correspondentes ao 75 polígonos. O código prevê sub-polígonos, pois é bem normal que eles existam. Não sei o que você pretende com esta operação, mas te confesso que fico pensando que ela não a melhor solução para "bordas". Sendo assim, acrescento outro trecho de código para unir todos os polígonos e feito isso, exportar uma borda única. ### une todos os poligonos e então exporta uma borda unica gpclibPermit() ## para permitir unionSpatialPolygons borda.spol <- (unionSpatialPolygons(ctba, ID=rep(1, length(ctba)), avoidGEOS=T)) borda.df <- data.frame(borda.spol@polygons[[1]]@Polygons[[1]]@coords) ######################## Estou usando R version 2.15.2 (2012-10-26) - Platform: x86_64-w64-mingw32/x64 (64-bit) Espero que ajude, -- ================================================ É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

2 1

extrair bordas de arquivos .shp
by Samuel Carvalho 04 Nov '12

04 Nov '12

Bom dia pessoal, Alguem poderia ajudar com a seguinta dúvida? Estou usando o pacote sp pra visualização e manipulação de dados espaciais no R. Um dos arquivos que estou trabalhando é do tipo .shp Qual procedimento para extrair as bordas dos polígonos deste arquivo que não seja tão lusitânica quanto o seguinte código: dados <- readShapePoly('dados.shp') length(dados@polygons) #supondo que o resultado deste comando seja == 2 length(dados@polygons[[1]]@Polygons) #e supondo que o resultado aqui seja == 10 length(dados@polygons[[2]]@Polygons) #e supondo que o resultado aqui seja == 10 Entao meu codigo para extrair as coordenadas ditas como bordas ou limites seria dados@polygons[[1]]@Polygons[[1]])@coords dados@polygons[[1]]@Polygons[[2]])@coords . . . dados@polygons[[1]]@Polygons[[10]])@coords e estes mesmos comandos para o poligono 2 o que totaliza 20 linhas de comando Talvez existisse algo como dados@borders mas não existe pois é possível verificar o que se pode extrair do objeto dados que é um objeto da classe "SpatialPolygonsDataFrame" pelo comando getSlots('SpatialPolygonsDataFrame') Obrigado Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

4 3

Fwd: Eclipse Juno
by Leandro Marino 03 Nov '12

03 Nov '12

*Izidro, não estou conseguindo acessar de jeito nenhum.* * * *Tentei também no juno e o erro se repete. Será que o servidor está fora do ar?* 2012/11/3 <izidro3.1416(a)gmail.com> > Olá Leandro > > Eu uso o eclipse Juno com o repositório > http://download.walware.de/eclipse-3.8 > > Izidro >

1 0

Eclipse Juno
by Leandro Marino 03 Nov '12

03 Nov '12

*Caros,* * * *formatei meu notebook e estou instalando alguns programas que considero adequados para o uso no dia a dia. Um deles é o StatET com o Eclipse.* * * *Baixei o eclipse e notei que é a versão 4.2 e ainda não há nada disponível no StatET. Baixei a versão indigo e tentei configurar e não ob**tive sucesso está fornecendo o seguinte erro.* * * * Unable to read repository at http://download.walware.de/eclipse-3.8/content.xml. Unable to read repository at http://download.walware.de/eclipse-3.8/content.xml. Cannot assign requested address: connect alguem sabe como posso resolver isto? *

1 0

Especificação de Pesos em glm(family=poisson)
by Diego Pujoni 02 Nov '12

02 Nov '12

Olá Pessoal, possuo uma variável preditora que é o número de espécies exóticas de peixes encontrados no lago i. Entretanto, para cada lago, eu possuo um esforço amostral diferente, isto é, existem lagos que foram amostrados apenas uma vez e lagos que foram amostrados 10 vezes. De forma que eu possuo uma incerteza maior para o lago que só foi amostrado uma única vez. Gostaria de incluir estas incertezas no ajuste e tentei utilizar o argumento "weights". Entretanto a escala dos meus pesos influencia muito o desvio residual. Existe uma forma correta para especificar estes pesos? Obrigado Sp = c(2,2,3,2,2,4,0,0,3,3,0,1,5,2,3,2,0,2,3) Peso = c(8,8,10,8,8,10,8,8,10,10,1,1,2,1,1,1,2,1,1) Var1 = c(1.48,1.65,2.50,2.20,1.48,2.29,2.83,3.59,2.29,1.65,2.70,4.94, 3.96,15.30,1.01,1.01,0.51,1.79,1.79) Var2 = c(1.48,0.00,0.00,0.00,0.00,0.00,2.83,3.59,0.00,0.00,14.99,4.94, 0.00,15.30,0.00,0.00,0.51,0.00,0.00) GLM = glm(Sp~Var1+Var2-1,family=poisson) summary(GLM) GLM1 = glm(Sp~Var1+Var2-1,weights = Peso, family=poisson) summary(GLM1) Diego PJ

2 1

Re: [R-br] selecionar intervalo de datas
by Eder Comunello 02 Nov '12

02 Nov '12

Exemplo de operação com datas com seleção por subset ############ dados<-data.frame( data=c("01/01/2004","01/02/2004","01/03/2004","01/04/2004","01/05/2004","01/06/2004"), tmin=c(17.82966,18.05029,22.47399,28.20188,13.46510,18.82721), tmax=c(44.23301,21.11526,22.79882,35.31366,35.84991,39.14286)); dados dados$data2 <- as.Date(dados$data,"%d/%m/%Y"); dados length(dados$tmax); mean(dados$tmax) data.cut <- as.Date("01/02/2004","%d/%m/%Y") novo <- subset(dados,data2>data.cut); novo length(novo$tmax); mean(novo$tmax) ############## Até mais, -- ================================================ É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

2 2

Re: [R-br] somar valores em intervalo pré-definido
by Eder Comunello 02 Nov '12

02 Nov '12

Lucas e colegas, bom dia! Tem que operar com as datas no formato interno do R, só transformando pra visualizar. Defini as datas de corte nas vars data.ini e data.fim... Alterei todos os valores de p1$num para '1', então se a soma estiver correta corresponderá ao número de eventos selecionados. Aqui deu certo! ;D #################### #Criando intervalo de datas Data<-seq(from=as.Date('01/01/1980',format='%d/%m/%Y'),to=as.Date('31/12/1980',format='%d/%m/%Y'),by='day') #transformando em data frame p1 <- data.frame(Data) p1$num <- 1 str(p1) data.ini <- as.Date('11/11/1980',format='%d/%m/%Y') data.fim <- as.Date('11/12/1980',format='%d/%m/%Y') data.dif <- data.fim-data.ini+1; data.dif #somando valores em um intervalo definido sum(p1$num[p1$Data >= data.ini & p1$Data <= data.fim]) #################### -- ================================================ É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

2 1

somar valores em intervalo pré-definido
by luc.souza＠usp.br 02 Nov '12

02 Nov '12

Bom dia. Estou precisando somar os valores de uma coluna em um intervalo pre-definido por mim. Por exemplo: tenho alguns valores referentes a algumas datas. Preciso somar alguns valores no período compreendio entre 11/11/1980 e 11/12/1980, conforme exemplo abaixo. Tentei fazer algumas coisas (abaixo) mas aparece o erro (tb abaixo). Só para lembrar, utilizo a data no formato americano e quando a faço a transformação, ela fica como character. Obriagdo pela ajuda Lucas #Criando intervalo de datas Data<-seq(from =as.Date('01/01/1980',format='%d/%m/%Y'),to=as.Date('31/12/1980',format='%d/%m/%Y'),by='day') #transformando em data frame p1 <- data.frame(Data) #colocando a data no formato americano que eu utilizo p1$Data <- as.Date(p1$Data, "%d/%m/%Y") p1$Data <- strftime(p1$Data, "%m/%d/%Y") #criando tabela de valores p1$num <- 1:366 #somando valores em um intervalo definido soma <- with(p1, sum(p1[p1$Data >= 11/11/1980 & p1$Data <= 11/12/1980, p1$num])) Error in `[.data.frame`(p1, p1$Data >= 11/11/1980 & p1$Data <= 11/12/1980, : undefined columns selected sum(which(p1$num[p1$Data >= '11/11/1980'] & p1$num[p1$Data <= '11/12/1980'])) erro: longer object length is not a multiple of shorter object length Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

1 0

selecionar intervalo de datas
by luc.souza＠usp.br 01 Nov '12

01 Nov '12

Boa tarde Possuo uma coluna com datas (m/d/aaaa), cujo formato é character. Preciso selecionar um intervalo entre duas datas. Já tentei subset, mas não funcionou. Alguma solução? obrigado Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

3 4

Problemas com ajuste de modelo ARIMAX
by Fatima do Nascimento Silva 01 Nov '12

01 Nov '12

Prezados, estou tentando ajustar um modelo ARIMAX no R, ao enviar as seguintes linhas de comando, aparece a seguinte mensagem de erro: "Erro em optim(init[mask], armaCSS, method = "BFGS", hessian = FALSE,: valor não finito fornecido por optim #linhas de comando (dados=read.table("H:vazao.txt",header=T,dec=","))#dados para ajuste do modelo attach(dados) (dados1=read.table("H:A.txt",header=F,dec=","))#data frame com as covariáveis #dimensão: [132,3] attach(dados1) library(forecast)#carregando pacote #Ajustando o modelo (mdl=arimax(dados, order = c(0, 2, 1), include.mean = T, xreg = dados1)) Atenciosamente, Fátima

2 4