It is a lot more complicated using manifold methods (esp with nearest neighbors).
The dimRed package can compute both of those (and more) but I think that there is no LLE projection for new data
For IsoMap, there is a recipe step that can be used to produce features for new data.
Relatedly, the GH version of embed has a recipe step for UMAP that projects for new data (supervised and unsupervised).
One last thing... for this small sample size, you might want to stick to simple linear rotational methods like PCA. I would be worried about how well these nonlinear projections would generalize with 39 samples.