Displaying 20 results from an estimated 2000 matches similar to: "calculating ratios from all combinations"
2011 Aug 18
4
splitting sample names
I have a data frame like this
xx<-data.frame(cbind(Sample=c('Ctrl_6h','1+0_6h','1+200_6h','1+5k_6h','Ctrl_5K_6h','ConA_6h'),
IFIT1=c(24,25,24.7,24.5,24.2,24.8)))
grep('[[:digit:]]h',xx$Sample)
yy<-xx$Sample
strsplit(yy,"_")
I have to extract the time information separated by '_' in the sample
2011 Nov 07
2
help with formula for clogit
I would like to know if clogit function can be used as below
clogit(group~., data=dataframe)
When I try to use in above format it takes a long time, I would appreciate
some pointers to get multiple combinations tested.
set.seed(100)
d=data.frame(x=rnorm(20)+5,
x1=rnorm(20)+5,
x2=rnorm(20)+5,
x3=rnorm(20)+5,
x4=rnorm(20)+5,
x5=rnorm(20)+5,
x6=rnorm(20)+5,
x7=rnorm(20)+5,
x8=rnorm(20)+5,
2011 Oct 27
2
help with paste
Hello gurus,
I have some variables, and i am creating combinations for analysis in the
end i need these variables to be displayed like "LEPTIN+SAA+PTH". currently
i am using loop to perform this. I would appreciate any pointers to do it
without the loop.
> mols=c("LEPTIN","SAA","PTH","sEGFR")
> samples=mols[1:3]
> samples
[1]
2011 May 19
3
A better way to do this
Hello gurus,
I have a dataframe containing two groups viz., 'control' and 'case', each of
these groups contains longitudinal data for 100 subjects. I have to plot all
these subjects on a single chart and then put a regression line for each of
the group for all the subjects. I have written a function to do the chart
grpcharts<-function (dat, group,group2,molecule,cutoff){
2010 Oct 20
3
Plot creates a straigth line
Hello all, I am using 'plot' to create standard curves for elisa data. when I
use 'plot' with type='b' i get the points connected with lines and one
straigth line from the lowest datapoint to the highest data point. how can i
avoid/remove it from the figure. i am using R2.9.1, below is the example of
the data.
od<-c(10, 8, 6,4,2,1, 10.5,7.8,6.4,3.8,2.1,0.95)
2011 Apr 13
1
error for ttest
Hello all,
I have arranged my data as per Dennis's suggestion in this post
http://www.mail-archive.com/r-help at r-project.org/msg107156.html.
the posted code works fine but when I try to apply it to my data, i get ">
u2 <- ddply(xxm, .(plateid, cytokine), as.data.frame.function(f))
Error in t.test.formula(conc ~ Self_T1D, data = df, na.rm = T) :
grouping factor must have
2011 Dec 20
5
Help with code
hello gurus,
i have a data frame like this
HTN HTN_FDR Dyslipidemia CAD t1d_ptype[1:25]
1 Y Y Y T1D
2 T1D
3 Ctrl_FDR
4 T1D
5 Y Ctrl
6 Ctrl
7
2007 Mar 01
1
problem with throwing lines out of matrix
Dear all,
again i have a problem with special case of dropping lines out of a
matrix.
i tried the following, where throwout is a vector of length 18 with
simple id values that should be compared to any
for (k in 1:length(throwout))
{
mymatrix=matrix(mymatrix[-(forfact[k]),],ncol=4)
}
this didnt work and i tried the following to find the error:
mymatrix[,1][mymatrix[,1]==throwout[7]]
2011 Oct 10
4
correlation matrix
Hello Gurus
I have two correlation matrices 'xa' and 'xb'
set.seed(100)
d=cbind(x=rnorm(20)+1,
x1=rnorm(20)+1,
x2=rnorm(20)+1)
d1=cbind(x=rnorm(20)+2,
x1=rnorm(20)+2,
x2=rnorm(20)+2)
xa=cor(d,use='complete')
xb=cor(d1,use='complete')
I want to combine these two to get a third matrix which should have half
values from 'xa' and half
2011 May 04
1
merging multiple columns from two dataframes
Hello,
I have data in a dataframe with 139104 rows which is multiple of 96x1449. i
have a phenotype file which contains the phenotype information for the 96
samples. the snp name is repeated 1449X96 samples. I haveto merge the two
dataframes based on sid and sen. this is how my two dataframes look like
dat<-data.frame(snpname=rep(letters[1:12],12),sid=rep(1:12,each=12),
2012 Jan 29
1
height of plots
Hello R gurus,
I have to create 12 plots, I have been using the following script, which
leaves a large white space between two plot. I would appreciate if someone
can suggest an alternative to reduce the white space.
par(mar=c(3,3,.5,.5))
split.screen(c(6,2)) # split display into two screens
for (i in 1:12)
{
if (i<11)
{
screen(i)
plot(1:10,xaxt='n', xlab='',
2011 Jan 08
4
Help How to create DSL for conditional validations
Hi All,
I am new to Ruby and ROR
I were trying to create small DSL for conditional validations
valid_with_cond :bypass_validation do
if self.addresses > 3
errors[:base] << "Can not have more than 3 addresses".
end
end
By this I wanted to create array of method and call them all in custom
validation method.
this above code I wanted to do
attr_accessor
2009 Apr 09
1
[Bridge] Out of memory problem
Hi, I'm using linux 2.6.21.5 and our kernel is freeze.
The problem is, if I create a Software bridge using $brctl command. and
add two interfaces say, eth0.0 and eth0.1 using
$brctl addbr br-lan
$brctl addif br-lan eth0.0
$brctl addif br-lan eth0.1
and when i send traffic from a host connected to one port to host
connected at other end, soon all the memory is dried up and and kernel
2009 Aug 28
1
extracting pvalues from ttest
Hello list,
I have a similar issue as this post
http://tolstoy.newcastle.edu.au/R/e6/help/09/04/11438.html#options2 and I
used the suggestion provided by Jorge with modifications to my data
do.call(c,lapply(your_list_with_the_t_tests,function(x) x$p.value))
but I am getting the following error after excuting the code
B<-by(eo,eo$PlateID, function(.sub) t.test(mcp1~Self_T1D,data=.sub,
2015 Apr 23
2
[LLVMdev] Get precise line/column debug info from LLVM IR
I am trying to locate instructions in an LLVM Pass by line and column
number (reported by an third-party tool) to instrument them. To achieve
this, I am compiling my source files with `clang -g -O0 -emit-llvm` and
looking for the information in the metadata using this code:
const DebugLoc &location = instruction->getDebugLoc();
// location.getLine()
// location.getCol()
2009 Mar 07
2
ttest in R
Dear list,
i am a biologist who needs to do some ttest between disease and non disease,
sex, genotype and the serum levels of proteins on a large number of
individuals. i have been using excel for a long time but it is very tedious
and time consuming. i am posting the data below and ask your help in
generating a code to get this analysis done in R. thanks
gender disease genotype data
M N CC
2008 Aug 23
3
Postfix and Dovecot SASL authentication error
ear All,
I configured Postfix 2.3.3-2 on CentOS 5.2 using SASL Dovecot
authentication put I am facing the problem to receive mails . My postconf -n
result is:
[root at mails ~]# postconf -n
alias_database = hash:/etc/aliases
alias_maps = hash:/etc/aliases
command_directory = /usr/sbin
config_directory = /etc/postfix
daemon_directory = /usr/libexec/postfix
debug_peer_level = 2
2003 Feb 27
2
epoch time conversion in R
I have a data file where each entry is indexed by the time in seconds since
epoch (e.g. 1046315697). Is there an easy way to convert this time value into
a more friendly time (such as Month-Year) when plotting it?
I searched through the manual, mailing lists, and functions like as.POSIXct
and strptime, but didn't find what I need.
Thanks,
Sharad.
2010 Sep 16
1
virConnectOpen question.
Hi,
I am using virConnectOpen(const char *uri) to connect to the local
hypervisor. When uri is set to "lxc:///", an error gets logged
in /var/log/messages about failed connection. But I do not see any error
for "xen:///". The machine does not have these drivers installed, and I am
using libvirt 0.6.3. My questions are -
1. Why I see error in /var/log/messages for container
2011 Aug 11
2
virInitialize dumps core intermittently.
Hi,
I have a host with 6 VMs set to autostart. When this host reboots, it
starts libvirtd, pegasus cimom and libvirt-cim. It also restarts the VMs.
Intermittently we noticed a core during start up that points to the last
call made by libvirt-cim to "virInitialize". I am thinking that it is a
timing issue where virInitialize is called while libvirt is trying to
autostart the VMs.