Chang ,JX

latex hw tempelate

peacelovekenzie — Tue, 10 Mar 2009 11:32:26 -0400

documentclass[a4paper,12pt,fleqn]{article} %extreport usepackage[top=3cm,bottom=1cm,left=2cm,right=2cm]{geometry} usepackage{fancyhdr,color,xcolor,dsfont} usepackage{listings,ulem,booktabs,graphics,lastpage} usepackage{bm,amsmath,amssymb,amsthm,latexsym} usepackage{graphicx,fouriernc,hyperref} %fouriernc pagestyle{fancy} fancyhf{} lhead{Macro theory IV:HW1} chead{} rhead{ d96627003 Chang,JunXiang quad thepage/pageref{LastPage} } %usepackage{pifont} fouriernc mathpazo concmath(it can’t utilize pdf) mathptmx defbe{begin{enumerate}} % Begin Enumerate defee{end{enumerate}} % End Enumerate defen{item} % ENtry (item) defbi{begin{itemize}} % Begin Itemize defei{end{itemize}} % End Itemize def[#1]{begin{align*}#1end{align*}} DeclareSymbolFont{letters}{OML}{ccm}{m}{it} DeclareMathAlphabet{mathcal}{OMS}{cmsy}{m}{n} renewcommand{labelitemi}{$blacktriangleright$} newcommand{nf}{textrm} newcommand{dsp}{displaystyle} newcommand{intff}{int^{infty}_{0}} newcommand{ert}{e^{-R(t)}} newcommand{ev}{mathcal{V}} newcommand{dpa}{partial} newcommand{lw}{leftarrow} newcommand{llw}{Leftarrow} newcommand{rw}{rightarrow} newcommand{rrw}{Rightarrow} newcommand{ham}{mathcal{H}} newcommand{lam}{pounds} newcommand{grl}{lambda} newcommand{gra}{alpha} newcommand{grb}{beta} newcommand{grd}{delta} newcommand{grt}{theta} newcommand{limu}{lim_{trw infty}} newcommand{cov}{mathsf{Cov}} newcommand{ex}{mathds{E}} newcommand{bx}{pmb{X}} newcommand{bgrb}{pmb{beta}} newcommand{bgrt}{pmb{theta}} newcommand{bA}{pmb{A}} newcommand{bu}{pmb{u}} newcommand{bve}{pmb{varepsilon}} newcommand{lp}{mathds{L}} newcommand{I}{nf{I}} newcommand{iid}{stackrel{mathrm{iid}}{sim}} linespread{1.3} pagecolor{white} color{black} %=========================here start=======================% begin{document}

Speed test

peacelovekenzie — Wed, 31 Dec 2008 02:24:43 -0500

/*Gauss*/ m=ones(200000,1); ret=0; starttime=date; for i (1,20,1); ret=ret+sumc(sumr(m)); endfor; endtime=date; print etstr( ethsec(starttime,endtime) ); /*OX*/ #include main(){ decl tim,i,m,ret,k; m = ones(200000,10); ret =0; tim = timer(); for (i=0;i<20;++i){ ret=ret+sumc(sumr(m)); } print(”time lapsed: “, timespan(tim, timer()), “n”); } #R # File name: loop_rowSums.r my.loop <- 20 m.dim <- list(nrow = 200000, ncol = 10) m <- matrix(1, nrow = m.dim$nrow, ncol = m.dim$ncol) ret <- 0 start <- Sys.time() for(k in 1 : my.loop){ ret <- ret + sum(rowSums(m)) } Sys.time() - start ========================== Speed: OX(0.18) Gauss(0.06) R(0.66) http://www.math.ncu.edu.tw/~chenwc/R_note/index.php?item=loop

mySQL

peacelovekenzie — Wed, 24 Dec 2008 04:51:12 -0500

/* IMPORT DATA CSV*/ PROC IMPORT DATAFILE=”C:DATASTATADASE17HOSP17.CSV” OUT=MYDATA DBMS=CSV REPLACE; GETNAMES=YES; RUN; /*OUTPUT DATA CSV*/ PROC EXPORT DATA=NEW OUTFILE=’c:DataNEW.csv’; /*CREATE NEW TABLE ORDER BY*/ PROC SQL; CREATE TABLE NEW AS SELECT* FROM MYDATA,IDX WHERE PRSN INV ORDER BY TAG; QUIT; /*USE CASE WHEN TO GREATE NEW VARIABLES*/ PROC SQL; CREATE TABLE NEW AS SELECT * , CASE WHEN PRSN LT 10 THEN ‘ZOO’ WHEN PRSN GE 10 THEN ‘BZOO’ END AS PD, CASE WHEN PRSN BETWEEN 1 AND 20 THEN ‘GOOD’ WHEN PRSN BETWEEN 21 AND 30 THEN ‘GOODT’ ELSE ‘NOT_BAD’ END AS DPI FROM MYDATA; QUIT; /*USE IN / NOT IN*/ PROC SQL; CREATE TABLE NEW AS SELECT * FROM MYDATA WHERE PRSN NOT IN (8,11,14,18,24,27,28,44,47,50, 52,56,57,61,62,64,65,78,80,91, 100,109,110,114,117, 124,126,132, 133,136,138, 143,145,148, 150,155,157, 160,164,166,167, 170,173,175,176,178,179, 180,181,183,184,185,186) ; QUIT; /*IN STATA*/ drop if inlist(prsn,8,11,186) /*R*/ indx=c(8,11,14,18,24,27,28,44,47,50, 52,56,57,61,62,64,65,78,80,91, 100,109,110,114,117, 124,126,132, 133,136,138, 143,145,148, 150,155,157, 160,164,166,167, 170,173,175,176,178,179, 180,181,183,184,185,186) two=one[!one$PRSN%in%indx,] /* USE EXCEPT IF TWO TABLES HAVE SAME COLS*/ DATA INDX; INPUT PRSN @@; CARDS; 2 3 4 5 7 8 10 12 14 15 19 31 84 92 ; QUIT; PROC SQL; CREATE TABLE NEW AS SELECT PRSN FROM MYDATA EXCEPT SELECT PRSN FROM INDX; QUIT; PROC PRINT DATA=NEW(OBS=20); VAR PRSN; RUN;

ce que vs ce qui

peacelovekenzie — Mon, 11 Aug 2008 02:37:33 -0400

que / ce que are used when they are the object complement in the relative clause as well as when they stand on behalf of inanimate (Je ne sais pas ce que le garçon a fait. Je connais un secret que je ne révèlerai jamais) qui/ ce qui are used when they are the subject in the relative clause (Je ne sais pas ce qui se passe; j’ai repris du gateau qui est si bon). Qui can also be used on behalf of persons when it is an object complement in the relative clause (je ne sais plus qui je dois croire). ===================================== “ce que” is used when a NOUN follows it “ce qui” is used when a VERB follows it for example: “Je ne sais pas ce que le garçon a fait” (le garçon being a noun) “Je ne sais pas ce qui se passe” (se passer being a verb)

heter job search in life

peacelovekenzie — Thu, 21 Aug 2008 23:37:20 -0400

last night, i read the article, described one’s feeling about job search ,in ppt/job board. the content is that he had no idea about what to do in his first time in searching job as he just graduated from national mba, or put it clearly, he didn’t find his posses as well as real interest/career. the first job is not so good, bz its job content is miscellaneous, not very specific, then he quitted the job. in second time, he didn’t hursh to the job, more patient about it, look on behalf of the more desirable job opportunities. this time, he had more experience about what he really want. at last, he got the satisfactory job in marketing dep. in economics, we utilize job search theory to analyse how long / the prob one find his job. but the story tells some (intelligent) switches on behalf of better position, as well as some switches for randomly. that is, the switching process is hetergeneous. i have no idea whether it is taken into account in literature. if not , it would be good start pt:P zero-tolerance last nigit, i did a little experiment. i put some BTL in the bicycle basket. i hope there's still two BTL, however theory tells us someone will throw down trash into ur basket in high possibility. this morning, i found someone did it.

Le président de la Fed voit l’inflation ralentir

peacelovekenzie — Fri, 22 Aug 2008 00:28:33 -0400

Le président de la Fed voit l’inflation ralentir le contexte économique est actuellement l’un des plus difficiles, a souligné ce vendredi Ben Bernanke, dans un discours lors du traditionnel symposium annuel de la Réserve fédérale. En revanche, la baisse des cours des matières premières et la récente stabilisation du dollar ont été jugés “encourageantes”. Le président de la Fed Ben Bernanke a estimé ce vendredi que la crise financière qui a débuté il y a un an ne s’était pas affaiblie et qu’elle commençait même à toucher le reste de l’économie américaine, créant “l’un des plus difficiles” contextes jamais vus. “Même si le fonctionnement de certains marchés s’est un peu amélioré, la tempête financière (…) n’a pas encore diminué, et ses effets sur l’ensemble de l’économie deviennent apparents sous la forme d’un ralentissement de l’activité économique et une hausse du chômage”, a estimé le président de la Réserve fédérale américaine (Fed). http://www.lesechos.fr/info/marches/300287578-le-president-de-la-fed-voit-l-inflation-ralentir.htm ========================== ralentir:slow down soulinger:underline lors:at the time of En revanche:on the other hand la basisse:fall matières premières:raw material affaiblir:weaken toucher :touch Même si:even though le fonctionnement : mechanism amélioré:improve créer->créant:create vus:voir hausse du chômage: rise of unemployment

???? ????? ?????

peacelovekenzie — Wed, 27 Aug 2008 02:24:14 -0400

???????? ????????? ????? ????? / ????? ????? / Lyrics ????? ????????? ?? ? ??? ??? ?? ??? ??????? ?????? ?? ???? ?????? ????? ?? ???????? ?? ??????? ??????????: “?? ?????? ????” ???? ????? ????? ??? ????????? ??????? ?????? ???????? ????? ????? ????? ? ??????? ??? ?? ??????? ???? ??????? ????????? ?????? ???? ????? ????? ???????? ? ???? ?????? ???????? ???-?? ?????? ???? ?????? ?????? ????? ?????? ???? ???? ? ????????? ????, ??? ??? ???? ?? ???????: “?? ?????? ????” ???? ????? ????? ??? ????????? ??????? ?????? ???????? ????? ????? ????? ? ??????? ??? ?? ??????? ???? ??????? ????????? ?????? ???? ????? ????? ???????? ? ???? ?????? ???? ??????? ???? ?? ??? ???????? ????? ????? ????? ? ??????? ??? ?? ??????? ???? ??????? ????????? ?????? ???? ????? ????? ???????? ? ???? ?????? ?? ? ???? ??????

??????

peacelovekenzie — Sun, 07 Sep 2008 00:48:43 -0400

???? ???????????????????????????????????????????????????????????????????????? ?????????????????????????? “??” ??????yu????????????“???????”??????????? ???????????????????????????????? ?????????, ????????????“???????????”??????“???????????”3????????, ????????????????????????????????????????????????????????????????????????? ?????????????????????“???????”?????????“???”??????“???????”???? “????”??“?”??????“c?????????”???“????????”?????“?????????”?????“??????? ??”???????“?????”????????????“?”???“???”????????“??????”??????????????? ?????????????????“?????”????????????“c”??“????”????????????????? “????????????????”??????“????????”?“????????????”?????????“?????”???“?? ?”, ?????????“????????????????”, ????????????????????????????????????????????????????? ??????????????? ??????????????????????????????“??”?“???”???“may”??, ?????????“???????”?????“?????????”???????????????????“??”????“???????”? ???October???????“???????”????“??????????”?????????? “???”???“??????”??????“novenber”???????“????????”????????“????????”?“??? ????” ??????????????????

mata poisson

peacelovekenzie — Wed, 10 Sep 2008 19:08:25 -0400

From Colin Cameron http://cameron.econ.ucdavis.edu/stata/stata.html * This program does * (1) Poisson regression using Stata MATA matrix commands * (2) Poisson regression using Stata command poisson as check * using number of bids data * Invoke MATA mata // Create y vector as well as X matrix from Stat data set using st_view command st_view(y=., ., “numbids”) st_view(X=., ., (”leglrest”, “whtknght”, “bidprem”, “cons”)) Xnames = (”leglrest”, “whtknght”, “bidprem”, “cons”) // used to label output // Obtain starting values from ols regression of lny on x b = invsym(X’X)*X’(ln(y+(y:==0)*0.1)) n = rows(X) // Newton-Raphson iterations to computer Poisson MLE // Loop until parameter modification less than 10^-16 cha = 1 iter = 1 do { mu = exp(X*b) sco = (X’(y-mu))/n hes = -(X’(X:*mu))/n bold = b b = bold + invsym(-hes)*sco cha = (bold-b)’(bold-b)/(b’b) iter = iter + 1 } while (cha > 1e-16) iter bmle = b // Compute standard errors using A^-1 * B * A^-1 k = cols(X) gradmatrix = X:*(y-mu) Vmle = invsym(-n*hes)*(gradmatrix’gradmatrix)*invsym(-n*hes)*(n/(n-k)) semle = sqrt(diagonal(Vmle)) tmle = bmle:/semle pmle = 2*ttail(n-k, abs(tmle)) // Print results (bmle, semle, tmle, pmle) // Easier to read output includes regressor names (Xnames’, strofreal(bmle), strofreal(semle), strofreal(tmle), strofreal(pmle)) end

svg in R

peacelovekenzie — Sat, 13 Sep 2008 21:22:50 -0400

library(”RSVGTipsDevice”) library(”SemiPar”) data(swiss) setwd(’c:/data/’) devSVGTips(”swiss.svg”, height=7, width=10, toolTipMode=1, title=”Swiss catholic as well as agriculture plot”) par(font.lab=8,font.main=9,cex=1.5) plot(swiss$Agriculture~swiss$Catholic,ylim=c(0,130), type=”n”, xlab=”catholic”,ylab=”Agri”, main=”catholic as well as agriculture”) grid() cnames<-rownames(swiss) swiss.1<-cbind(cnames,swiss) for (i in seq(len=nrow(swiss.1))) { setSVGShapeToolTip(title=swiss.1[i,"cnames"], desc=paste(” agri=”, swiss.1[i,"Agriculture"],” ,catholic”, swiss.1[i,"Catholic"])) points(swiss.1[i,"Agriculture"], swiss.1[i,"Catholic"], pch=12, cex=2, col=ifelse(swiss.1[i,"Catholic"]>50,”blue”,”red”)) } dev.off()