É possível desenhar, por exemplo Elipses de confiança de 95% em torno de amostras do mesmo grupo nos resultados da função plotMDS em edgeR? Se sim, como?
É possível desenhar, por exemplo Elipses de confiança de 95% em torno de amostras do mesmo grupo nos resultados da função plotMDS em edgeR? Se sim, como?
Posso ver que haveria um processo de três etapas para fazer isso:
Mesclar contagens de todas as amostras no grupo e, em seguida, reamostrar as pseudo-réplicas a partir delas. Se x
for uma matriz com as amostras no grupo sendo colunas e os genes sendo linhas
s <- amostra (row.names (x), n = média (colSums (x)), probs = rowSums (x) / sum (rowSums (x)) stab <- tabela (s) s <- as.vector (stab) nomes (s) <- nomes (stab)
Faça isso milhares de vezes.
Você, então, desejaria projetá-los em seu espaço de plotagem MDS - estou menos claro sobre como você poderia fazer isso sem perturbar o próprio espaço.
Calcule o parâmetro para uma elipise que conteria 95% desses pontos. Novamente, não tenho certeza sobre como fazer isso, mas Acho que o pacote car
pode ser um bom lugar para começar a procurar.