About Hanyu

Hanyu · ‎10-12-2015

I want to construct a template for 3 by 2 graphs using proc greplay. I don't quite understand the syntax of PROC GREPLAY and TDEF statement. I have found a template of 2 by 3 here http://support.sas.com/kb/44/973.html. Could anyone help me with the template language?

Hanyu · ‎10-09-2015

proc iml; start corr_mat(corr,p);/*p is the number of variables which is the dimension of the objective p*p correlation matrix*/ /*corr is the vector containing all nontrivial correlations in order*/ w=p:2; v=(1||cusum(w)+1); do i=1 to ncol(v); s=v[i]; corr=insert(corr,1,0,s); end; corr_mat=sqrvech(corr); return(corr_mat); finish; quit; Thank you Rick! Would the above code also works? Suppose we know how many variables we want, say p=7 variables, and we have the correlation vectors. Do you think this solution is correct?

Hanyu · ‎09-16-2015

I have read the Dr. Rick's post about how to extract nontrivial correlation from a correlation matrix. http://blogs.sas.com/content/iml/2012/08/16/extract-the-lower-triangular-elements-of-a-matrix.html My question is that I want to construct a correlation matrix from a vector only containing nontrivial corrleations. The vector only contains those correlations that are not 1. But I want to construct a correlation matrix whose diagonal elements are all 1.

Hanyu · ‎07-09-2015

When I use the proc sgpanel procedure, I want to remove the panel frame surrounding the panelby variable. Does anybody know how? Here is an example from Sanjay@SAS. proc sgpanel data=sashelp.cars(where=(type ne 'Hybrid')); panelby origin / onepanel rows=1 novarname layout=columnlattice spacing=20 colheaderpos=bottom noborder; scatter x=Type y=mpg_city / jitter markerattrs=(size=5); rowaxis display=(noline) valueattrs=(size=7) grid; colaxis display=(noline nolabel) valueattrs=(size=7) fitpolicy=none; run; Notice the frame surrounding 'Asia' 'Europe'. I want to remove that frame. Thank you.

Hanyu · ‎04-24-2015

Thank you very much. I think read({name sex}) would also work.

Hanyu · ‎04-24-2015

Thank you Dr. Regarding my first question, I am now trying to read two variables into one matrix. I am a bit struggling with this operation. proc iml; start read(s); use y; read all var s into su; close y; finish; run read("a b"); quit; If I define s to be something like "a b" or "a,b", sas say the variable a b is not in the dataset. Could you help me?

Hanyu · ‎04-24-2015

Thank you Dr. May I ask one more question? I want to also create dataset according to arguments of a module. My attempt is the following proc iml; start store_est(a,b); estimate=T(1:30); std=T(1:30); display=estimate||std; create estimate_a_b from display; append from display; close estimate_a_b; finish; store module=store_est; run store_est(s1,s2); However, the code does not work. Is there a way to parametrize the dataset name in the module?

Hanyu · ‎04-24-2015

Can we define a module which has variable names as arguments and reads the variables into a matrix? I tried the following code and it does not work. Suppose the variable "a" is a variable in dataset "y" and I want to read a into matrix su. The reason I want to do this is that I have many variables in y and I want to repeat my estimation for different variable. Thank you in advance. proc iml; start read(s); use y; read all var{s} into su; close y; finish; run read(a); quit; ERROR: s is not in the scope of variables for the data set.

Hanyu · ‎02-25-2015

Dr. I have defined a module that returns a vector containing individual likelihood value at each observation. Then I want to compute the gradient using NLPFDD but SAS reports that the module must return a scalar.

Hanyu · ‎02-25-2015

Thank you Dr.

Hanyu · ‎02-25-2015

As is well known that quasi maximum likelihood estimation requires computing first derivative of the individual likelihood at each observation. Does anybody know how to implement it in SAS?d

Hanyu · ‎02-05-2015

When using Newton Raphson algo, we better provide first derivative to the routine. However, I can only provide some of the analytical formula for the derivatives and I want to approximate other derivatives by finite difference. Is this possible for NLPNRA routine?

Hanyu · ‎02-02-2015

Hi Dr. Wicklin, I have read your blog. You use the following statement mu = param[1]; to indicate that mu is a parameter. I need to estimate many parameters. Do I need to define each of them as you define mu?

Hanyu · ‎02-02-2015

I have a very large system that has more than 100 parameters. I need to write a module to compute the value of the likelihood function. The first step is to define the parameter vector. Usually parameter is defined as a=parm[1];b=parm[2]. But if I define each parameter like that, it would take too long. I want to follow some rule to define the parameter vector, e.g. according to the data of the company name(microsoft apple sas...). Could any body give me a hint?

Hanyu · ‎01-29-2015

I try to locate this negative definite matrix in the matrix R using loc function. I cannot find any matrix has these values. I cannot understand how is this possible.

Online Status	Offline
Date Last Visited	‎09-19-2024 01:33 AM

Re: Transfer status window

Re: Transfer status window

Re: Transfer status window

Re: Transfer status window

Transfer status window

Re: Simulation and further calculation

Re: Simulation and further calculation

Re: Simulation and further calculation

Simulation and further calculation

Re: PROC Severity and Weibull Distribution

Re: Transfer status window

Re: PROC Severity and Weibull Distribution

Re: PROC Severity and Weibull Distribution

Re: define a module to read different variable in one dataset

Re: plotting two time series that shows alternatively using proc gplot

Re: Transfer status window

Re: plotting two time series that shows alternatively using proc gplot

template for 3 by 2 graphs

Re: construct correlation matrix from correlation vectors

construct correlation matrix from correlation vectors

proc sgpanel removing panel frame of panelby variable

Re: define a module to read different variable in one dataset

Re: define a module to read different variable in one dataset

Re: define a module to read different variable in one dataset

define a module to read different variable in one dataset

Re: computing derivative numerically at each observation

Re: computing derivative numerically at each observation

computing derivative numerically at each observation

SAS IML NLPNRA rountine

Re: Defining a very long parameter vector

Defining a very long parameter vector

Re: multivariate GARCH model using SAS IML