search for: loglog

I am trying to figure out how to apply a loglog link to a binomial model (dichotomous response variable with far more zeros than ones). I am aware that there are several relevant posts on this list, but I am afraid I need a little more help. The two suggested approaches seem to be: 1) modify the make.link function in GLM, or 2) use the...

I would like to do an R glm() with family = binomial(link="loglog") Right now, the cloglog link exists, which is nice when the data have a heavy tail to the left. I have the opposite case and the loglog link is what I need. Can someone suggest how to add the loglog link onto glm()? It would be lovely to have it there by default, and it certainly makes sense...

Dear all, A looong time ago, Witold Eryk Wolski asked here why there wasn't a log="xy" parameter to the hist() function <URL:http://www.R-project.org/nocvs/mail/r-help/2001/0267.html>, and Prof. Ripley responded that a loglog histogram does not make much sense, and that one should use a better density estimate if one seeks to plot log density. I understand the point and I hope you will forgive my ignorance: What alternatives do I have in R? (In astronomy, which is full of power laws, loglog histograms are not infrequen...

Hi, I wish to write a new link function for a GLM. R's glm routine does not supply the "loglog" link. I modified the make.link function adding the code: }, loglog = { linkfun <- function(mu) -log(-log(mu)) linkinv <- function(eta) exp(-exp(-eta)) mu.eta <- function(eta) exp(-exp(-eta)-eta) valideta <- function(eta) all(eta != 0)...

On Tue, 2019-05-21 at 20:30 +0100, Sérgio Basto via samba wrote: > On Tue, 2019-05-21 at 20:53 +0200, Maurizio Caloro via samba wrote: > > Hello Dear Members > > > > Please i have here the following error if i do the samba_dnsupdate > > --verbose > > > > the problem are that i can't join any new machine to me Samba AD > > machine. > >

Help needed with lm function: Dear R's, Could anyone tell me how to replace the link function (probit logit, loglog etc.) in lm with an abitrary user-defined function? The task is to perform ML Estimation of betas for a dichotome target variable. Maybe there is already a package for this (I did not find one). Any hints or a code excerpt would be welcome! Thank you -Jeff jeff.pr2 (at) added-insight (dot) net

...about 1000 columns. The 300 factors all pass a preliminary test of association with the outcome. Solving this with cph from Design takes about 3 hours. I have created the fit with x=T, y=T to save the model data. I want to validate the PH assumption by calling survplot(fit, gender=NA, logt=TRUE, loglog=TRUE) for many of the factors (here gender is one column name). Just creating this one plot takes 40m. I'd be happy to sample from the fitted model to create these tests, or figure out another way to check assumptions in the model. Has anyone done something similar, or have other suggestions...

Hi to every body. I would like assistance on how to implement the log-log link function for binary response. Is there any package that implements it? Many thanks Endy [[alternative HTML version deleted]]

...ex treat; DATALINES; 0 0 1 0 0 1 1 1 ; run; /* Survival estimates for the above combinations */ proc phreg data = pudat2; model withtime*wcens(0) = sex treat /ties = efron; baseline out = surv2 survival = survival lower = slower upper = supper covariates = test2 /method = ch nomean cltype=loglog; run; /* Survival estimates at 1 year */ proc print data = surv2 noobs; where withtime = 364; run; I now would like to do the same thing but in a competing risks setting (cumulative incidence rather than Cox model version (crr from cmprsk in R) however there is no pre-set routine for this in SAS....

Dear R users, I know one way to see the code for a hidden function, say function_x, is using default.function_x (e.g. summary.default). But how can I see the code for imported packages that have no namespace (in this case Design)? Many Thanks Eleni

Hi list,I want to make winbind kerberos ticket refresh work but I couldn't do it with configuration below: ------ smb.conf ------security = ADS workgroup = MYDOMAINrealm = MYDOMAIN.ORG log file = /var/log/samba/%m.loglog level = 6enable core files = no idmap config * : backend = tdbidmap config * : range = 3000-7999idmap config MYDOMAIN : backend = rid idmap config MYDOMAIN : range = 10000-999999 dedicated keytab file = /etc/krb5.keytabkerberos method = secrets and keytab template shell = /bin/bash template homedir...

Dear All, I am trying to produce log-log plots in R and I was wondering if any of you have a 'template' for generating these with 'nice' labels and log-log grids? I know I can set up axes individually and use the intervals I want, however, I will be producing a large number of these plots and would not like to do this manually for each of them + I am very new to R and at the

...assuming symmetry or equidistant thresholds - random effects via the Laplace approximation and adaptive Gauss-Hermite quadrature in the location-part of the model. - a range of standard link functions - flexible link functions where an extra link function-parameter bridges the log-log, probit and c-loglog links (log-gamma), and cloglog and logistic links (Aranda-Ordaz) - a suite of optimizers including an efficient Newton scheme. - works for binomial observations (a special case of ordinal data). - a suite of methods including anova, addterm, dropterm, profile, confint, plot.profile, predict, in add...

Dear experts, I would like to know how to plot the log-minus-log plot for survival analysis (to check the proportional assumption) in R. Using the AML example. fit <- survfit(Surv(time, status) ~ x, data=aml) length(fit$surv) #20 as the length of fit$surv is shorter than aml$x and aml$time. I don't know how to plot. Thank you. Regards, CH -- CH Chan Research Assistant - KWH

...assuming symmetry or equidistant thresholds - random effects via the Laplace approximation and adaptive Gauss-Hermite quadrature in the location-part of the model. - a range of standard link functions - flexible link functions where an extra link function-parameter bridges the log-log, probit and c-loglog links (log-gamma), and cloglog and logistic links (Aranda-Ordaz) - a suite of optimizers including an efficient Newton scheme. - works for binomial observations (a special case of ordinal data). - a suite of methods including anova, addterm, dropterm, profile, confint, plot.profile, predict, in add...

...0.6345925 1.4218795 2.3178357 1.5454827 OUTPUT # (B) In vglm.fitter(x = x, y = y, w = w, offset = offset, Xm2 = Xm2, : convergence not obtained in 155 iterations. OUTPUT # (C) > coef(fit3, matrix=TRUE) log(difference) log(scale) loglog(shape) (Intercept) -0.7954072 -0.1283776 -0.03175329 OUTPUT # (D) > Coef(fit3)) difference scale shape 0.4513974 0.8795212 2.6346374 (IS IT THAT THESE GIVE THE SCALE AND SHAPE PARAMETER ESTIMATES OF FRECHET.?) OUTPUT # (E) >...

...ow does "winbind refresh tickets" work? > > Hi list,I want to make winbind kerberos ticket refresh work > but I couldn't do it with configuration below: > ------ smb.conf ------security = ADS > workgroup = MYDOMAINrealm = MYDOMAIN.ORG > log file = /var/log/samba/%m.loglog level = 6enable core files = no > idmap config * : backend = tdbidmap config * : range = > 3000-7999idmap config MYDOMAIN : backend = rid > idmap config MYDOMAIN : range = 10000-999999 > dedicated keytab file = /etc/krb5.keytabkerberos method = > secrets and keytab > template sh...

Hello All, ? Am replicating in R an analysis I did earlier using SAS. See this as a test of whether I'm ready to start using R in my day-to-day work. ? Just finished replicating a Kaplan Meier analysis. Everything seems to work out fine except for one thing. The 95% CI around my estimate for the median is substantially larger in R than in SAS. For example, in SAS I have a median of 3.29 with a

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

...ease find the code and data following. Problem appears in lines 37 - 42 of the code. I am using R.2.13.0 on WinXP. Regards, Fayez PowerLawGraphsAblineQn.R? library('reshape') library(lattice) library(igraph) # to use power.law.fit function library(latticeExtra) # to use panel.lmlineq in loglog xyplot File2Open = "C:\\Documents and Settings\\All Users\\Documents\\Desktop\\Fayez\\RPractice\\PowerLawGraphsAblineQnData.txt" DataTable = read.table(File2Open, header = TRUE, sep = "\t") md <- melt(DataTable, id.vars = c('Range', 'TheValue', 'TheNo'...