- R-br - DINF UFPR

Re: [R-br] Coeficiente de Variação! Modelos não lineares!
by Thiago De paula protásio 25 Abr '11

25 Abr '11

Ok Walmes... Compreendo o seu ponto de vista..... Na verdade estou procurando um maneira de comparar modelos polinomiais com modelos não lineares. Não sei se o R² pode ser utilizado para esse fim, pois na obtenção do R² de modelos não lineares utilizamos a SQ total corrigida e geralmente o R² fica maior se compararmos com o método dos mínimos quadrados ordinários. O que vc recomenda? 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

Coeficiente de Variação! Modelos não lineares!
by Thiago De paula protásio 25 Abr '11

25 Abr '11

Pessoal; Gostaria de perguntar a vcs se existe algum comando no R para determinarmos o coeficiente de variação de modelos não lineares. Para modelos lineares temos a função cv.model() do pacote agricolae.........mas não consegui encontrar um comando semelhante para os modelos não lineares. Qualquer ajuda é de grande valia! Muito obrigado! A disposição! 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

Fatorial com reptições perdidas
by Edson Lira 25 Abr '11

25 Abr '11

Boa Páscoa a todos, estou analizando um experimento fatorial como variáveis respostas altura, diâmetro do peito, número de galhos, mdida da copa da árvore, os fatores são tipo de plantio (puro e mista) e aração (arado e aarado). Cada medida tem 3 repetições, o problema é que em alguns casos, a árvore morreu, estou com um impasse, não sei se considero zero, fato que infuenciaria na média ou resposta perdida "missing" (NA). Fiz uma tentativa usando o NA, rodei o modelo com aov, verifiquei os pressupostos, fiz o desdobramento das interações, quando tentei rodar o tukey, o R me deu a mensagem dizendo que x e y não tem a mesma dimensão. Alguem poderia me sugerir o que fazer com a arvore morta? Caso tenha que dropar os NA's como proceder? Edson Lira Estatístico Manaus-Amazonas

2 4

Retornando médias de tabela para tabela no R
by Emmanuel Arnhold 25 Abr '11

25 Abr '11

Pessoal, se eu tenho um objeto no R (tabela lida de um read.table) onde a primeira coluna são fatores (tratamentos) e as demais os dados (numéricos). Qual a maneira mais fácil ou se existe uma função para retornar a média das culunas numéricas dentro de cada categoria do fator. Gostaria de saber se é possível retornar as médias em tabela semelhante a forma original, mas somente com as médias (sem as repetições de tratamentos) dentro de cada categoria de meu fator. Agradeço desde já qualquer dica. Emmanuel

3 3

Problema com sample()
by Alexandre dos Santos 25 Abr '11

25 Abr '11

Bom Dia Pessoal, Tenho um banco de dados com três variáveis parc (número da parcela), mont2(número de insetos), dead2 (mortalidade), com 108 dados, estou tentando fazer uma amostragem com a função sample(), para um n da amostra igual a 30, ate ai tudo bem, mas quando repito a função sample() para um número de vezes superior ao n da amostra começam a aparecer NA nos sorteio, porque isso poderia estar acontecendo, uma vez que estou utilizando a função com reposição, fiz: > nde<-read.table("Tout_NDE.txt",header=T)# Todos os dados > # > dead2<-nde$Mortalite[nde$Temps==2]-nde$Mortalite[nde$Temps==1]# Mortalidade no segundo mês > mont2<-nde$Monticules[nde$Temps==2]-nde$Monticules[nde$Temps==1]# Número de insetos no segundo mês > parc<-1:108 # Número da parcela no campo > mois2<-cbind(parc,mont2,dead2) > # > nparc=30# Tamanho da amostra independente a ser sorteada com sample () > Nsim=31# Numero de simulações da função sample > > RES=NULL > > for(i in 1:Nsim){ + + res=NULL + spl=sample(parc,nparc, replace = TRUE) + resspl=mois2[spl[1:30],] + my.mont=mean(resspl[,2][i]) + my.dead=mean(resspl[,3][i]) + result=rbind(res,c(my.mont,my.dead,i)) + #print(result) + RES=rbind(RES,result) + colnames(RES)=c("monticule","mortalite","simulation") + } > print(RES) monticule mortalite simulation [1,] 5 4.000000 1 [2,] 0 0.000000 2 [3,] 1 12.000000 3 [4,] 1 5.263158 4 [5,] 0 0.000000 5 [6,] 0 0.000000 6 [7,] 2 9.090909 7 [8,] 0 0.000000 8 [9,] 1 0.000000 9 [10,] 0 0.000000 10 [11,] 1 11.764706 11 [12,] 3 6.451613 12 [13,] 3 12.000000 13 [14,] 3 4.545455 14 [15,] 1 3.225806 15 [16,] 0 0.000000 16 [17,] 3 4.545455 17 [18,] 1 0.000000 18 [19,] 5 13.636364 19 [20,] 2 12.121212 20 [21,] 2 6.896552 21 [22,] 1 5.263158 22 [23,] 2 0.000000 23 [24,] 10 8.333333 24 [25,] 5 14.285714 25 [26,] 8 10.000000 26 [27,] 8 10.000000 27 [28,] 3 4.761905 28 [29,] 0 0.000000 29 [30,] 1 3.571429 30 [31,] NA NA 31 > 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 6

(sem assunto)
by Juraci Mendes Moreira 25 Abr '11

25 Abr '11

Estou tentando rodar esta rotina no R mas esta dando este erro e tambem ele não esta estimando b2, poderiam me ajudar? lrp.fit<-nls(cv ~ (b0 + b1*x + b2*I(x^2))*(x<=x0p) +(b0 + b1*x0p + b2*x0p^2)*(x>x0p), start=list(b0=10, b1=7, b2=15, x0p=5), trace=F) Erro em nls(cv ~ (b0 + b1 * x + b2 * I(x^2)) * (x <= x0p) + (b0 + b1 * : gradiente singular

2 1

dificuldade em ajustar um quadratica a um plato de resposta linear
by Juraci Mendes Moreira 25 Abr '11

25 Abr '11

Estou tentando rodar esta rotina no R mas esta dando este erro e tambem ele não esta estimando b2, poderiam me ajudar? lrp.fit<-nls(cv ~ (b0 + b1*x + b2*I(x^2))*(x<=x0p) +(b0 + b1*x0p + b2*x0p^2)*(x>x0p), start=list(b0=10, b1=7, b2=15, x0p=5), trace=F) Erro em nls(cv ~ (b0 + b1 * x + b2 * I(x^2)) * (x <= x0p) + (b0 + b1 * : gradiente singular

1 0

Re: [R-br] Fatorial com reptições perdidas
by Thiago De paula protásio 24 Abr '11

24 Abr '11

Ok Walmes........ Rodei aqui e deu td certo!!! Vou estudar os seus comandos com calma.......... Abraço! 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

1 0

Re: [R-br] Fatorial com reptições perdidas
by Thiago De paula protásio 24 Abr '11

24 Abr '11

Walmes; Eu não estou conseguindo reproduzir o exemplo que vc indicou na página do LEG! Dê uma olhada no erro que apareceu....... Obrigado!!! Este exemplo será muito útil, pois tb tenho que avaliar um experimento com repetições perdidas....... # lista com as matrizes dos contrastes que específicam o desdobramento A dentro de B > > c.X <- sapply(jj, simplify=FALSE, + function(j){ + sapply(1:nrow(p.contr), + function(i){ + c.contr <- contrast(m0, + list(A=p.contr[i,1], B=j), + list(A=p.contr[i,2], B=j)) + c.contr$X + })}) Erro em library(sandwich) : não há nenhum pacote chamado 'sandwich' > c.X Erro: objeto 'c.X' não encontrado > > #------------------------------------------------------------------------------------------ > # usando a glht para desdobrar, cada entrada da lista é um nível de B > > mc.contr <- sapply(c.X, simplify=FALSE, + function(X){ + summary(glht(m0, linfct=t(X))) + }) Erro em is.vector(X) : objeto 'c.X' não encontrado > mc.contr Erro: objeto 'mc.contr' não encontrado > 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

Axelrod's Tournaments
by Augusto Ribas 24 Abr '11

24 Abr '11

Alguem saberia dizer se existe uma função para o R que replica o Axelrod's Tournaments? Ele vem do artigo: Robert Axelrod & William D. Hamilton 1981, "The Evolution of Cooperation", Science 211: 1390–96 Naquele documentario da BBC Nice Guys Finish First falam bastante desse artigo. Eu vi um programa que replica o torneio: http://www2.econ.iastate.edu/tesfatsi/demos/axelrod/axelrodt.htm Mas dei uma procurada e não encontrei nenhuma função pronta no R. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org/ http://dhost.info/augusto/ Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Re: [R-br] [R_STAT] Seleção de Quadrados Médios
by Greice Laureano 23 Abr '11

23 Abr '11

Fernando, (agora mandei para a lista nova :P) O problema é que o comando aov() não aceita esse procedimento... Se eu faço uma Anova do fit ele da a seguinte mensagem de erro: > anova(fit) Erro em UseMethod("anova") : método não aplicável para 'anova' aplicado a um objeto de classe "c('aovlist', 'listof')" 2011/4/23 Fernando Henrique Toledo <fernandohtoledo(a)gmail.com> > > > Greice, bom dia. > > Acesse: > > anova(fit)[,3] > > Isso resultará no vetor dos quadrados médios! a 1ª coluna é dos graus de > liberdade, a segunda das somas de quadrados, 3ª dos quadrados médios ... > > PS: Mudamos de lista, mude p/ a "r-br", esta em breve será desativada! > > att, > FH > > 2011/4/23 greice.laureano <greice.laureano(a)gmail.com> > > >> >> Bom dia, >> >> Gostaria de ajuda para extrair os quadrados médios de uma Anova. >> O banco de dados está abaixo assim como o código para fazer a Anova. >> Muito Obrigado; >> >> Greice. >> >> Dados: >> >> sujeito aferidor tempo valor >> 1 1 1 -2 >> 1 1 2 0 >> 1 1 3 1 >> 1 2 1 2 >> 1 2 2 1 >> 1 2 3 1 >> 2 1 1 16 >> 2 1 2 16 >> 2 1 3 15 >> 2 2 1 12 >> 2 2 2 14 >> 2 2 3 13 >> 3 1 1 5 >> 3 1 2 6 >> 3 1 3 6 >> 3 2 1 4 >> 3 2 2 4 >> 3 2 3 4 >> 4 1 1 11 >> 4 1 2 10 >> 4 1 3 10 >> 4 2 1 9 >> 4 2 2 7 >> 4 2 3 8 >> 5 1 1 7 >> 5 1 2 8 >> 5 1 3 6 >> 5 2 1 5 >> 5 2 2 6 >> 5 2 3 6 >> 6 1 1 -7 >> 6 1 2 -8 >> 6 1 3 -8 >> 6 2 1 -9 >> 6 2 2 -10 >> 6 2 3 -9 >> 7 1 1 18 >> 7 1 2 19 >> 7 1 3 19 >> 7 2 1 17 >> 7 2 2 17 >> 7 2 3 17 >> 8 1 1 4 >> 8 1 2 5 >> 8 1 3 5 >> 8 2 1 5 >> 8 2 2 5 >> 8 2 3 5 >> 9 1 1 0 >> 9 1 2 -3 >> 9 1 3 -2 >> 9 2 1 -7 >> 9 2 2 -6 >> 9 2 3 -5 >> 10 1 1 0 >> 10 1 2 0 >> 10 1 3 -2 >> 10 2 1 1 >> 10 2 2 2 >> 10 2 3 1 >> 11 1 1 -3 >> 11 1 2 -2 >> 11 1 3 -2 >> 11 2 1 -4 >> 11 2 2 -3 >> 11 2 3 -3 >> 12 1 1 3 >> 12 1 2 -1 >> 12 1 3 1 >> 12 2 1 -1 >> 12 2 2 -2 >> 12 2 3 1 >> 13 1 1 7 >> 13 1 2 9 >> 13 1 3 9 >> 13 2 1 4 >> 13 2 2 4 >> 13 2 3 2 >> 14 1 1 -6 >> 14 1 2 -7 >> 14 1 3 -6 >> 14 2 1 -8 >> 14 2 2 -10 >> 14 2 3 -9 >> 15 1 1 1 >> 15 1 2 1 >> 15 1 3 0 >> 15 2 1 -2 >> 15 2 2 -2 >> 15 2 3 -3 >> 16 1 1 -13 >> 16 1 2 -14 >> 16 1 3 -14 >> 16 2 1 -12 >> 16 2 2 -12 >> 16 2 3 -12 >> 17 1 1 2 >> 17 1 2 1 >> 17 1 3 0 >> 17 2 1 -1 >> 17 2 2 0 >> 17 2 3 0 >> 18 1 1 4 >> 18 1 2 4 >> 18 1 3 3 >> 18 2 1 7 >> 18 2 2 6 >> 18 2 3 4 >> 19 1 1 -10 >> 19 1 2 -9 >> 19 1 3 -10 >> 19 2 1 -10 >> 19 2 2 -11 >> 19 2 3 -10 >> 20 1 1 8 >> 20 1 2 9 >> 20 1 3 8 >> 20 2 1 2 >> 20 2 2 8 >> 20 2 3 8 >> 21 1 1 7 >> 21 1 2 6 >> 21 1 3 7 >> 21 2 1 8 >> 21 2 2 7 >> 21 2 3 7 >> 22 1 1 -3 >> 22 1 2 -2 >> 22 1 3 -4 >> 22 2 1 -5 >> 22 2 2 -5 >> 22 2 3 -5 >> 23 1 1 -5 >> 23 1 2 -5 >> 23 1 3 -7 >> 23 2 1 -6 >> 23 2 2 -8 >> 23 2 3 -7 >> 24 1 1 5 >> 24 1 2 5 >> 24 1 3 5 >> 24 2 1 3 >> 24 2 2 4 >> 24 2 3 4 >> 25 1 1 0 >> 25 1 2 -1 >> 25 1 3 -1 >> 25 2 1 -4 >> 25 2 2 -3 >> 25 2 3 -4 >> 26 1 1 7 >> 26 1 2 6 >> 26 1 3 6 >> 26 2 1 4 >> 26 2 2 4 >> 26 2 3 4 >> 27 1 1 -8 >> 27 1 2 -8 >> 27 1 3 -8 >> 27 2 1 -10 >> 27 2 2 -11 >> 27 2 3 -10 >> 28 1 1 1 >> 28 1 2 1 >> 28 1 3 2 >> 28 2 1 1 >> 28 2 2 -1 >> 28 2 3 0 >> 29 1 1 -3 >> 29 1 2 -3 >> 29 1 3 -3 >> 29 2 1 -5 >> 29 2 2 -4 >> 29 2 3 -5 >> >> sj=as.factor(dados$sujeito) >> af=as.factor(dados$aferidor) >> tempo=as.factor(dados$tempo) >> >> fit <- aov(dados$valor ~ sj*af + sj:tempo:af + Error(sj/af)) >> summary(fit) >> >> > __._,_.___ > > <fernandohtoledo@gmail.com?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| através > de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| Responder > através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJyNWprYWVzBF9TAzk3N…>| Adicionar > um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmNDRncmJ1BF9TAzk3N…> > Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/16039;_ylc=X3oDMTM3N2Q4Y2Vi…>( > 2) > Atividade nos últimos dias: > > - Novos usuários<http://br.groups.yahoo.com/group/R_STAT/members;_ylc=X3oDMTJnNmk2bWJoBF9TAz…> > 2 > > Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmbDdmMWM3BF9TAzk3NDkwND…> > Não deixe de ir atrás da pessoa dos seus sonhos!<http://global.ard.yahoo.com/SIG=15me74fii/M=758712.14532720.14422460.129601…> > [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlc2RmczJyBF9TAzk3NDkwNDM1BGdycElkAz…> > Trocar para: SÃ³ Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>, > Resenha DiÃ¡ria<R_STAT-digest@yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair > do grupo <R_STAT-unsubscribe(a)yahoogrupos.com.br?subject=Sair+do+grupo> • Termos > de uso <http://br.yahoo.com/info/utos.html> > . > > __,_._,___ >

3 5

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 20 Abr '11

20 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

crash mgcv
by Marilia Sá Carvalho 20 Abr '11

20 Abr '11

C@ros, Estou rodando R version 2.13.0 (2011-04-13), Platform: i486-pc-linux-gnu (32-bit), UBUNTU LUCID. Ao ajustar um modelo gam, com um termo bivariado - s(long,lat) - o R simplesmente sai, sem pedir licença. puff! Alguem ja tentou rodar a mgcv com essa versao? O problema e' do R ou do Ubuntu? O spline univariado roda sem problema. Todas as libs estao aualizadas, inclusive as que o mgcv usa. (Fiz um update.packages geral). Alguma sugestao? Abracos, Marilia -- -- Marilia Sá Carvalho Programa de Computação Científica - FIOCRUZ Rio de Janeiro - Brasil Tel: +55 21 38361116/1100 Fax: +55 21 22705141

2 1

Função nestednodf do pacote vegan
by Jéssica Caroline Falcão 20 Abr '11

20 Abr '11

Olá pessoal, alguém já usou a função "nestednodf" do pacote vegan? eu consigo rodá-la mas não consigo obter os valores de "p" Obrigada desde já! -- ------------------------------------------------------------------ Jéssica Falcão Mestranda em Ecologia e Conservação da Biodiversidade Universidade Federal de Mato Grosso IB- Departamento de Ecologia e Botânica Lab. Interação Inseto-Planta / Mirmecologia Cuiabá, Mato Grosso, BRASIL Email: jecafalcao(a)gmail.com

2 1

Resumo de comparações múltiplas.
by Olympio Neto 20 Abr '11

20 Abr '11

Boa noite, pessoal. Gostaria de saber se alguém conhece algum pacote do R que faz aqueles resumos com letras em comparações múltipla (estudo usando Tukey (multcomp) com modelo estimado via lme). No pacote agricolae, existe e sempre o utilizo quando estou sob um modelo linear ordinário (aov), mas com lme estou tendo um trabalho danado! No pacote multcomp existe o cld, mas o resultado é super pobre. Alguém conhece algo neste sentido? Abs e muito obrigado. Olympio

2 2

Re: [R-br] regressão multivariada
by Thiago De paula protásio 19 Abr '11

19 Abr '11

Concordo com o Allaman... Na análise da correlação linear entre grupos (ou a combinação linear) é mais apropriado utilizar a análise de correlação canônica e nesse caso o pacote CCA é bem interessante. A função cc() desse pacote permite extrair as cargas canônicas e as cargas canônicas cruzadas além da correlação canônica propriamente dita e de alguns gráficos. A disposição!!! 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

1 0

Re: [R-br] Análise de Correlação Canônica!!! Índices de Redundância!!!
by Thiago De paula protásio 19 Abr '11

19 Abr '11

Allaman; Rodei a função aqui e funcionou perfeitamente..... É isso mesmo que eu estava precisando..... Se precisar de alguma coisa estou a disposição! Novamente obrigado pela ajuda e esforço para responder ao tópico! Abraço! 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

1 0

construção de grafico complexo...
by Ayuni Sena 19 Abr '11

19 Abr '11

Pessoal, lá vai mais uma mega duvida.... Eu preciso gerar um gráfico muito parecido com a figura disponivel no link: http://www.datafilehost.com/download-f1ba4ceb.html Existe apenas uma pequena diferença: No meu caso, o eixo Y é repesentado pela variável N o eixo X pela variavel D Foram ajustados os seguintes modelos: 1 - N~SSasymp(D) #Regressao assintotica do livro do Pinheiro e Bates (2000) - Mixed effect models 2 - D ~ b0*N^b1*hdc^b2 3 - V ~ b3*D^b4*hdc^b5*N^b6 A partir destes modelos, quero plotar no gráfico as isolinhas que correspondem a: 1 - 100%, 55%, 40% e 15% da curva N~SSasymp(Dq) - isso vai gerar apenas um deslocamento da curva para baixo, mantendo o mesmo comportamento. 2 - a partir do modelo 2, obtemos a seguinte equação: N = (D/(b0*hdc^b2))^(1/b1) e assim, quero plotar no grafico as isolinhas correspondente a hdc=5, hdc=10, hdc=15, hdc=20, e assim por diante... (existirá o mesmo paralelismo entre as isolinhas relatado no passo acima) 3 = a partir do modelo 3, obtemos a seguinte equação: N = (V/(b3*D^b4*hdc^b5))^(1/b6) e assim, quero plotar as isolinhas correspondentes a V=50, V=100, V=150, e assim por diante.... (existirá o mesmo paralelismo entre as isolinhas relatado no passo acima) os passos nº 1 e nº 2 eu consegui fazer com auxilio dos comandos "seq" e "lines", da seguinte forma: #passo 1: N100<-779.36+(26947.92)*(exp(-exp(-1.29)*D)) N55<-428.648+(14821.356)*(exp(-exp(-1.29)*D)) N40<-311.744+(10779.168)*(exp(-exp(-1.29)*D)) N15<-116.904+(4042.188)*(exp(-exp(-1.29)*D)) plot(D, N, xlab=" Dq (cm)", ylab="Número de plantas por hectare (N/ha)", ylim=c(0,5000)) lines(dq,n100,type="l", lty=1, lwd=1) lines(dq,n55,type="l", lty=2, lwd=1) lines(dq,n40,type="l", lty=1, lwd=1) lines(dq,n15,type="l", lty=2, lwd=1) #passo 2: b0=25.360699 b1=-0.303046 b2=0.51502 b3=0.0000377 b4=2.2286023 b5=0.7798658 b6=0.9907735 D<-seq(from=5.5, to=21, lenght.out=1061) hdc<-seq(from=7, to=31,5, lenght.out=1061) N1<-(D/(b0*15^b2))^(1/b1) N2<-(D/(b0*20^b2))^(1/b1) #etc etc etc plot(DQ, NHA, ylab="Densidade de plantas por hectare (N/ha)", xlab="Dq (cm)") lines(D, N1,type="l", lty=1, lwd=1) lines(D, N2,type="l", lty=1, lwd=1) #etc etc etc Com auxílio desta rotina eu nao consegui fazer o passo 3.... Alguem teria uma sugestão melhor???????????????????????????? Agradeço a atenção de todos! Ayuni

1 0

Re: [R-br] Análise de Correlação Canônica!!! Índices de Redundância!!!
by Thiago De paula protásio 19 Abr '11

19 Abr '11

Allaman; Obrigado pela ajuda!! Sua namorada irá entender o esforço... Quanto ao R² canônico se eu não estiver errado.....é o seguinte: A quantidade de variância explicada, ou seja, o percentual de variância na variável estatística canônica dependente que pode ser explicada pela variável estatística canônica independente, e vice-versa, pode ser determinada elevando-se ao quadrado a correlação canônica (R² canônico). Obrigado novamente por mandar o pacote e a função....... A disposição! 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

1 0

livro geoestatística
by Mauro Sznelwar 19 Abr '11

19 Abr '11

Alguém conhece o livro Aplied Statistical Data Analysis with R? Sabe se tem algo mais ou diferente, que o Model basead Geostatistics do prof. Paulo Justiniano? Mauro

4 4

Re: [R-br] Análise de Correlação Canônica!!! Índices de Redundância!!!
by Ivan Bezerra Allaman 19 Abr '11

19 Abr '11

Quase apanhei da mulher tentando te responder este tópico! Bom, é o seguinte: Depois de umas duas horas na internet pesquisando sobre tal assunto, descobri o seguinte: Segundo o livro do Hair, o índice de redundância é calculado como a carga quadrada média vezes o R² canônico. No livro do gente fina do Hair e nem na net, eu não encontrei como se calcula o R² canônico. Se alguém souber me diga por favor. No entanto, encontrei um pacote que tem uma função o qual se extrai este índice de redundância. O pacote é o "yacca". Utilizando a base de dados do livro do Hair para conferir é claro se estava certo. Logo temos: setwd("/home/ivan/Área de Trabalho") daa <- read.table("da.txt", h=T) summary(daa) install.packages("yacca") library(yacca) X <- daa[,2:8] Y <- daa[,9:10] mod <- cca(X,Y) mod$xvrd#lhe retorna os índices de redundância mod$yvrd#lhe retorna os índices de redundância Olhe o pdf do pacote para maiores detalhes. Abraço. 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

plotar dados longitudinais
by Ayuni Sena 19 Abr '11

19 Abr '11

Boa Noite Pessoal Tenho um conjunto de dados longitudinais desbalanceados, onde cada amostra foi medida duas ou tres vezes (dependendo da amostra) com intervalos de 1 a 2 anos entre as medições. Eu gostaria de plotar em um gráfico X1~X2, de modo que os pontos que correspondem a mesma unidade amostral, porem em tempos diferentes, possam ser ligados, unidos um ao outro, mostrando o que aconteceu com a unidade amostral de um ano pro outro (aumentou X1 em função do aumento de X2, por exemplo). Lembrando que nem X1 nem X2 correspondem a variavel tempo. Quero fazer algo semelhante a figura que esta disponivel no seguinte link: http://www.datafilehost.com/download-e7370438.html Alguem sabe como fazer isso no R? Abraços Ayuni

1 0

Linux ubuntu
by Luciana Melo 19 Abr '11

19 Abr '11

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

1 0

Análise de Correlação Canônica!!! Índices de Redundância!!!
by Thiago De paula protásio 19 Abr '11

19 Abr '11

Prezados membros da R-Br; Novamente peço a ajuda de vcs...... Estou precisando de uma função para calcular o índice de redundância da análise de correlação canônica! Agradeço a qualquer ajuda!!! Bom feriado a todos!!!! 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

1 0

Otimização (paralelização) de aplicações no R
by Junior Beleti 19 Abr '11

19 Abr '11

Realizo minha pesquisa na área da programação paralela e pretendo estudar e desenvolver modelos de programação paralela para apoiar as aplicações no R. Realizei algumas buscas e encontrei o resumo do artigo: Neil Shah, Guruprasad Kora, Paul Breimyer, Yekaterina Shpanskaya. pR: Enabling Automatic Parallelization of Data-Parallel Tasks and Interfacing Parallel Computing Libraries in R with Application to Fusion Reaction Simulations. Link: http://www.r-project.org/conferences/useR-2010/abstracts/Shah+Kora+Breimyer… Gostaria de saber se alguém tem acesso ou sabe como eu posso acessar esse artigo ou outros que possam ser úteis em minha pesquisa. Obrigado, -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

3 2

Dúvida expression( ) - Resolvido
by Alexandre dos Santos 19 Abr '11

19 Abr '11

Muito obrigado Benilton, Abraço, Alexandre -----Original Message----- From: r-br-bounces(a)listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho Sent: terça-feira, 19 de abril de 2011 12:53 To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Dúvida expression( ) text(150,3,expression(italic("y(h)")==4.91+6.895~"["~1-italic("exp")(italic( "h")/134.90^2)~"]"),cex=0.75) 2011/4/19 Alexandre dos Santos <alexandresantosbr(a)yahoo.com.br>: > Boa tarde Pessoal, > > Esta tentando colocar em um gráfico a equação representada por: > > > > text(150,3,expression(italic("y(h)")==4.91+6.895[1-italic("exp")(italic("h") /134.90^2)]),cex=0.75) > > > > Porém, na realidade estou tentando exibir o símbolo de cochetes[] que faz > parte da equação, no entando, a função expression() interpreta como > > sobreescrito, alguem já passou por essa situação e 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 > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

1 0

Dúvida expression( )
by Alexandre dos Santos 19 Abr '11

19 Abr '11

Boa tarde Pessoal, Esta tentando colocar em um gráfico a equação representada por: text(150,3,expression(italic("y(h)")==4.91+6.895[1-italic("exp")(italic("h") /134.90^2)]),cex=0.75) Porém, na realidade estou tentando exibir o símbolo de cochetes[] que faz parte da equação, no entando, a função expression() interpreta como sobreescrito, alguem já passou por essa situação e 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

2 2

Dúvida expression( )
by Alexandre dos Santos 19 Abr '11

19 Abr '11

Boa tarde Pessoal, Esta tentando colocar em um gráfico a equação representada por: text(150,3,expression(italic("y(h)")==4.91+6.895[1-italic("exp")(italic("h") /134.90^2)]),cex=0.75) Porém, na realidade estou tentando exibir o símbolo de cochetes[] que faz parte da equação, no entando, a função expression() interpreta como sobreescrito, alguem já passou por essa situação e 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

1 0

Enc: BINOMIAL NEGATIVA TRUNCADA
by Eduardo Uliana 19 Abr '11

19 Abr '11

----- Mensagem encaminhada ---- De: Eduardo Uliana <eduardouliana(a)yahoo.com> Para: r-br-request(a)listas.c3sl.ufpr.br Enviadas: Quinta-feira, 14 de Abril de 2011 11:59:13 Assunto: BINOMIAL NEGATIVA TRUNCADA Bom Dia! Estou trabalhando com a distribuição binomial negativa truncada e gostaria de saber se tem como estimar os parâmetros e as probabilidades desta distribuição no R. Aguardo a resposta dos colegas. Obrigado. Eduardo Morgan Uliana Engenheiro Ambiental - Mestrando em Produção Vegetal(CCA-UFES)

3 4

Pacote "Revolution"
by Marcelo Rodrigo 19 Abr '11

19 Abr '11

Bom dia, alguém já usou o pacote "Revolution" (http://www.revolutionanalytics.com/) disponível nos repositórios do Ubuntu? Há realmente ganho? Abraços, -- Marcelo R. P. Ferreira Professor Assistente I Deparamento de Estatística/UFPB Tel.: +55 83 3216-7075/7785/7565 (ramal 210) http://www.de.ufpb.br/~marcelo

2 2

regressão multivariada
by Ayuni Sena 18 Abr '11

18 Abr '11

Boa Noite Colegas! Alguem saberia indicar uma função/pacote para analise de regressao multivariada? Agradeço a atenção Ayuni Sena

2 3

[Dúvida] Remover acentos de String
by Pedro Rafael 18 Abr '11

18 Abr '11

Pessoal como remover os acentos de um vetor de caracteres bem como caracteres especiais como Ç? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 1

Interação significativa em regressão logística
by Marcelo Cardoso mello 18 Abr '11

18 Abr '11

Prezados, tenho um modelo logístico do tipo: Prenhez (+/_) -> (Método de inseminação (0,1); Época do ano (0,1); interação) A interação foi significativa. Então, avaliei os métodos em cada época, ou seja, fiz dois modelos mais simples, um para cada época: Prenhez (+/_) -> (Método de inseminação (0,1)) (considerando a época 1 e teve efeito do método) Prenhez (+/_) -> (Método de inseminação (0,1)) (considerando a época 0 e nãp teve efeito do método) O procedimento é esse mesmo? Fiz via glm. Tem uma maneira diferente de fazer essa avaliação? Tenho mesmo que partir para modelos separados dentro de cada época? Agradeço a opinião dos colegas. Marcelo

3 7

Arvore Geradora Minima
by Emerson Cotta Bodevan 18 Abr '11

18 Abr '11

Prezados, existe algum pacote no R para construção de árvores geradoras mínimas? Abs, Emerson.

2 2

Fwd: [Dúvida] Limpando memória RAM
by Pedro Rafael 18 Abr '11

18 Abr '11

Realmente Jakson apesar de tá usando o WIndows que não tem um gerenciamento de memória que preste eu fiz o teste e não vi grandes melhoras. Consegui otimizar reescrevendo o programa de uma forma mais eficiente. Pedro ---------- Mensagem encaminhada ---------- De: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> Data: 18 de abril de 2011 11:01 Assunto: Re: [R-br] [Dúvida] Limpando memória RAM Para: "Jakson A. Aquino [via R-br]" < ml-node+3454967-637477867-223914(a)n4.nabble.com> Realmente Jakson apesar de tá usando o WIndows que não tem um gerenciamento de memória que preste eu fiz o teste e não vi grandes melhoras. Consegui otimizar reescrevendo o programa de uma forma mais eficiente. Pedro Em 17 de abril de 2011 01:13, Jakson A. Aquino [via R-br] < ml-node+3454967-637477867-223914(a)n4.nabble.com> escreveu: On Sat, Apr 16, 2011 at 03:03:49PM -0300, Henrique Dallazuanna wrote: > > remova o objeto com o rm e após execute gc(reset = TRUE) > > Eu não conhecia a opção reset da função gc(), mas num teste que fiz > aqui ela não parece ter feito muita diferença. Carreguei o banco de > dados do World Values Survey com read.dta (o arquivo está disponível > em três formatos, sendo um deles o do STATA). O comando top do Linux > mostrou que o R estava ocupando em torno de 2,7 GB da RAM. > > Depois dos comandos > > rm(list=ls()) > gc(reset = TRUE) > > a memória ocupada continuou um pouco acima de 1 GB. > > Existe alguma forma mais eficiente do que gc() de mandar o R > reorganizar (desfragmentar?) a memória? > > Atualmente, tenho feito o seguinte em situações em que preciso dessa > reorganização da memória: > > 1) Salvo a área de trabalho: save.image() > 2) Fecho o R: quit(save = no) > 3) Abro o R novamente e carrego a área de trabalho salva. > > [É claro que poderia ter usado quit(save = "yes"), mas prefiro > escolher um nome para o RData.] > > -- > Jakson > _______________________________________________ > R-br mailing list > [hidden email]<http://user/SendEmail.jtp?type=node&node=3454967&i=0&by-user=t> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Duvida-Limpando-memoria-RAM-tp345434… > To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

Movimentar título nos eixos x y de um gráfico
by Marcelo Claro de Souza 18 Abr '11

18 Abr '11

Além de utilizar xlab="nome da variável" e ylab="nome da variável" para determinar o título dos eixos, como faço para movimentar os títulos para mais perto ou mais longe da figura? Muito obrigado.

3 3

Re: [R-br] Movimentar título nos eixos x y de um gráfico
by Marcelo Claro de Souza 18 Abr '11

18 Abr '11

José, estou utilizando o script abaixo: cor=as.data.frame(cor) simbolos=c(19,19,19,19,19,19,19,2,2,2,2,2,2,2,3,3,3,3,3,3,3) N.P.cor=lm(N~P,data=cor) require(graphics) X11() jpeg('NPcor.jpeg') plot(N~P,data=cor,pch=simbolos,ylab=expression(N~(g~Kg^{-1})),xlab=expression(P~(g~Kg^{-1}))xlim=c(0,1.6),ylim=c(0,35)) abline(N.P.cor) legend("topleft",c("Ep","Es","Et"),pch=c(19,2,3)) dev.off() Eu gostaria de mudar o título do eixo y p mais próximo da figura (p o lado direito) pq os valores sobrescritos estão ficando muito colados na borda da figura enão estão aparecendo direito. Segue uma figura em anexo marcando o problema. Muito obrigado. --- Em seg, 18/4/11, José Cláudio Faria <joseclaudio.faria(a)gmail.com> escreveu: De: José Cláudio Faria <joseclaudio.faria(a)gmail.com> Assunto: Re: [R-br] Movimentar título nos eixos x y de um gráfico Para: r-br(a)listas.c3sl.ufpr.br Data: Segunda-feira, 18 de Abril de 2011, 8:44 Marcelo, Isso irá depender muito do ambiete gráfico que vc está usando (graphics, grig, rgl) e do pacote/função (graphics, lattice, ggplot2). Se você mandar um script de exemplo a ajuda poderá ser facilitada e mais certeira. Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

2 1

[Dúvida] Limpando memória RAM
by Pedro Rafael 18 Abr '11

18 Abr '11

Pessoal eu tenho duas base de dados de 600 MB totalizando 1200 MB e carro ela no R mas não necessito da base toda carregada na memória RAM. As variáveis que me interessam eu jogo em um vetor e com o comando rm() eu removo os dois data.frames. Esse comando retira os objetos completamente da memória RAM ou deixa lixo? Existe uma forma mais eficiente? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

7 10

Link para grupo no Site do R
by Augusto Ribas 18 Abr '11

18 Abr '11

Ola pessoal. Estava navegando pelo site do R (http://www.r-project.org/) e no menu lateral tem um link uma documantação com os user groups do mundo. Ai ja tinha o Brasil la mas cadastrada a lista R-Stat. Eu tomei a liberdade de atualizar para a lista r-br e coloquei o Leg como responsavel la. Espero que esteja tudo bem, se alguem quiser dar uma conferida la: http://rwiki.sciviews.org/doku.php?id=rugs:r_user_groups Assim, quem vier do site oficial do R ja vai cair nessa lista e nao no R-stat. Julguei nao necessario deixar o link para a lista R-stat la ja que ela vai ser desativa. Bem é isso Abraços -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org/ http://dhost.info/augusto/ Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Re: [R-br] Superficie de resposta publicavel no R?!
by José Cláudio Faria 17 Abr '11

17 Abr '11

Pessoal, Parece que o anexo passou nos servidores! Na parte da mensagem aqui para mim aparece assim: -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: srp.r Tipo: application/x-wine-extension-r Tamanho: 2317 bytes Descrição: não disponível URL:<http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110417/cfe9351f/atta…> Basta fazer o download do link URL acima e mudar o nome de attachement.bin para srp.r. Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Superficie de resposta publicavel no R?!
by José Cláudio Faria 17 Abr '11

17 Abr '11

Pessoal, Sempre que penso em fazer uma superfície de resposta publicável no R me dá um frio na barriga... A alguns anos fiz (com a ajuda do J.Fox e D.Murdoch) uma função (usando o rgl) muito flexível (disponível em: http://rwiki.sciviews.org/doku.php?id=graph_gallery:new-graphics&s=plotlm3d <http://rwiki.sciviews.org/doku.php?id=graph_gallery:new-graphics&s=plotlm3d> O problema dela é que labels do eixos não admitem expressões matemáticas, geralmente necessárias para expressar a escala da variável (mm^2, kg h^-1, etc). A limitação não é da função, mas do rgl. Dá para contornar editando o figura gerada ou inserindo bitmaps usando os recursos do rgl, mas ambas alternativas trabalhooosas. Pois bem, estou diante da necessidade de refazer várias superfícies de resposta, feitas no passado (9 anos) com o Statiscs for Windows (nunca mais usei, mas ainda sinto saudades de alguns recursos...). Com um frio na barriga e ralando no fim de semana, decidi tentar simplificar esta tarefa para o presente e futuro. Com a ajuda de alguma coisa que achei nas listas de discussão, a juda do Walmes e do Henrique (vcs não descansam no domingo né? ;-)) rascunhei umas mal traçadas linhas no R, que embora não perfeita, está resolvendo o básico. Decidi socializar, pois é assunto recorrente nas listas de discussão, e não fiz só. Possivelmente seja útil para outros usuários, melhor ainda se receber alguém fizer melhoramentos. No momento estou sem tempo, mas se tiver alguém interessado, podemos pensar em juntar a "plolm3d" e a atual "sr" em um pacote e publicar na CRAN. Penso que seria de muita utilidade! A função "sr", assim como um pequeno script demo estão em anexo (vou testar para ver se a lista recebe, se não for possível envio outra mensagem com o script). Abs, (obrigado novamente so Walmes e Henrique pela disposição de ajudar) -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0