[R] nlminb problem, convergence error code = 1

2012-10-07 Thread zugi young
Hi, I was trying to do a multi-level modeling of the data using the nlme package. I encounter the following message: > addRandomSlopeRange<- lme(FAN.range ~ CHN.range + Age + Block.T, data=data.block, random = ~CHN.range|CHI, method = "ML", na.action=na.omit) Error in lme.formula(FAN.range ~ CHN.

[R] plotting with a symbol on every nth point

2011-10-21 Thread zugi young
Hi, I would like to produce a plot with a symbol on every nth point in a time series data, like the one in the following: http://www.phon.ucl.ac.uk/home/yi/ProsodyPro/EnglishFocus.png x <- seq(-100,1000,25) plot(x,type="l") Could someone help me out with the above example? Thanks

[R] reporting ANOVA for nested models

2011-09-07 Thread zugi young
I have the following results for an ANOVA comparing two nested models. I wasn't sure how I am supposed to report this result in the area of psychology. Specifically, am I supposed to report the DF's or just the F ratio? I could manually calculate the degrees of freedoms, but there must be a reason

[R] breakpoints in rqss()

2010-05-24 Thread zugi young
Dear list, I used rqss() in quantreg package for a piecewise linear regression. Can someone tell me how to find the x values corresponding to the breakpoints and the slopes for the phases before and after the breakpoints? I searched the list and gather that there is another package "segmented" th