similar to: round, unique and factor

Hi I'm sorry about the triviality of my problem. I have a vector (v) of three columns (logA, logB, id). I want to compute (and plot) the correlation between logA and logB for different thresholds of id (e.g. >30, etc). So I tried: for(i in 1:100){ points(cor(v$logA[v$id>i], v$logB[v$id>i], use="complete.obs"), i)) } (i created a plot object already) but it comes with

Dear all, Suppose I have a linear mixed-effects model (from the package nlme) with nested random effects (see below); how would I present the results from the random effects part in a publication? Specifically, I?d like to know: (1) What is the total variance of the random effects at each level? (2) How can I test the significance of the variance components? (3) Is there something like an

Is there an R-language interface to the R API C-language functions Rf_logspace_add() and Rf_logspace_sub()? I don't see one but I may not looking under the right name. Various packages have functions which do that same sort of thing (log(exp(x)+exp(y)) and log(exp(x)-exp(y)) without unnecessary floating point errors). They have names like matrixStats::logSumExp(lx, na.rm=FALSE, ...)

Dear R users, according the help(log), the function log2(x) should give the natural logarithm of x. I expect in case of x=2 to to get 0.6931, however, R gives me 1 as a result. Similar, logb(2,2) gives 1 again. I'm wondering if I have missed something ? Yours Frank -- Frank Mattes, MD e-mail: f.mattes at ucl.ac.uk Department of Virology fax 0044(0)207 8302854 Royal Free Hospital

I stumbled onto this working on an update to coxph. The last 6 lines below are the question, the rest create a test data set. tmt585% R R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-unknown-linux-gnu (64-bit) # Lines of code from survival/tests/singtest.R > library(survival) Loading required package: splines

Hi I discovered a similar problem early this week and a colleage on irc was told it was a copy and paste error. in man syslogd, the following can be found: Sucker rod def. -- 3/4, 7/8 or 1in. hardened steel rod, male threaded on each end. Primary use in the oil industry in Western North Dakota and other locations to pump 'suck' oil from oil wells. Sec- ondary

Hi, I have just taken over the administration of our Samba Fileserver. Unfortantly, my colleage has not done his homework and sadly forgot keep the Samba version up-to-date. So, the first step for me would be to update from our current version 3.0.23c (SuSe rpms) to the most recent. Now I worry about some changes that could cause problems with my current configuration file. As I can not shutdown

Dear R-users, I am using R 2.13.0 and rms 3.3-0 , but can not reproduce figure 7.8 of the handouts *Regression Modeling Strategies* ( http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/course2.pdf) by the following code. Could any one help me figure out how to solve this? setwd('C:/Rharrell') require(rms) load('data/counties.sav') older <- counties$age6574 + counties$age75

Dear All, According to the Sturges rule, the number of classes of a histogram is the closest integer to 1 + logb(n,base=2) where n is the number of observations. The function hist(), by default, uses the Sturges rule. However, the code x <- 1:200 hist(x) produces a histogram with 10 classes and not 9 classes as determined by the Sturges rule. What am I missing? Thanks in advance, Paul

hello, I want to iterate through a list of names and use each element as an argument in a function. For instance: > a = c('one','two','three') > data= c() > for(elem in a){data=cbind(elem = 2,data)} > data elem elem elem [1,] 2 2 2 instead I want 'elem' to be substituted by the string in the list. Doing it by hand would be: > data =

Hi, I'd like to report something here, so it will not happen to others. We moved all disabled users in our samba AD to a dedicated folder in ADUC, which we called 'disabled'. A little while after we did that, our network started 'falling apart'. Some things still worked, others did not. I could for example no longer start ADUC, some users could not logon or map drives, etc,

[ Not sure why, but the first time I sent this it never seemed to go through; apologies if you're seeing this twice ... ] I have some fully functional code that I'm guessing can be done better/quicker with some savvy R vector tricks; any help to make this run a bit faster would be greatly appreciated; I'm particularly stuck on how to calculate using "row-wise" vectors

Yes, Your totaly correct, but readability get higer prio for me then not showing defaults im my setup. Just bit easier to read and it does not harm. This is also why my configs have the lines like this : # Keep no in production, set yes when debugging, this slows down your samba. Readability is also a big thing, and often forgoten/misplaced. Greetz, Lousi > -----Oorspronkelijk

Hi, I can't figure out why it is not possible to compute an exact p-value in cor.test if there are ties between values in one of the arrays like below: cor.test(c(1,2,2), c(5,6,7), method="k", alternative="two.sided"). Perhaps this is due to my lack of understanding of what is ment by p-value in this case. To me it seems reasonable that the p-value above should be the

I've seen a number of corrupted-index problems pop up on our test-dovecot servers (used by a few people at the office), and I'm wondering why corrupt indices don't just automatically get deleted ? Something like: if (!MAIL_INDEX_IS_IN_MEMORY(index)) { unlink(index->filepath); } in src/lib-index/mail-index.c:mail_index_set_error() comes to mind. It

Hi all, I'm a PhD student at sydney uni and am trying to run a non linear mixed model program to obtain estimates of parameters describing dairy cow lactation curves. At present, I have been able to get the data to converge using the S plus (S plus 2000) nlme function. However, when I put the same data into R (R 1.9.0), add in the nlme package and run the code, it does not converge by the

What is wrong in the following nested if-else statements: if (Condition_1) { # begin IF_1 statement_1 statement_2 statement_3 if (Condition_2) { # begin IF_2 a<- a +1 } # end IF_2 statement_4 statement_5 statement_6 statement_7 if (Condition_3) {

Hi, we have quite a few Asterisk machines running and try to keep them on a current version of the Asterisk 11 branch. But since we upgraded to 11.14.0 a couple weeks ago, we have to restart the Asterisk process every week because the load gets too high and our monitoring complains. Those machines are doing only SIP-to-SIP call relay, the dialplan is quite complex, transcoding is done only on a

I have a MLE task that for a small number of parameters finishes in a reasonable amount of time, but for my "real" case (with 17 parameters to be estimated) either takes far too long (over a day), or fails with "computationally singular" errors. So a) are there any parallel implementations of optim() (in R or otherwise) and b) how can I make my function more robust?

Well, im never trying to change your optionion Rowland, i wont dare.. > Just what has an init system got to do with mounting a filesystem ?? Compaired to fstab, i get a better result with my nfs mounts. Before, somethimes, with the fstab setup, when logging in on ssh on an automounted userhome dir, with NVSv4 kerberized. Did fail for now and then. Now with the systemd, i didnt have any