similar to: May I ask you a question about matrix population models?

Dear R user: I am studying the allele data of two populations. the following is the data: a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a11 a12 a13 a14 a15 a16 a17 pop1 0.0217 0.0000 0.0109 0.0435 0.0435 0.0000 0.0109 0.0543 0.1739 0.0761 0.1413 0.1522 0.1087 0.0870 0.0435 0.0217 0.0109 pop2 0.0213 0.0213 0.0000 0.0000 0.0000 0.0426 0.1702 0.2128 0.1596 0.1809 0.0957 0.0745 0.0106

I've been trying to learn about dynamically loading subroutines in R. Ultimately I would like to use this approach to implement some simulations. For the moment though, I've run into a problem compiling c++ code and loading it into R. The 'C/C++' and 'R'code below works fine when compiled and linked with gcc, but gives an error when the library is compiled and linked with

Hello, I am trying to do some simulation using the rmetasim package and I've run to this problem. --beginning of error msg-- Error in "[<-"(`*tmp*`, slice[l, ], slice[l, ], value = c(0.200000002980232, : number of items to replace is not a multiple of replacement length --end of error msg-- Here is the script I used. --script starts here-- ## load 'rmetasim'

hello,everyone. I have one question: example 1 > x=numeric(0) > y=5 > print(x+y) numeric(0) example 2 > x=numeric(1) > y=5 > print(x+y) [1] 5 why the print(x+y) is numeric(0) at the first example, but the result is 0 at the second example? __________________________________________________ 佈伵伝仮伱佲伔佈G佊伿佅佷仯伃佒佇伖侜伒佢佉伝伨侙佄佫伬伂伝侙佊伿伡侢伾仹伻伵伋伂伌侒佊伿佅佷

hai, I want to perform a simulation like this仯伜 Suppose that I have one population ,it's size is 500, is composed of x,y and z. The probability of x, y and z is respectively is 0.3, 0.5, 0.2. I wan to simulate a new same size population based ratio of x, y and z, how can I get and assess the number of x, y and z. ? luan Key Laboratory for Sustainable Utilization of Marine Fisheries

> On 26 Aug 2017, at 16:39, niharika singhal <niharikasinghal1990 at gmail.com> wrote: > > Hi, > > Thanks for your mail, and time > > It is not working for some arguments, when mean value is like >6. > > > case > > mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. > 1134463,0.2074027) > > rv

please i've a big problem. i've to do a econometric-quantitative methods assignment about the canadian lynx, the problem is that i really i don't know how to use r and how to apply all the steps. I begun the time plot, ACF and PACF but i'm not able to decide what is the correct model of ARIMA, Holt-winter, ecc to forecast the next 20 years of canadian lynx's cyle... if someone

Hello, <http://tolstoy.newcastle.edu.au/R/help/04/07/0312.html#0313qlink1> Could anybody tell me how to generate discrete Laplacian distribution? I need to sample uma discretised Laplacian density like this: J( g -> g´) ~ exp (-lambda | g´ - g |) g in {0,…, gmax} Thanks, Nicolette [[alternative HTML version deleted]]

I want to th current folder is the specific folder when R start.It is very boring to modify the current directory manually. How can I do it?

I want to speed up computations (involving matrices) by writing some C-code to be loaded. In the C-code, I need to invert matrices etc. As I've understood the "writing R extensions" doc, I can use use #include <R_ext/Linpack.h> in my .c-file and get access to linpack's facilities within my C-code. Is that correct? If so, can anyone point me to examples on how this is

Hello, I am having really hard time finding a good article about simulating genotypes of cases and controls at a disease locus using R. if you guys can point me or guide me where i can find more information, it will be helpful. thanks, Claire -- View this message in context: http://www.nabble.com/genotypes-simulation-tp17065607p17065607.html Sent from the R help mailing list archive at

Hello to everyone! I have built CRAN task view for genetics. For now I have not submit it to CRAN yet and it can be accessible from: http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.html http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.ctv I have not submitted it to CRAN, since I would like first some opinion about it. Genetics is really so broad field that I belive one person

I realized that the script I included in my original post (see below) was written for an older version of rmetasim and will not work with the current version. The only change that need be made to get my script to run with the latest release of rmetasim is to change the command 'simulate.landscape()' to 'sim.landscape()'. My apologies to anyone who tried to use my script and

Hi, Thanks for your mail, and time It is not working for some arguments, when mean value is like >6. case mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. 1134463,0.2074027) rv <-UnivarMixingDistribution(Norm(486.4255, 53.24133), Norm(664.0713, 3.674773), Norm(669.0484, 4.101381),

I'm trying to use a contributed package (rmetasim) to generate simulated genetic datasets. My scripts work fine when I run them on a Sun workstation running Solaris 7 and when I run them on a ~4 year old laptop PC that I have. However, when I run them on my new laptop (Dell Latitude D410 purchased in August 2005), the simulation will run for a short (variable) period of time, then R

---------- Forwarded message ---------- From: niharika singhal <niharikasinghal1990 at gmail.com> Date: Sun, Aug 27, 2017 at 11:57 AM Subject: Re: Find maxima of a function To: "David Winsemius [via R]" <ml+s789695n4745009h56 at n4.nabble.com>, "Ismail SEZEN [via R]" <ml+s789695n4744993h60 at n4.nabble.com>, Ulrik Stervbo <ulrik.stervbo at gmail.com>

Hi R community, I am new bird to R and moved recently from SAS. I am no means expert on either but very curious learner. So your help crucial for me to learn R. I have already got positive expression. I was trying to fit a mixed model in animal experiment but stuck at simple point. The following similar example is from SAS mixed model pp 212. # data genetic_evaluation <-

CRAN (and crantastic) updates this week New packages ------------ * bdoc (1.0) Michael Anderson http://crantastic.org/packages/bdoc This package contains a function that will classify DNA barcodes as well as a few test and reference data sets. * bdsmatrix (1.0) Terry Therneau http://crantastic.org/packages/bdsmatrix This is a special case of sparse matrices, used by coxme and

Hello all - I'm having a bit of a problem that I'm sure is being caused by my missing some trivial detail. But I haven't been able to find it, and I'm not even sure how I would construct the search to find relevant info in the archives. Here is the setup. I have Samba 3.0.20a running as a PDC against a LDAP back end. For the most part everything works fine. Users that are

Hola LLVMers, I was curious about the state of stack alignment on x86. I noticed there are a few bugs outstanding on the issue. I recently added some code which had the effect of throwing an extra function parameter on our stack at runtime, a 4 byte pointer. Esp is now not 16-byte aligned, so instructions like unpcklps xmm1, dword ptr [eps] cause grief. My AllocaInstr instructions are