Abril 2011 - R-br - DINF UFPR

Sweave Mostrar Codigo R
by eder＠leg.ufpr.br 23 Abr '11

23 Abr '11

Pessoal, Preciso que no documento final seja mantido o codigo R bem como os resultados processados por ele, por exemplo o arquivo .Rnw abaixo: \documentclass[a4paper]{article} \usepackage{Sweave} \title{Estimativa via verossimilhança} \author{Éder David Borges da Silva} \begin{document} \maketitle \section{Distribuição} <<echo=TRUE>>= ###---------------------------------------------------------------------------### dbetabinom1 <- function(x, size, prob, theta, log=FALSE) { alpha <- prob / theta beta <- (1 - prob) / theta dens <- lchoose(size,x) + lbeta(x+alpha,size-x+beta) - lbeta(alpha,beta) if (!log) dens <- exp(dens) return(dens) } rbetabinom1 <- function(n, size, prob, theta) { # n é o número de elementos da amostra # size, prob e theta os parâmetros da densidade r <- function(x) sample(0:x[1], size=1, replace=TRUE, prob=dbetabinom1(0:x[1], size=x[1], prob=x[2], theta=x[3], log=FALSE)) if (length(n)>1) n <- length(n) return(new("numeric", sapply(data.frame(rbind(size, prob, theta, rep(1,n))), r))) } ###---------------------------------------------------------------------------### @ \end{document} Em minha cabeça eu lembrava que o echo=TRUE fazia isso, ou seja mantinha o codigo, porem isso não esta acontecendo, quando compilo não aparece nada do codigo R. Onde será que estou errando? Att

3 3

Removendo pontos dentro de um grid
by Alexandre dos Santos 22 Abr '11

22 Abr '11

Bom dia pessoal, Estou criando uma malha de pontos da seguinte forma: p0<- c(15,15)#Ponto inicial malha10<- expand.grid(xp=seq(p0[1], by=10, length.out=20), yp=seq(p0[2], by=10, length.out=20))#Criando um grid de pontos a cada 10 metros plot(malha10)# Visualizando a malha de postos Agora, gostaria das coordenadas dos pontos a cada 20m, à partir de p0, mas gostaria de manipular o objeto já criado malha10, sem ter que fazer um novo objeto expand.grid(... by=20), na verdade, gostaria de criar um novo objeto malha20, mas fazendo malha20=malha10 e após remover as coordenadas dos pontos do objeto malha10 que não satisfaçam a condição , tentei um outer() para descobrir as distancias entre todos os pontos e após impor uma condição =20, mas não estou conseguindo não, alguem poderia dar um help, Obrigado, Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29

3 5

Transformar image{graphics} em plot{graphics}
by Augusto Ribas 22 Abr '11

22 Abr '11

Ola. Estou com uma duvidade aqui. Vamos supor que eu tenha uma matriz assim com 0 e 1(de presença e ausencia) exemplo<-matrix(rep(0,100),nrow=10) while(sum(exemplo)/length(exemplo)<0.1){ x<-sample(x=1:10,1,replace =T) y<-sample(x=1:10,1,replace =T) exemplo[x,y]<-1 } exemplo e represente ela graficamente com o comando image: image(1:10,1:10,exemplo,col=(c(0,1)),xlab="",ylab="") Mas eu gostaria de pegar essa matriz de exemplo e fazer um grafico com o comando plot. Então eu preciso transformar uma matriz de 0 e 1 em coordenadas x e y pra usar com o comando plot para por exemplo preencher com pontinhos o seguinte plot: Mais ou menos assim: Pense nesse grafico: plot(0,0,xlab="",ylab="",xlim=c(0,10),ylim=c(0,10),type="n") abline(h=0.5:9.5,v=0.5:9.5) Em cada local da matriz que tiver um 1 eu queria que tivesse um pontinho no grafico. Ai eu uso o comando points para representar a matriz dentro daquele plot com pontinhos. Mas eu nao estou conseguindo visualizar como transformar uma matriz de 0 e 1 em coordenadas de pontos X e Y. Eu queria arranja um jeito de transformar qualquer matriz de presença e ausencia (0 e 1) em coordenadas para plotar pontinhos com o comando plot. Será que alguem tem uma ideia simples que resolvera todos meus problemas e que eu não consigo visualizar? Bom feriado a todos. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org/ http://dhost.info/augusto/ Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Re: [R-br] Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

Alexandre!! Não sei se entendi direito, mais vai uma primeira tentativa. p0<- c(15,15)#Ponto inicial c0 <- seq(10,100,by=10) malha <- list() for(i in 1:length(c0)){ malha[[i]] <- expand.grid(xp=seq(p0[1], by=c0[[i]], length.out=20), yp=seq(p0[2], by=c0[[i]], length.out=20))#Criando um grid de pontos a cada 10 metros } malha[[2]] plot(malha[[2]]) plot(malha[[3]]) 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

Re: [R-br] Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

Sacola 2!! 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

Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

sacola 1!! 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

Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

Alexandre! Não sei se entendi direito, mais vai uma primeira tentativa aí. p0<- c(15,15)#Ponto inicial c0 <- seq(10,100,by=10) malha <- list() for(i in 1:length(c0)){ malha[[i]] <- expand.grid(xp=seq(p0[1], by=c0[[i]], length.out=20), yp=seq(p0[2], by=c0[[i]], length.out=20))#Criando um grid de pontos a cada 10 metros } malha[[2]] plot(malha[[2]]) plot(malha[[3]]) 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

Atualização do R no Linux
by Andre Oliveira 21 Abr '11

21 Abr '11

Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave? André Oliveira Souza

5 11

É possível aumentar a precisão numerica do R?
by Carolina Pena 21 Abr '11

21 Abr '11

Boa tarde pessoal, Estou trabalhando com operações que exigem uma precisão numerica muito grande. Por exemplo: (1.5692524257254742661/(10^41)) + 6.125125476214521152) No entanto, o R está gerando resultados inconsistentes. Alguém sabe se é possível aumentar a precisão do R? Se ninguém souber como fazer no R, mas souber de algum pacote que posso utilizar em C++ também seria muito útil. Obrigada, Carolina

3 2

Enc: Linux ubuntu
by Luciana Melo 21 Abr '11

21 Abr '11

Prezados, Não recebi minha mensagem , porque ? Grata, Luciana --- Em ter, 19/4/11, Luciana Melo <melo_lucavalcanti(a)yahoo.com.br> escreveu: De: Luciana Melo <melo_lucavalcanti(a)yahoo.com.br> Assunto: Linux ubuntu Para: r-br(a)listas.c3sl.ufpr.br Data: Terça-feira, 19 de Abril de 2011, 18:03 Prezados, Preciso instalar o R-project num sistema linux ubuntu e depois de vários erros na tentativa o sistema me indicou o R Kward. Pergunta: é a mesma coisa? Grata, Luciana Melo

2 1