About DaLack

DaLack · ‎07-14-2016

Thanks for your help. I thought about one legend but this solution is a good compromise.

DaLack · ‎07-13-2016

I want to plot the results of two quantile regressions. The output of my quantile regression was a grouped variable "quantile" plus the corresponding values. I tried to create a graph with the following code in SAS 9.4: ods graphics on / labelmax=42100 antialiasmax=42100; proc sgplot data=data1; scatter x=Var1 y=Var2; series x=Var1 y=Var3 / group=quantile legendlabel="Gew" name="G"; series x=Var1 y=Var4 / group=quantile lineattrs=(pattern=dot) legendlabel="OhneGew" name="OG"; keylegend "G" "OG" / location=inside position=topleft title="Quantile" titleattrs=(weight=bold) across=3; run; ods graphics off; The resulting plot is attached. My aim was to label the percentiles in the legend. The first row in the legend should begin with "Gew" and should belong to the first three percentiles (solid line). The second one should start with "OhneGew" and should belong to the other three percentiles (dotted line). Therefore I used the keylegend statement along with the legendlabel statement. But as a result, I only get a legend without the required legendlabels "Gew" and "OhneGew". How can I modify the legend in order to get the legendlabels? Many thanks in advance.

DaLack · ‎07-13-2016

Thanks for your efforts. My sgplot looks like the attached one. The corresponding code for this example is: ods graphics on / labelmax=42100 antialiasmax=42100; proc sgplot data=data1; scatter x=Var1 y=Var2; series x=Var1 y=Var3 / group=quantile legendlabel="Gew" name="G"; series x=Var1 y=Var4 / group=quantile legendlabel="OhneGew" name="OG"; keylegend "G" "OG" / location=inside position=topleft title="Quantile" titleattrs=(weight=bold) across=3; run; ods graphics off; I want to plot the results of two quantile regressions. I a got grouped variable called "quantile" with the 5th, 50th and 95th percentiles. My aim is to label the two outcomes with "Gew" for the three percentiles of the first regression and with "OhneGew" for the three percentiles of the second regression. Therefore I used the keylegend statement along with the legendlabel statement. But as a result, I only get a legend without the description "Gew" and "OhneGew". With respect to the attached sgplot, I want to modify the given legend. The first row of the legend should begin with "Gew" and the second one with "OhneGew". Many thanks for your help!

DaLack · ‎07-13-2016

I am trying to add legendlabels to grouped variables using proc sgplot in SAS 9.4. Therefore I considered the keylegend statement as mentioned below. ods graphics on / labelmax=42100 antialiasmax=42100; proc sgplot data=data1; title "QuantReg"; xaxis label="Var1" type=discrete; yaxis label="Var2"; scatter x=Var1 y=Var2 / markerattrs=(size=3 color=grey); series x=Var1 y=Var3 / group=quantile legendlabel="Gew" name="G"; series x=Var1 y=Var4 / group=quantile lineattrs=(pattern=dot) legendlabel="OhneGew" name="OG"; keylegend "G" "OG" / location=inside position=topleft title="Quantile" titleattrs=(weight=bold) across=3; run; title; ods graphics off; But I only get the legend without a description of the quantiles. Is there another way in order to add labels to the grouped variables? Thanks in advance.

DaLack · ‎07-08-2016

Thanks for the help!

DaLack · ‎07-08-2016

I have found the error. The problem is, that my data set contains around 7k observations. Therefore I have to add "plot(maxpoints=7000)" to the proc statement in order to get fit plots. And this results in missing confidence bands. I have already tried to increase the number but the intervals are still not visible. Is there a solution to this problem? Thanks in advance.

DaLack · ‎07-07-2016

I am using a nonparametric quantile regression with quadratic B-splines. Now I want to add the confidence bands to the fit plot. But I couldn't find an implementation in SAS 9.4. I have already checked the user's guide but got only a solution in SAS 9.2. Therefore I tried the following: ods graphics on; proc quantreg data=data1 algorithm=simplex ci=resampling; effect Var2=spline(Var1 / Basis=Bspline Details degree=2 knotmethod=list(14)); model Var3 = Var2 / plot=fitplot(showlimits) quantile= 0.05 0.5 0.95 seed=1268; output out=data2 pred=p / columnwise; run; ods graphics off; But it doesn't seem to work in SAS 9.4. I always get the fit plot but without any confidence limits. I did the same without a spline effect but again no additional bands. Is there another way in order to display the confidence intervals in SAS 9.4? Many thanks for the answer.

DaLack · ‎07-06-2016

Many thanks for the profound answer. I will consider your comments.

DaLack · ‎07-02-2016

I am using the effect statement in SAS 9.4 in order to estimate a nonparametric quantile regression with quadratic B-splines. I had chosen 2 internal knots and therefore got 5 basis functions. But I am struggling with the interpretation of the output. How can I read the table containing the parameter estimates? What is the meaning of the values? Why do I get an additional row with zeros? My output looks similar to the second table given on this page: http://support.sas.com/kb/57/682.html In this example, there is also an extra row for the parameter of "sWeight7" which contains zeros. Many thanks for your answer.

DaLack · ‎06-03-2016

Many thanks for your support.

DaLack · ‎06-01-2016

I am interested in estimating a distribution given several quantiles (0.05, 0.5, 0.95) and if necessary the mean. The shape should be somehow normal. It exists a useful package called "‘rriskDistributions" in R, explained here https://cran.r-project.org/web/packages/rriskDistributions/rriskDistributions.pdf I am looking for an analogue in SAS 9.4. Does anybody have an idea? Thanks in advance.

DaLack · ‎05-25-2016

Thanks for your answer. Yes, I have read this article and was wondering whether it is already available in SAS. But I agree, programming seems to be the only option to perform it...

DaLack · ‎05-25-2016

I am interested in estimating a mixture of quantile regressions similar to the mixture of linear regressions using proc fmm. I have checked the one-component-case in SAS (proc quantreg) but can't find any solution for the case of a mixture. A mixed model is not possible due to the unkown group membership of the observations. Thanks in advance.

DaLack · ‎05-17-2016

Okay, thanks. I think the problem is my SAS program. I have run your proc freq and the variable named component has only one value. I will try it with another version.

DaLack · ‎05-17-2016

But how can I label mixture components, if the partial option is somehow ignored and the predicted components aren't coincident with my given component membership? Thanks for your help.

Online Status	Offline
Date Last Visited	‎08-09-2016 04:42 PM

Re: How is it possible to add legendlabels to grouped variables using ...

How is it possible to add legendlabels to grouped variables using prog...

Re: How is it possible to add legendlabels to grouped variables using ...

How is it possible to add legendlabels to grouped variables using prog...

Re: How can I add confidence limits to a fit plot using proc quantreg ...

Re: How can I add confidence limits to a fit plot using proc quantreg ...

How can I add confidence limits to a fit plot using proc quantreg in S...

Re: How can I interpret the parameter estimates of the B-spline basis ...

How can I interpret the parameter estimates of the B-spline basis func...

Re: How to estimate a distribution based on quantiles?

Re: How is it possible to add legendlabels to grouped variables using ...

How is it possible to add legendlabels to grouped variables using prog...

Re: How is it possible to add legendlabels to grouped variables using ...

How is it possible to add legendlabels to grouped variables using prog...

Re: How can I add confidence limits to a fit plot using proc quantreg ...

Re: How can I add confidence limits to a fit plot using proc quantreg ...

How can I add confidence limits to a fit plot using proc quantreg in S...

Re: How can I interpret the parameter estimates of the B-spline basis ...

How can I interpret the parameter estimates of the B-spline basis func...

Re: How to estimate a distribution based on quantiles?

How to estimate a distribution based on quantiles?

Re: How to estimate a mixture of quantile regressions in SAS 9.4?

How to estimate a mixture of quantile regressions in SAS 9.4?

Re: How can I get the predicted components in the Iris data set using ...

Re: How can I use the option partial using proc fmm?