Displaying 20 results from an estimated 9000 matches similar to: "writing string values from a matrix to a file without enclosing quotes"
2010 May 20
1
Geneland error on unix: Error in MCMC(........ :, unused argument(s) (ploidy = 2, genotypes = geno)
I am receiving the above error ( full r session output below) the
script runs OK in windows. and "genotypes" and "ploidy" are both
correct arguments
any suggestions would be most welcome
Nevil Amos
MERG/ACB
Monash University School of Biological Sciences
> library(Geneland)
Loading required package: RandomFields
Loading required package: fields
Loading required
2010 May 29
1
warning In fun(...) : no DISPLAY variable so Tk is not available
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
2012 Jan 21
2
replacing "+" in string
I am trying to replace "+" in a string with another character
I am getting odd results using sub and gsub
> X<-"one + two"
> gsub("+","plus",X)
[1] "plusoplusnpluseplus plus+plus plustpluswplusoplus"
> sub("+","plus",X)
[1] "plusone + two"
> X<-"one ~ two"
it seems to work fine with other
2010 May 04
7
How to replace all <NA> values in a data.frame with another ( not 0) value
I need to replace <NA> occurrences in multiple columns in a data.frame
with "000/000"
how do I achieve this?
Thanks
Nevil Amos
2012 Apr 10
7
How to remove $ (Dollar sign) from string
How do I remove a "$" character from a string sub() and gsub() with "$" or
"\$" as pattern do not work.
> sub("$","","ABC$DEF")
[1] "ABC$DEF"
> sub("\$","","ABC$DEF")
Error: '\$' is an unrecognized escape in character string starting "\$"
>
2011 Oct 16
2
How to plot CI's (llim ulim) on ecodist mgram
I would like to put confidence intervals on a mantel corellogram
they are already calculated in the pmgram object but I am unsure how I
get the x value in order to plot them?
package(ecodist)
X<-1:100
Y<-rnorm(1:100)
Z<-rnorm(1:100)
XY<-dist(data.frame(X,Y))
YX<-dist(data.frame(Y,X))
my.mgram<-mgram(XY,XZ)
plot(my.mgram)
print(my.mgram)
> print(my.mgram)
$mgram
2012 May 20
5
removeing only rows/columns with "na" value from square ( symmetrical ) matrix.
I have some square matrices with na values in corresponding rows and
columns.
M<-matrix(1:2,10,10)
M[6,1:2]<-NA
M[10,9]<-NA
M<-as.matrix(as.dist(M))
print (M)
1 2 3 4 5 6 7 8 9 10
1 0 2 1 2 1 NA 1 2 1 2
2 2 0 1 2 1 NA 1 2 1 2
3 1 1 0 2 1 2 1 2 1 2
4 2 2 2 0 1 2 1 2 1 2
5 1 1 1 1 0 2 1 2 1 2
6 NA NA 2 2 2 0 1 2 1 2
7 1 1 1 1 1 1 0 2 1 2
8
2013 Aug 27
1
return the name of source.
Is there a fuction that will allow me to retrun the filename for a script
from within that script.
fir instance
If I have a script "myscript.r":
FileName<-unknown.fucntion()
print(FileName)
and run it
source("myscript.r")
will return
"myscript.r"
Thanks
Nevil Amos
[[alternative HTML version deleted]]
2012 Jan 11
2
Vegan(ordistep) error: Error in if (aod[1, 5] <= Pin) { : missing value where TRUE/FALSE needed
I am getting the following erro rmessage in ordistep. I have a number of
similarly structured datasets using ordistep in a loop, and the message
only occurs for some of the datasets.
I cannot include a reproducible sample - the specific datasets where this
is occur ing are fairly large and there are several pcnm's in the rhs of
the formula.
thanks for any pointers that may allow me to
2020 Oct 31
2
raster::levels() not working in packaged function.
Apologies, I cannot see how to make a rero for this issue.
I have a function that uses levels(r) tor return the RAT of a raster "r"
when the function is sourced from a script
source(".\R\function.r")
it works fine.
when the function is built into a package and sourced from there
library(mypackage) using the same script file to make the package
levels(r)[[1]]
the same line
2010 Oct 12
2
repeatability/intraclass with nested levels
I have a spectrophotometric dataset with repeated measures of a value
at 200 wavelengths for each of 150 individuals.
I would like to use the repeated samples to at each wavelength to look
at measurement/observer error, compared to difference between
individuals error
I have looked at doing this with icc{irr} or using an anova approach,
but I am unclear how to acheive this given that there
2011 Dec 03
1
partial mantel tests in ecodist with intential NA values.
I would like to perform partial mantel tests on only within group values, with "between group" values assigned to NA.
This is possible in package ncf partial.mantel.test, however this sues a different permutation to that used in ecodist.ecodist will not accept data with NA values, returning a "matrix is not square error.
is it possible to perform this test in ecodist?
many thanks
2010 Dec 17
1
Change legend position in s.value {ade4}
How do I change the postion of the legend in s.value {ade4} from the
defaul , bootom left?
thanks
Nevil Amos
2011 Sep 12
1
how to get xlab and ylab in bold?
A very basic query
This code plots OK the axis values are in bold but the axis labels are
not. how do I get them in bold too?
thanks
Nevil Amos
plot(c(1,1),xlim=c(0,450),ylim=c(0.7,1.4),xlab="Distance (cells) from
edge of grid",ylab="Resistance distance",
type="l",col="white",lwd=2,font=2,family='sans')
2010 Jun 01
1
Printing page nmbers in multi-page pdf
Is it possible to print page numbers in pdf() with multiple pages?
thanks
Nevil Amos
2010 Apr 29
1
using get and paste in a loop to return objects for object names listed a strings
I am trying to create a heap of boxplots, by looping though a series of
factors and variables in a large data.frame suing paste to constrcut the
facto and response names from the colnames
I thought I could do this using get()
however it is not working what am I doing wrong?
thanks
Nevil Amos
> sp.codes=levels(data.all$CODE_LETTERS)
>
> for(spp in sp.codes) {
>
>
>
2010 Sep 21
1
partial dbRDA or CCA with two distance objects in Vegan.
I am trying to use the cca/rda/capscale functions in vegan to analyse
genetic distance data ( provided as a dist object calculated using
dist.genpop in package adegenet) with geographic distance partialled out
( provided as a distance object using dist function in veganthis method
is attempting to follow the method used by Geffen et al 2004 as
suggested by Legendre and . FORTIN (2010).
I
2004 Oct 25
5
building a package under windows
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
2004 Oct 22
1
building an R package : where and how should my fortran library be loaded ?
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
2010 Apr 29
5
reduce size of pdf
is there a way to reduce the size of pdf files in R: ?
compression?
lower dpi ?
or some other option?