similar to: pseudo-R2 or GOF for regression trees?

Dear gentlemen, Can you explain me why the following happens (any OS I think, and even on 64 bits)? > sum(1000:1205)^2 [1] 51581223225 > sum(1000:1205)*sum(1000:1205) [1] NA Warning message: In sum(1000:1205) * sum(1000:1205) : NAs produced by integer overflow Best, Pierre -- Pierre Lafaye de Micheaux Adresse courrier: D?partement de Math?matiques et Statistique Universit? de

Hello List, I would like to orthonormalize vectors contained in a matrix X taking into account row weights (matrix diagonal D). ie, I want to obtain Z=XA with t(Z)%*%D%*%Z=diag(1) I can do the Gram-Schmidt orthogonalization with subsequent weighted regressions. I know that in the case of uniform weights, qr can do the trick. I wonder if there is a way to do it in the case of non uniform

Hello list, Sorry, these questions are not directly linked to R. If I consider an indefinte real matrix, I would like to know if the symmetry of the matrix is sufficient to say that their eigenvectors are real ? And what is the conditions to ensure that eigenvectors are real in the case of an asymmetric matrix (if some conditions exist)? Thanks in Advance, St?phane DRAY

Hello, suppose I have a list with matrices: a=list(x1=matrix(rnorm(10),5,2),x2=matrix(rnorm(10),5,2),x3=matrix(rnorm(10),5,2)) I want to compute for all combination of xi and xj (x1,x2 x1,x3 and x2,x3) a value. This value is given for the pair x1,x2 by trace(x1%*%t(x1)%*%x2%*%t(x2)) / trace(x1%*%t(x1))*trace(x2%*%t(x2)) I know that product matrices t(xi)%*%xi can be obtained by:

? Hallo, which function can i use to do (baseline) logistic regression + goodness of fit tests? so far i found: # logistic on binary data lrm combined with resid(model,'gof') # logistic on binary data glm with no gof-test # baseline logit on binary data

Hi, I'm not a statistician so sorry for possible trivial questions ... I want to perform a GoF test on sample data against several distribution (like Extreme Value, Phase Type, Pareto, ...). Since I suspect a long-tailed behaviour on data I want to use Anderson-Darling (AD) GoF test because it's well known it's more sensible to tail data. Looking at R packages the only AD test is

Hi, I'm not a statistician so sorry for possible trivial questions ... I want to perform a GoF test on sample data against several distribution (like Extreme Value, Phase Type, Pareto, ...). Since I suspect a long-tailed behaviour on data I want to use Anderson-Darling (AD) GoF test because it's well known it's more sensible to tail data. Looking at R packages the only AD test is

Hi Folks, Hopefully this question has enough R and not too much stats to be appropriate for this list. Based on,* Hosmer et al. 2013. Logistic regression for matched case-control studies. Applied Logistic Regression *(eqtn. 7.8)*, *I am assessing GOF of conditional (or matched) logistic regression models with the *standardized Pearson residuals*. The authors define ?large? as delta chi-squared

Dear R-users, I want to perform a One-Sample parametric bootstrapped Kolmogorov-Smirnov GoF test (note package "Matching" provides "ks.boot" which is a 2-sample non-parametric bootstrapped K-S version). So I wrote this code: ---[R Code] --- ks.test.bootnp <- function( x, dist, ..., alternative=c("two.sided", "less", "greater"), B = 1000 ) {

All, This may be a bit off topic so feel free to flame me ... my defence is that I am using R. I have data with case counts per family. I arrange the data in a simple table of frequency classes (e.g. how many families with 0 cases, how many with 1 case, &c.) and then GoF to Poisson and negative binomial. I treat each family as a natural sampling unit but families are of different size. I can

Hi all, I need to test gof of 3-dimensional t copula for my trivariate observed data set. So I used the command t.cop <- tCopula(c(0.785,0.283,0.613),dim=3,dispstr="un",df=6,df.fixed = TRUE) where c(0.785,0.283,0.613) is the correlation pattern of my data with 0.785 pearson correlation between variable 1-2, 0.283 correlation between 1-3 and 0.613 is the correlation between variable

Hello list, I have two vectors x and x2: x=runif(10) x2=runif(10) and one vectors with their names : my.names=c("x","x2") I would like to cbind these two vectors using their names contained in the vector my.names. I can create a string with comma ncomma=paste(my.names,collapse=",") and now, I just need a function to transform this string into a adequate

More information on the problem. A couple of days ago I posted a request for help for this error msg terre{testgb}73: smbpasswd Old SMB password: New SMB password: Retype new SMB password: /opt/samba/bin/smbpasswd: machine 127.0.0.1 rejected the session setup. Error was : ERRSRV - ERRbadpw (Bad password - name/password pair in a Tree Connect or Session Setup are invalid.). This happens only

I stumbled across this link and was wondering if this has any relation to the R project. http://www.fas.umontreal.ca/BIOL/Casgrain/en/labo/R/index.html -jennifer ----------------------------------- J. Steinbachs, Ph.D. Computational Biologist http://compbiology.org ----------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

Hello, Several machine cannot joint the cluster after a reboot. Here is a example of the error : caserta kernel: CMAN: Join request from como.dmz.lexum.pri rejected, node ID 3 already in use by ancona.dmz.lexum.pri All cluster.conf are the same (checked with md5sum) I googled this but did find nothing Any ideas? tx

I get these messages when compiling on a Sun Sparc 5 running Solaris 2.6: Compiling ipc.c ipc.c: In fonction `api_RNetServerEnum': ipc.c:1194: warning: passing arg 4 of `qsort' from incompatible pointer type Compiling nmbd_incomingrequests.c nmbd_incomingrequests.c: In function `process_node_status_request': nmbd_incomingrequests.c:380: warning: passing arg 4 of `qsort' from

No. But it has not been posted either. You got that message because you sent your message to the wrong address. You should have sent it to r-help at r-project.org You had probably sent it to r-help-request at r-project.org which would have had the effect that the server would have tried to interpret the contents of you message as commands (e.g. to unsubscribe, change your subscription

https://bugs.freedesktop.org/show_bug.cgi?id=40747 Summary: The new nouveau kernel module fails to use my monitor's native resolution Product: xorg Version: unspecified Platform: x86 (IA32) OS/Version: Linux (All) Status: NEW Severity: normal Priority: medium Component:

I don't need 3D acceleration but I'm interested in the nouveau driver because I hope it'll be able to properly set the mode of my card: currently nv is unable to make my card display 1600x1200 via DVI (it works OK via VGA but the quality is bit less stunning). So I compiled and installed it, but now I encountered a problem: it seems that (contrary to both the nv and the nvidia

Hello, I have a question that is not directly related to R ... but I try to do it in R ;-) : I would like to generate a matrix Q satisfying (for a given Z, X and W) the two following conditions: t(Q)%*%Q=Z (1) XQ=W (2) where: Q is m rows and r columns X is p rows and m columns D is p rows and r columns C is r rows and r columns with m>p,r e.g: m=6, p=2 r=3