similar to: nls with plinear and function on RHS

I am trying to understand what rbinom function does. Here is some sample code. Are both the invocations of bfunc effectively doing the same or I am missing the point? Thanks, Pieter bfunc <- function(n1,p1,sims) { c<-rbinom(sims,n1,p1) c } a=c() b=c() p1=.5 for (i in 1:10000){ a[i]=bfunc(30,p1,1) } b=bfunc(30,p1,10000)

I am looking for some help at removing low-frequency components from a signal, through Moving Average on a sliding window. I understand thiis is a smoothing procedure that I never done in my life before .. sigh. I searched R archives and found "rollmean", "MovingAverages {TTR}", "SymmetricMA". None of the above mantioned functions seems to accept the smoothing

Hi all, I have a very large binary vector, I wish to calculate the number of 1's over sliding windows. this is my very slow function slide<-function(seq,window){ n<-length(seq)-window tot<-c() tot[1]<-sum(seq[1:window]) for (i in 2:n) { tot[i]<- tot[i-1]-seq[i-1]+seq[i] } return(tot) } this works well for for reasonably sized vectors. Does

Hi All. I've run into a problem with the plinear algorithm in nls that is confusing me. Assume the following reaction time data over 15 trials for a single unit. Trials are coded from 0-14 so that the intercept represents reaction time in the first trial. trl RT 0 1132.0 1 630.5 2 1371.5 3 704.0 4 488.5 5 575.5 6 613.0 7 824.5 8 509.0 9

Dear sirs It seems like there is a bug in `profile.nls` with `algorithm = "plinear"` when a matrix is supplied on the right hand side. Here is the bug and a potential fix ##### # example where profile.nls does not work with `plinear` but does with # `default` require(graphics) set.seed(1) DNase1 <- subset(DNase, Run == 1) x <- rnorm(nrow(DNase1)) f1 <- nls(density ~ b1/(1 +

I was wondering if there is way to place constraints upon the "plinear" algorithm of nls, or rather is there a manner in which this can be achieved because nls does not allow this to be done. I only want to place constraints on one of the nonlinear parameters, a, such that it is between 0 and 1. I have attempted to use a=pnorm(a*) , but then the fitting procedure becomes

I want to fit a nonlinear model of the form: Y=A+B*X1+C*X2+log(X3/(X3+D)) I think that the best way is to use the plinear algorithm, but I don't know how to specify the formula in the nls function. I've tried: Y~cbind(rep(1,times=length(Y)),X1,X2,log(X3/(X3+D))) But this fits the model: Y=A+B*X1+C*X2+D1*log(X3/(X3+D)) How can I specify the formula correctly? Thanks, M

Hello, i have a problem with the function nls(). This are my data in "k": V1 V2 [1,] 0 0.367 [2,] 85 0.296 [3,] 122 0.260 [4,] 192 0.244 [5,] 275 0.175 [6,] 421 0.140 [7,] 603 0.093 [8,] 831 0.068 [9,] 1140 0.043 With the nls()-function i want to fit following formula whereas a,b, and c are variables: y~1/(a*x^2+b*x+c) With the standardalgorithm

hello all i havnt had a chance to read through the references provided for the "nls" function (since the libraries are closed now). can anyone shed some light on how the "plinear" algorithm works? also, how are the fitted values obtained? also, WHAT DOES THE ".lin" below REPRESENT? thanking you in advance ###################################### i have a quick

Hello sir, You have interpreted correctly that clustering will be done by generating the ring around the Document(i.e. the basic idea of LSI). But it is not like increasing the radius and the next shell will be another cluster, Rather it would pick one document (based on relevance score) and form a ring around it to cluster the document, then from the remaining documents(not in the cluster but

On 24 June 2014 01:55, Philip Reames <listmail at philipreames.com> wrote: > > On 06/19/2014 09:35 AM, Yi Kong wrote: >> >> Hi all, >> >> I'm currently working on implementing ACLE extensions for ARM. There >> are some memory barrier intrinsics, i.e.__dsb and __isb that require >> the compiler not to reorder instructions around their

HI all, Somebody had tried to do a rebuild of RedHat Storage Server 2.0 and RHEV 3.2 hypervisor with CentOS? Mock is an option to accomplish this? Thanks.

Dear all, I was wondering if it is possible in TableGen to either: 1. Selectively define an instruction depending on an SDNode's properties, e.g. if the SDNode is not commutative. 2. Override/ignore the TableGen error given when a commutative node has an immediate on the LHS. My case comes from trying to define a generic ALU operation multiclass for my target, which includes a

Paul.White@dsto.defence.gov.au wrote: > > > Hello bug people, > > > > Version: > > platform = i686-pc-linux-gnu > > arch = i686 > > os = linux-gnu > > system = i686, linux-gnu > > status = > > major = 1 > > minor = 3.1 > > year = 2001 > > month = 08 > > day = 31 > > language = R > > > >

(1) There are no lists here! You have numeric vectors, not lists. (2) The normal vector-recycling rules apply. It does 1 == 4, 2 == 5, 3 == 4, 4 == 5, 5 == 4, .... Try > 1:10 == 5:4 [1] FALSE FALSE FALSE TRUE TRUE FALSE FALSE FALSE FALSE FALSE (3) For what you appear to want > x %in% 4:5 [1] FALSE FALSE FALSE TRUE TRUE FALSE FALSE FALSE FALSE FALSE is the way to get it.

Feature request: I want to use update.formula to subtract an intercept (or other) term from a formula with a dot on the RHS. However, as this causes an error, I propose a patch below. Thus, I want: > update.formula(y ~ ., ~ . -1) [1] y ~ . - 1 Instead I get this error: Error in terms.formula(tmp, simplify = TRUE) : '.' in formula and no 'data' argument While the error

I'm using the following sequence to plot a scatter plot to PostScript. Those familiar with the Iris LDA example in MASS will recognise what I'm at. > postscript("hulda.eps", horizontal=FALSE, onefile=TRUE, height=6, width=6, pointsize=8, paper="special") > plot(hu.ld, type = "n", xlab= "first linear discriminant", ylab="second linear

Dear Colleagues, I would like to only inspect rules that contain a certain label substring on the rhs. In this special case the item labels are built like this: <itemtype>_<itemvalue> e.g. "Artikelgruppe_E0815" what I want to do is only show rules where "Artikelgruppe" is contained in the rhs - has anybody an idea how this could work? Sincerely

Hi Folks, I'm trying to cut my data inside the summaryBy function. Perhaps formulas don't work that way? I'd like to avoid adding another column if possible, but if I have to, I have to. Any ideas? Thanks, Allie require(doBy) df = dataframe(a <- rnorm(100), b <-rnorm(100)) summaryBy(a ~ cut(b,c(-100,-1,1,100)), data=df) # preferred solution, but it throws an

Does anyone know if there is a way to specify regressors dynamically rather than explicitly? More specifically, I have a data set in "long format" that details a number of individuals and their responses to a question (which can be positive, negative, or no answer). Each individual answers as many questions as they want, so there are a different number of rows per individual. For each