search for: outx

...I would like to calculate the Goodman-Kruskal gamma for the predicted classes obtained from an ordinal regression model using lrm in the Design package. I couldn't find a way to get gamma for predicted values in Design so have found previous positings suggesting to use : Rcorr.cens(x, S outx = TRUE) in the Hmisc package My question is, will this work for predicted vs observed factors? I.e. x = predicted class and S = observed class? Or is there a better way to obtain this? I used the maximum individual probability for each observation to determine the predicted class. Any he...

...normalmixEM (package:mixtools) to fit them in finite mixture narmal distributions,but get some mistakes.I don't know which steps make the error. I have used the following functions before >x<-read.table("data") >log.x<-log10(x$V1) >log.x<-sort(log.x) >outx = normalmixEM(log.x,arbmean = TRUE, arbvar = TRUE, k=2,epsilon = 1e-03,maxit = 10000, verb = FALSE) The output is given as: Need new starting values due to singularity... Need new starting values due to singularity... Need new starting values due to singularity... error in normalmixEM...

...rical variable* where 1=worst group and 5=best group - I presume that the c index would be calculated similarly but this time there would be many ties in the predictor (as regards predicted survival probability) - hence if I wanted to count all ties in such a case I would keep the default argument outx=FALSE? Does anyone have a clear reference which gives the formula used to generate the concordance index (with worked examples)? Many thanks for your help on these interpretations Kind Regards, Kim

...inal scale; ranked) for each plant species and relev?. I found a lot of r-packages dealing with cluster analysis, but none of them is able to calculate a distance measure for ranked data. Podani recommends the use of Goodman and Kruskals' Gamma for the distance. I found the function rcorr.cens (outx=true) of the Hmisc package which should do it. What I don't understand is how to define the input vectors x, y with my vegetation dataset. The other thing how I can use the output of rcorr.cens for a distance measure in the cluster analysis (e.g. in vegan or amap). Any help would be greatly app...

.... Now I am using function, rcorr.cens http://lib.stat.cmu.edu/S/Harrell/help/Hmisc/html/rcorr.cens.html This function involves use of Surv. Now here is a section of my syntax: > time<-data$ovsrecod > x1<-data$RMY.GROUPS > death<-data$death > rcorr.cens(x1,Surv(time,death),outx=FALSE) (2) As you can see, I have entered Surv(time,death)...this works (and complies with the example given in R for rcorr.cens) and all seems to be well...however, bearing in mind that in (1) we have: Surv(time,time2,event, type=, origin = 0) ...how does R know that 'death' in *my*...

Dear R-help, This is an example in the {Hmisc} manual under rcorr.cens function: > set.seed(1) > x <- round(rnorm(200)) > y <- rnorm(200) > round(rcorr.cens(x, y, outx=F),4) C Index Dxy S.D. n missing uncensored Relevant Pairs Concordant Uncertain 0.4831 -0.0338 0.0462 200.0000 0.0000 200.0000 39800.0000 19228.0000 0.0000 That S.D. confuses me...

...ot;setting owner of %s to %d:root" file uid; + Unix.chown file uid 0 (* Get the local user that libvirt uses to run qemu when we are * running as root. This is returned as an optional string @@ -205,8 +199,8 @@ let error_if_no_ssh_agent () = (* Create the directory containing inX and outX sockets. *) let create_v2v_directory () = let d = Mkdtemp.temp_dir "v2v." in - chown_for_libvirt_rhbz_1045069 d; On_exit.rm_rf d; + chown_for_libvirt_rhbz_1045069 d; d (* Wait for a file to appear until a timeout. *)

Dear all, I am interested to evaluate reclassification using net reclassification improvement and Integrated Discrimination Index IDI after survival analysis (Cox proportional hazards using stcox). I search a R package or a R code that specifically addresses the categorical NRI for time-to-event data in the presence of censored observation and, if possible, at different follow-up time points. I

...'t this need parens around the 'then' clause? > (* Get the local user that libvirt uses to run qemu when we are > * running as root. This is returned as an optional string > @@ -205,8 +199,8 @@ let error_if_no_ssh_agent () = > (* Create the directory containing inX and outX sockets. *) > let create_v2v_directory () = > let d = Mkdtemp.temp_dir "v2v." in > - chown_for_libvirt_rhbz_1045069 d; > On_exit.rm_rf d; > + chown_for_libvirt_rhbz_1045069 d; > d Yes, as you explain in the commit message this hunk is needed (and dependent)...

The Goodman/Kruskal gamma is a nice descriptive rank-order correlation statistic, often used in psychology. It is nice because it is easy to understand. It takes all pairs of values of each variable and asks whether they are congruent (S+ is the number in the same order for both variables) or discordant (S-, opposite ranking). The statistic is (S+ - S-)/(S+ + S-). It is like tau except for the

...cing "uid" out of scope). > >> (* Get the local user that libvirt uses to run qemu when we are >> * running as root. This is returned as an optional string >> @@ -205,8 +199,8 @@ let error_if_no_ssh_agent () = >> (* Create the directory containing inX and outX sockets. *) >> let create_v2v_directory () = >> let d = Mkdtemp.temp_dir "v2v." in >> - chown_for_libvirt_rhbz_1045069 d; >> On_exit.rm_rf d; >> + chown_for_libvirt_rhbz_1045069 d; >> d > > Yes, as you explain in the commit message th...

I have 10 SIP trunks that I'd really like to round-robin load balance. Currently I have a macro that switches between available lines, but there really must be a function in Asterisk to do this on its own. So my question is just that, are there any easy ways for Asterisk to either balance between SIP trunks or even just a built in function to find the next available SIP trunk. I think using

Dear list, I have a problem on calculating the standard error of Goodman-Kurskal's gamma using delta method. I exactly follow the method and forumla described in Problem 3.27 of Alan Agresti's Categorical Data Analysis (2nd edition). The data I used is also from the job satisfaction vs. income example from that book. job <- matrix(c(1, 3, 10, 6, 2, 3, 10, 7, 1, 6, 14, 12, 0, 1, 9,

...ope). > > > > >> (* Get the local user that libvirt uses to run qemu when we are > >> * running as root. This is returned as an optional string > >> @@ -205,8 +199,8 @@ let error_if_no_ssh_agent () = > >> (* Create the directory containing inX and outX sockets. *) > >> let create_v2v_directory () = > >> let d = Mkdtemp.temp_dir "v2v." in > >> - chown_for_libvirt_rhbz_1045069 d; > >> On_exit.rm_rf d; > >> + chown_for_libvirt_rhbz_1045069 d; > >> d > > > > Yes,...

...t;>> >>>> (* Get the local user that libvirt uses to run qemu when we are >>>> * running as root. This is returned as an optional string >>>> @@ -205,8 +199,8 @@ let error_if_no_ssh_agent () = >>>> (* Create the directory containing inX and outX sockets. *) >>>> let create_v2v_directory () = >>>> let d = Mkdtemp.temp_dir "v2v." in >>>> - chown_for_libvirt_rhbz_1045069 d; >>>> On_exit.rm_rf d; >>>> + chown_for_libvirt_rhbz_1045069 d; >>>> d >&gt...

Hi, Here is my patch queue I accumulated over quite a long time. Patches 1-6 are bugfixes, and rest is mostly RFC. Comments are welcome. Best regards, Maxim Levitsky

Bugzilla: https://bugzilla.redhat.com/show_bug.cgi?id=2182024 v1 was here: <https://listman.redhat.com/archives/libguestfs/2023-June/031910.html>. Make any "chown_for_libvirt_rhbz_1045069" failure a hard one; that way the user gets to see more direct and more uniform error messages (namely that we couldn't connect to libvirt). Document those connection problems, and a simple

...t nouveau_channel *evo = dev_priv->evo; int ret; - NV_DEBUG(dev, "\n"); + NV_DEBUG_KMS(dev, "\n"); ret = RING_SPACE(evo, 2 + (update ? 2 : 0)); if (ret) { @@ -193,7 +193,7 @@ nv50_crtc_set_scale(struct nouveau_crtc *nv_crtc, int scaling_mode, bool update) uint32_t outX, outY, horiz, vert; int ret; - NV_DEBUG(dev, "\n"); + NV_DEBUG_KMS(dev, "\n"); switch (scaling_mode) { case DRM_MODE_SCALE_NONE: @@ -301,7 +301,7 @@ nv50_crtc_destroy(struct drm_crtc *crtc) struct drm_device *dev = crtc->dev; struct nouveau_crtc *nv_crtc = nouve...