- R-br - DINF UFPR

Re: [R-br] microdados censo IBGE 2010
by Roney Fraga 14 Ago '12

14 Ago '12

Manoel, O Marcos F Silva criou algunas funções maravilhosas para ajudar a trablhar com microdados no R. No endereço: https://sites.google.com/site/marcosfs2006/textos baixar o arquivo Usando R para Importar Microdados.pdf e as funções no link .zip. Com base no censo de 2000, onde apenas a tabela para o estado de São Paulo tem aproximadamente 3,5 GB quando em formato .csv, acho que será um pouco complicado trabalhar com esse volume de dados diretamente no R. Para lidar melhor com isso optei por colocar os arquivos para cada estado como uma tabela do banco de dados PostgreSQL, isso resulta em um bom ganho de tempo, até porque as consultas básicas podem ser feitas no próprio PostgreSQL ou via algum pacote como RpgSQL ou RPostgreSQL. Atenciosamente Roney

3 2

Mapa de Bairros do Rio de Janeiro
by Sandro Sperandei 14 Ago '12

14 Ago '12

Srs, Estou tentando encontrar um mapa que tenha os bairros do Rio de Janeiro em um formato utilizável pelo maptools e que seja o mais atual possível. Os que poderiam existir no endereço do IBGE (http://www.ibge.gov.br/home/download/geociencias.shtm) estão vindo com a mensagem “Final inesperado do arquivo” quando o arquivo .zip é aberto. Alguém pode me ajudar? Att Sandro Sperandei

2 3

Cálculo entre todos os pares de distâncias em um objeto espacial
by ASANTOS 14 Ago '12

14 Ago '12

Boa noite pessoal, Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz: ### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662) ##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente .... lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda ##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long)) Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso, Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

4 4

Teste - mudança de email , desconsiderar
by marcelo claro de souza 14 Ago '12

14 Ago '12

1 0

Encontrar valores duplicados por linha em um data.frame
by Guilherme Gerhardt Mazzochini 14 Ago '12

14 Ago '12

Olá pessoal, Eu possuo um banco de dados com com quase 900 mil linhas e 23 colunas. Segue um exemplo simplificado abaixo: set.seed(321) v<-data.frame(rbind(sample(1:3,10,replace = T), sample(1:3,10,replace = T), sample(1:3,10,replace = T), sample(1:3,10,replace = T), rep(1,10))) > v X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 1 3 1 2 2 3 2 1 2 3 2 2 1 1 2 3 2 3 3 1 1 1 3 1 1 2 2 2 3 2 1 3 1 4 2 3 3 2 3 1 3 3 2 1 5 1 1 1 1 1 1 1 1 1 1 Para cada linha eu quero determinar qual a coluna do primeiro valor que não é duplicado do primeiro elemento da linha (coluna 1) e salvar em um vetor o numero da coluna do primeiro elemento não duplicado. Assim, > which(v[1,]!=v[1,1])[1] [1] 2 Ou seja, a segunda coluna possui o primeiro elemento diferente do elemento na primeira coluna. Eu tentei fazer um loop para fazer isso para todas as linhas, mas devido ao grande número de linhas está demorando muito tempo. x<-numeric(nrow(v)) for(i in 1:nrow(v)){ x[i]<-which(v[i,]!=v[i,1])[1] } > x [1] 2 3 3 2 NA # o valor NA é porque na quinta linha todos os valores são iguais. Alguém tem idéia de como eu posso fazer isso de forma mais rápida?? Obrigado antecipadamente, Guilherme

2 1

Re: [R-br] Duvida no TAPPLY
by Ivan Bezerra Allaman 14 Ago '12

14 Ago '12

?rowSums, ?colSums \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] Eliminar outliers de uma planilha.
by Ivan Bezerra Allaman 14 Ago '12

14 Ago '12

Uma alternativa seria: remove_outliers <- function(x, na.rm = TRUE, ...) { qnt <- quantile(x, probs=c(.25, .75), na.rm = na.rm, ...) H <- 2.5 * IQR(x, na.rm = na.rm) y <- x y[x < (qnt[1] - H)] <- NA y[x > (qnt[2] + H)] <- NA y } (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}

5 7

microdados censo IBGE 2010
by Manoel Nascimento 14 Ago '12

14 Ago '12

Olá! Sou novo na lista e no uso do R, ainda estou me adaptando e entendendo como funciona, e já apareceram várias dúvidas. Estou precisando ler os microdados do Censo IBGE 2010. Pelo que vi, os pacotes "IBGEPesq" e "dicionariosIBGE" ainda não foram atualizados para nos fornecer os dicionários para o censo, então as diversas funções que eles fornecem ainda não servem para o censo. Por isso, fiz o seguinte (podem dar risada que é coisa de novato mesmo): 1) Fiz um script para cada conjunto de dados, que é mais ou menos o seguinte (mudando alguns dados para cada conjunto de dados -- domicílios, mortalidade, pessoas, migrações): initial.dir<-getwd() setwd("/diretorio/local/IBGE/MicrodadosCenso2010") sink("DomiciliosCenso2010.out") DomicilioAC2010 <- read.fwf("/diretorio/local/IBGE/MicrodadosCenso2010/AC/Dados/Domicilios/Amostra_Domicilios_12.txt", widths = c(2, 5, 13, 8, 3, 1, 2, 3, 2, 1, 2, 2, 1, 6, 4, 1, 2, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 7, 5, 6, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)) [A mesma sintaxe é repetida para os microdados de domicílio de cada estado.] save.image("~/Documentos/cidades/IBGE/MicrodadosCenso2010/Censo2010Domicilios.RData") setwd(initial.dir) Feito deste modo, o script demora um pouco para terminar, mas funciona direitinho. Todas as tabelas são construídas de acordo com as informações que forneci. 2) Tentei fazer este script funcionar usando o "read.fwf" com a variável "col.names" para inserir os cabeçalhos de cada variável, mas ele dá um erro. (Terei que refazer mais tarde o script com esta variável para colar aqui o erro que dá.) Por isto, consegui construir as tabelas, mas elas estão sem cabeçalhos. Pergunto: há algum modo de inserir os cabeçalhos nas tabelas sem ser pela variável "col.names" do comando "read.fwf"? 3) Alguns números constantes nas tabelas são códigos, a serem substituídos por referências específicas. (Exs.: na coluna "Unidades da Federação", para a região Norte, 11=Rondônia, 12=Acre, 13=Amazonas, 14=Roraima, 15=Pará, 16=Amapá, 17=Tocantins; cada município tem um código específico; etc..) A única forma que tenho de saber se as tabelas que construí estão corretas é substituir os números pelas referências. Há algum meio simples de fazer isto? Obrigado pela atenção e paciência. Atenciosamente, Manoel Nascimento

1 0

gráfico de barras
by Vitor Aguiar 13 Ago '12

13 Ago '12

Boa tarde, gostaria de saber como fazer um simples gráfico de barras com o dataframe abaixo, onde as colunas 2 e 3 seriam barras consecutivas, com o nome na coluna 1 para esse grupo de 2 barras. Names NAgel NAcap 1 PentaE 0.14970140 0.09954859 2 D18S51 0.06891474 0.05968163 3 D21S11 0.13054984 0.11116655 4 TH01 0.26187191 0.07833215 5 D3S1358 0.06763797 0.05561891 6 FGA 0.08232084 0.08612497 Muito Obrigado, Vitor

3 4

Re: [R-br] Operação analitica de matrizes.
by Ivan Bezerra Allaman 13 Ago '12

13 Ago '12

Talvez "Ryacas " package? \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] Operação analitica de matrizes.
by Ivan Bezerra Allaman 13 Ago '12

13 Ago '12

Uma googlada iria lhe mostrar o seguinte link: "matrix algebra r" http://personality-project.org/r/sem.appendix.1.pdf (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 3

Operação analitica de matrizes.
by Julia Maria 13 Ago '12

13 Ago '12

Bom dia, será que alguem me dá uma luz? Estava querendo resolver algumas operações analiticamente no R, envolvendo Matrizes. Por exemplo, multiplicar uma matriz numerica 2x2 por uma outra 2x2 de elementos, a, b, c, d. O R é capaz de fazer isso? Obrigada

1 0

Regressão
by kaue veras 13 Ago '12

13 Ago '12

Bom dia caros colegas, Poderiam me dar uma ajuda? Eu gostaria de realizar dois tipos de regressões, linear e múltipla, podem me orientar como posso realiza-las? Através de algum pacote ou função? Acredito que não seja muito complicado, até já realizei uma vez, mas não estou me lembrando como fiz. Desde já agradeço. Atenciosamente, Kauê P. Veras da Cunha Estatística / 8º Período - UERJ MSN: kaueveras(a)hotmail.com E-mail: kaueveras(a)gmail.com Cel: 8254-9601

3 3

Re: [R-br] abrir mapa shp
by andrebvs 12 Ago '12

12 Ago '12

1 0

Eliminar linhas repetidas em um data.frame
by andrebvs 10 Ago '12

10 Ago '12

4 7

erro histograma
by Vanúcia Schumacher 10 Ago '12

10 Ago '12

Allaman e Benilton Carvalho,Problema resolvido, muito obrigado pelas contribuições. O problema era na verdade simples, os valores dos dados não respondiam ao valor passado para o script em relação ao eixo xsegue correção: hist(dados[[1]], seq(0,115,5),prob=TRUE, xlab="chuva (mm/dia)",ylab="frequência", main="", cex.lab=1.6, cex.sub=3,cex.axis=1.9,cex.main=6) Vanúcia SchumacherCurso de graduação em meteorologia - UFPELBolsista do Programa de Educação Tutorial - PET

1 0

erro histograma
by Vanúcia Schumacher 09 Ago '12

09 Ago '12

Olá Pessoal, Gostaria de uma ajuda com um erro que aparece ao tentar rodar um histograma no RErro: some 'x' not counted; maybe 'breaks' do not span range of 'x' Parece que há um problema com o tamanho dos meus dadosestou usando um vetor com mais de 1000 dados Alguém poderia ajudar? Valeu, hist(dados[[1]], seq(0, 20, 0.5), prob=TRUE, xlab="chuva (mm/dia)",ylab="frequência", main="", cex.lab=1.6, cex.sub=3,cex.axis=3,cex.main=6) Vanúcia SchumacherCurso de graduação em meteorologia - UFPELBolsista do Programa de Educação Tutorial - PET

2 1

Re: [R-br] erro histograma
by Ivan Bezerra Allaman 09 Ago '12

09 Ago '12

Imprima o seu objeto dados[[1]] e veja se realmente é o que vc quer! Nos envie uma parte dos seus dados ou CMR para ser possível ajudá-la. (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}

1 0

Formar grupos de linhas variáveis com linha fixa...
by andrebvs 09 Ago '12

09 Ago '12

1 0

Re: [R-br] Erro nls: matriz gradiente singular com estimativas de parâmetros iniciais
by Ivan Bezerra Allaman 09 Ago '12

09 Ago '12

Se você indicou exatamente os valores não deveria dar erro!! Perceba que primeiro vc nos informa o seguinte modelo: y = K*x(^a)*w(^b)/(z^c) e depois na função você informa outro modelo: 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)) Provavelmente o erro deve ser este!! (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}

4 7

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

2 2

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

Re: [R-br] Library Bayesthres
by Ivan Bezerra Allaman 07 Ago '12

07 Ago '12

Trabalha com medidas repetidas cuja a variável é do tipo multinomial ordinal? Abraço! \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}

2 1

Library Bayesthres
by Fabio Mathias Corrêa 06 Ago '12

06 Ago '12

Boa noite pessoal! Acabei de submeter ao CRAN a library Bayesthres. A library possui 3 algoritmos distintos para análises de dados categorizados, com 3 ou mais categorias. É uma library para modelos mistos, com inserção da fórmula idêntica a função lmer (library lme4). Ela possui 3 algoritmos implementados, sendo que você pode escolher entre a distribuição t-Student ou Gaussiana para a distribuição da variável latente. A versão deverá estar disponível nos repositórios da CRAN amanhã, mas para quem quiser dar uma olhada, pode ir no site: ftp://cran.r-project.org/incoming/ Esta não é a versão final, já estou trabalhando em uma versão com diretivas OpenMP e/ou matrizes sparsas para melhorar o processo de amostragem. Críticas, sugestões ou elogios será bem aceitos. 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

Ifelse para duas variáveis simultânemante
by Edson Lira 06 Ago '12

06 Ago '12

Boa tarde, tenho dificuldades em programação, tenho essa rotina abaixo, quero chamar de irmão aqueles nomes que se repetem em nome e nome1, caso contrário quero chama-lo de amigos. Estou tentando usar o ifelse para testar para duas variáveis simultâneamente, não descobri qual o meu erro, alguém tem alguma sugestão? nome<-c("joshua","jose","maria","maria","edinho","raimundo","raimundo","maria","jose") nome1<-c("jose","jose","maria","maria","edson","raimundo","raimundo","maria","jose") idade<-c(18,28,19,29,18,25,19,25,20) dados<-data.frame(nome,nome1,idade) dados$brother<-ifelse(dados$nome="jose"& ifelse(dados$nome1="jose", "irmãos","amigo")) [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 2

Contagem de combinações em um conjunto de dados
by Luciano F. Sgarbi 06 Ago '12

06 Ago '12

#Prezados, #Estou com um problema, e talvez alguém saiba a resposta. #Vejam o exemplo: #--------------------------------------------------------- mat<-matrix(sample(1:2,200,r=T),100,2) plot(mat,pch=16) #--------------------------------------------------------- #Nesse caso eu tenho 100 observações, mas o número de combinações(ou pontos) é apenas 4. #Alguém sabe se existe algum comando para contar esse número de combinações? #Desde já muito obrigado. #Abraços, #Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

3 4

Dúvidas sobre uma função
by Matheus Mendes 03 Ago '12

03 Ago '12

Bom dia... Estou com um trabalho a ser feito e gostaria de solicitar ajuda. Tenho uma planilha com a seguinte estrutura: Primeiro Segundo Terceiro Quarto 1 2 4 4 3 6 3 1 7 4 8 9 1 2 2 9 Preciso identificar quais as linhas que recebem os números 1, 2, 3 e 4, não importando a ordem, ou seja, se é primeiro, segundo, ... Estou com dúvidas para fazer uma função para este fim... Grato Matheus Henrique Silveira Mendes Engenheiro Agrônomo - UFLA Mestre em Genética e Melhoramento de Plantas UFLA Doutorando em Genética e Melhoramento de Plantas UFLA (35)91243078 (TIM) (35)98871518 (VIVO)

5 4

Combinar colunas de um banco de dados...
by andrebvs 03 Ago '12

03 Ago '12

5 16

Simulação Gama e weibull
by Simone toni ruiz corrêa 03 Ago '12

03 Ago '12

Boa noite, Gostaria de, primeiramente, parabenizar pela iniciativa desta lista!!! Sou "novata" no R e estou trabalhando com simulação de temperatura média do ar e radiação solar usando as *distribuições gama e weibull*. As simulações estão ok, mas para descrição em Material e métodos de minha Tese, necessito saber *qual a metodologia usada pelo R nestas simulações*. Por exemplo, se o R utiliza a simulação inversa ou direta e se há artigo referencia neste assunto. Neste mesmo assunto, quando verifiquei a *aderencia* (teste de kolmogorov-smirnov)* *da dist. weibull aos dados de radiação solar (série diária), muitos valores apresentaram a mensagem *optimization failled*. Alguém saberia me explicar o que acontece? Agradeço antecipadamente, Simone

4 6

write.csv
by Edwin Isensee 03 Ago '12

03 Ago '12

Prezados, Leio um arquivo texto com o comando: pingdata = read.csv("/Users/edwin/Desktop/PING_DATA_20120306000000.csv",header=TRUE, sep="|") Porém, quando eu escrevo o data frame para um arquivo com o comando: write.csv(pingdata, file="/Users/edwin/Desktop/xxPING_DATA_20120307000000.csv", append=FALSE, quote=FALSE,row.names=FALSE,col.names=FALSE) na primeira linha do arquivo o R coloca V1. Como faço para evitar a criação desta primeira linha com "V1"? Att, Edwin

2 1

Re: [R-br] Inserir Equação e Ajuste de Dados
by Ivan Bezerra Allaman 03 Ago '12

03 Ago '12

http://lmgtfy.com/?q=modelos+n%C3%A3o-lineares+%22r-br%22 \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

Inserir Equação e Ajuste de Dados
by Gustavo Dias Azevedo 03 Ago '12

03 Ago '12

Prezados, Estou trabalhando com ajuste de dados para uma equação e gostaria de saber se é possível, utilizando o R, saber se meus dados se ajustam a um modelo não-linear de equação e como eu faço isso. Obrigado,

1 0

Join my network on LinkedIn
by Julio Araujo via LinkedIn 03 Ago '12

03 Ago '12

LinkedIn ------------ Julio Araujo requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Julio Araujo http://www.linkedin.com/e/-33g73l-h5ev5rfp-4y/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Julio Araujo http://www.linkedin.com/e/-33g73l-h5ev5rfp-4y/rso/197235503/pw-P/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-h5ev5rfp-4y/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Valores ajustados de um modelo ARIMA(0,2,1)
by Fatima do Nascimento Silva 02 Ago '12

02 Ago '12

Olá pessoal, Bom dia! Fiz um ajuste de um modelo ARIMA(0,2,1) para uma determinada série temporal. Agora estou precisando obter o vetor de valores que foram ajustados pelo modelo, para que eu possa criar um gráfico da série observada versus valores ajustados. O comando que utilizei para ajustar o modelo foi: (mdl=arima(x,order=c(0,2,1),include.mean=F)) Alguém saberia que comando utilizo no software R para obter o vetor de valores que foram ajustados pelo modelo? Obrigada. Atenciosamente, Fátima Nascimento Programa de Pós-Graduação em Ciência e Engª de Petróleo-PPGCEP/UFRN Fone: (84)9936-9725

4 5

Re: [R-br] barplot + lines com dois eixos Y
by Marcelo Claro de Souza 02 Ago '12

02 Ago '12

Pode crer..hehe Eu tb não me atentei a esse detalhe. Vou mandar um cc p a lista :-) Abraço Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil ________________________________ De: Augusto Ribas <ribas.aca(a)gmail.com> Para: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> Enviadas: Quarta-feira, 1 de Agosto de 2012 21:09 Assunto: Re: [R-br] barplot + lines com dois eixos Y Nossa agora que vi, a gente tava trocando e-mail pessoais. Tava respondendo e nem vi. Era bom deixa as resposta la na lista, pra aparecer no nabble, ai se alguem procurar via google ou sei la ja acha como fazer. Eu sempre acho muitas soluções assim :) Em 31 de julho de 2012 23:09, Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> escreveu: > Olá Augusto. > Agradeço de coração todas as dicas que vc me deu. > Finalmente deu certo a figura que eu queria fazer. > > Abraço. > > Marcelo Claro de Souza > Biologist, PhD student in Plant Biology > Institute of Bioscience - UNESP, Brazil > > > ________________________________ > De: Augusto Ribas <ribas.aca(a)gmail.com> > Para: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> > Enviadas: Terça-feira, 31 de Julho de 2012 20:26 > > Assunto: Re: [R-br] barplot + lines com dois eixos Y > > x=c(seq(1:12)) > rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras > tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas > tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas > tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas > > > #o que vc vai ter que fazer é mudar as margens > #vc faz isso com um argumento dentro do par, o argumento mar > #basicamente é a distancia das coisas assim: c(bottom, left, top, right) > #os valores padrao são: c(5, 4, 4, 2) + 0.1., so subir o 2 pra mais > linhas pra caber o texto > #da ?par e da uma olhada na documentação > > par(mar=c(5, 4, 4, 5) + 0.1) > barras<-barplot(rainfall,axes=F,ylab="Nome do eixo esquerdo aqui") > axis(2) > axis(1,col="red") > > par(new=TRUE,mar=c(5, 4, 4, 5) + 0.1) > plot(tmed~barras,type="b",ylim=c(5,35),xlim=c(0.2,14.4),frame.plot=F,axes=F,ylab=NA) > #note que eu coloquei ylab=NA aqui em cima pq se vc escrever algo aqui > vai sobrepor > #com o que ja deixou de leganda no barplot() > points(tmin~barras,type="b",lty=2) > points(tmax~barras,type="b",lty=2) > axis(4) > #vc vai usar o comando mtext, que é pra escrever texto fora da area do plot. > #de ?mtext e veja as muitas possibilidades > #basicamente eu falei o lado no side, e a distancia do plot, line. > mtext(text="Nome do eixo direito aqui", side = 4,line = 3) > axis(1) > > #Mas lembre-se que existem mil jeitos de fazer as mesma coisa > #Talvez essa não seja a melhor forma, mas espero que possa resolver > seus problemas > > Em 31 de julho de 2012 08:50, Marcelo Claro de Souza > <marcelo_claro(a)yahoo.com.br> escreveu: >> Augusto, >> Muito obrigado. Vc me ajudou muito. >> Uma última dúvida. >> Para alterar a legenda do eixo Y do lado esquerdo eu uso o comando ylab, >> mas >> o mesmo não funciona para o segundo eixo Y. Qual comando devo usar para >> inserir a legenda do lado direito? >> Abraço. >> >> Marcelo Claro de Souza >> Biologist, PhD student in Plant Biology >> Institute of Bioscience - UNESP, Brazil >> >> >> ________________________________ >> De: Augusto Ribas <ribas.aca(a)gmail.com> >> Para: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> >> Enviadas: Segunda-feira, 30 de Julho de 2012 21:07 >> >> Assunto: Re: [R-br] barplot + lines com dois eixos Y >> >> So mais uma coisa, note que os eixos vão ficar bagunçados: >> >> barras<-barplot(rainfall,axes=F) >> axis(2) >> #vou colocar aqui em vermelho o eixo x do grafico de barras pra gente >> ver como fica em relação ao outro >> axis(1,col="red") >> par(new=TRUE) >> >> plot(tmed~barras,type="b",ylim=c(5,35),frame.plot=F,axes=F) >> points(tmin~barras,type="b",lty=2) >> points(tmax~barras,type="b",lty=2) >> axis(4) >> axis(1) >> >> A bolinha era pra ficar no centro da barra, mas não tava ficando. >> Isso por causa da escala de x, que do plot ficou diferente do barblot >> pra dar certo, bolinhas no meio da barrinha. >> dai tem que arrumar o limite do eixo x no plot: >> >> #denovo >> barras<-barplot(rainfall,axes=F) >> axis(2) >> #denovo vermelho pra gente ver >> axis(1,col="red") >> par(new=TRUE) >> >> >> >> plot(tmed~barras,type="b",ylim=c(5,35),xlim=c(0.2,14.4),frame.plot=F,axes=F) >> points(tmin~barras,type="b",lty=2) >> points(tmax~barras,type="b",lty=2) >> axis(4) >> #agora nao ficou bagunçado >> axis(1) >> >> #isso por causa do jeito que o barplot constroi o eixo X é diferente do >> plot >> #mas eu não entendo perfeitamente a diferença pra explicar >> #espero estar mais ajudando que atrapalhando agora. >> >> >> Em 30 de julho de 2012 19:28, Augusto Ribas <ribas.aca(a)gmail.com> >> escreveu: >>> Humm >>> >>> Veja se isso te ajuda: >>> >>> #Vai faze o grafico de barras >>> barplot(rainfall,axes=F) >>> #o eixo dele >>> axis(2) >>> >>> #vai recomeça o plot, na mesma janela >>> par(new=TRUE) >>> >>> #um plot normal de linha >>> plot(tmed~barras,type="b",ylim=c(5,35),frame.plot=F,axes=F) >>> #as outras 2 linhas >>> points(tmin~barras,type="b",lty=2) >>> points(tmax~barras,type="b",lty=2) >>> #vai fazer o eixo do outro lado (direito) >>> axis(4) >>> >>> #note que eu fiz os graficos sem eixos, e os adicionei a parte com o >>> comando axis, >>> #com o axis vc pode faze o eixo da forma que preferir. >>> #ai da pra controla cada grafico independente quanto a cores, tamanhos >>> eixos tudo >>> #se quiser mudar a escala, de uma olhada no comando ?scale(), mas acho >>> que talvez a solução acima seja mais >>> #o que vc deseja >>> >>> Em 30 de julho de 2012 17:27, Marcelo Claro de Souza >>> <marcelo_claro(a)yahoo.com.br> escreveu: >>>> Olá Augusto, >>>> Muito obrigado pela ajuda. >>>> Eu preciso que os dois eixos Y sejam independentes (escalas diferentes), >>>> de >>>> modo que as linhas fiquem sobre as barras. >>>> Tem como fazer isso? >>>> Muito obrigado. >>>> Abraço. >>>> >>>> >>>> Marcelo Claro de Souza >>>> Biologist, PhD student in Plant Biology >>>> Institute of Bioscience - UNESP, Brazil >>>> >>>> >>>> ________________________________ >>>> De: Augusto Ribas <ribas.aca(a)gmail.com> >>>> Para: r-br(a)listas.c3sl.ufpr.br; Marcelo Claro de Souza >>>> <marcelo_claro(a)yahoo.com.br> >>>> Enviadas: Segunda-feira, 30 de Julho de 2012 17:40 >>>> Assunto: Re: [R-br] barplot + lines com dois eixos Y >>>> >>>> Bom, de inicio acredito que vai ser algo assim: >>>> >>>> barras<-barplot(rainfall,ylim=c(0,35)) >>>> points(tmed~barras,type="l") >>>> points(tmin~barras,type="l",lty=2) >>>> points(tmax~barras,type="l",lty=2) >>>> axis(4) >>>> >>>> #vc vai salvar a posição das barras no vertor barras >>>> #ai depois note que vc esta fazendo as linhas >>>> # e por ultimo use o comando axes pra fazer o segundo ou mais eixos, >>>> se quiser algo no eixo X tb >>>> #leia a documentação pra deixa do jeito que gostar ?axis >>>> >>>> >>>> #agora não sei se entra nas suas necessidades mas >>>> #de uma olhada nesse pacote e nos graficos que da pra fazer com ele >>>> #é exatamente o tipo de dado que vc ilustrou >>>> #por exemplo: >>>> >>>> install.packages("climatol") >>>> require(climatol) >>>> >>>> data(datcli) >>>> #exemplo dele >>>> datcli >>>> #seus dados >>>> t(data.frame(rainfall,tmax,tmin,tmed)) >>>> >>>> #exemplo do comando >>>> diagwl(datcli,est="Example station",alt=100,per="1961-90",mlab="en") >>>> #seus dados de exemplo >>>> diagwl(t(data.frame(rainfall,tmax,tmin,tmed)),est="Example >>>> station",alt=100,per="1961-90",mlab="en") >>>> >>>> >>>> Em 30 de julho de 2012 15:18, Marcelo Claro de Souza >>>> <marcelo_claro(a)yahoo.com.br> escreveu: >>>>> Preciso plotar um gráfico de variáveis climáticas sendo o eixo Y >>>>> esquerdo >>>>> referente à precipitação e o eixo Y direito referente às temperaturas >>>>> min, >>>>> med e max. >>>>> A precipitação deve ser apresentada em barras e as temperaturas em >>>>> linhas. >>>>> Como posso fazer isso? >>>>> >>>>> x=c(seq(1:12)) >>>>> >>>>> >>>>> >>>>> rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras >>>>> >>>>> >>>>> >>>>> tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas >>>>> >>>>> >>>>> >>>>> tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas >>>>> >>>>> >>>>> >>>>> tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas >>>>> >>>>> Muito obrigado pela ajuda. >>>>> Abraço. >>>>> >>>>> Marcelo Claro de Souza >>>>> Biologist, PhD student in Plant Biology >>>>> Institute of Bioscience - UNESP, Brazil >>>>> >>>>> >>>>> _______________________________________________ >>>>> 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. >>>> >>>> >>>> >>>> -- >>>> Grato >>>> Augusto C. A. Ribas >>>> >>>> Site Pessoal: http://augustoribas.heliohost.org >>>> Lattes: http://lattes.cnpq.br/7355685961127056 >>>> >>>> >>> >>> >>> >>> -- >>> Grato >>> Augusto C. A. Ribas >>> >>> Site Pessoal: http://augustoribas.heliohost.org >>> Lattes: http://lattes.cnpq.br/7355685961127056 >> >> >> >> -- >> Grato >> Augusto C. A. Ribas >> >> Site Pessoal: http://augustoribas.heliohost.org >> Lattes: http://lattes.cnpq.br/7355685961127056 >> >> > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://augustoribas.heliohost.org > Lattes: http://lattes.cnpq.br/7355685961127056 > > -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Agregar fatores
by Sérgio Henrique almeida da silva ju 01 Ago '12

01 Ago '12

Qual a melhor forma de agregar fatores? Eu tenho um banco com as variáveis cidade, fumo (1 = sim, 2 = não) e diversas outras, porém eu gostaria de criar um outro banco só com a cidade e a frequência de quem fuma ou não. Exemplo cidade sim não 1 aracaju 163 1843 2 belem 185 1824 3 belo horizonte 280 1727 4 boa vista 191 1828 5 brasilia 230 1780 6 campo grande 228 1782 Eu conseguir fazer de forma não muito elegante: tabela <- aggregate(dados$fumante, by = list(dados$CIDADE),FUN = function(x) c(table(x))) cidade <- data.frame(tabela) Como fazer de forma mais elegante? Outra pergunta, gostaria de criar dento desse banco (cidade) uma variável que é o % de sim, como fazer? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 2

barplot + lines com dois eixos Y
by Marcelo Claro de Souza 30 Jul '12

30 Jul '12

Preciso plotar um gráfico de variáveis climáticas sendo o eixo Y esquerdo referente à precipitação e o eixo Y direito referente às temperaturas min, med e max. A precipitação deve ser apresentada em barras e as temperaturas em linhas. Como posso fazer isso? x=c(seq(1:12)) rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas Muito obrigado pela ajuda. Abraço. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Importar dados EpiData
by Pedro Emmanuel Alvarenga Americano do Brasil 30 Jul '12

30 Jul '12

Amigos de R, Há um software de entrada de dados que gosto muito chamado EpidData. ( www.epidata.dk) O Epidata utilizava até o ano passado o mesmo formato de dados que o antigo EpiInfo 6.4 em DOS. Os arquivos .rec. Então, para importar dados basta: library('foreign') dados <- read.epiinfo(...) Acontece que desde o ano passado o EpiData 3.1 esta sendo abandonado e substituído por um novo EpiData que está tentando se adequar a boas praticas clínicas com acesso hierarquizado aos dados entre outras coisas. O novo formato de arquivos são versoes XML do velho .rec... os epx e epz Alguem ja conseguiu importar dados nesse formato a partir do R? Isso está causando algumas dores de cabeça, pois não estou conseguindo migrar do velho para o novo EpiData devido a essa limitação! Cheguei mandar uma mensagem numa lista da turma do core do R mas o moderador rejeita as minhas mensagens. :-( Abraço, 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

2 2

Encontrar caminho de um arquivo no computador
by Luis Iván Ortiz Valencia 27 Jul '12

27 Jul '12

Prezados Gostaria de saber se é possível no R encontrar o caminho de um arquivo específico no computador. Por exemplo, existe um arquivo indicadores.pdf na pasta C:\ivan\R, existe uma função no R que responda com o caminho "C:\ivan\R" dando como entrada indicadores.pdf? atte IVAN

5 9

[Dúvida] Salvando data.frame em DB
by Pedro Rafael 27 Jul '12

27 Jul '12

Pessoal estou usando o comando dbWriteTable para salvar em um banco de dados PostgreSQL. O problema que mesmo colocando o argumento row.names = FALSE ele está salvando o banco com um ID. Só que eu gostaria que ele não criasse o ID pois ele já está criado nas tabelas do BD. Alguem sabe como resolver esse problema? [ ], Pedro Rafael Diniz Marinho.

1 0

regressão polinomial
by Tito Conte 26 Jul '12

26 Jul '12

Pessoal tenho o seginte conjunto de dados: f=c(-26.0000,-25.9862,-25.9343,-25.8822,-25.8433,-25.8054,-25.7948,-25.7872,-25.7668,-25.7284,-25.7015,-25.6282,-25.4612,-25.3016,-25.2564,-25.2412,-25.2412,-25.2232,-25.0869,-25.0000,-25.0000,-24.9856,-24.9397,-24.8976,-24.8533,-24.6587,-24.6373,-24.5740,-24.5406,-24.4379,-24.3934,-24.3628,-24.3277,-24.3042,-24.2972,-24.2973,-24.3020,-24.3068,-24.3068,-24.2956,-24.2920) w=c(-45.8529,-45.8302,-45.7575,-45.6680,-45.6107,-45.5050,-45.4574,-45.4058,-45.3538,-45.2913,-45.2634,-45.2264,-45.1667,-45.0719,-45.0213,-45.0000,-45.0000,-44.9748,-44.8726,-44.8405,-44.8405,-44.8351,-44.8259,-44.7972,-44.7764,-44.6867,-44.6674,-44.5912,-44.5587,-44.4927,-44.4515,-44.4115,-44.3314,-44.2664,-44.2171,-44.1729,-44.1288,-44.0820,-44.0509,-44.0118,-44.0000) se vocês plotarem eles, obteram uma curva, gostaria de obter a função desta curva ou uma regressão polinomial de grau elevado dela já tentei usar a função lm e não deu um bom resultado, quando chego em graus elevados e calculo utilizando os coeficientes obtidos os valores dão 10 vezes maiores (não sei bem porque) mas quando ploto ela aparece bonitinha no lugar que deveria aprecer em outro caso tentei na raça mesmo por multiplicação matricial mas também deu errado pois os valores da multiplicação da inversa são muito pequenos alguém pode me dar uma sugestão ou uma luz do que eu estou fazendo de errado? Tito Conte

4 14