Programa R para Biologia da Conservação
Gostaria de reagir a esta mensagem? Crie uma conta em poucos cliques ou inicie sessão para continuar.

Problemas com Gráfico Média e CI - plotCI

2 participantes

Ir para baixo

Problemas com  Gráfico Média e CI - plotCI Empty Problemas com Gráfico Média e CI - plotCI

Mensagem por souzaalves Seg Jun 04, 2018 10:48 pm

Pessoal,

estou tentando criar um gráfico de médias com o CI, e após todo script pronto, ele aparece um erro que não consigo, de forma alguma, resolver.
Como são poucos os dados, foi rodado um bootstrap antes e o gráfico gerado sobre ele. Segue todo o passo-a-passo.

Obrigadão!!!!

library(boot)
library(plotrix)

own_mean<-function(x,y) {
mean1 = mean(x[y])
}

a<-rnorm(105,1.6)
b<-rnorm(105,2.7)
c<-rnorm(105,2.3)
d<-rnorm(105,1.Cool
e<-rnorm(105,4.0)
f<-rnorm(105,2.Cool
g<-rnorm(163,0.4)
h<-rnorm(163,1.1)
i<-rnorm(163,2.2)
j<-rnorm(163,2.7)
k<-rnorm(163,1.2)
l<-rnorm(163,1.0)

### Mean bootstrapping and 99%CI generation ###
boot_a<-boot(a, own_mean, R=9999)
boot_b<-boot(b, own_mean, R=9999)
boot_c<-boot(c, own_mean, R=9999)
boot_d<-boot(d, own_mean, R=9999)
boot_e<-boot(e, own_mean, R=9999)
boot_f<-boot(f, own_mean, R=9999)
boot_g<-boot(g, own_mean, R=9999)
boot_h<-boot(h, own_mean, R=9999)
boot_i<-boot(i, own_mean, R=9999)
boot_j<-boot(j, own_mean, R=9999)
boot_k<-boot(k, own_mean, R=9999)
boot_l<-boot(l, own_mean, R=9999)


ci_a<- boot.ci(boot_a,conf = 0.99)
ci_b<- boot.ci(boot_b,conf = 0.99)
ci_c<- boot.ci(boot_c,conf = 0.99)
ci_d<- boot.ci(boot_d,conf = 0.99)
ci_e<- boot.ci(boot_e,conf = 0.99)
ci_f<- boot.ci(boot_f,conf = 0.99)
ci_g<- boot.ci(boot_g,conf = 0.99)
ci_h<- boot.ci(boot_h,conf = 0.99)
ci_i<- boot.ci(boot_i,conf = 0.99)
ci_j<- boot.ci(boot_j,conf = 0.99)
ci_k<- boot.ci(boot_k,conf = 0.99)
ci_l<- boot.ci(boot_l,conf = 0.99)

#Esse é o script para o plotCI

plotCI(c(6,3), c((mean(a),mean(g)), (mean(b),mean(h)), (mean(c),mean(i)), (mean(d),mean(j)), (mean(e),mean(k)), (mean(f),mean(l))), ui= c((ci_a[[4]][3],ci_g[[4]][3]), (ci_b[[4]][3], ci_h[[4]][3]), (ci_c[[4]][3], ci_i[[4]][3]), (ci_d[[4]][3], ci_j[[4]][3]), (ci_e[[4]][3], ci_k[[4]][3]), (ci_f[[4]][3], ci_l[[4]][3])), li=c((ci_a[[4]][2],ci_g[[4]][2]), (ci_b[[4]][2], ci_h[[4]][2]), (ci_c[[4]][2], ci_i[[4]][2]), (ci_d[[4]][2], ci_j[[4]][2]), (ci_e[[4]][2], ci_k[[4]][2]), (ci_f[[4]][2], ci_l[[4]][2])), err = "y", xaxt = "n" ,xlab="Quadrats",ylab="PIA (mean±99%CI)", gap= T,pt.bg=par("bg"),pch=16,xlim=c(1,12),ylim=c(0,12),main="T. pusillus")

Error: unexpected ',' in "plotCI(c(6,3), c((mean(a),"
souzaalves
souzaalves

Mensagens : 7
Data de inscrição : 04/10/2017

Ir para o topo Ir para baixo

Problemas com  Gráfico Média e CI - plotCI Empty Re: Problemas com Gráfico Média e CI - plotCI

Mensagem por Valeria Andrade Ter Jun 05, 2018 10:31 am

Ola SousaAlves,
Nao consigo ajuda-lo!!Tambem to "suando" para fazer um grafico de medias e IC..assim ..mas com Jacknife.....
Tomare que o Professor Marcus, explique sua "nossa" duvida no webnario #7 ou por aqui mesmo..
Se eu conseguir algum sucesso nesse grafico, envio imediatamente pra o forum o script!
Boa sorte!
Valeria

Valeria Andrade

Mensagens : 94
Data de inscrição : 12/03/2018

Ir para o topo Ir para baixo

Ir para o topo

- Tópicos semelhantes

 
Permissões neste sub-fórum
Não podes responder a tópicos