[R] Error in Quantile function

irishmhw ni.weir.ni at gmail.com
Fri Nov 14 15:33:13 CET 2008


That worked the charm, the line that you suggested age me a bad atomic number
error but I just assigned the spread$bootdataframe.age to a new vector and
it works great.  Thanks you have been a great help.  Best in all to you and
yours. 



David Winsemius wrote:
> 
> 
> On Nov 13, 2008, at 7:56 PM, irishmhw wrote:
> 
>>
>> This is strange the first function that I am doing this for is  
>> saying in the
>> typeof() is a double for some reason this is coming up as a list.
> 
> Well, dataframes _are_ lists.
> 
>>
>>> str(spread)
>> 'data.frame':   18 obs. of  1 variable:
>> $ bootdataframe.age: num  1.46e-09 1.67e-09 1.46e-09 1.45e-09  
>> 1.37e-09 ...
> 
> I'm pretty sure you cannot send an entire dataframe to quantile.  
> Pretty sure you need to give a vector argument. So what happens when  
> you try:
> 
>   exp_CIs <- quantile( spread$bootdataframe.age,  
> probs=c(0.01,0.05,0.95,0.99) )
> 
> -- 
> David Winsemius
> Heritage Labs
> 
>>
>>> typeof(spread)
>> [1] "list"
>>> class(spread)
>> [1] "data.frame"
>>>
>>
>>
>> David Winsemius wrote:
>>>
>>> Please show us what you get when you execute:
>>>
>>> str(spead)
>>> typeof(spread)
>>> class(spread)
>>>
>>> It is doubtful that either of your theories is correct. It's more
>>> likely that spread is a more complex object than you realize. We need
>>> to know what spread looks like to the R interpreter in the  
>>> environment
>>> where it is being called.
>>>
>>> -- 
>>> David Winsemius
>>> Heritage Labs
>>>
>>>
>>> On Nov 13, 2008, at 3:29 PM, irishmhw wrote:
>>>
>>>>
>>>> If anyone can assist with this problem you have my great thanks:
>>>>
>>>> I am trying to establish and plot confidence intervals on a
>>>> bootstrapped
>>>> function.  I have a more complicated function that has no problems
>>>> with
>>>> determining the confidence intervals using the quantile command.
>>>> This is
>>>> outside the bootstrap portion of the code that is working fine it is
>>>> just
>>>> determining everything for the more complicated model, not even sure
>>>> what
>>>> the error is telling me.
>>>>
>>>>>  spread <- expntl.dr(plot_k_list,pdose)
>>>>>  exp_CIs <- quantile(spread,probs=c(0.01,0.05,0.95,0.99))
>>>> Error in order(list(bootdataframe.age = c(1.56064428047387e-09,
>>>> 1.63862234892065e-09,  :
>>>> unimplemented type 'list' in 'orderVector1'
>>>>>
>>>>
>>>> I am not sure if this is problem with the quantile command trying to
>>>> order
>>>> the data or if the values in the spread are too close together.
>>>>
>>>> Thank you for any assistance you can offer.
>>>>
>>>> Mark
>>>> --
>>>
>>> ______________________________________________
>>> R-help at 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.
>>>
>>>
>>
>> -- 
>> View this message in context:
>> http://www.nabble.com/Error-in-Quantile-function-tp20488954p20492856.html
>> Sent from the R help mailing list archive at Nabble.com.
>>
>> ______________________________________________
>> R-help at 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.
> 
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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.
> 
> 

-- 
View this message in context: http://www.nabble.com/Error-in-Quantile-function-tp20488954p20501578.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list