Try this on your real data: > #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 Area Proportions > MultiLotBldgs..$Prop <- ave(MultiLotBldgs..$Area, MultiLotBldgs..$Bldgid, + FUN=function(x) x / sum(x)) > > # find not too small > notTooSmall <- !((MultiLotBldgs..$Area <= 45) | ((MultiLotBldgs..$Area > 45) & + (MultiLotBldgs..$Prop < 0.05))) > > MultiLotBldgs2.. <- MultiLotBldgs..[notTooSmall,] > # print out results > MultiLotBldgs2.. Bldgid Maplot Area Prop 2 1000 20001 170 0.8095238 3 1001 30000 50 1.0000000 4 1002 30001 100 1.0000000 5 1003 40000 100 0.9532888 > > >
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.
