Kabeli,
You seem to be doing your best to avoid working
your way through some introductory documentation
like An Introduction to R, which is sitting
right there on your computer.
So let's try to take it slowly, one step at a time.
#1. generate a vector of values to work with.
(forget the matrix() bit and the data.frame() bit
for now and forget about the sampling package
for now)
x <- rnorm(100, mean=50000, sd=5000)
#2. in order to select a random sample of values
from x, we'll randomly select 20 numbers from
the set {1, 2, 3, ..., 99, 100}. Those 20 numbers
will tell us which x-values get into our sample.
idx <- sample(100, 20, replace = FALSE)
Here's what I got (you'll get different numbers):
> idx
[1] 91 51 97 50 67 31 29 75 89 78 63 77
[13] 38 85 16 74 53 40 54 2
#3. now select x[91], x[51], etc as the elements
of x to be your sample:
xsample <- x[idx]
#4. now do whatever you want with xsample;
sum(xsample)
mean(xsample)
sd(xsample)
#5. in order to use something like the function
strata() in pkg:sampling, you will first have
to understand how to read its help page and
understand what's in the object returned by strata().
I do hope that this helps at least somewhat,
Peter Ehlers
KABELI MEFANE wrote:
Dear Mr Winsemius
I am sorry to have offended any of you by the mistakes i made. The package i loaded is sampling and there was an unwanted comma between size c(20, )and the bracket. What i wanted was to calculate the sum of H in a sample not in the original dataframe. If i do
sum(H) i get the sum value of H's in the original dataframe.
Load package sampling
H <- matrix(rnorm(100, mean=50000, sd=5000))
sampleframe=data.frame(type=c(rep("H",100)),value=c(H))
sampleframe
sum(H)
str=strata(sampleframe,c("type"),size=c(20), method="srswor")
sample.strat<-getdata(sampleframe,str)
sample.strat
Thanks for the input. Once again sorry for wasting your time.
Best Regards
--- On Fri, 2/10/09, David Winsemius <dwinsem...@comcast.net> wrote:
From: David Winsemius <dwinsem...@comcast.net>
Subject: Re: [R] How to use Subpopulation data?
To: "KABELI MEFANE" <kabelimef...@yahoo.co.uk>
Cc: R-help@r-project.org
Date: Friday, 2 October, 2009, 3:38 PM
On Oct 1, 2009, at 6:06 AM, KABELI MEFANE wrote:
Dear Helpers
I have a sample frame and i have sampled from it using three methods and now i
want to calculate the statistics but i only get the population parameters.
H <- matrix(rnorm(100, mean=50000, sd=5000))
sampleframe=data.frame(type=c(rep("H",100)),value=c(H))
sampleframe
str=strata(sampleframe,c("type"),size=c(20,), method="srswor")
sample.strat<-getdata(sampleframe,str)
sample.strat
If you want the number of rows in sample.strat then length(H) is the wrong
approach since that is the original (unsampled) object.
length(H)
i get:
length(H)
[1] 100
Desire to get:
length(H)
[1] 20
I cannot tell what packages you have loaded and strata is not in the sampling package
which I guessed (wrongly) was where you were getting "getdata". When you post
code you should precede that code with calls that load any non-base packages.
In later posting you ask for ways to calculate "the sum" but you do not say
what it is that you want the sum of.... . Our abilities to read minds is extremely
limited.
--David Winsemius
[[alternative HTML version deleted]]
------------------------------------------------------------------------
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
