ggbiplot | A biplot based on ggplot2 | Data Visualization library

I'm currently running principal component analysis. For the interpretation I want to create a profile (pattern) plot to visualize the correlation between each principal component and the original variables. Is anyone familiar with a package or code to create this in R? I'm using the prcomp() function in R.

See examples:

https://canadianaudiologist.ca/predicting-speech-perception-from-the-audiogram-and-vice-versa/ https://blogs.sas.com/content/iml/2019/11/04/interpret-graphs-principal-components.html

This is similar data to my db:

I need some help to modify my R script that could assign different color to each group. I always need to specify color name/code manually [e.g., scale_fill_manual(values = c("grey", "red", "yellow", "green"))].

Complete R script is here

I need some help to get PC3 and PC4 for my data.

I have a PCA plot that ive been working on a while now (I am not very good at R but this is teaching me a lot, just trying to make this one plot). I am now at the stage where the plot looks how we want it, now we just want to draw ellipses around the replicates. (see plot below)

Ideally, the ellipses would have the border colour the exact same colour and be filled in with a mostly transparent version of that same colour (see example plot below)

Many people have suggested using ggbiplot, which i tried but i found the manual very sparse in detail and i didnt really understand how to control the plot. To generate the plot i used the following code:

I use ggbiplot regularly and can control every aspect of the plot produced using ggplot2 tools, since it inherits from ggplot2...

In ggplot2, I usually control the number of columns in the legend with a line of the form:

How can I added the Species classifier next to the centroid in this PCA plot.

Code:

I am struggling with changing the loading labels of biplot. I used prcomp function to run PCA and ggbiplot. I've tried to change it with loading.label.label = c("a","b","c","d"), but it didn't work. I have attached my code and plot below.

I would like to plot the distribution of multiple columns of my data set. It has over 820.000 rows and 18 columns. I want to plot all columns except the columns with the dummy variables. I have already been able to create a graphic. But I want to have the values of the x-axis on the y-axis because these are the column values and I want to display their distribution for each column.