Agosto 2012 - R-br

teste de proporção
by Edson Lira 09 Ago '12

09 Ago '12

O que está errado na rotina abaixo? fvc04<-c(2236,2167,2061,2049,1903) #fluxo de veiculos DETRAN/Coroado 2004 fvc12<-c(2192,2214,2089,1947,1810) #fluxo de veiculos DETRAN/Coroado 2012 prop.test(fvc04,fvc12) movcd1<-c(1689,1661,1599,1585,1536) #fluxo de veiculos Coroado/DETRAN 2004 movcd2<-c(1671,1680,1759,1708,1603)#fluxo de veiculos Coroado/DETRAN 2004 prop.test(movcd1,movcd2) [ ]'s. Edson Lira Estatístico Manaus-Amazonas

1 0

Re: [R-br] Como estimar R2 de modelo linear a partir de previsões em novos dados?
by Ivan Bezerra Allaman 09 Ago '12

09 Ago '12

Talvez a ajuda seja mais efetiva aqui. https://groups.google.com/forum/?fromgroups#!forum/MedStats \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Re: [R-br] Como estimar R2 de modelo linear a partir de previsões em novos dados?
by Ivan Bezerra Allaman 09 Ago '12

09 Ago '12

Bom dia Pedro! Tu concordas comigo que a partir do momento que ajustastes um modelo por meio de uma amostra, o modelo irá fazer uma estimativa da variável resposta independentemente de qual amostra você utilize (desde que os pontos estejam dentro do intervalo no qual o modelo foi ajustado é claro) com aquela precisão no qual foi construído o modelo. Se vc utilizar um modelo ajustado com a amostra w1 e depois usar o modelo para fazer estimativas com a amostra w2, estas estimativas foram estimadas com a precisão dada no primeiro ajuste. Se você realmente, quer avaliar o ajuste do modelo feito com a amostra w1 e depois com a amostra w2 é simples, basta ajustar um modelo feito com a amostra w1 e depois ajustar outro modelo feito com a amostra w2 e comparar os R2, embora não vejo muito sentido nisso, pois é claro, que os valores serão diferentes sempre, pois é aquela velha história, se retirarmos 'n' amostras de uma população e retirarmos de cada amostra a média, estas médias serão diferentes obviamente pelo simples processo de amostragem. Abraço! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 4

Correlação de pearson
by Marcelo Claro de Souza 08 Ago '12

08 Ago '12

Estou rodando minhas correlações de Pearson segundo modelo abaixo: a=c(1,2,3,4,5,6) b=c(10,13,7,6,8,9) c=c(12,15,16,22,30,15) d=c(5,6,7,8,9,10) x=cbind(a,b,c,d) x=data.frame(x) cor(x,method='pearson') É possível determinar o número de casas decimais que será exibido na matriz de correlações? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 3

bibliog para estudo, preferencialmente em R.
by Cleber N.Borges 08 Ago '12

08 Ago '12

Prezados, Boa Tarde, Peço por gentileza indicações de Referencias Bibliograficas para 02 tópicos de estatística, que são: 1) Intervalo de confiança de uma prevalência; 2) Análise não paramétrica para respostas subjetivas/discretas. Acredito que o primeiro seja referente à Epidemologia, porém nunca estudei nada a respeito. Caso alguém puder confirmar que o termo 'Prevalência' seja usual, ficaria grato. Eu sou químico de formação, porém tenho uns 30 livros de estatística e nunca me deparei com esses 02 tópicos. Agradeço antecipadamente por qualquer ajuda. Obrigado. Cleber N.Borges

1 0

Erro nls: matriz gradiente singular com estimativas de parâmetros iniciais
by Gustavo Dias Azevedo 08 Ago '12

08 Ago '12

Prezados, Estou tentando ajustar dados a um modelo de cinética química do tipo y = K*x(^a)*w(^b)/(z^c) x, w e z são as concentrações dos meus reagentes e produtos (dados) K, a, b e c são os parâmetros que eu preciso determinar no modelo. Utilizei o seguinte comando dados <- read.table("Teste2.txt",header=TRUE) modelo <- nls(y~K*(x^a)*(w^b)*(z^c),data=dados,start=list(K=1,a=0.33,b=0.66,c=-1.33)) Porém, ao executar a última linha, recebo a seguinte mensagem: *matriz gradiente singular com estimativas de parâmetros iniciais* Já procurei alguns tópicos e algumas pessoas tiveram o mesmo problema em ajustes diferentes, e encontrei diversas "soluções" diferentes, como usar o argumento alg="plinear", trocar os chutes iniciais etc, porém nenhuma delas resolveu meu problema. Alguém sabe como me ajudar? Devo trocar de comando ou fazer alguma manipulação no modelo? Obrigado!

2 2

Como estimar R2 de modelo linear a partir de previsões em novos dados?
by Pedro Emmanuel Alvarenga Americano do Brasil 08 Ago '12

08 Ago '12

Amigosr de R, Ajustei um modelo linear com parte dos dados que tenho. Depois fiz uma previsão desse modelo com outra parte dos dados que tenho. O que gostaria de fazer é estimar o R2 e o R2 ajustado ou outra medida de ajuste geral do modelo para as previsões do segundo conjunto de dados para comparar com o ajuste do primeiro e ter uma ideia se é muito diferente. Procurei um pouco por ai mas não encontrei uma função que estimasse. Alguem tem alguma sugestão? library('rms') lmfit6 <- ols(desfecho ~ FOISini+idade+lung+intercorrencia+VM2+disphagyaSev+degen+rcs(numSessions)+degen * rcs(numSessions),w1,x=T,y=T) lmfit6 Linear Regression Model ols(formula = desfecho ~ FOISini + idade + lung + intercorrencia + VM2 + disphagyaSev + degen + rcs(numSessions) + degen * rcs(numSessions), data = w1, x = T, y = T, penalty = p.bic$penalty) Model Likelihood Discrimination Ratio Test Indexes Obs 372 LR chi2 403.24 R2 0.662 sigma 1.5341 d.f. 14.96145 R2 adj 0.648 d.f. 356.0385 Pr(> chi2) 0.0000 g 2.279 Residuals Outcome Min 1Q Median 3Q Max -4.06353 -1.09756 0.09737 1.09347 4.32595 predict(lmfit6,w2) # previsoes do segundo conjunto de dados mas não as medidas de discriminação. Abraço forte, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

Linhas com transparência
by Leandro Marino 08 Ago '12

08 Ago '12

*Caros,* * * *certa vez vi aqui na lista um exemplo de gráfico de dispersão aonde os pontos tinham alguma transparência de forma que quanto mais sobrepostos ficavam, mas escura aquela região ficava. * * * *Gostaria de fazer o mesmo para um gráfico de linhas. Irei plotar multiplas linhas em um gráfico e gostaria que elas tivessem esta transparência. Alguém pode me ajudar a fazer algo deste tipo? Não tenho nem idéia de qual foi o pacote gráfico utilizado para o exemplo dos pontos transparentes.* * * *Grato por qualquer ajuda!* * * *Leandro Marino*

2 2

Bayesthres para Bayesthresh
by Fabio Mathias Corrêa 07 Ago '12

07 Ago '12

Caros colegas, Por sugestão do Prof. Brian Ripley a library Bayesthres teve seu nome alterado para Bayesthresh. Valeu!!! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076

1 0

Eliminar outliers de uma planilha.
by Marcelo Claro de Souza 07 Ago '12

07 Ago '12

Após plotar uns boxplots observei a presença de alguns outliers. É possível gerar uma nova planilha sem os outliers ou tenho de remover um a um manualmente? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1