search for: sstar

I try to code the ULS factor analysis descrbied in ftp://ftp.spss.com/pub/spss/statistics/spss/algorithms/ factor.pdf # see PP5-6 factanal.fit.uls <- function(cmat, factors, start=NULL, lower = 0.005, control = NULL, ...) { FAfn <- function(Psi, S, q) { Sstar <- S - diag(Psi) E <- eigen(Sstar, symmetric = TRUE, only.values = TRUE) e <- E$values[-(1:q)] e <- sum(e^2/2) e } FAgr <- function(Psi, S, q) { Sstar <- S - diag(Psi) E <- eigen(Sstar, symmetric = TRUE) L <- E$vectors[, -(1:q), drop = FALS...

...nif(1,a-da,a+da) logalpha = sum(dcauchy(y,location=fn(x,astar,bstar),scale=s,log=T) - dcauchy(y,location=fn(x,a,b),scale=s,log=T)) logu = log(runif(1,0,1)) acc = (logu < logalpha) b = acc*bstar + (1-acc)*b a = acc*astar + (1-acc)*a list(b=b,a=a,acc=acc) } samples = function(x,y,a,b,s,ds){ sstar = runif(1,s-ds,s+ds) while(sstar <= 0){ sstar = runif(1,s-ds,s+ds) } logalpha = sum( dcauchy(y,location=fn(x,a,b),scale=sstar,log=T) - dcauchy(y,location=fn(x,a,b),scale=s,log=T)) - log(sstar/s) logu = log(runif(1,0,1)) acc = (logu < logalpha) s = acc*sstar + (1-acc)*s list(s=s,acc...

...- pchisq(fit$STATISTIC, dof, lower.tail = FALSE) } fit$n.obs <- n.obs fit$call <- cl fit } factanal.fit.mle <- function(cmat, factors, start=NULL, lower = 0.005, control = NULL, ...) { FAout <- function(Psi, S, q) { sc <- diag(1/sqrt(Psi)) Sstar <- sc %*% S %*% sc E <- eigen(Sstar, symmetric = TRUE) L <- E$vectors[, 1:q, drop = FALSE] load <- L %*% diag(sqrt(pmax(E$values[1:q] - 1, 0)), q) diag(sqrt(Psi)) %*% load } FAfn <- function(Psi, S, q) { sc <- diag(1/sqrt(Psi))...

...- pchisq(fit$STATISTIC, dof, lower.tail = FALSE) } fit$n.obs <- n.obs fit$call <- cl fit } factanal.fit.mle <- function(cmat, factors, start=NULL, lower = 0.005, control = NULL, ...) { FAout <- function(Psi, S, q) { sc <- diag(1/sqrt(Psi)) Sstar <- sc %*% S %*% sc E <- eigen(Sstar, symmetric = TRUE) L <- E$vectors[, 1:q, drop = FALSE] load <- L %*% diag(sqrt(pmax(E$values[1:q] - 1, 0)), q) diag(sqrt(Psi)) %*% load } FAfn <- function(Psi, S, q) { sc <- diag(1/sqrt(Psi))...

...g the function grad() and I have programmed some simulated studies. I have testes with 72,576 different combinations of the initial values and now I got an error in only 4 situations. The error comes from the function grad() and I don't understand why. This is the output: Error in grad.default(Sstar, x = t, dist1 = T1pdist, dist2 = T2pdist, : function returns NA at 3.90625e-071.01927622979798e-087.7105723770202e-075.27096223871154e-087.28540377612885e-071.25230637383194e-076.56019362616806e-072.21329924168263e-075.59920075831737e-073.324709613353e-074.487790386647e-070.0001000016964206937....

...ed, result 58 (WERR_BAD_NET_RESP) 10623 consecutive failure(s). Last success @ Thu Mar 8 14:34:14 2018 AEDT Any thoughts on why this DC is now not replicating properly? Any thoughts on how to remedy this? -- Cheers, David Minard. Ph: 0247 360 155 Fax: 0247 360 770 ITDS - ACE - SSTaRS Western Sydney University Building Y - Penrith Campus (Kingswood) Locked bag 1797 Penrith NSW 2751 [Sometimes waking up just isn't worth the insult of the day to come.] -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

...I guess it boils down to: How would a large corporation have one source of username/passwords, and multiple separate areas where their users's attributes, policies, group membership etc are managed separately? -- Cheers, David Minard. Ph: 0247 360 155 Fax: 0247 360 770 ITDS - ACE - SSTaRS Western Sydney University Building Y - Penrith Campus (Kingswood) Locked bag 1797 Penrith NSW 2751 [Sometimes waking up just isn't worth the insult of the day to come.] -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

...t-level=2012_R2 ERROR: Forest function level can't be higher than the domain function level(s). Please raise it/them first! Am I missing something? os: Centos 6 Samba 4.7 (will go to 4.8 soon) Self compiled -- Cheers, David Minard. Ph: 0247 360 155 Fax: 0247 360 770 ITDS - ACE - SSTaRS Western Sydney University Building Y - Penrith Campus (Kingswood) Locked bag 1797 Penrith NSW 2751 [Sometimes waking up just isn't worth the insult of the day to come.] -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

...t;> >> Any thoughts on why this DC is now not replicating properly? Any >> thoughts on how to remedy this? >> >> >> -- >> >> Cheers, >> David Minard. >> Ph: 0247 360 155 >> Fax: 0247 360 770 >> >> ITDS - ACE - SSTaRS >> Western Sydney University >> Building Y - Penrith Campus (Kingswood) >> Locked bag 1797 >> Penrith NSW 2751 >> >> [Sometimes waking up just isn't worth the insult of the day to come.] >> >> -- >> This message has been scanned for viruses...

...cn=MS-SMS-Roaming-Boundary-Range. Error code = 8203. ?<11-13-2018 13:06:09> Failed to extend the Active Directory schema. Your Active Directory does not allow schema updates, please find details in "C:\ExtADSch.log". -- Cheers, David Minard. Ph: 0247 360 155 ITDS - ACE - SSTaRS Western Sydney University Building Y - Penrith Campus (Kingswood) Locked bag 1797 Penrith NSW 2751 For IT assistance contact the Service Desk - Ext 5111 or W: MyIT.westernsydney.edu.au; E: ITServiceDesk at westernsydney.edu.au; Chat: MyIT.westernsydney.edu.au, login and click on "Chat&quot...

...Any thoughts on why this DC is now not replicating > properly? Any > > thoughts on how to remedy this? > > > > > > -- > > > > Cheers, > > David Minard. > > Ph: 0247 360 155 > > Fax: 0247 360 770 > > > > ITDS - ACE - SSTaRS > > Western Sydney University > > Building Y - Penrith Campus (Kingswood) > > Locked bag 1797 > > Penrith NSW 2751 > > > > [Sometimes waking up just isn't worth the insult of the day > to come.] > > > > -- > > This message has been sc...

...success @ Thu Mar 8 14:34:14 2018 AEDT > > > Any thoughts on why this DC is now not replicating properly? Any > thoughts on how to remedy this? > > > -- > > Cheers, > David Minard. > Ph: 0247 360 155 > Fax: 0247 360 770 > > ITDS - ACE - SSTaRS > Western Sydney University > Building Y - Penrith Campus (Kingswood) > Locked bag 1797 > Penrith NSW 2751 > > [Sometimes waking up just isn't worth the insult of the day to come.] > > -- > This message has been scanned for viruses and > dangerous content by Mai...