Я использую приведенный ниже код для создания графика PCA в Rstudio:
col <- c(rep("orange2", 9), rep("red", 16), rep("forestgreen", 11))
pca <- prcomp(t(counts))
plot(pca$x[,1:2], xlab="", ylab="", col=col, pch=20, cex=2, type="n") +
text(pca$x[,1:2], sample, col=col) +
lam <- pca$sdev
v <- lam^2
v <- v/sum(v)
k1 <- round(100*v[1], 2)
k2 <- round(100*v[2], 2)
title(main=paste("All samples,", k1+k2, "% of variance explained"), xlab=paste(k1, "% of
variance explained"), ylab=paste(k2, "% of variance explained")) +
legend(-80, -25, c("PEDIS:2", "PEDIS:3", "PEDIS:4"), text.col=c("orange2", "red",
"forestgreen"), cex=0.8)
Последовательный запуск этих линий работает нормально, чтобы создать график для экспорта в виде pdf или png, но я надеялся запустить этот код для генерации объекта данных, чтобы я мог объединить эту фигуру с другим. Как мне go сделать это? Я попытался добавить example <- plot... в начало строки 3, но я получаю эту ошибку:
Ошибка в заголовке (main = paste ("All samples,", k1 + k2, "% of объяснение дисперсии "),: non-Numberri c аргумент бинарного оператора