Я пытаюсь построить анализ основных компонентов, используяprcomp а такжеggbiplot, Я получаю значения данных за пределами круга устройства и не смог изменить их масштаб до вызоваprcomp таким образом, что я могу ограничить данные единичным кругом.

data(wine)
require(ggbiplot)
wine.pca=prcomp(wine[,1:3],scale.=TRUE)
ggbiplot(wine.pca,obs.scale = 1, 
         var.scale=1,groups=wine.class,ellipse=TRUE,circle=TRUE)

Я пытался масштабировать, вычитая среднее и деление на стандартное отклонение, прежде чем звонитьprcomp:

wine2=wine[,1:3]
mean=apply(wine2,2,mean)
sd=apply(wine2,2,mean)
for(i in 1:ncol(wine2)){
  wine2[,i]=(wine2[,i]-mean[i])/sd[i]
}
wine2.pca=prcomp(wine2,scale.=TRUE)
ggbiplot(wine2.pca,obs.scale=1, 
         var.scale=1,groups=wine.class,ellipse=TRUE,circle=TRUE)

ggbiplot Пакет установлен следующим образом:

require(devtools)
install_github('ggbiplot','vqv')

Вывод любого фрагмента кода:

Согласно приведенному ниже комментарию @Brian Hanson, я добавляю дополнительное изображение, отражающее вывод, который я пытаюсь получить.