You might want to put some debugging statements to see what the values are. As you say the example seems to work, so look at the data:
> MultiLotBldgs2.. <- MultiLotBldgs..[Is$NotTooSmall.X, ] > # some debugging statements > str(MultiLotBldgs2..) 'data.frame': 4 obs. of 3 variables: $ Bldgid: num 1000 1001 1002 1003 $ Maplot: num 20001 30000 30001 40000 $ Area : num 170 50 100 100 > str(Is) List of 1 $ NotTooSmall.X: Named logi [1:6] FALSE TRUE TRUE TRUE TRUE FALSE ..- attr(*, "names")= chr [1:6] "10001" "10002" "1001" "1002" ... > print(sum(Is$NotTooSmall)) [1] 4 > > On Wed, Feb 17, 2010 at 6:58 PM, ROLL Josh F <jr...@lcog.org> wrote: > Sorry Just a generic list > > Is<-list() > > forgot to add that from my actual code > ------------------------------ > *From:* jim holtman [mailto:jholt...@gmail.com] > *Sent:* Wednesday, February 17, 2010 3:58 PM > *To:* ROLL Josh F > *Cc:* r-help@r-project.org > *Subject:* Re: [R] Procedure not working for actual data > > Your example does not work since "Is" is not defined. What is it > supposed to be? > > On Wed, Feb 17, 2010 at 6:34 PM, LCOG1 <jr...@lcog.org> wrote: > >> >> Hello all, >> I have what i feel is a unique situation which may not be resolved with >> this inquiry. I have constructed the below data set so that i may give an >> example of what im doing. The example works perfectly and i have no >> issues >> with it. My problem arises with my actual data, which includes another 11 >> columns of data (used in later analysis) and a total of about 7000 >> cases(rows). i mention the dimensions of the actual data because im >> wondering if my below process would encounter problems with more data. >> To be sure the problem occurs in the last step. Is$NotTooSmall gives me >> a >> binary output that is then put back in MultiLotBldgs.. (as shown in the >> example) to return the cases i want to keep. >> In my actual data the binary designation is correct but when >> MultiLotBldgs2.. returns it doesnt remove the cases that are False in >> Is$NotTooSmall. Like i said my sample data works fine but my actual >> implementation does not. Any suggestions? I know this is not easy to >> answer without seeing the problem but this is the best i can do without >> sending you all of my data. >> >> Cheers, >> JR >> >> >> >> >> #Sample data >> Bldgid<-c(1000,1000,1001,1002,1003,1003) >> Maplot<-c(20000,20001,30000,30001,40000,40001) >> Area<-c(40,170,50,100,100,4.9) >> #Construct Sample dataframe >> MultiLotBldgs..<-data.frame(Bldgid,Maplot,Area) >> #Get Building Areas >> MultiLotBldgArea.X <- unlist(tapply(MultiLotBldgs..$Area, >> MultiLotBldgs..$Bldgid, >> function(x) x)) >> >> # Calculate the proportion of the total building area in each piece of the >> building >> MultiLotBldgProp.X <- unlist(tapply(MultiLotBldgs..$Area, >> MultiLotBldgs..$Bldgid, >> function(x) x/sum(x))) >> >> #Identify buildings that should be considered for joining >> Is$NotTooSmall.X <- !(((MultiLotBldgArea.X <= 45) | >> ((MultiLotBldgArea.X > 45) & >> (MultiLotBldgProp.X >> < 0.05)))) >> >> MultiLotBldgs2.. <- MultiLotBldgs..[Is$NotTooSmall.X, ] >> >> -- >> View this message in context: >> http://n4.nabble.com/Procedure-not-working-for-actual-data-tp1559492p1559492.html >> Sent from the R help mailing list archive at Nabble.com. >> >> ______________________________________________ >> 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<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 that you are trying to solve? > -- Jim Holtman Cincinnati, OH +1 513 646 9390 What is the problem that you are trying to solve? [[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.
