population vector in pcadapt

[abilska]

Hello, we would like to analyze SNPs under selection across different populations (in our case, we have six), and the sample sizes (belove) provided are just examples. Unfortunately, we keep encountering an error, and it seems that the term related to different populations (pop_vector) has been deleted from the program.
Could you please help us understand what we might be doing wrong? Is this population selection option available in another feature/version?

Best regards,

################
Load required libraries
library(pcadapt)

# Load your SNP data (ensure it's in a suitable format)
# Assuming the file 'snp_data.csv' is structured with samples as rows and SNPs as columns
snp_data <- read.csv("snp_data.csv", row.names = 1)

# Create a population vector based on the specified sizes
pop_vector <- c(rep(1, 20),  # Population 1
                rep(2, 30),  # Population 2
                rep(3, 10),  # Population 3
                rep(4, 40),  # Population 4
                rep(5, 8),   # Population 5
                rep(6, 11))  # Population 6

# Run pcadapt analysis
pcadapt_result <- pcadapt(snp_data, K = 6)  # K is the number of populations

# Plot the results
plot(pcadapt_result)

# Extract SNPs under selection based on a defined threshold
# Set an appropriate threshold; this example uses a placeholder
threshold <- 0.05  # Adjust this value based on your criteria
snp_selection <- which(pcadapt_result$loadings[,1] > threshold)

# Print SNPs under selection
print(snp_selection)

# If you want to visualize the loadings
plot(pcadapt_result, type = "scores")

[privefl]

What is the error you get exactly?