?subset
subset(data, column == " E")
On Nov 26, 2007 3:07 PM, Katherine Jones <kajones at connect.carleton.ca>
wrote:> Hi,
>
> I am having trouble with a large dataset I am importing from SPSS.
> The problem is I have to merge two datasets (which seems to be
> working OK) then select rows based on attributes. I have a column
> with either blank cells, B or E entered. I want to select all rows
> with E. I have other columns with numerical data which I will then do
> analyses on.
> data[column==" E"] does not work. I use " E" not
"E", because levels
> (column) returns " " " " " B" " E".
>
> Any help on what I am doing wrong is much appreciated. I'm getting
> quite stressed as I have 10 files with approx 100,000 records in each
> to analyse so manipulating data becomes a pain.
>
> Here is the code below, not sure it makes much sense without seeing
> the dataset:-
>
>
> chaff<-read.spss("/Users/Kat/Desktop/papers in progress/btopaper/
> edited BTO data/fatnewchaff.sav", to.data.frame=TRUE)
> chafffat<-read.spss("/Users/Kat/Desktop/papers in
progress/btopaper/
> edited BTO data/fatmethods.sav")
> chaffmerge2<-merge(chaff, chafffat, by.x=c("RINGNO",
"FAT",
> "FATMTD"), by.y=c("RINGNO", "FAT",
"FATMTD"), all=T)
> attach(chaffmerge2)
> chaffhabfactor<-factor(chaffmerge2$HYBRID_A)
> levels(chaffhabfactor)
> Echaff<-chaffmerge2[FATMTD==" E",]
> attach(Echaff)
> names(Echaff)
> plotmeans(Echaff$FAT~Echaff$HYBRID_A)
> chaffFat<-factor(Echaff$FAT)
> levels(chaffFat)
> chaffzeros<-table(chaffFat, Echaff$HYBRID_A)
> chaffzeros
>
> ****
> chaffFat 1 2 3 4 5
> 0 261 354 345 1003 235
> 1 38 23 17 6 2
> 2 19 0 4 2 0
> 3 7 0 1 0 1
> 4 2 0 0 0 0
> 5 145 34 123 100 60
> 8 0 0 0 0 0
> 10 202 141 248 279 101
> 15 73 12 79 51 9
> 20 84 60 64 133 19
> 25 14 6 20 22 3
> 30 30 25 22 54 13
> 35 3 0 7 4 4
> 40 7 10 2 12 5
> 45 2 0 3 1 0
> 50 1 0 0 2 1
> 60 0 1 0 1 1
> ****
> The 1,2,3,4,5, values of chaffFat above correspond to "B" which
> should have been removed!!!!
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.
>
--
Jim Holtman
Cincinnati, OH
+1 513 646 9390
What is the problem you are trying to solve?