The following function shows one possible approach to solve this problem . I assume that the dataset contains 1000x1 double values as a table.
function res= f(dataset, Vs)
y=dataset.var; %var is name of the column containing elements
n=length(y);
res= sqrt(sum((y-Vs).^2)./n);
end
