- R-br - DINF UFPR

Teste
by Fabiana Goncalves Barbosa 07 Dez '14

07 Dez '14

Tenho os seguintes dados indep inter 288 1887 187 2195 377 111 191 154 165 379 350 156 189 153 100 250 173 186 176 122 431 108 249 1841 213 426 314 252 122 184 635 138 168 169 151 166 154 117 109 141 100 313 1030 125 350 116 329 726 253 420 218 245 123 111 108 135 313 124 275 118 192 100 197 246 222 247 103 567 557 558 299 211 102 184 116 180 119 609 259 155 151 634 343 472 299 454 242 381 167 249 162 248 145 124 126 169 125 151 157 156 208 122 161 106 139 104 101 434 NA 190 NA 171 NA 122 NA 117 NA 2016 NA 349 NA 309 NA 271 NA 255 NA 259 NA 239 NA 158 NA 199 NA 200 NA 168 NA 1368 NA 622 NA 530 NA 373 NA 271 NA 121 NA 188 NA 123 NA 194 NA 125 NA 1811 NA 403 NA 322 NA 229 NA 260 NA 215 NA 151 NA 129 NA 135 NA 104 NA 543 NA 460 NA 286 NA 265 NA 248 NA 164 NA 130 NA 125 NA 110 NA 192 NA 122 NA 103 NA 361 NA 314 NA 262 NA 193 NA 155 NA 513 NA 301 NA 243 NA 177 NA 155 NA 188 NA 193 NA 122 NA 108 NA 109 NA 191 NA 324 NA 163 Que teste uso para verificar se existe diferença entre as duas amostras?

3 4

extracao da Longitude correta para dados de SST
by Humberto 06 Dez '14

06 Dez '14

Boa noite pessoal, Estou extraindo de um arquivo ncdf os valores de latitude (lat) longitude( lon) e temperatura da superfície da agua do mar( SST), entretanto a longitude (lon) ela está na projeção 0 a 360 graus. Minha pergunta é como faço para converter ela para -180 a 180 graus? Segue o script que estou utilizando juntamente com o link dos dados. Desde já agradeço a ajuda Humberto lon lat data sst Lon como deveria ser 299.5 10.5 11/30/1981 27.85 -60.000 300.5 10.5 11/30/1981 27.96 -59.500 301.5 10.5 11/30/1981 28.02 -58.500 302.5 10.5 11/30/1981 28.02 -57.500 303.5 10.5 11/30/1981 27.97 -56.500 304.5 10.5 11/30/1981 27.84 -55.500 305.5 10.5 11/30/1981 27.81 -54.500 306.5 10.5 11/30/1981 27.78 -53.500 307.5 10.5 11/30/1981 27.75 -52.500 Link para acesso dos dados https://www.dropbox.com/s/xnab2g365mduykx/SST_81_14.nc?dl=0 Script: require(ncdf) #'ncdf' versao 1.6.6 require(reshape) require(chron) ###sst nc <- open.ncdf("SST_81_14.nc") nc lat <- get.var.ncdf(nc, "lat") lon <- get.var.ncdf(nc, "lon") tunits <- att.get.ncdf(nc,"time","units") # print tunits tunits$value # criar tempo "real" Times time<- get.var.ncdf(nc,"time") Time<-chron(time,origin=c(month=1, day=1, year=1900)) ## array da variavel sst wm <- get.var.ncdf(nc, "sst") str(wm) ## atribuindo dimnames conforme lon, lat e tempo dimnames(wm)[[1]] <- as.character(lon) dimnames(wm)[[2]] <- as.character(lat) dimnames(wm)[[3]] <- as.character(Time) str(wm) ## melt para reestruturação dos dados mwm <- melt(wm) head(mwm)

6 6

RCurl: passando parâmetros e executando função javascript
by Éder Comunello 04 Dez '14

04 Dez '14

Colegas, boa tarde! Estou precisando de ajuda com relação ao pacote *RCurl*. Tomando por exemplo o endereço < http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001>, ao acessá-lo via *browser *pode-se notar dois botões que executam trechos de código em *javascript *e geram links para o *download *de arquivos de dados: * Botão [Arquivo Access] remete ao comando "javascript:criarArq(2254001,2)" * Botão [Arquivo Texto] remete ao comando "javascript:criarArq(2254001,1)" Estou tentando usar o *RCurl *para passar diretamente os parâmetros para criação do *link *para o "Arquivo Texto". Isto feito, posteriormente pretendo usar o *link *gerado para *download* dos dados. Contudo, não estou tendo sucesso e ficarei muito grato por qualquer ajuda. Abaixo segue o código que tentei aplicar: ### <BEGIN> ### library(RCurl) myURL <- 'http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001' ### URL resultante da operação via browser: ### myURL2 <- ' http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001&CriaArq=true&TipoArq=1 ' myCurl <- getCurlHandle(cookiefile = "myCookies.txt", useragent="Mozilla/5.0", verbose = T, autoreferer = T, followlocation = T) myParams <- list(Codigo=2254001, CriaArq='true', TipoArq=1) ### Codigo: ID da estação ### CriaArq: true (executa função?), ### TipoArq: 1 (tipo texto) e 2 (Access) postForm(myURL, .params=myParams, curl=myCurl) getURL(myURL, curl=myCurl) ### Operação falha, retornando página original! ### Em caso de sucesso embutiria no código um link no formato: ### <a href="ARQ/A20131008-164720-14/CHUVAS.ZIP"> ### <END> ### -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

4 3

Anisotropia, verificação gráfica
by Felinto COSTA 04 Dez '14

04 Dez '14

Prezados. Para o traçado da elipse ("rose diagram") tal como sugerido por Isaaks e Srivastava devo usar os alcances com cor=0.05 (practical ranges) para as diversas direções estudadas ou os alcances (phi) dos patamares (sigmasq) calculados por variofit ou likft (GeoR)? Antecipadamente agradeço a atenção dos colegas. FCosta

2 2

Formatação de várias variáveis
by Alessandro Corrêa 03 Dez '14

03 Dez '14

Prezados colegas, Estou com um banco de dados (x) no seguinte formato sexo idade tempo Q1 Q2 Q3 Q4 Q5 Q6 ...Q63 1 44 24.0 4 3 3 3 2 4 0 62 28.0 1 5 2 2 3 1 1 52 6.0 4 5 4 2 5 4 0 51 27.0 4 4 4 3 3 3 1 32 3.5 1 3 2 1 4 3 0 33 5.0 1 5 1 1 1 1 As variáveis de Q1 a Q6 são variáveis na escala de Likert1 = Discordo totalmente2 = Discordo 3 = Neutro 4 = Concordo 5 = Concordo totalmente Gostaria de saber se é possível criar uma rotina ou se há uma função capaz deformatar as variáveis de Q1 a Q63 em fatores com as etiquetas referentes a cada nível. Se eu fosse fazer uma a uma, repetiria 63 vezes o seguinte comando para cada variável: x$Q1<-factor(x$Q1, label=c("Discordo totalmente", "Discordo", "Neutro", "Concordo", "Concordo totalmente"), levels=1:5) Acredito que haja uma maneira mais prática. Obrigado desde já. Alessandro

4 7

Numero de dias acima de um valor
by Jefferson Ferreira-Ferreira 02 Dez '14

02 Dez '14

Olá amigos!! Estou com uma questão que não estou sabendo como resolver. Eu tenho uma série de dados de nivel d'água (3 vezes por dia). O que eu quero saber é em qual data começa a enchente. Como o nível sofre algumas "pequenas enchentes" que de fato não representam o real início da enchente, o que eu preciso é saber em qual data começa o seguinte evento: nível > 2cm por 45 dias consecutivos. Isso porque se o nível permanecer por menos de 45 dias consecutivos acima de 2cm eu não estou considerando isso como o início da enchente, mas como uma subida intermitente da água. Outro complicador é que preciso saber essa informação para cada um dos 10 equipamento (datalogger), discriminado na culona que chamei de logger. Ou seja, a data (em dia juliano) do início da enchente para cada logger. Em anexo envio os dados de apenas um dos meus aparelhos (logger == 1) Alguma ideia?

4 8

salvando codigos carregados
by Hélio Gallo Rocha 02 Dez '14

02 Dez '14

Bom dia a todos da lista. Vendo alguns posts vi um sobre semente, onde o listeiro Alejandro onde salva uma simulação e reinicia a partir de um ponto com .Random.seed Haveria algo parecido em se salvar relatórios e dados gerados a partir de códigos rodados numa seção e abrir a partir daquele ponto em outra seção? Agradeço desde já. -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Re: [R-br] Digest R-br, volume 48, assunto 3
by Fabiana Goncalves Barbosa 02 Dez '14

02 Dez '14

Olá Cesar A variável "total" é uma variável de contagem. Em 2 de dezembro de 2014 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Teste (David Feitosa) > 2. Re: Teste (Cesar Rabak) > 3. Re: Teste (Cesar Rabak) > 4. Re: Teste (Cesar Rabak) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 1 Dec 2014 20:18:17 -0300 > From: David Feitosa <contato(a)davidfeitosa.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > < > CAL1rLto7Dn-1HFEmxx70vh56D6n96YBdR-+1oTAyhfxf8Pop9Q(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Diferença em que sentido? > > > > Atenciosamente, > > David Feitosa > > (\_(\ > (=°;°) > (("")("") > > Em 1 de dezembro de 2014 20:12, Fabiana Goncalves Barbosa < > fabibarbos(a)gmail.com> escreveu: > > > Tenho os seguintes dados > > > > indep inter > > 288 1887 > > 187 2195 > > 377 111 > > 191 154 > > 165 379 > > 350 156 > > 189 153 > > 100 250 > > 173 186 > > 176 122 > > 431 108 > > 249 1841 > > 213 426 > > 314 252 > > 122 184 > > 635 138 > > 168 169 > > 151 166 > > 154 117 > > 109 141 > > 100 313 > > 1030 125 > > 350 116 > > 329 726 > > 253 420 > > 218 245 > > 123 111 > > 108 135 > > 313 124 > > 275 118 > > 192 100 > > 197 246 > > 222 247 > > 103 567 > > 557 558 > > 299 211 > > 102 184 > > 116 180 > > 119 609 > > 259 155 > > 151 634 > > 343 472 > > 299 454 > > 242 381 > > 167 249 > > 162 248 > > 145 124 > > 126 169 > > 125 151 > > 157 156 > > 208 122 > > 161 106 > > 139 104 > > 101 434 > > NA 190 > > NA 171 > > NA 122 > > NA 117 > > NA 2016 > > NA 349 > > NA 309 > > NA 271 > > NA 255 > > NA 259 > > NA 239 > > NA 158 > > NA 199 > > NA 200 > > NA 168 > > NA 1368 > > NA 622 > > NA 530 > > NA 373 > > NA 271 > > NA 121 > > NA 188 > > NA 123 > > NA 194 > > NA 125 > > NA 1811 > > NA 403 > > NA 322 > > NA 229 > > NA 260 > > NA 215 > > NA 151 > > NA 129 > > NA 135 > > NA 104 > > NA 543 > > NA 460 > > NA 286 > > NA 265 > > NA 248 > > NA 164 > > NA 130 > > NA 125 > > NA 110 > > NA 192 > > NA 122 > > NA 103 > > NA 361 > > NA 314 > > NA 262 > > NA 193 > > NA 155 > > NA 513 > > NA 301 > > NA 243 > > NA 177 > > NA 155 > > NA 188 > > NA 193 > > NA 122 > > NA 108 > > NA 109 > > NA 191 > > NA 324 > > NA 163 > > > > Que teste uso para verificar se existe diferença entre as duas amostras? > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/98b0d5a3/atta… > > > > ------------------------------ > > Message: 2 > Date: Mon, 1 Dec 2014 21:20:40 -0200 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > <CAKrF98=_+f4DZ1RwUKn8GxHhHgu6nzdfst6cFUwB-= > L8h0osdg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Fabiana, > > O teste estatístico mais apropriado para seus dados depende de uma série de > premissas que se põe a respeito de como eles são produzidos. > > A variável "total" é uma variável de contagem ou é uma variável de outra > natureza? > > Comentando as respostas que você já obteve, afirmo que com os dados [que > você mostrou] sem nenhuma transformação um teste t (ou mesmo o F como > propõe o Fernando) seriam uma "loteria" pois esses testes são para detetar > diferenças nas médias e essa estatística é muito afetada por valores > extremos (cuja distribuição dos seus dados mostra que ela é assimétrica). > > O inst que tenha a amostra com o maior valor "levaria" vantagem pois o > máximo é 7,5 vezes a média enquanto o menor valor é apenas 0,34 vezes ela. > > Recomendo que você faça uma análise exploratória dos seus dados e > identifique face ao processo que os geram qual seria a melhor forma de > comprovar a hipótese experimental que você que verificar com o teste > estatístico. > > meus 0,019999. > > -- > Cesar Rabak > > > 2014-12-01 16:03 GMT-02:00 Fernando Antonio de souza < > nandodesouza(a)gmail.com > >: > > > Teste t ou F. > > Em 01/12/2014 13:14, "Luciano F. Sgarbi" <luciano.f.sgarbi(a)gmail.com> > > escreveu: > > > > Fabiana. > >> Não seria o caso de um teste t? > >> Abraço > >> Em 01/12/2014 11:27, "Fabiana Goncalves Barbosa" <fabibarbos(a)gmail.com> > >> escreveu: > >> > >>> Olá Pessoal > >>> > >>> Tenho os seguintes dados: inst = variavel independente e total = > >>> variavel dependente. Sendo que, 54 valores são para "indep" e 119 para > >>> "inter". > >>> > >>> Quero verificar se há diferença significativa no "total" nos dois > grupos > >>> "indep" e "inter". Que teste usar? > >>> > >>> Segue os dados: > >>> > >>> inst total indep 288 indep 187 indep 377 indep 191 indep 165 > >>> indep 350 indep 189 indep 100 indep 173 indep 176 indep 431 indep > >>> 249 indep 213 indep 314 indep 122 indep 635 indep 168 indep 151 > >>> indep 154 indep 109 indep 100 indep 1030 indep 350 indep 329 > indep > >>> 253 indep 218 indep 123 indep 108 indep 313 indep 275 indep 192 > >>> indep 197 indep 222 indep 103 indep 557 indep 299 indep 102 indep > >>> 116 indep 119 indep 259 indep 151 indep 343 indep 299 indep 242 > >>> indep 167 indep 162 indep 145 indep 126 indep 125 indep 157 indep > >>> 208 indep 161 indep 139 indep 101 inter 1887 inter 2195 inter 111 > >>> inter 154 inter 379 inter 156 inter 153 inter 250 inter 186 inter > >>> 122 inter 108 inter 1841 inter 426 inter 252 inter 184 inter 138 > >>> inter 169 inter 166 inter 117 inter 141 inter 313 inter 125 inter > >>> 116 inter 726 inter 420 inter 245 inter 111 inter 135 inter 124 > >>> inter 118 inter 100 inter 246 inter 247 inter 567 inter 558 inter > >>> 211 inter 184 inter 180 inter 609 inter 155 inter 634 inter 472 > >>> inter 454 inter 381 inter 249 inter 248 inter 124 inter 169 inter > >>> 151 inter 156 inter 122 inter 106 inter 104 inter 434 inter 190 > >>> inter 171 inter 122 inter 117 inter 2016 inter 349 inter 309 > inter > >>> 271 inter 255 inter 259 inter 239 inter 158 inter 199 inter 200 > >>> inter 168 inter 1368 inter 622 inter 530 inter 373 inter 271 > inter > >>> 121 inter 188 inter 123 inter 194 inter 125 inter 1811 inter 403 > >>> inter 322 inter 229 inter 260 inter 215 inter 151 inter 129 inter > >>> 135 inter 104 inter 543 inter 460 inter 286 inter 265 inter 248 > >>> inter 164 inter 130 inter 125 inter 110 inter 192 inter 122 inter > >>> 103 inter 361 inter 314 inter 262 inter 193 inter 155 inter 513 > >>> inter 301 inter 243 inter 177 inter 155 inter 188 inter 193 inter > >>> 122 inter 108 inter 109 inter 191 inter 324 inter 163 > >>> Grata > >>> > >>> > >>> > >>> > >>> _______________________________________________ > >>> 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/f257fc3c/atta… > > > > ------------------------------ > > Message: 3 > Date: Mon, 1 Dec 2014 21:28:48 -0200 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > <CAKrF98nABg0== > h+mHEX14i9xTWecm8kTAdhj+qC+-QaNResJtw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > David, > > Via de regra considera-se mais "seguro" usar testes bicaudais para > verificações entre tratamentos para a resposta a uma variável numérica > (como parece ser o caso aqui) por isso o "sentido" da diferença não deveria > ser uma preocupação do pesquisador (a menos que haja uma razão dada pela > natureza do problema muito grande para descartar um dos lados da cauda). > > []s > > > 2014-12-01 21:18 GMT-02:00 David Feitosa <contato(a)davidfeitosa.com>: > > > Diferença em que sentido? > > > > > > > > Atenciosamente, > > > > David Feitosa > > > > (\_(\ > > (=°;°) > > (("")("") > > > > Em 1 de dezembro de 2014 20:12, Fabiana Goncalves Barbosa < > > fabibarbos(a)gmail.com> escreveu: > > > >> Tenho os seguintes dados > >> > >> indep inter > >> 288 1887 > >> 187 2195 > >> 377 111 > >> 191 154 > >> 165 379 > >> 350 156 > >> 189 153 > >> 100 250 > >> 173 186 > >> 176 122 > >> 431 108 > >> 249 1841 > >> 213 426 > >> 314 252 > >> 122 184 > >> 635 138 > >> 168 169 > >> 151 166 > >> 154 117 > >> 109 141 > >> 100 313 > >> 1030 125 > >> 350 116 > >> 329 726 > >> 253 420 > >> 218 245 > >> 123 111 > >> 108 135 > >> 313 124 > >> 275 118 > >> 192 100 > >> 197 246 > >> 222 247 > >> 103 567 > >> 557 558 > >> 299 211 > >> 102 184 > >> 116 180 > >> 119 609 > >> 259 155 > >> 151 634 > >> 343 472 > >> 299 454 > >> 242 381 > >> 167 249 > >> 162 248 > >> 145 124 > >> 126 169 > >> 125 151 > >> 157 156 > >> 208 122 > >> 161 106 > >> 139 104 > >> 101 434 > >> NA 190 > >> NA 171 > >> NA 122 > >> NA 117 > >> NA 2016 > >> NA 349 > >> NA 309 > >> NA 271 > >> NA 255 > >> NA 259 > >> NA 239 > >> NA 158 > >> NA 199 > >> NA 200 > >> NA 168 > >> NA 1368 > >> NA 622 > >> NA 530 > >> NA 373 > >> NA 271 > >> NA 121 > >> NA 188 > >> NA 123 > >> NA 194 > >> NA 125 > >> NA 1811 > >> NA 403 > >> NA 322 > >> NA 229 > >> NA 260 > >> NA 215 > >> NA 151 > >> NA 129 > >> NA 135 > >> NA 104 > >> NA 543 > >> NA 460 > >> NA 286 > >> NA 265 > >> NA 248 > >> NA 164 > >> NA 130 > >> NA 125 > >> NA 110 > >> NA 192 > >> NA 122 > >> NA 103 > >> NA 361 > >> NA 314 > >> NA 262 > >> NA 193 > >> NA 155 > >> NA 513 > >> NA 301 > >> NA 243 > >> NA 177 > >> NA 155 > >> NA 188 > >> NA 193 > >> NA 122 > >> NA 108 > >> NA 109 > >> NA 191 > >> NA 324 > >> NA 163 > >> > >> Que teste uso para verificar se existe diferença entre as duas amostras? > >> > >> _______________________________________________ > >> 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/d73c49f0/atta… > > > > ------------------------------ > > Message: 4 > Date: Mon, 1 Dec 2014 21:29:15 -0200 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > < > CAKrF98nTx_GZraqEdnUuv+ZS0RL2gA0H9Hmz5yZ8B3FdE18iyg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Fabiana, > > Os dados desse post são os mesmos do post original, apenas a forma de > organizá-los está diferente. A teoria para testá-los é a mesma. > > []s > > > 2014-12-01 21:12 GMT-02:00 Fabiana Goncalves Barbosa <fabibarbos(a)gmail.com > >: > > > Tenho os seguintes dados > > > > indep inter > > 288 1887 > > 187 2195 > > 377 111 > > 191 154 > > 165 379 > > 350 156 > > 189 153 > > 100 250 > > 173 186 > > 176 122 > > 431 108 > > 249 1841 > > 213 426 > > 314 252 > > 122 184 > > 635 138 > > 168 169 > > 151 166 > > 154 117 > > 109 141 > > 100 313 > > 1030 125 > > 350 116 > > 329 726 > > 253 420 > > 218 245 > > 123 111 > > 108 135 > > 313 124 > > 275 118 > > 192 100 > > 197 246 > > 222 247 > > 103 567 > > 557 558 > > 299 211 > > 102 184 > > 116 180 > > 119 609 > > 259 155 > > 151 634 > > 343 472 > > 299 454 > > 242 381 > > 167 249 > > 162 248 > > 145 124 > > 126 169 > > 125 151 > > 157 156 > > 208 122 > > 161 106 > > 139 104 > > 101 434 > > NA 190 > > NA 171 > > NA 122 > > NA 117 > > NA 2016 > > NA 349 > > NA 309 > > NA 271 > > NA 255 > > NA 259 > > NA 239 > > NA 158 > > NA 199 > > NA 200 > > NA 168 > > NA 1368 > > NA 622 > > NA 530 > > NA 373 > > NA 271 > > NA 121 > > NA 188 > > NA 123 > > NA 194 > > NA 125 > > NA 1811 > > NA 403 > > NA 322 > > NA 229 > > NA 260 > > NA 215 > > NA 151 > > NA 129 > > NA 135 > > NA 104 > > NA 543 > > NA 460 > > NA 286 > > NA 265 > > NA 248 > > NA 164 > > NA 130 > > NA 125 > > NA 110 > > NA 192 > > NA 122 > > NA 103 > > NA 361 > > NA 314 > > NA 262 > > NA 193 > > NA 155 > > NA 513 > > NA 301 > > NA 243 > > NA 177 > > NA 155 > > NA 188 > > NA 193 > > NA 122 > > NA 108 > > NA 109 > > NA 191 > > NA 324 > > NA 163 > > > > Que teste uso para verificar se existe diferença entre as duas amostras? > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/8116ae2c/atta… > > > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 48, assunto 3 > **************************************** >

1 0

Teste
by Fabiana Goncalves Barbosa 01 Dez '14

01 Dez '14

Olá Pessoal Tenho os seguintes dados: inst = variavel independente e total = variavel dependente. Sendo que, 54 valores são para "indep" e 119 para "inter". Quero verificar se há diferença significativa no "total" nos dois grupos "indep" e "inter". Que teste usar? Segue os dados: inst total indep 288 indep 187 indep 377 indep 191 indep 165 indep 350 indep 189 indep 100 indep 173 indep 176 indep 431 indep 249 indep 213 indep 314 indep 122 indep 635 indep 168 indep 151 indep 154 indep 109 indep 100 indep 1030 indep 350 indep 329 indep 253 indep 218 indep 123 indep 108 indep 313 indep 275 indep 192 indep 197 indep 222 indep 103 indep 557 indep 299 indep 102 indep 116 indep 119 indep 259 indep 151 indep 343 indep 299 indep 242 indep 167 indep 162 indep 145 indep 126 indep 125 indep 157 indep 208 indep 161 indep 139 indep 101 inter 1887 inter 2195 inter 111 inter 154 inter 379 inter 156 inter 153 inter 250 inter 186 inter 122 inter 108 inter 1841 inter 426 inter 252 inter 184 inter 138 inter 169 inter 166 inter 117 inter 141 inter 313 inter 125 inter 116 inter 726 inter 420 inter 245 inter 111 inter 135 inter 124 inter 118 inter 100 inter 246 inter 247 inter 567 inter 558 inter 211 inter 184 inter 180 inter 609 inter 155 inter 634 inter 472 inter 454 inter 381 inter 249 inter 248 inter 124 inter 169 inter 151 inter 156 inter 122 inter 106 inter 104 inter 434 inter 190 inter 171 inter 122 inter 117 inter 2016 inter 349 inter 309 inter 271 inter 255 inter 259 inter 239 inter 158 inter 199 inter 200 inter 168 inter 1368 inter 622 inter 530 inter 373 inter 271 inter 121 inter 188 inter 123 inter 194 inter 125 inter 1811 inter 403 inter 322 inter 229 inter 260 inter 215 inter 151 inter 129 inter 135 inter 104 inter 543 inter 460 inter 286 inter 265 inter 248 inter 164 inter 130 inter 125 inter 110 inter 192 inter 122 inter 103 inter 361 inter 314 inter 262 inter 193 inter 155 inter 513 inter 301 inter 243 inter 177 inter 155 inter 188 inter 193 inter 122 inter 108 inter 109 inter 191 inter 324 inter 163 Grata

4 3

Raiz Quadrada do Erro Médio Quadrático e Teste de Médias
by Helder Gramacho 30 Nov '14

30 Nov '14

Olá Pessoal, boa noite! A seguir tem as explicações do problema e o código para reproduzir os resultados, procurei detalhar bem para facilitar a vida de quem se dispõe a ajudar. No trabalho disponível no link abaixo: http://www.lsie.unb.br/rbc/index.php/rbc/article/download/433/428 os autores validam o resultado de uma interpolação utilizando pontos de controle e calculando a Raiz Quadrada do Erro Médio Quadrático (REMQ) dado pela seguinte equação: [image: Imagem inline 2] Onde VMLTC é o Valor do Metro Linear de Testada Corrigida. Repeti o procedimento para 4 diferentes superfícies de interpolação e calculei o REMQ para cada uma como mostrado abaixo: #Download do arquivo no dropbox links <- c("https://www.dropbox.com/s/1grlssvxpihtcmt/valid.csv") tokens <- gsub("^.*/s/","",dirname(links)) fileNames <- basename(links) newLinks <- file.path("http://dl.dropbox.com/s", tokens, fileNames); newLinks for (a in newLinks) { tryCatch(download.file(a, dest=basename(a), mode='wb'), error=function(...) print("Falha no download!"))} valid<-read.table(file="valid.csv",sep=",",header=T,dec=".") # Função que calcula o REMQ (ve=valor estimado vr=valor real) remq<-function(ve,vr) { r<-ve-vr r2<-r^2 soma<-sum(r2) n<-length(r2) remq<-sqrt((1/(n-1))*soma) result<-c("Raiz do Erro Médio Quadrático"=remq) return(result) } # Calculando o REMQ remq.orig<-remq(valid[,2],valid[,1]) remq.mod.1.1<-remq(valid[,3],valid[,1]) remq.mod.2.1<-remq(valid[,4],valid[,1]) remq.mod.2.1.sar<-remq(valid[,5],valid[,1]) resumo.remq<-rbind(orig=remq.orig,mod.1.1=remq.mod.1.1,mod.2.1=remq.mod.2.1,mod.sar=remq.mod.2.1.sar);resumo.remq Agora gostaria de saber se as diferenças entre os REMQ são significativas estatisticamente e é aqui que está a minha dúvida, gostaria que avaliassem se o procedimento abaixo está correto, utilizo as diferenças entre os pontos de controle e os pontos obtidos nas superfícies para realizar a análise de variância e em seguida o Teste de Tukey. # Preparando para Teste de Médias m1<-valid[,2]-valid[,1] m2<-valid[,3]-valid[,1] m3<-valid[,4]-valid[,1] m4<-valid[,5]-valid[,1] dados.mod<-cbind(m1,m2,m3,m4);dados.mod dados.org <-data.frame(dif.mod=c(dados.mod[,1],dados.mod[,2],dados.mod[,3],dados.mod[,4]),Tipo=factor(c(rep("m1",20),rep("m2",20),rep("m3",20),rep("m4",20))),observ=factor(rep(1:20,4))) # Realizando a análise de variância tapply(dados.org$dif.mod,dados.org$Tipo,mean) ajuste<-aov(dados.org$dif.mod~dados.org$Tipo+dados.org$observ) summary(ajuste) posthoc<-TukeyHSD(x=ajuste,'dados.org$Tipo',conf.level=0.95);posthoc E então? Está correto utilizar as diferenças entre os pontos de controle e os pontos obtidos na superfície para realizar o teste de médias? Observe que os REMQ são bem diferentes, mas no teste médias as diferenças não são significativas, é aí que está a dúvida. Não sei se da forma como realizei o teste de médias está realmente avaliando as diferenças entre os REMQ. Desde já agradeço toda ajuda, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 2