Julho 2013 - R-br - DINF UFPR

comando ftable
by geovane barbosa 02 Jul '13

02 Jul '13

Olá pessoa estou executando o seguinte comando. ftable(mytable), como eu faço para ele me retornar os percentuais? abraços

6 8

Comparações múltiplas em modelos lineares generalizados
by Gilenio Borges Fernandes 02 Jul '13

02 Jul '13

Prezados(as) Comparações múltiplas (inferência simultânea) feitas através de funções paramétricas com o pacote multcomp do R são válidas para modelos lineares ou para o preditor linear dos modelos lineares generalizados. No caso de modelo de Poisson (ligacao log), modelo logit, etc., os resultados (conclusões) nos preditores lineares podem ser estendidos para a resposta média? Veja o exemplo seguinte do modelo logit com 4 tratamentos e uma covariável onde fiz os intervalos simultâneos na resposta média: Trat = rep(gl(4, 4, label = c('T1','T2', 'T3', 'T4')),3) y=c(1,0,1,1,0,1,0,0,1,1,0,1,1,0,0,0,0,0,0,1,1,1,0,0,0,1,1,1,1,0,0,1,0,0,0,1,1,1,0,0,1,1,1, 1,0,0,1,0) y_ant=c(0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,1,1,0,1,0,0,1,1,1,1,0,0,0,0,0,0,1,1,1,0,0,1, 1,1,1,0,0,1,0) d1=data.frame(Trat,y_ant,y) attach(d1) gmod1 <- glm(y ~ Trat+y_ant, data = d1,family = binomial()) summary(gmod1) # # Comparacao de pares de tratamentos (no preditor linear) library(multcomp) comp1=glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE)) summary(comp1, test = adjusted("Westfall")) # Wastfall (1997) # # Teste e IC's para diferencas de pares de tratamentos no preditor linear gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) summary(gmod1_ci) gmod1_ci$confint # # Comparacao de pares de tratamentos (nas prob condicionais ajustadas) gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) gmod1_ci$confint <- apply(gmod1_ci$confint, 1, binomial()$linkinv) gmod1_ci$confint # gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) gmod1_ci$confint <- apply(gmod1_ci$confint, 2, binomial()$linkinv) gmod1_ci$confint Desde já Grato, -- Gilenio Borges Fernandes Professor Associado III Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

2 1

Pacote(mclust)
by Carlos Pereira 02 Jul '13

02 Jul '13

estou tendo problema com a função densityMclust(data, ...) do pacote(mclust)o objeto data pode ser vetor ou matriz, porem quando uso um vetor diferente do exemplo do pacote dar erro. Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

2 1

Dúvida sobre stepwise
by Victor Eduardo 02 Jul '13

02 Jul '13

Tem como fazer stepwise com teste de significância? Parece que no R, o stepwise faz pelo critério do AIC

2 2

trabalhar com médias
by alanarocha＠sapo.pt 02 Jul '13

02 Jul '13

Olá, eu tenho estes dados de 4 grupos: A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 112 277 317 1153 126 208 454 1199 210 475 2794 257 444 2081 246 493 1576 264 599 901 589 871 583 722 599 1110 548 1412 454 732 382 601 533 1130 335 642 598 682 448 1035 355 1203 462 1136 500 1109 354 805 369 912 547 1454 842 2384 760 675 810 1215 1565 1585 942 827 1015 631 647 709 2192 1392 1028 1201 892 1287 1579 1071 669 1103 964 926 1501 979 1178 625 678 1054 855 983 1076 876 que codigos eu posso fazer no R para responder a estas perguntas: 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. Ana

2 1

trabalhar com médias
by alanarocha＠sapo.pt 02 Jul '13

02 Jul '13

Olá, eu tenho estes dados de 4 grupos: A. pthi <150 B. pthi 150-300 C. pthi 300-600 D. >600 112 277 317 1153 126 208 454 1199 210 475 2794 257 444 2081 246 493 1576 264 599 901 589 871 583 722 599 1110 548 1412 454 732 382 601 533 1130 335 642 598 682 448 1035 355 1203 462 1136 500 1109 354 805 369 912 547 1454 842 2384 760 675 810 1215 1565 1585 942 827 1015 631 647 709 2192 1392 1028 1201 892 1287 1579 1071 669 1103 964 926 1501 979 1178 625 678 1054 855 983 1076 876 que codigos eu posso fazer no R para responder a estas perguntas: 7. Obter a Dose média em cada grupo e o gráfico da sua evolução ao longo do tempo. Ana

1 0

diferenca em horas
by Edson Lira 02 Jul '13

02 Jul '13

Caros amigos da Lista, tenho um banco com as variáveis abaixo: Dia_da_Semana Dia_do_Mês Entrada Saída Semana No link abaixo está um banco teste para vocês analisarem. http://www1.datafilehost.com/d/320c3ab5 Quero calcular a diferença de tempo entre a entrada e a saída. Para poder analisar tempo médio por semana. Não estou conseguindo implementar nos pacotes chron e lubridate. Alguém poderia me ajudar? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

Instalação rmysql ubuntu
by Leandro Marino 01 Jul '13

01 Jul '13

Caros, executei a instalação do Rmysql no ubuntu sem sucesso. já executei o primeiro passo das instruções mas nao sei como fazer o segundo passo... alguma sugestão? Configuration error: could not find the MySQL installation include and/or library directories. Manually specify the location of the MySQL libraries and the header files and re-run R CMD INSTALL. INSTRUCTIONS: 1. Define and export the 2 shell variables PKG_CPPFLAGS and PKG_LIBS to include the directory for header files (*.h) and libraries, for example (using Bourne shell syntax): export PKG_CPPFLAGS="-I<MySQL-include-dir>" export PKG_LIBS="-L<MySQL-lib-dir> -lmysqlclient" Re-run the R INSTALL command: R CMD INSTALL RMySQL_<version>.tar.gz 2. Alternatively, you may pass the configure arguments --with-mysql-dir=<base-dir> (distribution directory) or --with-mysql-inc=<base-inc> (where MySQL header files reside) --with-mysql-lib=<base-lib> (where MySQL libraries reside) in the call to R INSTALL --configure-args='...' R CMD INSTALL --configure-args='--with-mysql-dir=DIR' RMySQL_<version>.tar.gz ERROR: configuration failed for package ‘RMySQL’ * removing ‘/usr/local/lib/R/site-library/RMySQL’ The downloaded source packages are in ‘/tmp/RtmpIdWQKu/downloaded_packages’ Warning message: In install.packages("RMySQL") : installation of package ‘RMySQL’ had non-zero exit status a instalação do my sql se deu com o comando: *sudo apt-get install mysql-server apache2 libapache2-mod-php5 php5 php5-mysql phpmyadmin* * * * * * * <http://www.showmetech.com.br> *LEANDRO MARINO* | Showmetech Estatístico | Fotógrafo Mobile: 55 21 9845-7707 Email: contato(a)leandromarino.com.br Profiles SMT: facebook <http://www.facebook.com/showmetech> | twitter<http://www.twitter.com/showmetech> | Google+ <https://plus.google.com/115440786884851522659> <http://www.leandromarino.com.br>

2 2

interaction plot
by Natalia Martins 01 Jul '13

01 Jul '13

Prezados, bom dia! Estou fazendo uns graficos utilizando a função interaction. plot. O grafico sai em ordem alfabética no eixo x, como eu poderia proceder para que saísse em outra ordem? interaction.plot(CP,Clone,Altura, xlab = 'Corpo de prova', ylab = 'Altura', col=c("blue","red"),lty = 1) O brigada, Natália da Silva Martins

2 2

carregar arquivos
by alanarocha＠sapo.pt 01 Jul '13

01 Jul '13

olá a todos, alguém me sabe esplicar qual a diferença de eu fazer: load ("Cinacalset_DosetH3DosetC3.RData") ou load (".RData") ? e porque no primeiro caso dá o erro abaixo? como resolver? Error in readChar(con, 5L, useBytes = TRUE) : cannot open the connection In addition: Warning message: In readChar(con, 5L, useBytes = TRUE) : cannot open compressed file 'Cinacalset_DosetH3DosetC3.RData', probable reason 'No such file or directory' Cumprimentos Ana

2 1