Displaying 20 results from an estimated 3000 matches similar to: "group data"
2010 Nov 29
2
FW: how to use by() ?
Thank you for the suggestion, Bill. The result is not quite what I would like. Here's sample code for you or anyone else who may be interested:
Al1 = c('A','C','C','C')
Al2 = c('G','G','G','T')
Freq1 = c(0.0078,0.0567,0.9434,0.9908)
MAF = c(0.0078,0.0567,0.0566,0.0092)
m1 = data.frame(Al1=Al1,
2009 Sep 18
1
Incorrect Dimension
I am new in R and i am having trouble here. I?ve already searched in the list
but hasn?t helped
When i run this script above i get the message "Error in gen[j, i] :
incorrect number of dimensions". However gen is 1000x200 (ind x loc) and so
is g
could anybody help me
for (i in 1 : loc) { #loc=200
for (j in 1 : ind) { #ind=1000
g1 <= function ( gen ) matrix ( if (gen[j, i]
2010 Nov 29
3
how to use by() ?
Hello, All!
How might one accomplish this using the by() function?
m1 is a data frame.
# populate column "m1$major_allele"
for ( i in 1:length(m1$major_allele)) {
if ( m1$Freq1[i] == m1$MAF[i]){
m1$major_allele[i] = m1$Al1[i]
}
else{
m1$major_allele[i] = m1$Al2[i]
}
}
Jim
[[alternative HTML version deleted]]
2008 Mar 07
3
Combine two columns
Is there a way to combine two columns within a data frame?
Example data:
id snp AL1 AL2
1500 30 A B
1510 30 A A
1520 30 A B
This is what I would like:
indv snp AL1AL2
1500 30 AB
1510 30 AA
1520 30 AB
Any help is greatly appreciated.
Alysta
2009 Jun 03
1
Would like to add this to example for plotmath. Can you help?
Greetings:
I would like comments on this example and after fixing it up, I need
help from someone who has access to insert this in R's help page for
plotmath.
I uploaded a drawing
http://pj.freefaculty.org/R/Normal-2009.pdf
that is created by the following code
http://pj.freefaculty.org/R/Normal1_2009_plotmathExample.R
This will be a good addition to the plotmath help page/example.
2010 Dec 09
4
lapply getting names of the list
Hello All,
I have a toy dataframe like this. It has 8 columns separated by tab.
Name SampleID Al1 Al2 X Y R Th
rs191191 A1 A B 0.999 0.09 0.78 0.090
abc928291 A1 B J 0.3838 0.3839 0.028 0.888
abcnab A1 H K 0.3939 0.939 0.3939 0.77
rx82922 B1 J K 0.3838 0.393
2006 Mar 25
1
Trouble phrasing an R command that will run the model i need (ANOVA, nes
Hi John,
Thank-you very much for replying, - in answer to your questions:
*"As I understand what you wrote, all factors but linelabel are crossed,
with linelabel nested within all the other factors?"
That's exactly it,
*" What is in the highest order cells? Single observations, or multiple
scores?"
I would say "single observations" - but i have to
2012 Jan 22
1
passing "..." arguments to (plot.ca)
I am hoping someone can give me a couple of pointers on how to pass
arguments using "..."--in this specific case in an attempt to plot an
object created by ca(), a tool for correspondence analysis.
#some illustrative code
library(ca)
set.seed(123)
dat1 <- data.frame(matrix(ceiling(runif(30, 1, 100)), nc=5)) #some fake data
ca1 <- ca(dat1) #the correspondence analysis
plot(ca1)
2013 May 07
4
create unique ID for each group
Hey All,
I have a dataset(dat1) like this:
ObsNumber ID Weight
1 0001 12
2 0001 13
3 0001 14
4 0002 16
5 0002 17
And another dataset(dat2) like this:
ID Height
0001 3.2
0001 2.6
0001
2013 Apr 12
2
split date and time
Hi R experts,
For example I have a dataset looks like this:
Number TimeStamp Value
1 1/1/2013 0:00 1
2 1/1/2013 0:01 2
3 1/1/2013 0:03 3
How can I split the "TimeStamp" Column into two and return a new table like
this:
Number Date Time Value
1 1/1/2013 0:00 1
2 1/1/2013 0:01 2
3 1/1/2013 0:03 3
Thank!
[[alternative HTML version
2010 Mar 09
3
Shade area under curve
I want to shade the area under the curve of the standard normal density.
Specifically color to the left of -2 and on. How might i go about doing
this?
Thanks
--
View this message in context: http://n4.nabble.com/Shade-area-under-curve-tp1586439p1586439.html
Sent from the R help mailing list archive at Nabble.com.
2007 Apr 23
1
problem with 3-way conferenicing
Hi,
I am trying to achieve 3-way conferencing taking hint from wiki link
http://www.voip-info.org/wiki/view/Asterisk+n-way+call+HOWTO
Here is the scenario:
1. user "ua1" calls user "ca1"
2. "ua1" then presses the feature code "*0" to redirect "ca1" to
conference room 300
3. "ua1" then dials the user "33"
4. user
2014 Apr 22
2
Re: TLS and intermediate CA
Thanks for the response.
My current chain is as follows:
caroot -> child-ca1 -> server cert
My cacert.pem file has both the caroot and the child-ca1 certs. I have
recompiled libvirt on my machine with some extra debug statements and
verified that both the caroot cert and the child-ca1 certs are being
loaded. But when I try to connect the caroot and child-ca1 certs only
appear under the
2008 Nov 20
2
Reformatting a table
Hi !
I am new to R. Can somebody help me in reformatting huge output files ,i.e, rearranging sets of columns in specific order.
For example: I have data for three compunds 1, 2 and 3
file1:
ID CA1 CA3 CA2 MA2 MA1 MA3
1 14 15 13 7 12 3
2 19 7 12 10 14 5
3 21 12 19 6 8 9
to
File 2:
ID CA1 CA2 CA3 MA1 MA2 MA3
1 14 13 15 12 7 3
2 19 12 7 14 10 5
3 21 19 12 8 6 9
or File3:
ID
2009 Sep 17
3
How to separate a function by 2 probabilities
Good Mourning,
I have a function to generate a matrix as I show part of it;
g[j,i]<-if (gen[j,i]==0) al1[i,1]+al1[i,1] else ...
However i would like that this function occurred with a probability P and
that another function (another formula to generate g matrix) with
probability P-1
That?s it, if P is .7, i would like that in 70% of the times (for random i
and j) the matrix g was generated
2012 May 10
0
Time series and stl in R: Error only univariate series are allowed
I am doing analysis on hourly precipitation on a file that is disorganized.
However, I managed to clean it up and store it in a dataframe (called CA1)
which takes the form as followed:
Station_ID Guage_Type Lat Long Date Time_Zone Time_Frame H0
H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H11
H12 H13 H14 H15 H16 H17 H18 H19 H20 H21 H22 H23
1
2005 Sep 30
0
p-value for non-linear variable in overdispersed glm()
Dear all,
I am fitting an nonlinear glm() using optim() by first minimising
glm(resp~ var1 + var2, family=binomial, data=data)$deviance
where var1= exp(-a1*dist1), and var2= exp(-a2*dist2), where a1 and a2 are
parameters and dist1 and dist2 are independent variables.
Next, I calculate the value of var1 (and var2) by plugging in the value of
al1 (and al2) that minimises deviance,
and fit
2008 Mar 25
2
Compare two data sets
I would like to compare two data sets saved as text files (example below) to determine if both sets are identical(or if dat2 is missing information that is included in dat1) and if they are not identical list what information is different between the two sets(ie output "a1", "a3" as the differing information). The overall purpose would be to remove "a1" and
2012 Jun 07
3
conditional statement to replace values in dataframe with NA
Hello and thanks for helping.
#some data
L3 <- LETTERS[1:3]
dat1 <- data.frame(cbind(x=1, y=rep(1:3,2), fac=sample(L3, 6, replace=TRUE)))
#When x==1 and y==1 I want to replace the 1 values with NA
#I can select the rows I want:
dat2<-subset(dat1,x==1 & y==1)
#replace the 1 with NA
dat2$x<-rep(NA,nrow(dat2)
dat2$y<-rep(NA,nrow(dat2)
#select the other rows and rbind
2013 Mar 27
9
conditional Dataframe filling
Hi everyone:
This may be trivial but I just have not been able to figure it out.
Imagine the following dataframe:
a b c d
TRUE TRUE TRUE TRUE
FALSE FALSE FALSE TRUE
FALSE TRUE FALSE FALSE
I would like to create a new dataframe, in which TRUE gets 0 but if
false then add 1 to the cell to the left. So the results for the
example above should be something like:
a b c