Preciso comparar meu valor vindo de um bootstrap de matriz, com o valor único original, mas não sei exatamente como fazer isso com meus dados. Aqui eu tenho os dados de exemplo
###DADOSEXEMPLO##
valores <- c(1, 5, 3, 8, 2, 9, 3, 2, 3 ,2, 5, 3, 6, 5, 1, 5, 3, 4, 2, 2, 5, 6, 5, 7, 1, 2,
8, 12, 5, 1, 5, 3, 6, 5, 9, 3,3,4,5,7,8,2,7,8,4,3,1,4,4 )
matrixnet<-matrix(valores , 7, 7)
matcor <- cor(matrixnet, method = "spearman")
matrixnetwork = graph.adjacency(matcor, mode="undirected",
weighted
= TRUE,
add.colnames=NULL, diag=FALSE)
plot(matrixnetwork)
Aqui calculo o valor único do gráfico, neste caso para a medida entre
centr_betw(matrixnetwork, directed = FALSE, normalized = TRUE)$res
Aqui aplico a permutação "alterar valor em vertes" em conjunto com a função específica que desejo aplicar, neste caso Betw
set.seed(123
bootcentr_Betw <- lapply(seq_len(Nperm), function(x){
randomnet <- rewire(matrixnetwork, with=each_edge(0.5)) #rewire vertices with constant probability
E(randomnet)$weight <- sample(E(matrixnetwork)$weight) #shuffle initial weights and assign
them randomly to edges
return(centr_betw(randomnet)$res)
})
Agora preciso comparar o valor original com o bootstrap do gráfico aleatório.
Abaixo tentei traçar o valor único original em relação ao bootstrap, replicando o código que obtive de outro lugar
# 4. Plote o valor observado em rela??o ? distribui??o de valores aleat?rios
niveis<-row.names(bootstrength )
for(k in niveis)
{
if(any(is.na(bootstrength [k,]) == TRUE))
{
print(c(k, "metrica tem NA"))
} else {
nome.arq<- paste("modularity",k,".png", sep="")
png(filename= nome.arq, res= 300, height= 15, width=21, unit="cm")
plot(density(bootstrength [k,]), main="Observed vs.
randomized",)
abline(v=obs[k], col="red", lwd=2, xlab="")
dev.off()
print(k)
nome.arq<- paste("Patefield_Null_mean_sd_",k,".txt", sep="")
write.table(cbind(mean(bootstrength [k,]),sd(bootstrength [k,])),
file=paste(nome.arq,sep=""),
sep=" ",row.names=TRUE,col.names=FALSE)
}
}
Confesso que não sei onde está o arquivo que ele gerou.
Depois disso tentei replicar as estatísticas que encontrei aqui no fórum, mas não funcionou.
#5. Estime o valor P (signific?ncia)
Estou tentando comparar um valor observado e sua permutação e gerar um valor p pra isso.
Eu fiz o teste t sugerido no código em uma resposta abaixo. Porém vi que algumas pessoas utilizam uma fórmula chama de teste de permutação. Que é Número de valores (estatísticas) que são maiores ou iguais ao valor observado e dividem pelo número de valores. No código, pval = soma(s) >= s0)/N; Alguém pode me ajudar a aplicar essa fórmula aos dados acima? Pois é diferente do teste que eu apliquei