Displaying 20 results from an estimated 3000 matches similar to: "How to separate a function by 2 probabilities"
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
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,
2013 Apr 11
2
group data
Hey,
I have a dataset and I want to identify the records by groups for further
use in ggplot.
Here is a sample data:
ID Value
AL1 1
AL2 2
CA1 3
CA4 4
I want to identify all the records that in the same state (AL1 AND A2),
group them as "AL", and do the same for CA1 and CA4. How can I have an
output like:
ID Value State
AL1 1 AL
AL2 2 AL
CA1 3 CA
CA4 4
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
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
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
2011 Sep 26
1
compute probabilities on a Bayesian Network
Deal R Users,
I'm trying to find out how can I compute probabilities on a Bayesian
Network using R. The Bayesian Network I modelled is shown at
http://www.dsr.inpe.br/~mello/1727/BNgrapmodel.png and I'd like to know
how can I proceed (commands on R) to answer questions like: (1) what is
the probability of S=T given C=T, L=T, R=F, H=F and D=F; or (2) what is
the probability of S=T
2009 Nov 25
5
How to sum only a few elements in a line
Hello,
I have a matrix with the numbers 0,1 and 9
I would like to write a function that could sum each line skiping everytime
a number 9 appears
for example
[0 1 0 1 1 9 1]
the sum would be 4.
However I cannot replace 9 by 0 otherwise after the sum is done I wouldn?t
be able to distiguish which ones were real zeros and which ones were nines
replaced by zero just to sum.
Thank you very much
--
2009 Sep 27
2
Count number of zeros in a collumn
I have a matrix 700x2000 which is sampled in each cycle from another matrix
788x2000 with the numbers 0,1 and 9
There is one specific collumn of this matrix, dart[,1977], that usually,
after the samplimg procedure has only 1 and 9 (because the zero frequency in
this collumn is low).
However, when this happens, I want to include an IF conditional in my code.
so basically what i wanted to do was to
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.
2009 Sep 23
1
BLUP with missing data
hello guys, I need to do a BLUP in the simplest model
y = Xm + Zg + e
however I have missing data in the analysis which I can?t consider as
0(zero). So I need to generate the matrix X'Z, Z'X and Z'Z step by step; I
can?t use
crossprod(x) #neither
X'X <- t(x)%*%x
because I should skip the elements with missing data in the matrix
I?ll try to be more clear,
supposing
a matrix x
2019 Nov 14
2
Dovecot Master User: Access user's mailbox without owner's password
Hi, I hope you're all right. I describe below the scenario where the
problem occurs.
I'm trying to activate a master user [1] to be able to access all the
boxes of all users by imap.
I have configured the dovecot-master-users [2] file with the
appropriate permissions.
When I try to access, for example with roundcube, through
user at mydomain.com*my_master_user at not-exist.com
I see in the
2009 Sep 22
2
chage the output file name
After I run my script I save all my output e.g.
write.table(d, file = "C:/PINEheight.txt)
write.table(v, file = "C:/PINEvolume.txt)
write.table(v, file = "C:/PINEdiameter.txt)
write.table(v, file = "C:/PINEdensity.txt)
and then I would like to run again with another tree and save again e.g
write.table(d, file = "C:/POPULUSheight.txt)
Is there anyhow I can define a
2009 Oct 15
1
Sampling procedure
I would like to divide a vector in 9 groups in a way that each number is
present in only one group.
In a vector of 783 I would like to divide in 9 different groups of 87
Example <- matrix(c(1:783),ncol = 1)
s1 <- as.matrix(sample(Example,87, re = FALSE))
Example <- Example[-s1]
s2 <- as.matrix(sample(Example,87, re = FALSE))
#however I don?t know how to remove the second group from
2010 Mar 26
2
Is there a faster way to do this?
Hi guys, I am still learning R, and not well familiar with all the apply
functions.
I am trying to find faster alternatives to replace the for cycle.
Is there a faster way to do the example below?
nm <- 1000
b <- matrix (rnorm (5000, 0, 1), nrow = 500, ncol = nm)
a <- matrix (0, nm, nm)
for (i in 1 : nm) {
for (j in 1 : nm) {
if ( j == i) {
next }
a[i, j] <- t (b [, i]) %*% b[, j]
}
}
2009 Oct 28
1
Is there a faster way to do it?
#Mdarts is a matrix 2343x788
#frequencia is a vector 2343x1
# 9 in Mdarts[fri,frj] stands for my missing values which i want to replace
by the value in the vector frequencia
Mdarts<-t(matrix(scan("C:/GWS/CNB/dartg.txt"),ncol=nindT,nrow=nm, byrow=T))
frequencia <- matrix(scan("C:/GWS/CNB/freq.txt"),ncol=1)
for (fri in 1:nindT){
for (frj in 1:nm){
Mdarts[fri,frj] <- if
2019 Aug 04
5
problems with authentication
Hi,
I set up a samba 4 in Debian 9.9 as a Domain member server, but
authentication is not working as follows:
root at srv-proxy:/etc/samba# wbinfo -a marcio at EMPRESA.COM.BR
Enter marcio at EMPRESA.COM.BR's password:
plaintext password authentication succeeded
Enter marcio at EMPRESA.COM.BR's password:
challenge/response password authentication failed
wbcAuthenticateUserEx(+marcio at
2012 Oct 19
2
PDC and BDCs : net rpc testjoin
People,
I have one PDC and a BDC on the matrix side and two BDCs on the branch
office.
I don't know if it is a problem. Anybody could help me?
PDC # net rpc testjoin
get_schannel_session_key: could not fetch trust account password for domain
'DOMAIN_NAME'
net_rpc_join_ok: failed to get schannel session key from server PDC for
domain DOMAIN_NAME. Error was
2015 Mar 04
4
Samba4+ClamAV
Senhores, bom dia!
Tem como instalar o Samba 4 com o antivirus ClamAV? Achei v?rios
artigos, mas n?o consegui fazer funcionar o ClamAV.
Desde, j? agrade?o a aten??o de todos.
Atenciosamente,
Marcio Fiorette