similar to: random sampling inside a dataset

Dear all, I was wondering how to read DNA sequences in R, is there a specific function and/or a specific package for that? Thank you very much in advance, Gian ** [[alternative HTML version deleted]]

Dear List, I will appreciate any advice regarding how to convert the following numbers [I got in return by taxondive()] in numeric integers without the e.g. 6.4836e+01 abbreviations. Thank you very much in advance, Gian > taxa_dive Species Delta Delta* Lambda+ Delta+ S Delta+ Nat1 5.0000e+00 6.4836e+01 9.5412e+01 6.7753e+02 8.7398e+01 436.99 Nat2

Hi Members, Do exist the possibility to delete a command line into the history? Fo example. If I' ve typed a line code that is wrong, Can I delete it from the hostory in reason to do not save it in the .Rhistory file? Thanks for helping, Gian [[alternative HTML version deleted]]

Hi all, I would love to plot my data with R. I have abundance and frequency of fungal taxonomic data that should be plotted in the same graph. In Microsoft Excel is that possible but the graphic result is, as always, very poor. Is there a function that may let me plot these data in R? I have a matrix made of two columns, on is the relative abundance and the other is the relative frequency for

Hello R community, I have a question regarding correlation and regression analysis. I have two variables, x and y. Both have a standard deviation of 1; thus, correlation and slope from the linear regression (which also must have an intercept of zero) are equal. I want to probe two particular questions: 1) Is the slope significantly different from zero? This should be easy with the lm

Hi all! I've got this error while running example(Glm) library("rms") > example(Glm) Glm> ## Dobson (1990) Page 93: Randomized Controlled Trial : Glm> counts <- c(18,17,15,20,10,20,25,13,12) Glm> outcome <- gl(3,1,9) Glm> treatment <- gl(3,3) Glm> f <- glm(counts ~ outcome + treatment, family=poisson()) Glm> f Call: glm(formula = counts ~

Hi ! I'd like to create a vector that has this kind of numeration 001 002 003 . . . 099 I have looked at format help page but couldn't get any hint on how to do it. Thanks Anna Anna Freni Sterrantino Ph.D Student Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO. [[alternative HTML version deleted]]

Dear all, I have a dataset representing several geographical points (locations), each one having a specific value after performing a PCA. Now I'm trying to understand how to separately represent on a geographic map (of Europe) the pattern of PC1,2, ecc.with colors (e.g.heatmap) I have to add that the dataset includes relatively few points. Therefore, I would also like to interpolate the values

Hello, I have a matrix with the distances among sites. And I have another matrix with the presence and absence of each species in each site. I would like to test the spatial autocorrelation among sites. I have tried to use the function gearymoran of the ade4 package, but error messages keep popping up. Do you know any function for me to test the spatial autocorrelation of my data? Thanks,

Hi, I am writing this in a wrong way, can someone please correct me? > A <- matrix() > length(A) <- 6 > dim(A) <- c(3,2) > colnames(A) <- c("X","Y") > A X Y [1,] NA NA [2,] NA NA [3,] NA NA > A$X Error in A$X : $ operator is invalid for atomic vectors > Thanks, cruz

Hi, I am very new to R and wanted to know if there is a package that, given very long nucleotide sequences, searches and identifies short (7-10nt) motifs.. I would like to look for enrichment of certain motifs in genomic sequences. I tried using MEME (not an R package, I know), but the online version only allows sequences up to MAX 60000 nucleotides, and that's too short for my needs..

Hello! I'm having some trouble  trying to replicate in R a Stata function  invbinomial(n,k,p)        Domain n:     1 to 1e+17        Domain k:     0 to n - 1        Domain p:     0 to 1 (exclusive)        Range:        0 to 1        Description:  returns the inverse of the cumulative binomial; i.e., it                          returns the probability of success on one trial such              

Hi Everyone, I was hoping someone could help me with the settings for Tinn-R. I see in the screen shots that it has syntax help, or something similar (tips on what functions, etc). I can not seem to get this to turn on in the program, and I am wondering if I have to set up a few options. I quickly read through the help and could not figure it out. Many thanks! - Brock P.S. It appears as

Hi, I'm a new user. I've been reading through the manual and looking at various examples but am still trying to make sense of the most efficient ways to handle matrices of data. If I have a 2D matrix of data, how do I get the mean, min, max value of each row? I see the "mean" function on a matrix will give me averages by row, but min and max give me the value for the entire

Hi, Why this function doesn't work? function (x) { if (is.factor(x)) { if (!is.ordered(x)) { warning("La mediana non si puo' calcolare!!!") return(NA) } me <- median(unclass(x)) if (me - floor(me) != 0) { warning("Mediana indeterminata") return(NA) } else

Hi everybody, I am trying to fit a GAM model without intercept using library mgcv. However, the result has nothing to do with the observed data. In fact the predicted points are far from the predicted points obtained from the model with intercept. For example: #First I generate some simulated data: library(mgcv) x<-seq(0,10,length=100) y<-x^2+rnorm(100) #then I fit a gam model with

Hi all, I have some data x, which are actualy consisted of numerical enties. But the class of this matrix is set to be "factor" by someone else. I used "class(x)", it turns out to be "factor". So I can not calculate them. How can I turn them into numerical data so that I can apply math operations on them? Thanks a lot for your help. Selina [[alternative HTML

Hello Chris,  I 've found two other issues  with MAE and CRPS, giving warning  when running examples. I've the same issue on my data.  Hope that you could find some time to take a look here. Thank you Anna > library(ensembleBMA) Loading required package: chron > example(MAE) MAE>   data(ensBMAtest) MAE>   ensMemNames <-

Hi, The call waiting tone is very annoying (you hear nothing while it plays the beep). I need callwaiting because of the queues (the phone has to ring as soon as you hangup) but I want to remove the beep on my dahdi channels, how can I do? Thanks, Niccol? -- http://www.linuxsystems.it

Hi, I have a problem with asterisk 10.8, hints and fxs dahdi phones: if a remote peer and an fxs phone gets connected and the remote peer hangsup, then asterisk sends the "Idle" state to notify the watcher before you hangup the fxs phone! Such a way if the user forgets to hangup the fxs phone (which is a cordless for example) then the operators will keep sending calls to him