similar to: Predicting and Plotting "hypothetical" values of factors

Hi. Suppose I have a vector that I partition into disjoint, contiguous subvectors. For example, let v = c(1,4,2,6,7,5), partition it into three subvectors, v1 = v[1:3], v2 = v[4], v3 = v[5:6]. I want to find the maximum element of each subvector. In this example, max(v1) is 4, max(v2) is 6, max(v3) is 7. If I knew that the successive subvector maxima would never decrease, as in the example,

Hi Function sequence() repeatedly concatenates its output, and this is slow. It is possible to improve on the performance of sequence by defining myseq <- function(x){unlist(sapply(x,function(i){1:i}))} The following session compares the performance of myseq(), and sequence(), at least on my G5: > identical(sequence(1:50),myseq(1:50)) [1] TRUE > system.time(ignore <-

Hi Function sequence() repeatedly concatenates its output, and this is slow. It is possible to improve on the performance of sequence by defining myseq <- function(x){unlist(sapply(x,function(i){1:i}))} The following session compares the performance of myseq(), and sequence(), at least on my G5: > identical(sequence(1:50),myseq(1:50)) [1] TRUE > system.time(ignore <-

Hi, I have encountered this problem quite a few times and thought I would ask. Let's say that I have two endpoints, a and b, which are integers. If a <= b, I would like to get a:b, but if a > b, then numeric(0), for example: myseq(3,5) => 3:5 myseq(3,3) => 3 myseq(3,2) => numeric(0) The operator : just gives decreasing sequences in the latter case, and I could not coax seq

Hi, I wonder why factor scores produced by factanal are correlated, and I'd appreciate any hints from people that may help me to get a deeper understanding why that's the case. By the way: I'm a psychologist used to SPSS, so that question my sound a little silly to your ears. Here's my minimal example: *********************************************** v1 <-

Dear colleagues, following John Fox' advice in this article (http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-cox-regression.pdf), I'm trying to create a new data frame to examine the differential survival curves from a combination of covariates. These are derived from a Cox Proportional Hazards model I fit to data about the diffusion of a particular policy across American

Dear All, I am using the plotrix library to plot some matrices. I have a problem: some of my data are outliers, hence using a linear color scale does not work very well (you would see too many cells having a similar, indistinguishable color). See the code snipped at the end of the email. Plotting the logarithm of the data gets the job done, but my problem is that I would like to write in every

On 15 Jun 2004 at 13:52, Vito Muggeo wrote: > Dear Petr, > Probably I don't understand exactly what you are looking for. > > However your "plot(x,c(y,z))" suggests a broken-line model for the > response "c(y,x)" versus the variables x. Therefore you could estimate > a segmented model to obtain (different) slope (and breakpoint) > estimates. See

Dear listers, I am trying to fit a model using nlsList() using alternately a SSfol() selfstart function or its developped equivalent formulae. This preliminary trial works well mydata<-groupedData(Conc~Tps|Organ,data=mydata) mymod1<-nls(Conc~SSfol(Dose,Tps,lKe,lKa,lCl),data=mydata) as well as a developped form: mymod2<-nls(Conc~Dose * exp(lKe+lKa-lCl) *

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..

Hi there, I am looking for R-packages that can help me visualize properties on nucleotide sequences. I want to display sequences in the 1-100K base range as lines and plot features above and below those lines. Any ideas would be welcome. Thanks, Bernd [[alternative HTML version deleted]]

Dear R users: > append(1:5, 0:1, after=2) [1] 1 2 0 1 3 4 5 If I want to repeat the appended value every 2 like the following: [1] 1 2 0 1 3 4 0 1 5 How should I modify? Thank you for any help.

I am a new R user. As a test, I want to write a simple code that does the following simulation: 1. Randomly generate a number from a distribution, say, Poisson. Let's say that number is 3. 2. Randomly generate 3 numbers from another distribution, say, Normal. 3. Compute the sum of the numbers generated in step 2 and read it into a vector, V. 4. Repeat steps 1 through 3 for 100,000 times. 5.

Hi all, I would like to factorize the entries in a dataframe given some groupings. E.g: mydf = data.frame( a = rnorm(100,10), b = rnorm(100,10), c = rgamma(100, 1, scale=1)) group = hist(mydf$c, breaks="FD") group$breaks The idea is to create a factor "mydf$d" with levels corresponding to the ranges in group$breaks. There must be an easy way to do this that I

In my database the id for a field is 810. RoR thinks that it is 809. My log file (below) shows that just before the insert it grabs the next sequence value. I dont think it should do that. This causes problems further down the line when trying to access the Person object. Development Log file : "SQL (0.016000) select people_seq.nextval id from dual Person Create (0.015000) INSERT

Hi Peter, You didn't give a very specific example, but it seems to me that what you wish to do is not really complicated. I suppose you have created a table of sequences vs. say hyprophobicity, charge, etc..., something like... seq hydroph arom b0001 0.104762 0.000000 b0002 0.035122 0.065854 b0003 0.024193 0.070968 b0004 -0.096729 0.084112 b0005 -0.973469 0.091837 b0006

Hi all, I've read through much of the list archives and other various docs and am now more confused than ever regarding how samba interacts with external services. All of the docs seem to assume that a Windows AD server is part of the picture when ldap/krb5 are involved. What I'm trying to do is set up a samba fileserver that Windows clients can connect to via something like: "net

First of all, this isn't a troll or a flame, it's a serious question. (Well, OK, two.) Is *ANYONE* known to have *EVER* gotten WINE to work with any Xserver other than XFree86? (AccelX, MetroX?) When was the last time anyone *tested* it on anything other than XFree86? I ask this because, about five years or so ago, I was using the WINE that was current at the time on XFree86 with a

Two * servers: *a and *b. Outside call comes in *b, and is automatically routed to *a. Someone on a sip phone connected to *a then decides to transfer the call to someone on a sip phone connected to *b. The transfer works. At this point, is *a still in the converstation? Or is * smart enough to see where the data stream is going/coming from? Thanks for any help in advanced, and sorry if

Hello everybody. I know this is user forum - but I also know that developers "lurks" in here from time to time. So I would like to ask... what do you (developers) think is the biggest issue in wine development? Is it money? How much would wine development process speed up if Codewavers would start selling... let say, 10 times more licenses for Crossover then they sell now? Would it be a