search for: ses

Hi for the loop section runif needs curved brackets Try IAP <-NA for (i in 1:Sample.Size){ if (DataSet$SES[i]>0) { IAP[i] <- ifelse(runif(1)>0.75, 1, 0) # High SES, higher chance to be in Treatment # } else { IAP[i] <- ifelse(runif(1)<=0.25, 1, 0) # Low SES, lower chance to be in Treatment # } } # End loop # IAP IAP zjiaqi19880219 wrote > Hi, all, > &gt...

...ter on the coefficients and significance tests, and prompted my question. The code I used is shown below, with the results at the bottom. Thanks in advance! Mike head(mort) ## data clearly contains multiple entries for some of the dataframe ids no.id.aft <- aftreg(Surv(enter, exit, event) ~ ses, data = mort) ## Inital model id.aft <- aftreg(Surv(enter, exit, event) ~ ses, data = mort, id=id) ## overlapping intervals error mort.sort <- ## ensure records ordered mort[ order(mort$id, mort$enter),] ## remove overlap for (i in 2:nrow(mort.sort)){ if (mort.sort[i,'...

I'm trying to replicate a time-series cross-sectional analysis (countries over years) with SEs clustered by country. ?The original analysis was done in Stata 10 with: xtreg [DV] [IVs] fe cluster(country). Using plm() in R (cran.r-project.org/web/packages/plm/index.html), I've replicated the coefficients. I sought to estimate country-clustered SEs with vcovHC(), and tried a variety of op...

Hi all, Im a little bit confused concerning the effect() command, effects package. I have done several glm models with family=quasipoisson: model <-glm(Y~X+Q+Z,family=quasipoisson) and then used results.effects <-effect("X",model,se=TRUE) to get the "adjusted means". I am aware about the debate concerning adjusted means, but you guys just have to trust me - it

In a function, say foo.glm for glm objects I want to use the name of the object as a label for some output, but *only* if a glm object was passed as an argument, not a call to glm() producing that object. How can I distinguish these two cases? For example, I can use the following to get the name of the argument: foo.glm <- function(object) { oname <- as.character(sys.call())[2] oname } > indep <- glm(Freq ~ mental + ses, family = poisson, data = Mental) > foo.glm(indep) [1] "indep" But in foo.glm()...

Dear All, I am trying to use the function ses from the forecast package. >From its help I have : Usage: ses(y, h = 10, level = c(80, 95), fan = FALSE, initial = c("optimal", "simple"), alpha = NULL, lambda = NULL, biasadj = FALSE, x = y, ...) My query is that if I do not mention the initial value will its...

...hey gave the same results. That is, the codes ?y~z+s(x,by=z)? and ?y~z+s(x,by=z)-1? gave the same results. The following is my result: ########################################################################### ?anova(model1,model2,test="Chisq") Analysis of Deviance Table Model 1: FLBS ~ SES + s(FAFR, by = SES) + s(byear, by = SES) + s(FAFR, byear, by = SES) Model 2: FLBS ~ SES + s(FAFR, by = SES) + s(byear, by = SES) + s(FAFR, byear, by = SES) - 1 Resid. Df Resid. Dev Df Deviance P(>|Chi|) 1 1.2076e+03 1458.4 2 1.2076e+03 1...

Hi, I was trying to get error bars in my matplot. I looked at an earlier thread, and the sample code that I made is: #------------------ library(plotrix) mat1 <- matrix(sample(1:30,10),nrow=5,ncol=2) ses <- matrix(sample(1:3,10,replace=T),nrow=5,ncol=2) vect <- seq(20,100,20) rownames(mat1) <- rownames(ses) <- vect colnames(mat1) <- colnames(ses) <- letters[1:2] matplot(mat1,pch=c('x','o'),type = "b",lwd = 2,lty = c(1,2), col = c("green",&q...

The first one, i.e. "optimal"; check help for match.arg() for the idiom. -pd > On 2 Oct 2017, at 11:48 , Ashim Kapoor <ashimkapoor at gmail.com> wrote: > > Dear All, > > I am trying to use the function ses from the forecast package. > > From its help I have : > > Usage: > > ses(y, h = 10, level = c(80, 95), fan = FALSE, initial = c("optimal", > "simple"), alpha = NULL, lambda = NULL, biasadj = FALSE, x = y, ...) > > My query is that if I do...

...following problem: I have some code below, but first I'll describe the problem. When I compute predicted logits from a logistic regression, the standard errors of the predicted logits are way off (but the predicted logits are fine). Furthermore, the model logit coefficients have appropriate SEs. As a comparison, I ran the same model without the survey design; the predicted SEs come out fine. Here is example code (first no survey design model and predictions; then survey design model and predictions): > #MODEL COEF. ESTIMATES (NO SURVEY DESIGN) > model.l.nosvy <- glm(qn58~t8l,da...

So I'm curious how SAS JBOD arrays and linux MDraid as implemented in CentOS6, and SES (SCSI/SAS Enclosure Services) backplane controllers 'get along' and how much configuration is needed to get the warning lights to work properly. scenario: whitebox server with a SAS backplane or two, daisy chained on a SAS HBA (like an LSI Logic 2008), and disks organized as several...

What is 1. par.ests <- optimfit$par 2. fisher <- hessb(negloglik, par.ests, maxvalue=maxima); 3. varcov <- solve(fisher); 4. par.ses <- sqrt(diag(varcov)); Thanks a lot, fit.GEV <- function(maxima) { sigma0 <- sqrt((6. * var(maxima))/pi) mu0 <- mean(maxima) - 0.57722 * sigma0 xi0 <- 0.1 theta <- c(xi0, mu0, sigma0) #10/5/2007: removed assign() for maxima.nl #10/5/2007: passed additional parameter to avoid us...

Hi there, as a student sociology, I'm starting to learn about SEM. The course I follow is based on LISREL, but I want to use the SEM-package on R parallel to it. Using LISREL, I found it to be very usable to be able to see the total direct and total indirect effects (standardized and unstandardized) in the output. Can I create these effects using R? I know how to calculate them

...dly slow when trying to read the complete set. Does anybody has some pointers for me for improving the import or handling such large data sets? Thanks in advance! Uli read.DET.data <- function ( f ) { counter <- 1 spk.v <- c() imp.v <- c() score.v <- c() th.v <- c() ses.v <- c() rec.v <- c() type.v <- c() fid <- file( f ,"rb") tempi <- readBin(fid , integer(), size=1, signed=FALSE) while ( length(tempi) != 0) { spk.v[ counter ] <- tempi imp.v[ counter ] <- readBin(fid, integer(), size=1, signed=FALSE) score.v[ counter ] &...

I am unsure what is being returned, and what is supposed to be returned, when using 'predict' with "type='expected'" for an aftreg survival model. The code below first generates a weibull model, then uses predict to create a vector of the linear predictors, then attempts to create the 'expected' vector, which is empty. The final two steps in the code generate a lognormal model with the same data, and the same empty 'expected' vector. My expectation had been that 'expected'...

On 04/26/2017 04:22 AM, Gordon Messmer wrote: > On 04/25/2017 03:25 PM, Robert Moskowitz wrote: >> This made the same content as before that caused problems: > > I still don't understand, exactly. Are you seeing *new* problems > after installing a policy? What are the problems? > >> #!!!! The file '/var/lib/mysql/mysql.sock' is mislabeled on your system.

On Tue, Mar 14, 2017 at 02:46:19PM -0400, Ron Wheeler wrote: > https://docs.fedoraproject.org/en-US/Fedora/11/html/Security-Enhanced_Linux/sect-Security-Enhanced_Linux-Working_with_SELinux-Enabling_and_Disabling_SELinux.html > > If disabling Selinux solves your problem, then your problem may be related > to Selinux. > If it does not change yout problem, you may want to look

...[guest] context=google-in disallow=all allow=ulaw allow=g729 connection=jp_jabber jabber.conf [general] debug=yes ;autoprune=no autoregister=yes [jb_jabber] type=client serverhost=talk.google.com username=XXXXXXXXX at gmail.com/Talk secret=XXXXXXX port=5222 usetls=yes usesasl=yes ;status=Available statusmessage="Connected via Asterisk" ;timeout=100 ;keepalive=yes Extensions.conf [google-in] exten => s,1,NoOp(Call from GTalk) exten => s,n,Set(CallerID(Name)="From GoogleTalk") exten => s,n,Dial(SIP/1000) jabber show co...

The recently-left programmer did *something*, and he didn't know what, and the guy who picked it up is working with me to find out why /var/log/messages is getting flooded with Oct 26 11:01:06 <servername> kernel: type=1105 audit(1477494066.569:642430): pid=108551 uid=0 auid=4294967295 ses=4294967295 subj=system_u:system_r:unconfined_service_t:s0 msg='op=PAM:session_open grantors=pam_keyinit,pam_keyinit,pam_limits,pam_systemd,pam_unix,pam_krb5,pam_xauth acct="<user>" exe="/usr/bin/su" hostname=? addr=? terminal=? res=success' Oct 26 11:01:06 <serv...

...:s0 tcontext=unconfined_u:object_r:var_lib_t:s0 tclass=file > type=SYSCALL msg=audit(1489588773.253:1171): arch=c000003e syscall=2 success=no exit=-13 a0=aa5080 a1=80000 a2=1a4 a3=aa5080 items=0 ppid=1485 pid=3838 auid=4294967295 uid=0 gid=0 euid=0 suid=0 fsuid=0 egid=0 sgid=0 fsgid=0 tty=(none) ses=4294967295 comm="asterisk" exe="/usr/sbin/asterisk" subj=system_u:system_r:asterisk_t:s0 key=(null) > type=AVC msg=audit(1489588777.432:1172): avc: denied { getattr } for pid=3844 comm="asterisk" path="/var/lib/asterisk/astdb.sqlite3" dev="dm-0&quo...