Message-ID: <577cb4ea4b334b8ba7f847891eeb7a1d@SRVEXCHCM1301.precheza.cz>
Date: 2022-03-10T09:55:04Z
From: PIKAL Petr
Subject: NA's in segmented
In-Reply-To: <VI1PR06MB5069823FEB2F232B5EF356569F099@VI1PR06MB5069.eurprd06.prod.outlook.com>
Hi
1. Do not use HTML formated mail, your message is scrambled
2. With your code only you can get result all of us get error
> Mod = lm(as.numeric(data[,"meanMSD25"])~(as.numeric(data[,"Slice"])), weights=1*sqrt(as.numeric(data[,"Detec"])))
Error in data[, "meanMSD25"] :
object of type 'closure' is not subsettable
Try to provide reproducible example.
3. with lm it is not necessary to use data[,"meanMSD25"], you should provide your data object to data option in lm call (it is silly to call your data data - fortune("dog") applies here)
4. why do you use as.numeric(data[,"meanMSD25"]), shouldn't the column be numeric itself?
So my wild guess is that you read the data into R wrong way and instead of numeric they are character, which could be one source of your problem.
Cheers
Petr
> -----Original Message-----
> From: R-help <r-help-bounces at r-project.org> On Behalf Of M?lina Cointe
> Sent: Tuesday, March 8, 2022 9:45 PM
> To: r-help at r-project.org
> Subject: Re: [R] NA's in segmented
>
> Hi,
>
> I?m contacting you because I have some trouble with the slope() function of
> segmented. When I plot the predicted value everything seems to have worked
> well. But when I use slope(), the slope for the first segment is a line with 0
> and Nas? Also I can get a negative slope whereas on the graph it?s clearly a
> positive slope?
>
> Here are the lines I?m using :
>
> Mod = lm(as.numeric(data[,"meanMSD25"])~(as.numeric(data[,"Slice"])),
> weights=1*sqrt(as.numeric(data[,"Detec"])))
> x = (as.numeric(data[,"Slice"]))
> o <- segmented(Mod, seg.Z=~x, psi=NA, control=seg.control(display=FALSE,
> K=2))
>
> Thanks in advance for your help,
>
> Best regards,
>
> M?lina COINTE
>
>
>
> [[alternative HTML version deleted]]