Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
From | Janet Hill <janethill73@yahoo.co.uk> |
To | "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu> |
Subject | Re: st: RE; Summarise similar curves |
Date | Wed, 24 Apr 2013 11:56:13 +0100 (BST) |
Nick, Many thanks, as usual. I will investigate the interpolations. Janet ----- Original Message ----- From: Nick Cox <njcoxstata@gmail.com> To: "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu> Cc: Sent: Wednesday, 24 April 2013, 11:48 Subject: Re: st: RE; Summarise similar curves I can't judge from your description how well a restricted cubic spline will work. I am happy to regard -rcspline- (SSC) as a useful command for what it does. But the gaps in your data do not prohibit much if it's reasonable to interpolate to a regular grid. Check out -ipolate- (official) -cipolate- (SSC) -csipolate- (SSC) -nnipolate- (SSC) -pchipolate- (SSC) All these commands support applications to panel data. A loop foreach m in ci csi nni pchi { ssc desc `m'polate } will bring up capsule descriptions of the unofficial commands. Nick njcoxstata@gmail.com On 24 April 2013 11:34, Janet Hill <janethill73@yahoo.co.uk> wrote: > I have a dataset which should have consisted of a patient characteristic measured at 6 week intervals over a two year period. There are 400 patients in three groups and I intended to graphically summarise the important features of the groups using the principal component method described by Jones,MC and Rice,JA, 'Displaying the important features of large collections of similar curves, American Statistician,46(1992) 140-145'. Unfortunately the patients did not always turn up on the specified date so I do not have an equispaced dated grid thus precluding the use of this method. > > I am now considering using a restricted cubic spline, (rcspline from http://fmwww.bc.edu/RePEc/bocode/r) to show a smoothed curve for each group. > > My questions are a) is this a reasonable approach or b) is there anything else in the Stata armamentarium which is more suitable. I am using Stata 12. * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/