search for: gad

...n the original SAS code in case some people are also SAS users and might find this helpful or interesting.   Thanks,   Paul      R Code:   ################################## #### Chapter 6: One-Way ANOVA #### ##################################   #### Example 6.1: One-Way ANOVA using HAM-A Scores in GAD ####   connection <- textConnection(" 101 LO 21  104 LO 18 106 LO 19  110 LO  . 112 LO 28  116 LO 22 120 LO 30  121 LO 27 124 LO 28  125 LO 19 130 LO 23  136 LO 22 137 LO 20  141 LO 19 143 LO 26  148 LO 35 152 LO  .  103 HI 16 105 HI 21  109 HI 31 111 HI 25  113 HI 23 119 HI 25  123 HI 18...

I have Samba PDC and use tdbsam backend. I use the version 3.0.11, but recently I have exchanged it for the version 3.0.23c. Now I have a number of problem. Here is one of then. I wish to move one user "gad" from Domain Users (RID 513) to Domain Guests (RID 514). Now Primary group of user "gad" is "Domain Users": pdbedit -Lv gad -d0 Unix username: gad NT username: Account Flags: [U ] User SID: S-1-5-21-3738234612-4042114426-1695072666-3024...

...(a) [1] "coef" "sigma2" "var.coef" "mask" "loglik" "aic" "arma" "residuals" "call" "series" "code" "n.cond" "model" Thanks, Gad -- Gad Abraham Department of Mathematics and Statistics The University of Melbourne Parkville 3010, Victoria, Australia email: g.abraham at ms.unimelb.edu.au web: http://www.ms.unimelb.edu.au/~gabraham

...plots on a page, however the second plot replaces the first one (although it only takes up the upper half as it should): par(mfrow=c(2,1)) a<-sin(1:100) b<-cos(1:100) lag.plot(a) lag.plot(b) What's the trick to this? I'm using R 2.2.1 (2005-12-20 r36812) on Ubuntu Linux. Thanks, Gad -- Gad Abraham Department of Mathematics and Statistics University of Melbourne Victoria 3010, Australia email: g.abraham at ms.unimelb.edu.au web: http://www.ms.unimelb.edu.au/~gabraham

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn rev 40228 language R version.string R version 2.4.1 (2...

...39;m using survreg() from the survival package for parametric survival regression (modelling inter-arrival times of patients to a waiting list as exponentially distributed, with various regressors such as queue size and season). Does anyone know which algorithm survreg() uses for this? Thanks, Gad -- Gad Abraham Department of Mathematics and Statistics The University of Melbourne Parkville 3010, Victoria, Australia email: g.abraham at ms.unimelb.edu.au web: http://www.ms.unimelb.edu.au/~gabraham

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn rev 40228 language R version.string R version 2.4.1 (2...

...label "2006" in the middle of the x axis, with no months labelled, although there's plenty of room for labels. How do I get zoo to label the months too, without having to manually work out which day was the first day of each month? I'm using zoo 1.0-3 with R 2.2.1. Thanks, Gad -- Gad Abraham Department of Mathematics and Statistics University of Melbourne Parkville 3010, Victoria, Australia email: g.abraham at ms.unimelb.edu.au web: http://www.ms.unimelb.edu.au/~gabraham

Hi, everyone. GAD package analyses complex ANOVA models with any combination of orthogonal/nested and fixed/random factors, as described by Underwood (1997). There are two restrictions: (i) data must be balanced; (ii) fixed nested factors are not allowed. Homogeneity of variances is checked using Cochran's C tes...

...(c(100, 100))) diag(M) <- 1 This doesn't work: M <- new("dsyMatrix", uplo="U", x=0, Dim=as.integer(c(100, 100))) Error in validObject(.Object) : invalid class "dsyMatrix" object: length of x slot != prod(Dim) Is there an easier way of doing this? Thanks, Gad -- Gad Abraham PhD Student, Dept. CSSE and NICTA The University of Melbourne Parkville 3010, Victoria, Australia email: gabraham at csse.unimelb.edu.au web: http://www.csse.unimelb.edu.au/~gabraham

...R 2.7.1 Ubuntu packages from CRAN. The CRAN r-base-dev package depends on refblas3-dev or atlas3-base-dev, but octave3.0-headers depends on libblas-dev. The two blas packages seem to conflict, so r-base-dev and octave3.0-headers can't be installed together. Is this a known issue? Thanks, Gad $ sudo apt-get install octave3.0-headers Reading package lists... Done Building dependency tree Reading state information... Done The following extra packages will be installed: libblas-dev libfftw3-dev libhdf5-serial-dev liblapack-dev Suggested packages: libhdf5-doc The following packages w...

On 8/3/07 22:12, "Gad Abraham" <g.abraham at ms.unimelb.edu.au> wrote: > Laura Hill wrote: >> >> >> On 7/3/07 00:15, "Gad Abraham" <g.abraham at ms.unimelb.edu.au> wrote: >> >>>> On 6 Mar 2007, at 08:54, Laura Hill wrote: >>>> >>&g...

..., x2) : 'newdata' does not have the correct number of columns > dim(x2) [1] 5 20 > dim(p$rotation) [1] 20 5 All predict.prcomp really does is x2 %*% p$rotation, which is a perfectly legal operation that should yield a 5x5 matrix: > dim(x2 %*% p$rotation) [1] 5 5 Thanks, Gad -- Gad Abraham Dept. CSSE and NICTA The University of Melbourne Parkville 3010, Victoria, Australia email: gabraham at csse.unimelb.edu.au web: http://www.csse.unimelb.edu.au/~gabraham

...svd(x)$v) [1] 5 3 > length(svd(x)$d) [1] 3 > dim(svd(x, nv=5)$v) [1] 5 5 > length(svd(x, nv=5)$d) [1] 3 > Is there a way of getting more PCs and eigenvalues than rows? Is the eigen-decomposition of the covariance matrix really that numerically bad? (i.e., eigen(cov(x)) ) Thanks, Gad -- Gad Abraham Dept. CSSE and NICTA The University of Melbourne Parkville 3010, Victoria, Australia email: gabraham at csse.unimelb.edu.au web: http://www.csse.unimelb.edu.au/~gabraham

...an.r-project.org testing/ Packages [515kB] Fetched 1B in 3s (0B/s) W: Failed to fetch http://debian.cran.r-project.org/cran2deb/debian-amd64/testing/Packages.bz2 Hash Sum mismatch My /etc/apt/sources contains: deb http://debian.cran.r-project.org/cran2deb/debian-amd64 testing/ Any idea? Thanks, Gad -- Gad Abraham PhD Student, Dept. CSSE and NICTA The University of Melbourne Parkville 3010, Victoria, Australia email: gabraham at csse.unimelb.edu.au web: http://www.csse.unimelb.edu.au/~gabraham

...s a namespace but generates a warning: > library(test2) Attaching package: 'test2' The following object(s) are masked from package:stats : predict Is this an intended behaviour? If I ignore this masking warning, do I risk unintended consequences later down the track? Thanks, Gad > sessionInfo() R version 2.8.1 (2008-12-22) i386-apple-darwin8.11.1 locale: en_AU.UTF-8/en_AU.UTF-8/C/C/en_AU.UTF-8/en_AU.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] test2_0.1 -- Gad Abraham MEng Stude...

...39; cannot exceed nrow(x) = 20 > traceback() 6: .Call("La_chol2inv", x, size, PACKAGE = "base") 5: chol2inv(fit$qr$qr) 4: ols(y ~ x1 * x2, data = x) 3: eval.with.vis(expr, envir, enclos) 2: eval.with.vis(ei, envir) 1: source("t.R", echo = TRUE) Any ideas? Thanks, Gad > sessionInfo() R version 2.6.2 (2008-02-08) x86_64-pc-linux-gnu locale: LC_CTYPE=en_AU.UTF-8;LC_NUMERIC=C;LC_TIME=en_AU.UTF-8;LC_COLLATE=en_AU.UTF-8; LC_MONETARY=en_AU.UTF-8;LC_MESSAGES=en_AU.UTF-8;LC_PAPER=en_AU.UTF-8;LC_NAME=C; LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_AU.UTF-8;LC_IDEN...

...the predictions grow without bound. Setting the boundary knots to extreme values such as (0, 10000) doesn't help, because there are no intermediate knots to constrain the spline and you get a very large cubic curve between these boundaries and the observed sizes.) Any suggestions? Thanks, Gad -- Gad Abraham Department of Mathematics and Statistics The University of Melbourne Parkville 3010, Victoria, Australia email: g.abraham at ms.unimelb.edu.au web: http://www.ms.unimelb.edu.au/~gabraham

...stic D_boot is derived by applying the function to the bootstrap sample. Then optimism O = average(D_boot - D_orig), and the bias-corrected estimate is D_app - O. Can someone explain this difference, and especially why boot doesn't evaluate the frozen model again on the original data? Thanks, Gad -- Gad Abraham Department of Mathematics and Statistics The University of Melbourne Parkville 3010, Victoria, Australia email: g.abraham at ms.unimelb.edu.au web: http://www.ms.unimelb.edu.au/~gabraham

...3CVN/P Totalling 2631 Now we won so have 1000 points Cray should have (2631/2805)*1000 = 938 points They have been given 991 points! Why! Perhaps Peter can point out my error before I raise the issue with the adjudicator. Neil, G3RIR SPONSORED LINKS Craft hobby<http://groups.yahoo.com/gads?t=ms&k=Craft+hobby&w1=Craft+hobby&w2=Hobby+and+craft+supply&w3=Ham+radio&w4=De+montfort+university&c=4&s=88&.sig=bL0By1yBwB-v0XKMui6r0w> Hobby and craft supply<http://groups.yahoo.com/gads?t=ms&k=Hobby+and+craft+supply&w1=Craft+hobby&w2=Hobby+and+c...