similar to: Error in the inverse of a diagonal matrix?

dear R experts: ?I am writing my own little newey-west standard error function, with heteroskedasticity and arbitrary x period autocorrelation corrections. ?including my function in this post here may help others searching for something similar. it is working quite well, except on occasion, it complains that Error in solve.default(crossprod(x.na.omitted, x.na.omitted)) : system is

thank you, achim. I will try chol2inv. sandwich is a very nice package, but let me make some short suggestions. I am not a good econometrician, so I do not know what prewhitening is, and the vignette did not explain it. "?coeftest" did not work after I loaded the library. automatic bandwidth selection can be a good thing, but is not always. as to my own little function, I like the

Dear all, The background is that I'm trying to fix this bug in the geometry package: https://r-forge.r-project.org/tracker/index.php?func=detail&aid=1993&group_id=1149&atid=4552 Boiled down, the problem is that there exists at least one matrix X for which det(X) != 0 and for which solve(X) fails giving the error "system is computationally singular: reciprocal condition

Hi, Let?s consider aux1. I wanted to order its data considering its first column as the index > aux1[1:4,] V1 V2 V3 V4 V5 V6 1 0.08506724 150 1956 15.08 233.82 463.295 14 2 0.76205323 203 1922 15.35 218.74 463.295 14 3 0.86274890 195 1835 12.54 203.39 463.295 14 4 0.37940715 152 1567 24.73 190.85 463.295 14 > aux2 <-- aux1[order(aux1[1:4,1]) , ] > aux2

Hi all, is there any R function that loads a bmp file, or any other format, and converts its grid of pixels into a matrix or data frame? Thanks in advance. -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20000628/f3fc8b21/attachment.html

Dear all, I''ve just discovered the ''Reference Classes''. In a previous attempt ---a year ago--- to re-implement in a Object Oriented fashion the AMORE package using S4 classes I strongly felt the need of such capability. It''s great to have the Reference Classes now available. Along with the discovery of the Rcpp package, this new programming paradigm has boosted

Hello everyone, this is a very simple question but I could not find the answer. :-( Is there any function in R for calculating the rank of a matrix? (In the sense of the dimension of the vectorial space that spans) Thank you in advance! Manuel --- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi R-helpers, Can some one tell me how to train 'mynn' of this type?: mynn <- nnet(y ~ x1 + ..+ x8, data = lgist, size = 2, rang = 0.1, decay = 5e-4, maxit = 200) I assume that this nn is untrained, and to train I have to split the original data into train:test data set, do leave-one-out refitting to refine the weights (please straighten this up if I was wrong). I just don't know

Hi all, I wish to implement various landscape ecology measures through R, such as: box counting dimension; twist number statistics; contagion and lacuniarity indices; angular second moment; adjacency measures; dominance indices; etc ... Some of the measures can be applied to shape analysis and classification. Is anyone implementing any of these measures? If so I would like to contribute and

Hi, I would like to share my ZFS filesystem over the network and make it in addition fault tolerant. I am out after performance and fault tolerance, but I do not want to miss the advantages of deduplication, cloning and snapshoting offered by ZFS. I have read something about Lustre being integrated with ZFS, so that could be an option, right? Could I also use, for example, MooseFS? Thanks!

Hello, I found rather surprising the behaviour of POSIXct and POSIXlt classes when combined with min and tapply. The details can be deduced from the script below: ############# Start of the script #################### before <- Sys.time() Sys.sleep( 1 ) now1 <- now2 <- Sys.time() my.times <- c( before, now1, now2 ) class( my.times ) ## [1] "POSIXct"

Hi, Is there a way to do the following with lartc tools : I would like to limit any entering user to not use more than Xkb/mb to my website. The IPs they use are changing all the time so static IP limitation cannot be used. Is there a way doing so ? EC. _______________________________________________ LARTC mailing list / LARTC@mailman.ds9a.nl http://mailman.ds9a.nl/mailman/listinfo/lartc HOWTO:

Full_Name: Tadashi Kadowaki Version: 2.4.0 OS: Redhat Linux 9 Submission from: (NULL) (58.12.166.67) Doesn't nls function support subset? It seems not to work. And, there are no information in the online help. Has it sunk into oblivion?

Hi all, I''m in the throes of the REST conversion and I''m wondering: How important is it to retain backward compatibility? The language will clearly be consistent between the two, but it looks like it''s going to be a heckuva lot more complicated to keep compatibility for all network services (as in, for each of them, I''ll have to write a shell that

Hi everyone! I want to manage mail server resource part (like it can CGP) and with it I have one question. Is any way to limit overall max simultaneous connections to imap/pop3 server from one(each) host, except use iptables/ipfw and so on? Like a patch to dovecot or, maybe, it can be released in future versions? I know about mail_max_userip_connections in dovecot 1.1 but it can't limit

Dear R-users, i have a simple problem maybe, but i don't see the solution. i want to find the entry-wise closest element of an vector compared with another. ind1<-c(1,4,10) ind2<-c(3,5,11) for (i in length(ind2):1) { print(which.min(abs(ind1-ind2[i]))) } for ind2[3] it should be ind1[3] 10, for ind2[2] it should be ind1[2] 4 and for ind2[1] it should be ind1[1] 1. but with the

Hello all, I'm trying to use the apply function on a data frame, by applying a function that takes a one row data.frame as argument . Here's the example : myfun = function(x) paste(x$f1 , x$f2) df = data.frame(f1 = c(1,4,10),f2 = "hello") apply(df,1,myfun) ==> Does not work (I get "character(0)" ) Though : myfun(df[1,]) works, and myfun(df) works as well. So if

Dear R-experts, I have a data file df as the following: genename variable at 1hr variable at 4 hr variable at 10 hr gene1 gene2 . . . gene5000 I would like to have a graph with X-axis as the time point (1hr, 4hr, 10hr), and Y-axis as the value of the variable. So, basically, I want to do a line with the three different values at 1hr , 4hr and 10hr for all the 5000 genes. My purpose is

Dear all, I encounter some discrepancies when comparing the deviance of a weighted and unweigthed model with the AIC values. A general example (from 'nls'): DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) This is the unweighted fit, in the code of 'nls' one can see that 'nls' generates a vector