Março 2011 - R-br

Reamostragem
by Fabio Mathias Corrêa 17 Mar '11

17 Mar '11

----- Mensagem encaminhada ---- De: Paulo Justiniano Ribeiro Jr <paulojus(a)c3sl.ufpr.br> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Sexta-feira, 31 de Outubro de 2008 8:21:23 Assunto: Re: [R_STAT] Reamostragem Andre 1. obtenha as reamostras: reamostras <- array(sample( e, length(e)*1000, rep=T), dim=c(length( e)/2,2,1000) ) 2. faca os testes T nas 1000 reamostras TT <- apply(reamostras, 3, function(x) t.test(x[,1] , x[,2])) 3. Veja o resultado em uma das amostras: TT[[1]] 4. extraia uma quantidade que quiser como por exemplo os p-valores pvals <- sapply(TT, function(x) x$p.value) OBS: note que estao sendo usadas as opcoes default do teste t para comparacao de duas amostras dadas em t.test. (bilateral, nao pareado, var. diferentes, etc) Se quiser mudar acrescente argumentos desejados na chamada do apply(): TT <- apply(reamostras, 3, function(x) t.test(x[,1] , x[,2]), var.equal=TRUE) Paulo Justiniano Ribeiro Jr LEG (Laboratorio de Estatistica e Geoinformacao) Universidade Federal do Parana Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 3361 3573 Fax: (+55) 41 3361 3141 e-mail: paulojus AT ufpr br http://www.leg. ufpr.br/~ paulojus On Fri, 31 Oct 2008, André Oliveira Souza wrote: > ola; queria fazer uma amostragem e tirar as informações. > > Olha tenho duas mostras. > > a<-c(1,2,4,5, 3,4,5,6,6, 7,2) > b<-c(5,3,4,5, 3,4,5,3,4, 3,5) > > Aplico o teste t e tenho um P-valor Po. > > Eu gostaria de juntar estes dois vetores num mesmo e fazer 1000 > reamostragem neste vetor de tamanho do vetor a e com reposição. > > e<- c(5,3,4,5,3, 4,5,3,4,3, 5,1,2,4,5, 3,4,5,6,6, 7,2) > > Depois eu queria ver aplicar o teste t(nas amostras) e ver como se > comporta estes p valores. > > att > obrigado. > > __._,_.___ Mensagens neste tópico (4) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 10 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

1 0

Material de programação em R!!!
by Fabio Mathias Corrêa 17 Mar '11

17 Mar '11

----- Mensagem encaminhada ---- De: Edimeire Alexandra Pinto <economatistica(a)yahoo.com.br> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Segunda-feira, 27 de Outubro de 2008 16:37:07 Assunto: Re: [R_STAT] Material de programação em R!!! Rafael, achei mais este material. Está anexado. Abraços ________________________________ Novos endereços, o Yahoo! que você conhece. Crie um email novo com a sua cara @ymail.com ou @rocketmail. com.__._,_.___ Mensagens neste tópico (4) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 10 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

1 0

Medidas
by Fabio Mathias Corrêa 17 Mar '11

17 Mar '11

----- Mensagem encaminhada ---- De: André Oliveira Souza <andreolisouza(a)gmail.com> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Terça-feira, 21 de Outubro de 2008 12:19:04 Assunto: [R_STAT] Re: Medidas --- Em R_STAT@yahoogrupos. com.br, "jalmarfc" <jalmarfc@.. .> escreveu >olha se ajuda !!!! library(fBasics) x<-rnorm(4444) basicStats(x) André > Óla pessoal..existe algum comando que calcula todas estas medidas > (media, variancia, assimetria e a curtose). > Obrigado pela ayuda > __._,_.___ Mensagens neste tópico (8) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 7 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

1 0

Gráfico Probability Plot no R
by Fabio Mathias Corrêa 17 Mar '11

17 Mar '11

----- Mensagem encaminhada ---- De: Marcio Nicolau <marcio.nicolau(a)gmail.com> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Terça-feira, 21 de Outubro de 2008 9:30:41 Assunto: Re: [R_STAT] Gráfico Probability Plot no R Olá Daniela, Veja a função probplot do pacote e1071 Ex x <- rnorm(100, mean=5) probplot(x) p <- probplot(x, line=FALSE) lines(p, col="green", lty=2, lwd=2) ## Median and 90% confidence interval lines(p, h=.5, col="red", lwd=3, bend=TRUE) lines(p, h=c(.05, .95), col="red", lwd=2, lty=3, bend=TRUE) Abraços 2008/10/21 Fabio Mathias <fabio.ufla@yahoo. com.br> Tente a library (car), caso vc queira construir qq-plot. >comando. qq.plot(y) > > Fábio Mathias Corrêa >Estatística e Experimentação Agropecuária > > > > >----- Mensagem original ---- >De: Daniela Recchia <daniela_recchia@ yahoo.com. br> >Para: R_STAT@yahoogrupos. com.br; stat-math <stat-math@yahoogrup os.com.br> >Enviadas: Segunda-feira, 20 de Outubro de 2008 23:17:29 >Assunto: [R_STAT] Gráfico Probability Plot no R > > > >Boa noite! > >Estou procurando um comando que contrua o gráfico Probability Plot (aquele com a >região de confiança) para dados discretos. > >Já pesquisei alguns pacotes, mas não encontrei nada. Alguém teria alguma dica? > >Obrigada! >Daniela >____________ _________ _________ _________ _________ __ >Fale com seus amigos de graça com o novo Yahoo! Messenger >http://br.messenger .yahoo.com/ > >________________________________ Novos endereços, o Yahoo! que você conhece. Crie um email novo com a sua cara @ymail.com ou @rocketmail.com. -- Márcio Nicolau __._,_.___ Mensagens neste tópico (2) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 7 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

1 0

odfWeave
by Fabio Mathias Corrêa 17 Mar '11

17 Mar '11

----- Mensagem encaminhada ---- De: Paulo Justiniano Ribeiro Jr <paulojus(a)c3sl.ufpr.br> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Segunda-feira, 20 de Outubro de 2008 20:06:33 Assunto: Re: [R_STAT] Re: odfWeave Caros Vou me permitir palpitar um pouco embora nao tenha certeza se posso ajudar (já explico porque) 1. em anexo um arquivo odfWeave (ex.odt) do tipo que o usuario tem que editar e uma saida do comando: > odfWeave('ex. odt', 'exW.odt') que gera o arquivo exW.odt que é o relatorio com as aná;k8ises incluidas. 2. tem algums confusao nas mesagens aqui na lista: ou se usa o LaTeX + Sweave **OU ** o odfWeave !! Nao precisa saber LaTeX para usar o odfWeave. 3. **Desconfio* * que o problema dos acentos está na codificação de caracteres (confusoes entre UTF-8 e ISO-8859-1) Nao sei como o exW.odt que estou enviando vai aparecxer ai para vcs, mas aqui na minha maquina (Ubuntu Linux/ e tb numa Debian Linux) aparece tudo direitinho O windows (estou assumindo que as msg sao de usuarios desta coisa esquisita e complicadissima de usar chamada windows) até onde sei nao implementa uft-8 e **talvez** isto possa explicar o problema. Deve ter como acertar estas configurações --- mas aqui nao posso ajudar mais pq nao tenho e nao sei usar este sistema operacional. Mas eu veria alguma coisas assim tentando procurar (google) por coisas como: internacionalizatio n odfWeave accents odfWeave e similares Boa sorte P.S.: De forma mais geral: Melhor solucao na minha opiniao: instale o linux, passe a usar o LaTeX e Sweave (embora seja possivel usar o odfWeave no Linux, é claro) E ... aproveite para fazer melhor uso dos 400 a 800 reais embutidos no preco de computadores que vc compra.... Paulo Justiniano Ribeiro Jr LEG (Laboratorio de Estatistica e Geoinformacao) Universidade Federal do Parana Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 3361 3573 Fax: (+55) 41 3361 3141 e-mail: paulojus AT ufpr br http://www.leg. ufpr.br/~ paulojus On Mon, 20 Oct 2008, André Oliveira Souza wrote: > --- Em R_STAT@yahoogrupos. com.br, "Guilherme Moraes Ferraudo" > <guiferraudo@ ...> escreveu > > Mas deve ter uma forma de resolver isto porque estou sem saida > eu tenho um livro sobre o R e se eu tiver que tirar todos os > acentos ....não vai dar ......e ainda depois concertar no arquivo de > saida. Alias minhas figuras tambem estão muito grande e tenho mais > 80 figuras se eu for concertar uma por uma tambem vai levar muito > tempo. > att. > obrigado. > > André eu tenho o mesmo problema que você mas não sei resolver o > problema via > > código. > > Consegui resolver rodando todo o código sem acento e, depois, no > arquivo > > gerado (output.odt) consertei o texto. > > > > Estou esperando ajuda. rs > > > > Abraços > > > > >__._,_.___ Mensagens neste tópico (8) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 8 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

1 0

Elipsóides de Confiança
by Fabio Mathias Corrêa 17 Mar '11

17 Mar '11

----- Mensagem encaminhada ---- De: walmes zeviani <walmeszeviani(a)yahoo.com.br> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Segunda-feira, 20 de Outubro de 2008 19:08:34 Assunto: Res: [R_STAT] Elipsóides de Confiança Saudações Willian, Sobre o que você relatou eu consegui a abordagem que está apresentada na seguência de comandos a seguir. #====== Contruir contornos de confiança ======# #======= para dois modelos de regressão ======# library(ellipse) #para lidar com contornos data(mtcars) #conjunto de dados fit <- lm(mpg ~ am + vs , mtcars) #modelo 1 fit2 <- lm(mpg ~ drat + vs , mtcars) #modelo 2 # Plotando os contornos para o 'modelo 1' plot(ellipse( fit, which = c('am', 'vs'), level = 0.90), type = 'l', xlim=c(2,9), ylim=c(2,11) , xlab="drat e am", ylab="vs", col=2, lwd=2, main="Contornos de confiança") par(new=T) #permite a sobreposição de gráficos # Sobrepondo os contornos do 'modelo 2' plot(ellipse( fit2, which = c('drat', 'vs'), level = 0.90), type = 'l', xlim=c(2,9), ylim=c(2,11) , xlab="", ylab="", col=3, lwd=2, main="") # Adicionando informações úteis (acessórios) points(fit$coeffici ents['am' ], fit$coefficients[ 'vs']) points(fit2$ coefficients[ 'drat'], fit2$coefficients[ 'vs']) text(7,10,"mpg~ am+vs") text(3.5,3," mpg~drat+ vs") #=========== ========= ========= ========= =======# Tomando o cuidado de delarar o mesmo xlim ylim para os dois gráficos é o que importa. Procure delarar xlab ylab apenas no primeiro gráfico. Espero ter ajudado você. Uma ótima semana. Walmes. ----- Mensagem original ---- De: William Marciano <fwpm2004@yahoo. com.br> Para: R_STAT@yahoogrupos. com.br Enviadas: Segunda-feira, 20 de Outubro de 2008 16:54:34 Assunto: [R_STAT] Elipsóides de Confiança Por favor, gostaria de saber se alguém já construiu em um mesmo gráfico elipsóides de confiança para modelos de regressão padrão diferentes, já tentei com a função confidence.ellipse( ), mas não consegui colocar em um mesmo gráfico dois modelos que eu ajustei para compará-los. Desde já agradeço a atenção de todos! ____________ _________ _________ _________ _________ __ Fale com seus amigos de graça com o novo Yahoo! Messenger http://br.messenger .yahoo.com/ ____________ _________ _________ _________ _________ __ Fale com seus amigos de graça com o novo Yahoo! Messenger http://br.messenger .yahoo.com/ __._,_.___ Mensagens neste tópico (1) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 8 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

1 0

Nova Lista
by Andre Korzenowski 17 Mar '11

17 Mar '11

Parabéns pela iniciativa de tornar mais dinâmico o processo da lista. Ela já está plenamente funcional? E quanto ao histórico, conseguirá ser importado? Um abraço. -- Ms. André Luis Korzenowski Curriculum Lattes<http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4792588H3> - Página Web <http://www.korzenowski.com/>

1 0

Soluções interessantes
by Fabio Mathias Corrêa 17 Mar '11

17 Mar '11

Pessoal, Tenho guardado em uma pasta de email diversas soluções para os mais variados problemas. Logo estarei postando as soluções para a nova lista. Assim evitamos muitas perguntas cuja a solução já foi postada! Quem quiser fazer o mesmo. Coloque o nome da mensagem no assunto e no corpo da mensagem o script de ajuda! Valeu!!! Fábio Mathias Corrêa Departamento de Estatística Universidade Estadual de Santa Cruz Tel.: 73-3680-5076 Cel.: 73-9991-8155

1 0

Testando a bagaça 2!!!
by Fabio Mathias Corrêa 16 Mar '11

16 Mar '11

1 + 1 = 11 (p < 0.05)***** Valeu!!! Fábio Mathias Corrêa Departamento de Estatística Universidade Estadual de Santa Cruz Tel.: 73-3680-5076 Cel.: 73-9991-8155

1 0

Parabéns pela iniciativa
by Guilherme Moraes Ferraudo 16 Mar '11

16 Mar '11

Boa noite a todos, gostaria de externar os meus sinceros agradecimentos aos mantenedores e criadores dessa lista. Isto é mais um passo importante na consolidação do R no Brasil. Abraços -- Guilherme Moraes Ferraudo Jaboticabal/Campinas - SP http://lattes.cnpq.br/2096118558794430

2 1