Как создать попарно диаграмму рассеяния с данными разложения по сингулярным значениям - PullRequest
Новая
       3

Как создать попарно диаграмму рассеяния с данными разложения по сингулярным значениям

0 голосов
/ 04 марта 2020

За последние 14 часов я пытался создать попарную диаграмму рассеяния, учитывая данные, полученные из изображения. Данные следующие: 

library (imager)
library(graphics)
im <- load.image('Denny.jpg')
im
plot(im)
img=(im[,,,1]+im[,,,2]+im[,,,3])/3
img
img= apply(apply(img,1,rev))
image(img,col=grey(0:12/12))


    # See more examples at https://cran.r-project.org/web/packages/imager/vignettes/gettingstarted.html


#####################
library(graphics)
library(ggplot2)
library(GGally)

A=apply(im[,,,1],1,rev)
A
image(A)

### singular value decomposition 
svd.fit=svd(A)
svd.fit
plot(svd.fit$d)

## choose first 20 components
A.hat=svd.fit$u[,1:20]%*%diag(svd.fit$d[1:20])%*%t(svd.fit$v[,1:20])
image(A.hat)
A.hat

pairs(A.hat)
svd.fit=svd(A)

plot(svd.fit$d)

## choose first 20 components
A.hat=svd.fit$u[,1:20]%*%diag(svd.fit$d[1:20])%*%t(svd.fit$v[,1:20])
image(A.hat)
ggpairs(A.hat)

Любая помощь будет принята с благодарностью !!!

...