hello all, I am new to GEOquery and things are not very clear to me as my understanding is limited for the package. I need to know phenotypic data for GSE21653 (the example of one tutorial).in the example of tutorial there are 266 samples,but by this function:

dim(pData(gse[[1]]))

head(pData(gse[[1]])[, 1:3])

the result in R is:

dim(pData(gse[[1]]))

[1] 266 47

head(pData(gse[[1]])[, 1:3])

      title geo_accession                status

GSM540108 BC1 GSM540108 Public on May 05 2010

GSM540109 BC2 GSM540109 Public on May 05 2010

GSM540110 BC3 GSM540110 Public on May 05 2010

GSM540111 BC4 GSM540111 Public on May 05 2010

GSM540112 BC5 GSM540112 Public on May 05 2010

GSM540113 BC6 GSM540113 Public on May 05 2010

where are other samples?

another question, for another experiment how can I make this function?