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 <JRoll@lcog.org> wrote:

>  Sorry Just a generic list
>
> Is<-list()
>
> forgot to add that from my actual code
>  ------------------------------
> *From:* jim holtman [mailto:jholtman@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 <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?
>



-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem that you are trying to solve?

	[[alternative HTML version deleted]]

