similar to: cat in windows vs linux

Hello all, A colleague at work set me the challenge to convert some MATLAB code into R, to see which is faster. We'd seen that benchmark comparing MATLAB 6.5 to R1.90 (and others), and so I thought that I should be able to get roughly comparable speeds. The code has lots of multiplications of matrixes, transposes, and MATLAB's "repmat". I did the code conversion, and R was about

Hello all, A really noddy question for you all: I''m trying without success to do some subhypothesis testing. Using simple anova model, with a toy dataset from a book. I have four factors A,B,C,D, and wish to test mu_C = mu_D. This is what I have tried: > contrasts(infants$group,how.many=1) <- c(0,0,1,-1) > contrasts(infants$group) [,1] A 0 B 0 C 1

Full_Name: Rob Foxall Version: 2.2.0 OS: Windows XP Submission from: (NULL) (149.155.96.5) I was using multinom from nnet package, when I did something stupid -- I entered in an incorrect factor variable as response. This factor had only one level. Instead of R telling me not to be so dumb, it crashed, clicking on debug coming up with the message "An exception 'Unhandled Win32

Dear R-helpers, I have a model created by lda, and I would like to use this model to make predictions for new or old data. The catch is, I want to do this without using the "predict" function, i.e. only using information directly from the foo.lda object to create my posterior probabilities. In anticipation of likely responses, I will be brushing up my lda knowledge using the

Dear R-Help, I have written a function, which in simplified format can be represented as: temp.fn <- function(my.mean,my.sd){ Parameters <- list(mean = my.mean, sd = my.sd) curve(do.call(dnorm,c(list(x), Parameters)), from = my.mean-my.sd, to = my.mean+my.sd) } This works as I want it do. Wishing to immortalise this function into my very own package however, results in the following

I have a 'cat' function within a for loop. I would like it to print out the result everytime it goes through the for loop so that I can monitor the progress, but it only prints out (execute the cat function) once at the very end of the for loop. A simple example: for(i in 1:100000){ cat(i, '\n') } What should I do in order for the cat function to be executed everytime it goes

Hello, Probably a stupidly easy question, but I have done the following in order to make predictions from a fitted glm with new data: my.glm <- glm(lt96~so296[,1:17],family=binomial(link=logit)) p96 <- predict.glm(my.glm,newdata=so293[,1:17],type="response") but I always get the fitted linear predictors from the original model, ie there doesn't seem to be acknowledgement of

Hello All I am having a bit of trouble formatting my HTML with the desired number of digits after the decimal place. Am I doing something wrong/misunderstanding or is it a bug? Looking at the example supplied with ?HTML.data.frame: HTML(iris[1:2,1:2],nsmall=c(3,1),file="") Gives html output that includes the lines: </tr> <tr><td

using an example from r online help > state <- c("tas", "sa", "qld", "nsw", "nsw", "nt", "wa", "wa", "qld", "vic", "nsw", "vic", "qld", "qld", "sa", "tas", "sa", "nt", "wa", "vic",

Hi R guru coders I wrote a bit of code to add a new column onto a "topTable" dataframe. That is a list of genes processed using the limma package. I used a for loop but I kept feeling there was a better way using a more vector oriented approach. I looked at several commands such as "apply", "by" etc but could not find a good way to do it. I have this feeling there

All, Does anyone know if there is an option I can set to turn screen-buffered output on and off with the win32 rgui? (Apart from the point and click method). I am running some simulations where it is useful to watch output but it gets mildly tiresome having to manually switch things on and off via the gui. Thanks Toby. > version _ platform i386-pc-mingw32 arch

Dear All, It is not a typical R question (though I use R for this) but I thought someone will help me. For the list of P values, I have calculated FDR using p.adjust() in R (bioconductor). But my FDR values are same for all the P values. When do we get same FDR values? Does the smallest P values should less than 1/N? (where N is the number of P values) Thanks in advance. Kind regards, Ezhil

Hello, I have an expression for the y axis, which I wish to rotate by 90 degrees (so that this label has the same orientation as for the x axis). Thought the graphical parameter crt (or srt) would do the trick, but haven't had any joy so far. Any ideas? (Using Version 1.2.2, linux, redhat 7.0) Cheers, Rob. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Dear R-help, I have a character vector, some elements will be numeric, some not, and some even empty. E.g.: temp1 <- c("abcd"," 2 ","") I'm only interested in the numeric elements, the rest I can just throw away. It is easy enough to loop through the vector: temp <- try(eval(parse(text=temp1[1])), silent=TRUE); class(temp) # try-error temp <-

Since I only seem to get questions, and no feedback, from the Wiki page, I'll ask here. There seems to be no lack of opinions here... I have a working wakeup call system on my home * system. The architecture is something I'm not perfectly happy with, though. There are two AGI scripts, written in Perl, which handle (a) scheduling, confirming, and cancelling a wakeup call, and (b) the

Hi Jenny If you want a general solution I understand. However I just downloaded the file fine (as far as I can tell) so you are welcome to have a copy. I can email it to you if you want. I do not think your test for NA is valid. i.e if(test != "NA"){ } I think you should use if(is.na(test)){ } Or something similar. J --- John Seers Institute of Food Research Norwich

Hi R, I wanted to know how do we access the elements of a list. In particular, v=list(c(1,2,3,4,5),c(1,2,33,4,5),c(1,2,333,4,5),c(1,2,3333,4,5)) I want to access all the thirds items of the elements of the list. i.e., I want to access the elements, 3,33,333,3333. This can be done through sapply as: sapply(v,function(x) x[3]) But I need to access this without using

Hi All, My R code takes very long time to finish the processing. I want to see at what stage the script is running. So I wrote some output messages using cat. But instead of displaying the cat messages at different stages they are buffered and displayed in the end when entire processing is done. Can you please suggest how to stop this buffering or some alternative way to display messages Thank

Dear r-users, Suppose I have three datasets: Dataset-1: Date x y Jan-1,2005 120 230 Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21 Jan-7,2005 11 299 Mar-5,2005 200 311 Dataset-2: Date x y Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21

Hi all I'm looking for some help on winbind/idmap for a new host The debian version is new on this host : debian 10 buster so samba is > root at homedir10:~# samba --version > Version 4.9.5-Debian We want to use security=ads so we join this host to the domain No problem for windows clients : they can mount shares that are accessible to their primary unix group and secondary unix