similar to: what is the "NOT IN" operator

Good day, I'm ran aov () and lme() on a split-plot using a 64-bit machine. For aov() I don't see the values for ErrorA, F-value and p-value in the output. For lme(), output is different from results from a 32-bit. Please see codes used and corresponding output. Is my code wrong and/or not sufficient or is this a compatibility issue? ************** model1<-aov(Y~Main*Sub +

I'm running 10, 000 iterations each for the bridge and blasso. 3, 000 iterations roughly takes a week in a core-duo processor with 16GB RAM. I'll have access to a 6C processor machine and I came across the multicore package. Can I use multicore with the bridge and blasso function? That is, is parallel processing possible for the bridge and blasso functions by using the multicore

For various reasons, it suits our workplace to have a cygwin version of R. I am pretty sure that cygwin is still not a supported environment for R, but we have managed to compile R-2.5.1 under cygwin without too many dramas. Our procedure is described below. We still have a few problems compiling libraries without manually changing files from .so to .dll, but it seems ok. I was wondering

Hi! I have fitted two glms assuming a poisson distribution which are: fit1 <- glm(Aids ~ Year, data=aids, family=poisson()) fit2 <- glm(Aids ~ Year+I(Year^2), data=aids, family=poisson()) I am trying to work out how to represent the fitted regression curves of fit1 and fit2 on the one graph. I have tried: graphics.off() plot(Aids ~ Year, data = aids) line(glm(Aids ~ Year,

Dear all, I imagine that this is a trival question, but it has perplexed for most of the day. Any help would be greatly appreciated. Below is an example of what I'm trying to do. Essentially I want to produce all unique 1 x 1 combinations of certain columns from a dataframe, and join these to other columns from the same dataframe. I'm having problems with the nested loop as I can only

Hi all, How can I have a dynamic list for different combinations of grouping factors in the following example? Thanks. dat <- data.frame(x=rnorm(100), a=sample(letters[1:5], replace = T), b=sample(letters[1:5], replace = T), c=sample(letters[1:5], replace = T), d=sample(letters[1:5], replace = T)) ## define grouping

How about this. Its not as 'swish' as Doug's response, but it might be a little simpler to think about and adapt... # show the test vector > tmp [1] A B C B C A C B A A B Levels: A B C > # create a named vector to store the number of occurences of each level > index <- numeric(length=length(levels(tmp))) > names(index) <- levels(tmp) > index A B C [1] 0 0

Hi, I am using the scan function to read from a file. In general the syntax I am using is as follows: scan("myfile", list(a = "", b="", c="", d=""), sep=",") The file I am trying to read contains a #< in the 'd' field. When the file is read into an object, the d field comes as "" for that particular string. I

I have written a function to plot data which will be used for various different chemistries. A simplified version is: plot_data <- function(risk,levels,chem,sd2,measure){ plot(risk, levels,main=paste ("per", sd2, measure, "\n in usual", chem)) } The problem is with the title. This works fine if the variable "chem" is just text, but if it is an expression then

I have recently compiled and installed openssh-2.1.1p1 on a linux box. The login reporting does not seem to work properly. When logging into the box via ssh (protocol 1) utmp shows the user logged in and the tty properly, but the field for the login date/time and the field for originating host contain all NULLs. Is anyone else seeing this same behavior, or have I just done something really

Hello, my enviroment: All Servers are Ubuntun 16.04-18.04 SAMBA AD DC Server and several SAMABA DOMAIN MEMBER (connected via WINBIND). In ADDC I've created a group "restrictaccess" and added some users. Now when im typing "id <username>" on a Domain Member, for some users the group "restrictaccess" are listed for some not! For example: ON DC: #

Hello, I’m reaching out because our online radio station (sites.up.edu/kdup/) has stopped streaming in the last week - the website itself works fine but when anyone attempts to listen to the live music stream they receive a message that says “cannot connect to the server “icecast.up.edu.” Do we need to update the server, or is this a problem that someone else can take care of? Let me know if I

rsync does not sync the timestamp on symlink (Solaris 8). It is probablly due to the limitation of Unix implementation of symlink, but I would like to know why rsync/Unix does not do this, and what we can do about it. Is the conclusion that "rsync syncs everything except the timestamp on symlink"? Why do I need timestamp on symlink? Supposed something stopped working because something

Hi, please can you help me. When I add a legend to a boxplot it appears inside the plot area, how do you get it into the margins? I have already changed the parameters so there is space for it on the margin on the right hand side of the graph. Thanks, Emily [[alternative HTML version deleted]]

I'm trying to figure out if it's possible to use a character string as an equation, e.g: eqn1string <- "x^2 + x + 5" Then I want to tell R: 1) that eqn1string is actually an equation (even though it was stored as a character string), and 2) to apply the equation to a specified value of x (e.g. given x <- 6, what is the result of the equation). Thanks in advance for any

Hello List, ? I have a dataset consisting of strings that I want to split while saving the delimiter. ? Some example data: ?leucocyten + gramnegatieve staven +++ grampositieve staven ++? ?leucocyten ? grampositieve coccen +? ? I want to split the strings such that I get the following result: c(?leucocyten +?, ??gramnegatieve staven +++?, ??grampositieve staven ++?) c(?leucocyten ??, ?grampositieve

I thought replacing the spaces following instances of +++,++,+,- with "\n" and then reading with scan should succeed. Like Ivan Krylov I was fairly sure that you meant the minus sign to be "-" rather than "?", but perhaps your were using MS Word as an editor which is inconsistent with effective use of R. If so, learn to use a proper programming editor, and in any case

I have a data frame with three columns, one coded as a factor. I would like to separate my data out into separate data frames according to the factor level. Below is a simple example to illustrate. I can get R to return the data in the correct format but cannot work out how to get separate data frames. I am a newcommer to programming in R so am not sure what I am missing! Thanks, Emily

Hi, I'm new to R and I have a file with many rows of values. Each row contains a title and values for a contingency table e.g. row 1= title 8 0 37796 47 which is a table called 'title' with values 8 0 37796 47 I would like to know how I can import this using R and for each row calculate a p value using the fisher test. Using each p value I will do multiple a correction. I am

Dear all, I have many csv files whose contents I want to change a bit en masse. So far, I've written code that can change them in a for loop, like so: # Subset of files in the folder I want to change subset = "somestring" # Retrieve list of files to change filelist=list.files() filelist = filelist[grep(subset, filelist)] for(i in 1:length(filelist)){ setwd(readdir)