similar to: Great difference for piecewise linear function between R and SAS

Dear friends, I used the B-spline (degree=1) method to fit the piecewise linear function and the results are listed below. m.glm<-glm(mark~x+poly(elevation,2)+bs(distance,degree=1,knots=c(16.13,24)) +bs(y,degree=1,knots=c(-0.4357,-0.3202 )),family=binomial(logit),data=point) summary(m.glm) Coefficients: Estimate Std.

Hi! Something funny happens to me in lnme. I don't know if it's a bug or not. Here's a short sketch of the problem (you could skip the introduction and just go to the models): I'm counting the number of spider webs ( = response variable) along two rivers, with a two-level treatment: free-flowing and dammed. There are three field sites (random factor) along each river, and

I have a program which needs to compute squared Euclidean distance between two vectors million of times, which the Rprof shows is the bottleneck. I wondered if there is any faster way than my own simple function distance2 = function(x1, x2) { temp = x1-x2 sum(temp*temp) } I have searched the R-help archives and can not find anything except when the arguments are matrices. Thanks for any

Hi all, To follow up on an older thread, it was suggested that the following would produce confidence intervals for the estimated BLUPs from a linear mixed effect model: OrthoFem<-Orthodont[Orthodont$Sex=="Female",] fm1OrthF. <- lmer(distance~age+(age|Subject), data=OrthoFem) fm1.s <- coef(fm1OrthF.)$Subject fm1.s.var <- fm1OrthF. at bVar$Subject fm1.s0.s <-

> myfit1 <- gls(col1 ~ col2+col3+col4+col5+col6-1, data=data2, corr=corAR1( 0.3202), method='ML') Error: cannot allocate vector of size 199712 Kb if I get rid of the "corr=corAR1(0.3202)" option, it works okay... can anybody help me? thanks a lot! [[alternative HTML version deleted]]

Dear R team, I don't get what I think I should get when using anova.glm() with family="gaussian" -- please ignore this and forgive me if this turns out to be another example of a fundamental misunderstanding on my part (a highly likely event!) For example: S <- as.factor(rep(c(rep("m",2),rep("f",2)),2)) A <-

Hi, I want to fit a kernel density estimator by bkde of library KernSmooth. I need only the density value at the point 0. I do not understand the following behaviour: > q <- rnorm(100) > bkq <- bkde(q, bandwidth=0.11, gridsize=1, range.x=c(0,0)) Error in 0:L : NA/NaN argument > bkq <- bkde(q, bandwidth=0.11, gridsize=1, range.x=c(-1,1)) > bkq $x [1] -1 $y [1] NA > bkq

The freenum.org project wants to use your trunks! The freenum.org project is an ENUM parallel tree, which has as an eventual goal the distribution of ENUM numbering in nations or areas which due to political or other issues are not able to get secure, inexpensive, or functional ENUM capabilities. As a preliminary round, we're putting toll-free gateways into the system for various nations.

Dear list, I have been checking my package ('analogue') using R2.13-0-alpha (details of exact svn version appended below) and the R CMD check procedure is generating an error rebuilding a vignette in the package, which raises a NOTE in the check. The log printed to screen during check shows: .... * checking re-building of vignettes ... NOTE Error in re-building vignettes: ... Error in

Dear R-helpers, In my dataset I have two continuous variable (A and B) and one factor. I'm investigating the regression between the two variables usign the command lm(A ~ B, ...) but now I want to know the regression coefficient (r2) of A vs. B for every factors. I know that I can obtain this information with excel, but the factor have 68 levels...maybe [r] have a useful command. Thanks,

Dear List People Someone here could please send me some advice on changin' LDAP passwords? It seems that everything is working nice, except by dovecot: I can send emails using my new LDAP password trought postfix smtp authentication, I can login to my Intranet website and to the Mod-LDAP protected pages on my company. I just can't POP my mail messages, because Dovecot says

Dear R users, I came up with some simple functions to give me the standard betas and tolerance values from a predefined lm() model. I have been trying to insert the results from these functions into the coefficients matrix in a modified summary.lm function that I'm calling summary2 (I'd never edit the summary.lm function directly!!). I managed to get the results inserted into the

Dear list,I am quite new to the world of biostatistics and I encounter some issues in the precise understanding of the coxph function of the survival package.I have a set of survival data (patient who had (or died from) a breast cancer) I'd like to see which are the variables that might cause dead or not.When trying> summary(coxph(Surv(Time_to_distant_recurrence_yrs, !Distant_recurrence)~

Hi, is it, somehow, possible to distinguish between multiple simultaneous logins in MailLog plugin? For example, which login deleted which message. Is that p[lugin able to log e.g. session ID? azur

Cit?t Aki Tuomi <aki.tuomi at open-xchange.com>: > On 7.11.2018 16.13, azurit at pobox.sk wrote: >> Hi, >> >> is it, somehow, possible to distinguish between multiple simultaneous >> logins in MailLog plugin? For example, which login deleted which >> message. Is that p[lugin able to log e.g. session ID? >> >> azur >> >> >

On 23. Feb 2023, at 16.13, Leander Beernaert <leander.beernaert at proton.ch> wrote: > > Hey, > > We recently announced Gluon (https://github.com/ProtonMail/gluon/) our IMAP server library we are using in Proton Bridge(https://github.com/ProtonMail/proton-bridge). We would love to have it have it listed in the IMAP Server Compliancy Status wiki page

Hey Timo, Thanks for the quick turnaround, once we have the test results I'll contact you again. Should I also include instructions on how to run the a self contained server with a dummy backend so you can independently verify our results? Leander Beernaert Proton AG ------- Original Message ------- On Thursday, February 23rd, 2023 at 8:59 PM, Timo Sirainen <timo at sirainen.com>

Hello, I also have issues using Dovecot Backup. I am trying to Backup (or possibly Sync - one way) a single user IMAP account from a remote server, unknown type, to my own server. For that purpose, I defined the following parameters : imapc_host = <host name> imapc_user = <user name> imapc_password = <password> imapc_features = rfc822.size (<user name> has format:

On Sat, Jun 01, 2013 at 06:45:48AM +0200, Duncan Sands wrote: > > These results are very disappointing, I was hoping to see a big improvement > somewhere instead of no real improvement anywhere (except for gas_dyn) or a > regression (eg: mdbx). I think LLVM now has a reasonable array of fast-math > optimizations. I will try to find time to poke at gas_dyn and induct: since >

Something I can't quite place finger on here. Altogether too much Mafia, in the bulk email business generally, and I know Switzerland borders on Italy ... This sounds, (albeit vaguely,) altogether too much like the thieves I seem to have fallen amongst lately. Two stolen trucks, three stolen laptops, another one wrecked, three or four stolen cell phones, passwords GPG keys, city hall hookers