similar to: aggregating data

I am really new to R and would appreciate some help to sort out a problem with data extraction from a large file. I have entered the following command to filter data from a large data set called "upanddown" such that those records in upanddown corresponding to blockType=3'UTR are transferred to a new file called "new": new<-sqldf("select * from upanddown where

Greetings list, I am new to programming in R, and am using nnet() function for a project on neural networking. Firstly I wish to ask if there is any pdf explaining the algorithm nnet uses, which could tell me what the objects of the nnet class, like 'conn', 'nconn, 'nsunits', n and 'nunits' mean, and how weights are calculated. The package pdf has little or no

Hi, I am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

Hello, Thanks everyone for helping me with the previous queries. step 1: Here is the orginal data: short sample ProbeID Sample_1_D Sample_1_C Sample_2_D Sample_2_C 1 224588_at 2.425509867 11.34031409 11.46868531 11.75741478 step 2: i calculate the variance of the sample using this R code x<-1:20000 y<-2:141 data.matrix<-data.matrix(data[,y])#create data.matrix

I have a file with a data in columnar format like below: probeID rc_AI104113_at rc_AI178259_f_at rc_AI179134_i_at rc_AI179134_f_at rc_AI104113_at rc_AA819429_f_at How can I rewrite it in the format below: 'rc_AI104113_at', 'rc_AI178259_f_at', 'rc_AI179134_i_at', 'rc_AI179134_f_at', 'rc_AI104113_at', 'rc_AA819429_f_at' Is there any function to do

Here is my R Code x<-1:20000 y<-2:141 data.matrix<-data.matrix(data[,y])#create data.matrix variableprobe<-apply(data.matrix[x,],1,var) variableprobe #output variance across probesets hist(variableprobe) #displaying histogram of variableprobe write.table(cbind(data[1], Variance=apply(data[,y],1,var)),file='c://variance.csv') #export as a .csv file. Output in Excel all in 1

I am trying to use lme to fit a mixed effects model to get the same results as when using the following SAS code: proc mixed; class refseqid probeid probeno end; model expression=end logpgc / ddfm=satterth; random probeno probeid / subject=refseqid type=cs; lsmeans end / diff cl; run; There are 3 genes (refseqid) which is the large grouping factor, with 2 probeids nested within each refseqid,

Hello R-help, I currently have R-project 3.4.2 and R-studio 1.1.383 installed on some of our universities computer labs. Since the installation of Visual studio the default version of R installed has changed to the version installed by VS. Users can change the default R version: found in R studio global options but users need to do this every session. Is there a file or

> On Mar 15, 2018, at 10:24 PM, Adrian Friskin <adrian.friskin at qut.edu.au> wrote: > > Hello R-help, > I currently have R-project 3.4.2 and R-studio 1.1.383 installed on some of our universities computer labs. Since the installation of Visual studio the default version of R installed has changed to the version installed by VS. Users can change the

Hi Everyone, Thanks for all the help with the previous queries. Here is what i want to do. i have 20000 probesets-->calculate all the variance accross all the probesets-->filter out probesets that are low so now i ended up with only 10000. The 10000 is fine but when i export to excel, it is missing the probeID. Here are my code and examples. #########calculate the variance across the

G'Day, I want to access in a function called from tapply the current factor. In my example below, all I want to do is to write the current factor on each histogram. Needless to say my example does not work. I would be grateful for pointers in the right direction. Many thanks Bernie McConnell Sea Mammal Reserach Unit cc <- 1:10 ff <- rep(c("a","b"),5) pp<-

Hi, I'm having problems with the queue on a Tru64 box running Samba 2.2.4. When the queue is viewed from a client machine it displays a print job with fields equivalent to the first line the lpq command outputs. For example, lpq returns: # lpq -P itsgp_css1 crow.qut.edu.au: Tue May 14 09:23:54 2002: no entries Then if I get on a client machine (Win2k) and look at the queue there is a

Again, this is off topic, not about statistics or R, but I think of interest to many on this list. The title is: "So you got a null result. Will anyone publish it?" https://www.nature.com/articles/d41586-024-02383-9 Best to all, Bert

Hi, I was wondering if there is a paper or a technical report that documents the performance impact of the different optimizations passes on a some set of benchmarks. Is something like this available ? Best regards, Ahmed -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.llvm.org/pipermail/llvm-dev/attachments/20150619/c3c0a941/attachment.html>

Dear R users, This is a problem that has puzzled me for quite some time, and I wonder if anyone could offer me an insight to what is happening? I generated a multinomial matrix by using A<-rmultinom(100,1,c(0.4,0.3,0.2,0.1), and multiply this matrix by a set of values say B<-runif(n,100,50). I then calculate var(colSums(A*B)). I was able to work out the theoretical mean of

Hi everyone, I?m a graduate student in a Brazilian Intitute of Technology, and I?m doing some academic research regarding secure voice transmission over phone lines. One of our reserach goals is to implement a hardware vocoder, with low bit rates, and a preferably free algorithm, to be used in this secure voice system. Actually, there is a functional system using a proprietary AMBE

Hi, I install a guest Fedora inside the host Fedora. When I use the service xend status to check the status of xend. It shows xend is running. However when I use xm create xen1 to start the guest domain xen1, it failed. The corresponding file xen1 is in /etc/xen/xen1. The error msg is as follows: [root@phillis xen]# xm create xen1 Using config file "./xen1". Going to boot Fedora

I have a similar question, I know about --vsize, but I use R under Emacs using ESS. I know there is an easy way to call R from ESS with command line options but I don't remember how. I can't find help on this in ESS documentation or in R documentation. Could someone please remind me how this is done. ------------------------------------ | Robert Denham | |

Hello, I do not have an account on R Bugzilla, so I will post my bug report here. I want to report a very old bug in base R *source()* function. It relates to sourcing some R scripts in UTF-8 encoding on Windows machines. For some reason if the UTF-8 script is containing cyrillic letter *"?"*, the script execution is interrupted directly on this letter (btw the same scripts are sourcing

Dear R users, I am using the beadarray package. I am trying to upload raw bead-level data using these commands: ######################################################## library(beadarray) datadir <- ("C:/Computer_programs/R/beadarray/cecilia") targets = read.table("targets.txt", sep = "\t", header = TRUE, as.is = TRUE) BLData = readIllumina(arrayNames =NULL,