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 <jroll@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]]