similar to: Geneland error on unix: Error in MCMC(........ :, unused argument(s) (ploidy = 2, genotypes = geno)

I do the these passages: library(Geneland) set.seed(1) data <- simdata(nindiv=200, coord.lim=c(0,1,0,1) , number.nuclei=5 , allele.numbers=rep(10,20), IBD=FALSE, npop=2, give.tess.grid=FALSE) geno <- data$genotypes coord <- t(data$coord.indiv) path.mcmc <-

I am getting the above warning following loading of Geneland 3.1.5 on unix , while a simple plot sends output to the pdf file ( see attached code) no output results from Geneland functions, resulting in empty pdf files > library (Geneland) Loading required package: RandomFields Loading required package: fields Loading required package: spam Package 'spam' is loaded. Spam version

Hi, I have a package of my own which seems to work fine under linux. I want to make a compiled version for windows. (I work with windows 2000 and R 2.0) I have followed the steps described in the file readme.packages (in the top-level directory of the binary installation) and I had a trouble at this step (which aim is not documented): C:\Program Files\R\rw2000\src\gnuwin32>make libR.a

Hi , I'm currently trying to make available a few fortran subroutines and R functions (which make interface to these subroutines), as an R package. I'm doing it under linux with R 1.9.0 (but hope to do it for windows too). I have trouble for loading my fortran code. Here are the steps involved : 1) In R: # Create tree for package

I am trying to use the function 'filter.NA=TRUE' in Geneland. The function appears to be set on TRUE by default, as it appears as TRUE in the 'parameter.txt' file output and hence I do not need to enter the function per se (as it is an 'Unused argument otherwise') . Hence all my missing data (individuals that I have not yet scored at that specific loci) are scored as

Hi all, I would like to multithread that script, to detect structure from multilocus genetic data : >library(Geneland) > >geno = read.table("cot966gen_test.txt") #the file is show after >MCMC(geno.dip.codom = geno, varnpop=T, npopmax=20, spatial = F, nit=100000, thinnin=100, path.mcmc="./") >PostProcessChain(path.mcmc="./", nxdom=100, nydom=100,

Hello, anybody know the procedure to plot the geneland result in a map? thank you ********************************************************** Vincenzo Landi Post Doctorate student Animal genomic and breeding cell:0039/3395388713 Fax. 075-5857122 [[alternative HTML version deleted]]

I am trying to format data for the Geneland package I need to write a series of paired numbers and paired "000" to a plain ascii file without surrounding the numerals with quotes, the original data is in a matrix formatted paired strings, and written to file using write.table(x,"output filename",col.names=F,row.names=F) thus "77 79" "132 132" "000

I am trying to install R on a SLED 10.1 machine. R-base-2.7.0-7.1-i586.rpm fails with stas at linux-6b8s:~/RPMs> rpm -Uvh R-base-2.7.0-7.1.i586.rpm warning: R-base-2.7.0-7.1.i586.rpm: Header V3 DSA signature: NOKEY, key ID 14ec5930 error: Failed dependencies: libgfortran.so.1 is needed by R-base-2.7.0-7.1.i586 I tried to trick it into believing there's the library by setting up

Hello all, I have two matrices, pop and pop2, each the same number of rows and columns that I want to compare for equality. I am concerned about efficiency in this operation. I've tried a few things without success so far. Doing something simple like: if (pop==pop2) { cat('equal') } else { cat('NOT equal') } results in the warning: 1: In if (pop == pop2) { : the

Suddenly WINE have started to bring a Select Script File window whenever I try to start application with its parameters from a shortcut. Image:[Image: http://img523.imageshack.us/img523/4844/wineselectscript.jpg ] This is the full copypaste from a terminal when I try to execute the line that the shortcut would too: Code: epri at epri-ubu32-desk:~$ wine "/home/epri/.wine/drive_c/Program

Is there any way to echo comments from an R source file into an SWeave->LaTeX document? Example: # Npop is population total # Npoph0..Npoph2 are stratum totals # Npoph is vector of stratum totals Npop<-sum(to2000) Npoph0<-sum(to2000[bg==0]) Npoph1<-sum(to2000[bg==1]) Npoph2<-sum(to2000[bg==2]) Npoph<-c(Npoph0,Npoph1,Npoph2) In the final LaTeX document, I'd like the

Hello all, I have the following function call to create a matrix of POP_SIZE rows and fill it with bit strings of size LEN: pop=create_pop_2(POP_SIZE, LEN) I have 3 questions: (1) If I did keep_pop[1:POP_SIZE] == pop[1:POP_SIZE] to keep a copy of the original data structure before manipulating 'pop' potentially, would this make a deep copy or just shallow? Ie

Hi all, Thanks to all who offered advice on using F95 in R. Now I'm trying to compile a test package using gfortran, Linux 2.4.21 and R 2.1.0. I was able to successfully compile and use a test F95 routine by setting my environment variables as follows in bash: export PATH=~/bin/:$PATH export F77=gfortran export LD_LIBRARY_PATH=~/bin/irun/lib export GFORTRAN_STDIN_UNIT=-1 Now I'm

Hello there! I have to run Samba with ads under AIX 5.3. I have krb5 running, and openldap is running too. I can configure and make samba 3.0.2.0 with the following configure options: ./configure --prefix=/usr/local/samba --with-ads --with-winbind --with-included-popt --with-aio-support ./configure --prefix=/usr/local/samba --with-winbind --with-included-popt --with-aio-support If irun the

can somebody that has these hardwares(X101P and S400P) working on his asterisk system please assist.................. you can send the solution to austino@skannet.com..... the error message is what i have below. ---------- Forwarded message ---------- Date: Thu, 1 May 2003 21:01:21 +0100 (WAT) From: austino@skannet.com To: asterisk-user@lists.digium.com Subject: Error working with X101P and

Dear all I'm developing a grafical interface for Geneland that is presently one of the R packages available. My first question is about refreshing the interface while running a very long process. run<-function(){ ? ? ? ? tttry <- tktoplevel() ? ? ? ? tkwm.geometry(tttry, "+200+200") ? ? ? ? tkwm.title(tttry,"wait") ? ? ? ?

Hello, I am working with a matrix of multilocus genotypes for ~180 individual snail samples, with substantial missing data. I am trying to calculate the pairwise genetic distance between individuals using the stats package 'dist' function, using euclidean distance. I took a subset of this dataset (3 samples x 3 loci) to test how euclidean distance is calculated: 3x3 subset used

Hi, Just find a scaling factor that would make the two sets of data comparable. Here I divided the second row by 5 and did the same for the second axis. Charles-?douard F1 <- as.table(matrix(c(50,11,6,17,16,3,1,2237,611,403,240,280,0,0), 2,7)) barplot(F1, beside = TRUE, col = c("blue", "grey")) axis(2, at=c(0,10,20,30,40,50,60, labels=c(0,10,20,30,40,50,60))) axis(4, at =

Dear Charles-Edouard Thanks a lot. Yes indeed barplot sounds excellent. Unfortunately, the scale of the smaller axis is fixed, even If I am able to draw to axes. The idea is to expand the scale to the scale to the second axis for comparison. F1 <- as.table(matrix(c(50,11,6,17,16,3,1,2237,611,403,240,280,0,0), 2,7)) barplot(F1, beside = TRUE, col = c("blue", "grey"))