- R-br - DINF UFPR

R ler bases SAS
by Richard Santos 23 Nov '12

23 Nov '12

Um colega disse que "Especificamente com o formato sas7bdat tem o pacote sas7bdat. Eu já li bases do sas com esse pacote no R open source. Funciona! http://cran.r-project.org/web/packages/sas7bdat/index.html" Se ler o Censo Demográfico, que nem é um BIGDATA, eu já fico feliz. Rsrs. Lê? Rsrsrs. Abraços, Richard

1 0

Seminários dos Usuários de R no Brasil
by Ary Casado 23 Nov '12

23 Nov '12

Pelo que vejo na lista de usuários temos representantes de vários estados brasileiros e de várias instituições. Acredito que seja a hora de organizar os usuários e desenvolvedores em torno deste foco de interesse: a utilização e desenvolvimento do R tanto no meio científico como no âmbito da aplicabilidade prática. Não sei quantos interessados existem no Recife, mas acredito que já exista o interesse de iniciar encontros e seminários com o objetivo maior gerar um grupo de usuários e desenvolvedores. Arykerne Nascimento Engenheiro Químico - UFPE Especialista em Informática - UFPE Mestrando em Ciências Nucleares - UFPE

2 1

Movimento Valorização do profissional de Estatística, já!
by Richard Santos 23 Nov '12

23 Nov '12

Caros, Devido aos comentários recentes e conclusões, que tal, quem concordar, inclusive tendo em vista o ANO INTERNACIONAL DA ESTATÍSTICA, em 2013, lançarmos o movimento "valorização do profissional de Estatística já"? Esse movimento contemplaria: - Só repetindo, de novo, acolher e valorizar os outros profissionais que se especializam em alguma técnica estatística, mas não os confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos preciso. - Parabenizar as instituições que fazem concurso para Estatístico e vagas privadas que contemplam o estatístico, quando o trabalho for em sua maioria fazer estatística; - Motivar o trabalho em equipe, economista, com estatístico, com engenheiro, com TI, com cientista social, etc; - Motivar os estatísticos a fazerem o concurso de Analista Geral existentes atualmente no IPEA e no IBGE, e para nos fortalecer lá dentro, fazendo frente aos economistas; - Definir estatístico, na regulamentação: graduação, outro profissional que tenha pós em estatística? Matemática com pós em estatística? Economista com pós em estatística? Manter a definição de estatístico e técnico em estatística; - conscientizar nos cursos de Estatística básica, nos mini-cursos, etc; - Fazer marketing à sociedade, à sociedade civil, aos empregadores e intelectuais brasileiros com relação à Estatística e ao profissional de estatística; - Só repetindo, de novo, acolher e valorizar os outros profissionais que se especializam em alguma técnica estatística, mas não os confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos preciso. Porque se um dia no Brasil a estatística virar uma coisa que qualquer um pode fazer ou se virar totalmente "promíscua", eu só recomendo livros, ensino quais botões apertar, e não mudo de país não, eu vou fazer outra coisa, aqui no Brasil mesmo, rs. Abraços, Att. Richard Estatístico - CONRE 1

1 0

Comando install.packages(pacote)
by Richard Santos 23 Nov '12

23 Nov '12

Pessoal, Eu dei o comando install.packages('pacote'), e o R me deu as mensagens: --- Please select a CRAN mirror for use in this session --- Aviso: unable to access index for repository http://cran.fiocruz.br/src/contrib Aviso: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/src/contrib Aviso: unable to access index for repository http://cran.fiocruz.br/bin/windows/contrib/2.15 Aviso: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.15 Mensagens de aviso perdidas: 1: In open.connection(con, "r") : não foi possível conectar a 'cran.r-project.org' na porta 80. 2: package ‘pacote’ is not available (for R version 2.15.2) What the hell? Em 23/11/12, Richard Santos<jamesrichardsantos(a)gmail.com> escreveu: > Excelentes dicas, Flávio, obrigado pela atenção. Gostei do que foi > falado sobre a linguagem, e sobre a instalação de pacotes local. Vou > testar, realmente, para aprender o R aos poucos, também. Só uma > confirmação: essas opções do Revolution e dataframe, etc, são para ler > bases grandes também, né? > Grato. > Richard > > Em 22/11/12, Flavio Barros<flaviomargarito(a)gmail.com> escreveu: >> - sempre prefiro a programação, mas a programação do R é mais extensa >> do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no >> SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1; >> >> Os programas em R são bem compactos na verdade (em termos de sintaxe), >> mesmo se comparados a outras linguagens como python, ruby ou mesmo o >> matlab. Vide ( >> http://blog.revolutionanalytics.com/2012/11/which-programming-language-is-t… >> ) >> >> Quando a usar X4 = 1 ou dados$x4 = 1 no R é possível das duas formas, >> basta >> 'fixar' o conjunto de dados antes com o attach(dados). Particularmente >> eu >> prefiro escrever x <- 4 (padrão do google também! >> http://google-styleguide.googlecode.com/svn/trunk/google-r-style.html#assig… >> que diga-se de passagem é um usuário ativo do R.) >> >> - para grandes bases, tem-se algumas opções como Revolution, mas mesmo >> o Revolution parece que apresenta um problema que encontrei similar ao >> Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou >> formato. Eu não queria encontrar nenhuma limitação quanto a ler >> grandes bases; >> >> Na verdade isso não é necessário. No R existem algumas estruturas de >> dados >> como o matrix, o data.frame e o list. De fato, se você quiser ler os seus >> dados em uma matriz, todos os dados tem que ser do mesmo tipo, mas em um >> data.frame não! O formato XDF da revolution também não impõe esse tipo de >> limitação. Caso as variáveis tenham tamanhos diferentes ainda é possível >> ler como um list de forma que nem o tamanho das variáveis precisa ser >> igual. Mande um e-mail para o suporte da revolution que acredito que o >> problema foi sua implementação. Eles respondem rápido. Já recebei uma >> resposta no domingo depois de 20 minutos! >> >> - para cada procedimento que se quer fazer, pode ser necessário >> instalar um pacote diferente, eu tentei pedir para o administrador do >> meu trabalho que uma pasta minha fiique livre para eu instalar >> pacotes, mas não sei se isso vai ser feito. >> >> Acho que já foi falado aqui, mas não custa lembrar: não é necessário >> privilégio de administrador para instalar pacotes. Caso você não seja >> administrador basta dar um install.packages('pacote') que ele vai criar >> um >> biblioteca local de pacotes. Esse não é um problema. >> >> No mais, é normal para aprender qualquer software ter que gastar >> tempo, mas as questões que citei acima se quiserem podem considerar >> como gosto pessoal mesmo. Mas acho importante eu citar essas questões >> para os empregadores saberem, pelo menos. Pois creio que essas >> questões não permitem ainda dizer para todos que o R é uma opção fácil >> de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não >> que não seja ótimo. Apenas acho que ainda há argumentos para que uma >> empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS, >> quando for necessário. >> >> Na minha opinião existem casos e casos para migrar do SAS ou SPSS para o >> R. >> Uma questão importante é o treinamento: se a empresa já tem profissionais >> bem treinados usando SAS, muito código pronto que seria difícil migrar a >> melhor solução é continuar a usar essas ferramentas. Agora no caso de um >> consultor, uma startup ou uma empresa pública que vai optar por um pacote >> estatístico acredito que o R é sim uma melhor opção. >> >> O R é muito utilizado no Google por exemplo. Se é bom para o Google dá >> para >> usar em qualquer empresa. >> >> Não advogo todo mundo usar o R, mas a linguagem tem muitas >> características >> interessantes, é bem intuitiva para quem programa em outras linguagens >> (python, C, Java por ex) e tem custo muito menor que o SAS ou SPSS. E é >> open source. >> >> No seu caso Richard que acredito que deve ter uma boa experiência em SAS, >> trocar pode não ser a melhor opção para você. Se você tem acesso a uma >> licença do SAS ou do SPSS e tem experiência com essas ferramentas não >> vejo >> porque trocar. No entanto é interessante ir aprendendo o R em paralelo. >> >> Flavio Barros >> [image: Facebook] >> <http://www.facebook.com/flavio.barros.1650?ref=tn_tnmn> [image: >> LinkedIn] <http://www.linkedin.com/profile/view?id=61839390&trk=tab_pro> >> Contact me: [image: Google Talk] flaviomargarito(a)gmail.com >> “"Sing like no one's listening, love like you've never been hurt, dance >> like nobody's watching, and live like its heaven on earth." - Mark >> Twain<http://www.quotesdaddy.com/quote/1040448/mark-twain/sing-like-no-ones-liste…> >> ” Get this email app! >> <http://www.wisestamp.com/apps/quotes?utm_source=extension&utm_medium=email&…> >> >> Create your free signature: >> <http://r1.wisestamp.com/r/landing?promo=33&dest=http%3A%2F%2Fwww.wisestamp.…> >> CLICK >> HERE!<http://r1.wisestamp.com/r/landing?promo=33&dest=http%3A%2F%2Fwww.wisestamp.…> >> >> >> >> 2012/11/22 Richard Santos <jamesrichardsantos(a)gmail.com> >> >>> ** >>> >>> >>> Caros, >>> >>> Apenas gostaria de informar que eu, como estatístico, particularmente, >>> não vou divulgar para as empresas e instituições que o R é uma opção >>> com os mesmos níveis de facilidade do que o SAS e o SPSS. Eu não tenho >>> dúvidas de que alguém que fique expert no R consiga fazer tudo e ainda >>> mais rápido do que o SAS. >>> >>> Apenas que eu pessoalmente não vou optar pelo R agora, pelo menos por >>> enquanto, ainda que eu tenha que trazer notebook para o trabalho, por >>> alguns motivos que não tive adesão por questão de gosto pessoal: >>> >>> - sempre prefiro a programação, mas a programação do R é mais extensa >>> do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no >>> SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1; >>> >>> - para grandes bases, tem-se algumas opções como Revolution, mas mesmo >>> o Revolution parece que apresenta um problema que encontrei similar ao >>> Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou >>> formato. Eu não queria encontrar nenhuma limitação quanto a ler >>> grandes bases; >>> >>> - para cada procedimento que se quer fazer, pode ser necessário >>> instalar um pacote diferente, eu tentei pedir para o administrador do >>> meu trabalho que uma pasta minha fiique livre para eu instalar >>> pacotes, mas não sei se isso vai ser feito. >>> >>> No mais, é normal para aprender qualquer software ter que gastar >>> tempo, mas as questões que citei acima se quiserem podem considerar >>> como gosto pessoal mesmo. Mas acho importante eu citar essas questões >>> para os empregadores saberem, pelo menos. Pois creio que essas >>> questões não permitem ainda dizer para todos que o R é uma opção fácil >>> de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não >>> que não seja ótimo. Apenas acho que ainda há argumentos para que uma >>> empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS, >>> quando for necessário. >>> >>> Obrigado. >>> >>> Abraços, >>> Att. >>> Richard >>> >>> >> >

5 13

pacote fora do ar
by Rodrigo Plei 23 Nov '12

23 Nov '12

Prezados, O pacote "PBSmapping" foi retirado recentemente do CRAN. Alguém teria o arquivo zip do mesmo para compartilhar? Agradeço desde já. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777 Google Scholar: http://scholar.google.com.br/citations?user=zX_EezEAAAAJ&hl=pt-BR

2 2

Função Bartlett
by tiago souza marçal 23 Nov '12

23 Nov '12

Alguém sabe como selecionar o valor p após a realização do teste de Bartlett? Att. Tiago.

2 2

Re: [R-br] Criando Rotina
by Giovane Pretto Azambuja 23 Nov '12

23 Nov '12

Boa tarde, pessoal! Consegui realizar. Obrigado pela ajuda de vocês. Grande abraço Giovane P. Azambuja Enviado do Email do Windows *De:* Luis Iván Ortiz Valencia *Enviado:* 22 de novembro de 2012 12:33 *Para:* r-br(a)listas.c3sl.ufpr.br *Assunto:* Re: [R-br] Criando Rotina Como sao varios objetos, pode criar uma lista contendo esse objetos, depois pode chamar ( com o for ) um a um os elementos da lista(objetos) e fazer o summary. Ivan Em 22 de novembro de 2012 12:27, Alisson Lucrecio <alissonluc(a)yahoo.com.br>escreveu: Use o pacote plyr e ver do for. # exemplo para fazer uma aov list.dados <- lapply(data[,-c(1)], function(y){ dm <- data[,c(1)] y <- y dfm <- cbind(dm, y) return(list( dm = dm, y = y, dfm = dfm)) }) list.av <- llply(list.dados, function(list){ with(list, aov(y ~ trar, data=dfm)) }) list.anova <- lapply(list.av, anova) Alisson Lucrécio da Costa ------------------------------ *From:* Giovane Pretto Azambuja <gibuja(a)gmail.com> *To:* "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> *Sent:* Thursday, November 22, 2012 11:47 AM *Subject:* [R-br] Criando Rotina Olá, pessoal Sou um pouco novato em R, mas vamos tentar. Tenho de fazer a análise estatística, pela função summary(), de uma pancada de objetos. Gostaria de saber se existe uma maneira em que o R faça essas análises, e as armazene em um objeto diferente para cada análise, de maneira automática. Tipo uma rotina. Creio que consigo fazer isso com o comando for, porém como são vários objetos não me veio em mente como organiza-lo. OBS: Não precisa formular a rotina inteira, apenas dê uma luz. Agradeço quem puder ajudar Giovane P. Azambuja Enviado do Email do Windows _______________________________________________ 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. -- Luis Iván Ortiz Valencia Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ Estatístico Msc. Spatial Analyst Msc. _______________________________________________ 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

Path Analysis
by Malcon do Prado Costa 23 Nov '12

23 Nov '12

Caros colegas usuários do programa R, gostaria de saber se existe um pacote ou rotina para elaboração de uma análise de caminhos (Path Analysis) no R. Estou ajudando um amigo que trabalha com palmeiras do RS e ele quer fazer relações entre variáveis ambientais e biométricas, a partir da análise de caminhos. Os senhores conhecem algum pacote que trabalha com esta análise? Grato -- Malcon do Prado Costa "Aprenda como se você fosse viver para sempre. Viva como se você fosse morrer amanhã.(...) A vida é um todo indivisível.(...) É melhor que fale por nós a nossa vida, que as nossas palavras.(...) A minha fé, nas densas trevas, resplandece mais viva." (Mahatma Gandhi) -- Malcon do Prado Costa - APRENDIZ de jardineiro, mateiro e curandeiro. "Aprenda como se você fosse viver para sempre. Viva como se você fosse morrer amanhã.(...) A vida é um todo indivisível.(...) É melhor que fale por nós a nossa vida, que as nossas palavras.(...) A minha fé, nas densas trevas, resplandece mais viva." (Mahatma Gandhi)

5 4

uma ajuda
by alanarocha＠sapo.pt 23 Nov '12

23 Nov '12

bom dia a todos, oEu estou a fazer um trabalho onde tenho uma variavel condicional cm dois niveis "H" e "C" para relacionar com duas variaveis ou seja verificar se no H ha relaçao entre pthi e dose e em C ver se se ha relacao entre pthi e dose. eu fiz analise de regressão e uma anova que nao sei se será necessário. mostro meus dados e codigo. ao grafico que fiz nao cnsigo cola-lo, existe algum mil para onde possa envialo em anexo se for necessario. 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 pthi dose #'pthi ~ dose+grupo' indica: modele pthi como função estatística das variáveis dose e grupo (efeito aditivo dos modelos lineares); #a variável grupo entrou alterando o intercepto da regressão entre (pthi) e (dose). class(grupo) [1] "factor" class(pthi) [1] "numeric" class(dose) [1] "integer" modelo <- lm(pthi~dose+grupo) modelo Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 760.1821 0.1729 -159.5969 summary(modelo3) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -858.70 -322.92 -64.76 219.93 1939.70 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 760.18206 84.59910 8.986 7.13e-16 *** dose 0.17291 0.08286 2.087 0.0385 * grupoH -159.59687 92.60311 -1.723 0.0868 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 541.9 on 159 degrees of freedom Multiple R-squared: 0.0325, Adjusted R-squared: 0.02033 F-statistic: 2.671 on 2 and 159 DF, p-value: 0.07231 anova(modelo3) Analysis of Variance Table Response: pthi Df Sum Sq Mean Sq F value Pr(>F) dose 1 696390 696390 2.3711 0.12559 grupo 1 872380 872380 2.9703 0.08675 . Residuals 159 46698728 293703 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 x11() par(col="green") plot(pthi~dose+grupo,x1lab="dose",x2lab="grupo",ylab="pthi", main="grafico de dspersão") #Adicionar a reta da regressão no gráfico: abline(modelo3) #recta da regressão text("y=760.1821+0.1729dose-159.5969grupo, p=0.07231") ue preciso obter um grafico dedispersão onde se veja uma curva do H e outra do C mesmo não havend relação a minha onclusao O modelo entre pthi~dose entre grupos Ao relacionar pthi~dose sendo o grupo H ou C a estimativa da recta dos mínimos quadrados/regressão/relação foi Pthi estimado=760.1821+0.1729*dose-159.5969*grupo a tabela ANOVA revela cada pvalue superior a 5 e 10% o que traduz a ineficiência do modelo de regressão ou seja não há relação. Estes modelos de regressão não têm qualquer significância estatística. o que posso escrever sobre a anova para quem nao sabe o que é perceber melhor? Agradeço desde ja a vossa ajuda. Ana Rocha

1 0

uma ajuda
by alanarocha＠sapo.pt 23 Nov '12

23 Nov '12

bom dia a todos, oEu estou a fazer um trabalho onde tenho uma variavel condicional cm dois niveis "H" e "C" para relacionar com duas variaveis ou seja verificar se no H ha relaçao entre pthi e dose e em C ver se se ha relacao entre pthi e dose. eu fiz analise de regressão e uma anova que nao sei se será necessário. mostro meus dados e codigo. ao grafico que fiz nao cnsigo cola-lo, existe algum mil para onde possa envialo em anexo se for necessario. 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 pthi dose #'pthi ~ dose+grupo' indica: modele pthi como função estatística das variáveis dose e grupo (efeito aditivo dos modelos lineares); #a variável grupo entrou alterando o intercepto da regressão entre (pthi) e (dose). class(grupo) [1] "factor" class(pthi) [1] "numeric" class(dose) [1] "integer" modelo <- lm(pthi~dose+grupo) modelo Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 760.1821 0.1729 -159.5969 summary(modelo3) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -858.70 -322.92 -64.76 219.93 1939.70 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 760.18206 84.59910 8.986 7.13e-16 *** dose 0.17291 0.08286 2.087 0.0385 * grupoH -159.59687 92.60311 -1.723 0.0868 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 541.9 on 159 degrees of freedom Multiple R-squared: 0.0325, Adjusted R-squared: 0.02033 F-statistic: 2.671 on 2 and 159 DF, p-value: 0.07231 anova(modelo3) Analysis of Variance Table Response: pthi Df Sum Sq Mean Sq F value Pr(>F) dose 1 696390 696390 2.3711 0.12559 grupo 1 872380 872380 2.9703 0.08675 . Residuals 159 46698728 293703 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 x11() par(col="green") plot(pthi~dose+grupo,x1lab="dose",x2lab="grupo",ylab="pthi", main="grafico de dspersão") #Adicionar a reta da regressão no gráfico: abline(modelo3) #recta da regressão text("y=760.1821+0.1729dose-159.5969grupo, p=0.07231") ue preciso obter um grafico dedispersão onde se veja uma curva do H e outra do C mesmo não havend relação a minha onclusao O modelo entre pthi~dose entre grupos Ao relacionar pthi~dose sendo o grupo H ou C a estimativa da recta dos mínimos quadrados/regressão/relação foi Pthi estimado=760.1821+0.1729*dose-159.5969*grupo a tabela ANOVA revela cada pvalue superior a 5 e 10% o que traduz a ineficiência do modelo de regressão ou seja não há relação. Estes modelos de regressão não têm qualquer significância estatística. o que posso escrever sobre a anova para quem nao sabe o que é perceber melhor? Agradeço desde ja a vossa ajuda. Ana Rocha

1 0