- R-br - DINF UFPR

Conflito de resultados + plotagem
by Rodrigo Plei 19 Mai '11

19 Mai '11

Prezados, Utilizei a função "nls" do R para estimar as constantes de uma função potencial ajustada aos meus dados de comprimento e peso de lulas utilizando o método dos mínimos quadrados: m <- nls(y~a*x^b, data = ds, start = c(a=.1, b=2)) Obtive as seguintes constantes: a = 0.0003967 , b = 2.3933797 e o R2 = 0.9681 Ocorre que o mesmo ajuste no Excel 2007 resulta em a = 0.0002 , b = 2.6097 e R2 = 0.9819 O mesmo ajuste no Statistica resulta em constantes e R2 idênticas aos encontrados no R. Não deveria ter dado o mesmo resultado nos três programas? Não sei em qual resultado confiar. Também gostaria de saber como faço para plotar a curva dessa função potencial sobreposta a minha dispersão de dados (peso vs. comprimeto). Agradeço desde logo. 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

4 4

R studio - linux
by Marcelo Claro de Souza 19 Mai '11

19 Mai '11

Eu costumo abrir um terminal dentro da pasta onde estão meus arquivos e utilizo o comando write.csv2('nome do arquivo.csv') para salvar minhas planilhas dentro da própria pasta (não sei salvar os arquivos fora da pasta onde abro o terminal). Entretanto, ao atualizar o ubuntu para a versão 11.04 esse comando não está funcionando, o R roda o comando mas não salva o arquivo. Gostaria de saber como possa fazer isso pelo R studio como uma alternativa, ou alguma outra solução p o R clássico. Muito obrigado.

3 2

fatorial fracionado em R
by jonas bodini alonso 19 Mai '11

19 Mai '11

Bom dia aos membros da lista Alguem sabe uma função do R que trabalhe com experimentos fatoriais fracionados? Se sim como faço para fazer para escolher a fração de interesse? Cordialmente Jonas

2 2

Desdobramento da Interação! Função split()
by Thiago De paula protásio 19 Mai '11

19 Mai '11

Pessoal; Recorro novamente a vcs para pedir um auxílio. Estou querendo utilizar a função split para desdobrar o efeito da interação entre fatores.Contudo, não estou obtendo sucesso. Meu experimento é um fatorial completo 3 x 2 instalado em esquema fatorial, ou seja, tenho 3 clones plantados em 2 locais distintos. Alguém já trabalhou com essa função ou me aconselha a utilizar outra?Posso enviar alguns dados para elucidar o problema. Muito obrigado!! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Erro ao executar arquivo Rpatch.sh!
by Ivan Bezerra Allaman 18 Mai '11

18 Mai '11

Bom dia senhores!! Estou tentando atualizar a versão do R por meio do arquivo Rpatch.sh que o professor PJ disponibilizou em um tópico anterior. Está ocorrendo o seguinte erro: ivan-Aspire-5920 Área de Trabalho # chmod 777 Rpatch.sh ivan-Aspire-5920 Área de Trabalho # ./Rpatch.sh ./Rpatch.sh: linha 1: -------------------------------------------------------: comando não encontrado ============================================= ============================================= ARQUIVOS FONTE R - VERSAO 2-13 ============================================= ============================================= ============================================= ============================================= ARQUIVOS FONTES R - VERSAO 2-13 ============================================= ============================================= Nao instalando dependencias, o R ou alguns pacotes podem nao compilar dependendo do que voce tem instalado no seu sistema. Para instalar dependencias automaticamente utilize opcao --deps Veja opcoes com --help ============================================= ============================================= ARQUIVOS FONTE R-VERSAO 2-13 ============================================= ============================================= ./Rpatch.sh: linha 70: svn: comando não encontrado R version 2.11.1 (2010-05-31) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R é um software livre e vem sem GARANTIA ALGUMA. Você pode redistribuí-lo sob certas circunstâncias. Digite 'license()' ou 'licence()' para detalhes de distribuição. R é um projeto colaborativo com muitos contribuidores. Digite 'contributors()' para obter mais informações e 'citation()' para saber como citar o R ou pacotes do R em publicações. Digite 'demo()' para demonstrações, 'help()' para o sistema on-line de ajuda, ou 'help.start()' para abrir o sistema de ajuda em HTML no seu navegador. Digite 'q()' para sair do R. Vejam que não atualizou! Eu copiei exatamente o que estava no tópico. Se alguém puder me ajudar eu ficarei grato. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

3 3

Leitura de arquivos
by Emerson Cotta Bodevan 18 Mai '11

18 Mai '11

Prezados. Tenho um arquivo texto com 100 linhas, sendo que as 50 primeiras possuem três colunas e as 50 últimas, quatro colunas. Como faço para ler esse arquivo, separando-o em dois data-frames (Um com 3 colunas e o outro com 4 colunas)? Agradeço a ajuda, -- Emerson

3 4

combinar tipos diferentes de gráficos
by Rodrigo Plei 18 Mai '11

18 Mai '11

Prezados, Como combinar um gráfico de barras com um gráfico de dispersão? Ambas plotagens possuem os mesmos dados para o eixo x. Gostaria que as barras ficassem atreladas ao eixo y principal (esquerdo) e os pontos da dispersão atrelados ao eixo y secundário (direito). Agradeço desde logo. 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

2 2

Re: [R-br] Erro ao executar arquivo Rpatch.sh!
by Ivan Bezerra Allaman 18 Mai '11

18 Mai '11

Resolvi do jeito tradicional. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Usando AIC para seleção de modelo com o comando drop1(). Alguém poderia dar uma olhadinha :)
by Augusto Ribas 18 Mai '11

18 Mai '11

Ola colegas de lista. Eu estava lendo sobre seleção de modelos e fiquei na duvida se estou indo pelo caminho certo no R. Outra coisa, eu não achei exatamente nos livros que li o que eu deveria reportar em um artigo, ao longo destes passos (se estes estiverem fazendo sentido). Eu tentei fazer um exemplo com uma com uma uma variavel resposta e duas variaveis explicativas. Usei o dataset do anscombe que vem no R, usando o primeiro exemplo dele mais uma variavel ao acaso para exemplificar. Aqui vai um codigo de como estava pensando: #Abrindo os dados anscombe x2<-runif(11) #fazendo um plot para olhar panel.hist <- function(x, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(usr[1:2], 0, 1.5) ) h <- hist(x, plot = FALSE) breaks <- h$breaks; nB <- length(breaks) y <- h$counts; y <- y/max(y) rect(breaks[-nB], 0, breaks[-1], y, col="cyan", ...) } panel.cor <- function(x, y, digits=2, prefix="", cex.cor, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- abs(cor(x, y)) txt <- format(c(r, 0.123456789), digits=digits)[1] txt <- paste(prefix, txt, sep="") if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt) text(0.5, 0.5, txt, cex = cex.cor * r+1) } panel.lm = function(x, y, ...) { tmp <- lm(y ~ x, na.action = na.omit) abline(tmp) points(x, y, ...)} #legal, y1 é relacionado com x1 mas nao com x2 olhando esse grafico pairs(cbind(y1=anscombe$y1,x1=anscombe$x1,x2),diag.panel=panel.hist, lower.panel=panel.cor,upper.panel=panel.lm) #agora fazendo um modelo modelo1<-lm(anscombe$y1~anscombe$x1+x2) modelo1 #vendo premissas par(mfrow=c(2,2)) plot(modelo1) #vendo o quanto os parametros ajudam drop1(modelo1) #testando quem ajuda significativamente drop1(modelo1,test="F") #novo modelo que eu vou ficar ja que vi no drop1 que posso tirar o x2 modelo2<-lm(anscombe$y1~anscombe$x1) modelo2 #vendo premissas pro novo modelo par(mfrow=c(2,2)) plot(modelo2) #o que reportar num artigo de todo o processo? #todo esse processo pode ser inverso tb, começando de um modelo simples e indo pra um mais complexo add1(modelo2, ~anscombe$x1+x2,test="F") É isso. Agradeço a paciencia de te -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 1

coordenadas valor máximo num grafico xy
by jose coelho 18 Mai '11

18 Mai '11

Vivam , ainda sou aprendiz no R... mas gostaria de saber como calcular as coordenada do valor máximo de y em função de x para colocar num gráfico xy. Existe alguma função ? max(y) eu sei qual valor máximo do y mas também quero saber qual o valor de x no qual o y atinge o seu máximo. Muito obrigado jmmcoelho

3 3