[R] Odp: error in my selection

[Petr PIKAL]

Hi

Is it possible that there is a problem with attaching data? Look what 
search() tels you.

Regards
Petr


r-help-bounces at r-project.org napsal dne 15.01.2008 18:38:36:

> 
> Hi everybody,
> 
> Well, I really hoped that until now i know how to do a selection ..... 
but 
> obviously i don't. So risking your laughing at me - here you go.
> 
> I have a table with lots of vegetation data, and one column specifies in 
which
> region the data is. There are only 2 regions, region 1 and region 3. I 
want to
> select all data that belong to region 1 and after do some xtabs on it. 
If i do
> everything correctly i thought that i will have only the data from 
region one 
> in my result .... it seems not. The uniqueID - is the location name that 
has 
> incorporated into it the region as well, cat3_name is the category name 
and 
> Total is the actual numeric data i am interested in tabulating. There 
are few 
> other columns but for my immediate purpose they are not of interest now.
> 
> Code (and i hope it does not come scrambled):
> ----------------------
> 
> > x3 <- read.csv("unique_veg_id.csv")
> ### columns in my data:
> 
> > dimnames(x3)[[2]] [1] "LocatID"   "Region"    "Total"     "LocID" 
"UTMX"
> "UTMY"      "Class"     "cat3_name" "EventID"   "day_sampl" "uniqueID" > 
> 
> x3_reg1 <- x3[x3$Region == "JELA_reg1",]> > ma1 <- ftable(xtabs(Total ~ 
> uniqueID + cat3_name, data=x3_reg1))> ma1  
> cat3_name Anthropogenic Asteraceae - herb Bare ground Canopy tree Coarse 
woody
> debris Cypress knees Emergent aquatic Epiphyte Floating aquatics Grasses 
Herbs
> Ilex Midstory tree Palm Sapium sebiferum Sedge Shrub Terrestrial herb 
Vine 
> WateruniqueID  
> JELA_reg1_1_2_6554_2007Oct24_Class1                       0    0 
> 0           0                   0             0               25 0  !
>           0       5     0    0             0    0                0    55 
    0
> 10    0     0JELA_reg1_1_3_6985_2007Oct17_Class2                       0 
 
> 15           0           0                   0             0   15 
> 0                 0       0     5    0             0    0 0 
> 50     7                0   10     0.
> .
> .
> .
> 
> JELA_reg3_8_1_8015_2007Nov08_Class1                       0    0 
> 0           0                   0             0                0 0 
> 0       0     0    0             0    0                0     0     0   
> 0    0     0JELA_reg3_9_2_8216_2007Nov08_Class1                       0  
 
> 0           0           0                   0             0   0 
> 0                 0       0     0    0             0    0 0 
> 0     0                0    0     0> 
> 
> -----------------------
> 
> As you see instead of having only the uniqueID locations that belong 
only to 
> region 1, i have all locations for both regions, 1 and 3. If i save on 
disk 
> x3_reg1 and after read in back with read.csv and do the same xtabs on 
that i 
> do have only the locations belonging to region 1. Do i really need this 
extra 
> step? What i am doing wrong in my selection?
> 
> Thank you for your input,
> 
> Monica
> _________________________________________________________________
> 
> 
> GLM_CPC_VideoChat_distantfamily_012008
>    [[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.