Março 2014 - R-br

Re: [R-br] Teste Paralelismo no r
by Robert Iquiapaza 07 Mar '14

07 Mar '14

Geovane Assumindo que as regressões propostas nos gráficos (p8) fazem sentido e satisfazem os pressupostos do modelo de regressão linear poderia fazer o seguinte: library(car) escol='Superior/Pós graduação' model2=lm(Risco~0+(Idade*Quintis),subset(dados,Escolaridade3cat==escol)) summary(model2) # o teste para ‘Idade:Quintis5‘ na verdade está testando se os coeficientes angulares são iguais ( ou a diferença deles é zero) #para testar a diferença dos interceptos linearHypothesis(model2, "Quintis1 - Quintis5 =0") linearHypothesis(model2, "Idade:Quintis5 =0") Aqui mostrei os testes para uma escolaridade (escol), na verdade vc tem mais informação, teria que rodar um só modelo com todas as informações e fazer os testes nele: model1=lm(Risco~0+Idade*Escolaridade3cat*Quintis,data=dados) summary(model1) Boa sorte Robert From: Benilton Carvalho Sent: Thursday, March 06, 2014 2:15 PM To: r-br ; geovane barbosa Subject: Re: [R-br] Teste Paralelismo no r Geovane, Para realizar uma analise de dados, o estatistico precisa entender a (hi/e)storia dos dados. Por exemplo, sem saber do que se tratam os dados, e' dificil acreditar que a regressao linear simples que voce apresenta no seu exemplo seja adequada. A outra coisa: como/onde entra o fator "Quintis" existente no seu conjunto de dados? Enfim, supondo adequacidade em todos os sentidos (cabera' a vc verificar os pressupostos), se o interesse e' averiguar a identidade das regressoes nos tres grupos escolares (vou ignorar Quintis aqui, pois realmente nao sei o que e', nem como entra na sua analise), eu ajustaria os modelos com interacao (vis-a-vis Escolaridade) e compararia com o modelo simplificado (que contem apenas Idade): fit1 = lm(Risco ~ Idade*Escolaridade3cat, data=dados) fit0 = lm(Risco ~ Idade, data=dados) anova(fit0, fit1) benilton Em 6 de março de 2014 10:43, geovane barbosa <geovanecb(a)yahoo.com.br> escreveu: Olá pessoal , segue em anexo um exemplo reprodutível. Quero verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e banco num link. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") https://www.dropbox.com/s/ymyzh7p0xhmumda/bancoteste.csv Em Quinta-feira, 6 de Março de 2014 10:12, geovane barbosa <geovanecb(a)yahoo.com.br> escreveu: Olá Prof, segue em anexo um exemplo reprodutível. Como eu lhe disse anteriormente, queria verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e o banco de dados ficticio. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") Em Quarta-feira, 5 de Março de 2014 16:54, Benilton Carvalho <beniltoncarvalho(a)gmail.com> escreveu: A gente pede que exemplos reprodutíveis sejam dados para que possamos dar uma resposta definitiva sem a necessidade de idas e vindas. Como eu não entendi quais são os três cenários, fica a tarefa p vc: o caso a que vc se refere consiste em testar se as interações das variáveis de interesse com os referidos grupos são conjuntamente iguais a zero. b On Mar 5, 2014 11:00 AM, "geovane barbosa" <geovanecb(a)yahoo.com.br> wrote: Olá pessoal bom dia. Estou com uma questão que não consigo soluciona-la. Realize um regressão linear múltipla para comprar um grau de risco por idade e por escolaridade. Quero verificar se a inclinação proporcionada por essas retas é estatisticamente significativa, ou seja comparando a inclinação nos três cenários. bem como tambem seu intercepto. Como eu faço isso? Seria o test de Brant? Não consigo encontrar tal ferramenta no r. abraços _______________________________________________ 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 kappa por categorias
by branco.sk8＠ig.com.br 06 Mar '14

06 Mar '14

prezados, gentileza. como se faz para fazer o teste kappa por categorias (condicionar a outra variável?) (ex. testar concordância kappa entre dois teste em função do sexo, ou em função da idade etc) (OBS. já sei fazer o teste kappa, só não sei fazer o kappa por outras variáveis... obrigado

2 1

Re: [R-br] Teste Paralelismo no r
by geovane barbosa 06 Mar '14

06 Mar '14

Olá pessoal , segue em anexo um exemplo reprodutível. Quero verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e banco num link. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") https://www.dropbox.com/s/ymyzh7p0xhmumda/bancoteste.csv Em Quinta-feira, 6 de Março de 2014 10:12, geovane barbosa <geovanecb(a)yahoo.com.br> escreveu: Olá Prof, segue em anexo um exemplo reprodutível. Como eu lhe disse anteriormente, queria verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e o banco de dados ficticio. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") Em Quarta-feira, 5 de Março de 2014 16:54, Benilton Carvalho <beniltoncarvalho(a)gmail.com> escreveu: A gente pede que exemplos reprodutíveis sejam dados para que possamos dar uma resposta definitiva sem a necessidade de idas e vindas. Como eu não entendi quais são os três cenários, fica a tarefa p vc: o caso a que vc se refere consiste em testar se as interações das variáveis de interesse com os referidos grupos são conjuntamente iguais a zero. b On Mar 5, 2014 11:00 AM, "geovane barbosa" <geovanecb(a)yahoo.com.br> wrote: > >Olá pessoal bom dia. > > > > >Estou com uma questão que não consigo soluciona-la. > > > > >Realize um regressão linear múltipla para comprar um grau de risco por idade e por escolaridade. > > >Quero verificar se a inclinação proporcionada por essas retas é estatisticamente significativa, ou seja comparando a inclinação nos três cenários. bem como tambem seu intercepto. Como eu faço isso? Seria o test de Brant? Não consigo encontrar tal ferramenta no r. > > >abraços > > > > > >_______________________________________________ >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.

2 1

Coordenadas Geográficas no Pacote GeoR
by Helder Gramacho 06 Mar '14

06 Mar '14

Boa tarde, É possível trabalhar com coordenadas geográficas (em graus decimais) no pacote GeoR? Caso seja possível, para estimar o semivariograma utilizamos as distâncias calculadas em metros ou em graus decimais? A maioria dos trabalhos que tenho visto tem utilizado coordenadas planas UTM, e aí surgiu a dúvida: Quando a região estiver situada entre dois fusos UTM como proceder? Att, Hélder Santos *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 4

Teste Paralelismo no r
by geovane barbosa 05 Mar '14

05 Mar '14

Olá pessoal bom dia. Estou com uma questão que não consigo soluciona-la. Realize um regressão linear múltipla para comprar um grau de risco por idade e por escolaridade. Quero verificar se a inclinação proporcionada por essas retas é estatisticamente significativa, ou seja comparando a inclinação nos três cenários. bem como tambem seu intercepto. Como eu faço isso? Seria o test de Brant? Não consigo encontrar tal ferramenta no r. abraços

2 1

Colocar estatística Z na curva normal - TH
by Andre Oliveira 04 Mar '14

04 Mar '14

Pessoal, estou aqui tentando, mas sem sucesso colocar a estatística Ztab, Zcal, mas sem sucesso. media <- 170 stder <- 8 x <- seq(media-3*stder, media+3*stder, 0.01) y <- dnorm(x, media, stder) plot(x, y, type='l',axes=F,xlab="",ylab="") abline(h=0) hx <- seq(media+2*stder, media+3*stder, 0.01) hy <- dnorm(hx, media, stder) n <- length(hy) polygon(c(hx, rev(hx)), c(hy, rep(0, n)), col=2) # Tentei text(locator(n=1),"Zcal=3.53") text(locator(n=1),"Ztab=1,64") Mas, preciso aumentar a área abaixo do gráfico para plotar estes valores e não consegui! Alguém tem alguma sugestão? Outra situação é que Z_{cal} para ficar mais elegante não funcionou com locator. obrigado

3 4

Corpus não funciona dentro de um TermDocumentMatrix() - Package TM
by Fernando Gama 04 Mar '14

04 Mar '14

Galera, Estou tendo um problema em relação a leitura de um corpus, quando executo a função TermDocumentMatrix(). :: O trecho de erro do meu código: term.sample<-TermDocumentMatrix(corpus, control=list(wordLengths=c(4,100))) :: Quando executo este é o erro lançado: Error: inherits(doc, "TextDocument") is not TRUE Realizei alguns testes para ver se funcionava como: Dei o inspect(corpus) antes do (term.sample) ele retorna normalmente com os dados pós parsing. E também testei os dados nativos da biblioteca, no caso o "crude" no lugar do corpus, o que funcionou também. Por esse motivo creio eu, que o problema pode estar no corpus, mas imagino ao mesmo tempo, que se estivesse lá, teria retornado algum erro, ou não? Alguém tem alguma idéia? Já enfrentou esse tipo de problema? Abs, Fernando.

1 0

transformação de booxcox
by jorgeas＠ufcg.edu.br 04 Mar '14

04 Mar '14

Ola pessoal estou tetando utilizar a transformação de boxcox nestes dados e o resultado é uma reta alguem dpo Prof. Jorge Alves de Sousa, DSc. Professor Adjunto II Universidade Federal de Campina Grande - UFCG Centro de Educação e Saúde - CES Unidade Acadêmica de Educação - UAE Olho D'Agua da Bica s/n , Cuité-PB, Cep.: 58175-000 Telefone: (83) 3372-1900

2 1

RStudio - Shiny Server Pro
by Benilton Carvalho 03 Mar '14

03 Mar '14

Aos interessados: http://ow.ly/ubO8b

1 0

Ogiva - fdth
by Andre Oliveira 03 Mar '14

03 Mar '14

Prezados estou tentando produzir uma Ogiva aqui usei os códigos. library(fdth) x=c(45,51,50,58,50,44,46,57,42,41,60,58,41,50,54,60,52,46,52,51) tab=fdt(x,start=41,h=6,end=60) tab plot(tab, type= "fp") Porém o meu gráfico não é ligado nas frequências zero. Alguém tem alguma sugestão de como ligar nas frequências zero anterior e posterior as ultimas frequências do histograma? André Oliveira Souza

2 2