search for: 95110

...57.86592 9.755003 4.343534 4.280714 93610 93 1 53.36475 39.47247 4.381618 1.622119 1.159044 94110 94 0 23.47514 58.23936 10.789339 3.690415 3.805741 94610 94 1 53.34524 39.32675 4.602774 1.520247 1.204988 95110 95 0 23.76256 58.18757 9.583178 4.186825 4.279866 95610 95 1 53.14037 39.55376 4.313944 1.632805 1.359124 and want to get it in the form: age sex cat Int.Prev.Est 93 0 1 23.75482 93 0 2 57.86592 93 0 3 9.755003 93 0...

Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries of the matrix are: (cos\theta) exp(j\alpha); -(sin\theta)exp(j(\alpha-\Omega)); (sin\theta)exp(j(\beta+\Omega)); (cos\theta) exp(j\beta); where "j" is for complex. However, since for any two unitary matrices, their product should also be a unitary matrix. When I

Dear all, I have a population with three age-classes, at time t=0 the population is: n.zero <- c(1,0,0) I have a transition matrix A which denotes "fertility" and "survival": A <- matrix(c(0,1,5, 0.3,0,0, 0,0.5,0), ncol=3, byrow=TRUE) To obtain the population at t=1, I calculate: A %*% n.zero To obtain the population t=2, I calculate: A %*% (A %*% n.zero) ... and so

...I searched the bugs and new and didn't see anything regarding this. It appears to me to be a bug, can anyone confirm? -- ------------------------------------- Scott Machtmes, Sr. Systems Administrator Engineering Lab Services Infrastructure Lead Brocade 1745 Technology Drive San Jose, Ca. 95110 408-333-7203 scottm@brocade.com

Hi, I am a user of R computer language.Several days ago,I download R-2.1.0.From the help manual,I saw there is a function "mle" to find the maximum likelihood estimators.But when I try to use this function,there is no this function.Could you help me to solve this problem? I am looking forward to your answer. Thank you. Jun Liu

Hi, The question is as follows: has anyone coded panel data unit root tests with R? Even the "first generation" tests (see e.g. Levin & Lin 1993; Pesaran, & Smith & Im 1996; Maddala & Wu 1999) would be sufficient for my needs. To my understanding, these are rather easy to code, but as I have taken just my first steps in coding with R, existing code would save me

Hey, I am estiamting a partially linear model y=X\beta+f(\theta) where the f(\theta) is estiamted using wavelets. Has anyone heard of methods to test if the betas are significant or to address model fit? Thanks for any thoughts or comments. Elizabeth Lawson __________________________________________________ [[alternative HTML version deleted]]

Dear R-users, Is anyone aware of any function/package for generating a random vector from a joint distribution defined by an independent graph? Or I have to work it out myself? Thanks. Jinfang ------------------------------ Jinfang Wang, Associate Professor Chiba University, Japan

I often use simple nested random-effect models for interlaboratory data. The variance components are important things to know. Is there an R function or package that gets variance components from lm objects? Or can someone point me to a method of doing so?

...the .Internal qbinom code. However, that's beyond my current R skill level. Thanks for all your efforts to make R what it is today. Best Wishes, spencer graves -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at pdf.com www.pdf.com <http://www.pdf.com> Tel: 408-938-4420 Fax: 408-280-7915

Hello and sorry to disturb. I'm trying to plot the confidence intervals for the fixed effects of a lme. I want to obtain graphically, if it is possible, a bar with Estimate, upper and lower CI for each level of the factors. I know how to do for a lm model but for a lme one, I tried with plot(intervals(...)) and plot(ci(...)) from the gmodels package but it doesn't work well. Thanks

...ckage 'Matrix' required by 'lme4' could not be found What do you suggest? Spencer Graves p.s. I get the same result using several different (US) mirrors. -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at pdf.com www.pdf.com <http://www.pdf.com> Tel: 408-938-4420 Fax: 408-280-7915

...) } Qp.mean } set.seed(1) plot(simCus5(Qp0=4.5)) ################################################# Thanks for your time in reading this. Best Wishes, Spencer Graves Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at pdf.com www.pdf.com <http://www.pdf.com> Tel: 408-938-4420 Fax: 408-280-7915

Dear R users; Ive got two questions concerning nlme library 3.1-65 (running on R 2.2.0 / Win XP Pro). The first one is related to augPred function. Ive been working with a nonlinear mixed model with no problems so far. However, when the parameters of the model are specified in terms of some other covariates, say treatment (i.e. phi1~trt1+trt2, etc) the augPred function give me the following

Hi all, I encountered this problem with parameterization in R: I have two factors in a regression. how about if I want to set constraint so that for each factor, the sum of their coefficients equals to zero(instead of choosing a reference category)? for example, I have factor(variable) A(with three categories) and factor(variable) B(with 4 categories), and I want to parameterize so that the sum

...___________________ > On the road to retirement? Check out MSN Life Events for advice on how > to get there! http://lifeevents.msn.com/category.aspx?cid=Retirement > -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at pdf.com www.pdf.com <http://www.pdf.com> Tel: 408-938-4420 Fax: 408-280-7915

Dear Friends, I am seeking for any help on an error message in lme functions. I use mixed model to analyze a data with compound symmetric correlation structure. But I get an error message: "Error in corMatrix.corCompSymm(object) : NA/NaN/Inf in foreign function call (arg 1)". If I change the correlation structure to corAR1, then no error. I have no clue how to solve this problem.

All, I've coded a function and it works manually if I copy it line by line into R. However, when I try to "load" (copy and paste) the entire function into R, I get the following error after the listed line of code: + N.j.list = lapply(rej.hyp, length) Error: symbol print-name too long Does anyone you know what this error means? Strangely, when I copy the same line verbatim

...ps://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide! > http://www.R-project.org/posting-guide.html > > -- > Spencer > Graves, PhD > Senior Development Engineer > PDF Solutions, Inc. > 333 West San Carlos Street Suite 700 > San Jose, CA 95110, USA > > spencer.graves at pdf.com > www.pdf.com <http://www.pdf.com> > Tel: 408-938-4420 > Fax: 408-280-7915 > > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE...

Hello, I want to solve some two dimensional equation system with R. Some systems are not solvable analytically. Here is an example: (I) 1/n*sum{from_i=1_to_n}(Xi) = ln lambda + digamma(c) (II) mean(X) = x / lambda I want to find lambda and c, which R-function could do that task? Carsten [[alternative HTML version deleted]]