[R] strange error : isS4(x) in gamm function (mgcv package). Variable in data-frame not recognized???

[Joris Meys]

Follow up:

I finally succeeded to more or less reproduce the error. The origin
lied in the fact that I accidently loaded a function while being in
browser mode for debugging that function. So something went very much
wrong with the namespaces. Teaches me right...

Cheers
Joris

On Wed, Jul 28, 2010 at 2:27 PM, Joris Meys <jorismeys at gmail.com> wrote:
> Dear all,
>
> it gets even more weird. After restarting R, the code I used works
> just fine. The call is generated in a function that I debugged using
> browser(). Problem is solved, but I have no clue whatsoever how that
> error came about. It must have something to do with namespaces, but
> the origin is dark. I tried to regenerate the error, but didn't
> succeed.
>
> Somebody an idea as to where I have to look for a cause?
>
> Cheers
> Joris
>
> On Wed, Jul 28, 2010 at 1:16 PM, Joris Meys <jorismeys at gmail.com> wrote:
>> Dear all,
>>
>> I run a gamm with following call :
>>
>> result <- try(gamm(values~ s( VM )+s( RH )+s( TT )+s( PP
>> )+RF+weekend+s(day)+s(julday) ,correlation=corCAR1(form=~ day|month
>> ),data=tmp) )"
>>
>> with mgcv version 1.6.2
>>
>> No stress about the data, the error is not data-related. I get :
>>
>> Error in isS4(x) : object 'VM' not found
>>
>> What so? I did define the dataframe to be used, and the dataframe
>> contains a variable VM :
>>
>>> str(tmp)
>> 'data.frame':   4014 obs. of  12 variables:
>>  $ values : num  73.45 105.45 74.45 41.45 -4.55 ...
>>  $ dates  :Class 'Date'  num [1:4014] 9862 9863 9864 9865 9866 ...
>>  $ year   : num  -5.65 -5.65 -5.65 -5.65 -5.65 ...
>>  $ day    : num  -178 -177 -176 -175 -174 ...
>>  $ month  : Factor w/ 156 levels "1996-April","1996-August",..: 17 17
>> 17 17 17 17 17 17 17 17 ...
>>  $ julday : num  -2241 -2240 -2239 -2238 -2237 ...
>>  $ weekend: num  -0.289 -0.289 -0.289 0.711 0.711 ...
>>  $ VM     : num  0.139 -1.451 0.349 0.839 -0.611 ...
>>  $ RH     : num  55.2 61.4 59.8 64.1 60.7 ...
>>  $ TT     : num  -23.4 -23.6 -19.5 -16.1 -15.3 ...
>>  $ PP     : num  6.17 4.27 -4.93 -9.23 -2.63 ...
>>  $ RF     : Ord.factor w/ 3 levels "None"<"<2.5mm"<..: 1 1 1 1 1 1 1 1 1 1 ...
>>  - attr(*, "means")= num
>>
>> Any idea what I'm missing here?
>>
>> Cheers
>> Joris
>>
>>
>> --
>> Joris Meys
>> Statistical consultant
>>
>> Ghent University
>> Faculty of Bioscience Engineering
>> Department of Applied mathematics, biometrics and process control
>>
>> tel : +32 9 264 59 87
>> Joris.Meys at Ugent.be
>> -------------------------------
>> Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php
>>
>
>
>
> --
> Joris Meys
> Statistical consultant
>
> Ghent University
> Faculty of Bioscience Engineering
> Department of Applied mathematics, biometrics and process control
>
> tel : +32 9 264 59 87
> Joris.Meys at Ugent.be
> -------------------------------
> Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php
>



-- 
Joris Meys
Statistical consultant

Ghent University
Faculty of Bioscience Engineering
Department of Applied mathematics, biometrics and process control

tel : +32 9 264 59 87
Joris.Meys at Ugent.be
-------------------------------
Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php