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? [[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.
