## VMU on portrait value questionnaire ratings
fit_vmu <- vmu(PVQ40agg) ## fit 2D VMU
fit_vmu
plot(fit_vmu, cex = c(1, 0.7)) ## call biplot from stats
## VMU with fixed column coordinates (circular)
tuv <- matrix(0, nrow = 10, ncol = 2)
alpha <- -360/10
for (i in 1:10){
alpha <- alpha+360/10
tuv[i, 1]<- cos(alpha*pi/180)
tuv[i, 2] <- sin(alpha*pi/180)
}
fit_vmu2 <- vmu(PVQ40agg, col.coord = tuv) ## fit 2D circular VMU
fit_vmu2
plot(fit_vmu2, cex = c(1, 0.7))
Run the code above in your browser using DataLab