similar to: fit a deterministic function to observed data

Hi everyone, I'm trying to perform a bi exponential Fit with the package NLS. the plinear algorithm seems to be a good choice see: p<-3000 q<-1000 a<--0.03 b<--0.02 t<-seq(0:144);t y<-p*exp(a*t) + q*exp(b*t)+rnorm(t,sd=0.3*(p* exp(a*t) + q*exp(b*t))) fittA <- nls(y~cbind(exp(a*t), exp(b*t)), algorithm="plinear",start=list(a=-.1, b=-0.2), data=list(y=y, t=t),

Hi, In R, I know the method table(), which builds a contingency table of the counts y at each level for the factor x. But I would like to know what is the inverse method of table(), if it exists, to obtain the vector z, from the two vectors x and y? x <- (86, 90, 94, 98, 102, 106, 110, 114, 118, 122, 126, 130, 134, 138, 142, 146, 150, 154, 158, 162, 166, 170, 174) y <- c(2, 5, 10, 17,

Hi all, I would like to fit a gamma pdf to my data using the method of RSS (Residual Sum-of-Squares). Here are the data: x <- c(86, 90, 94, 98, 102, 106, 110, 114, 118, 122, 126, 130, 134, 138, 142, 146, 150, 154, 158, 162, 166, 170, 174) y <- c(2, 5, 10, 17, 26, 60, 94, 128, 137, 128, 77, 68, 65, 60, 51, 26, 17, 9, 5, 2, 3, 7, 3) I have typed the following code, using nls method:

I'm using a rather strange config, borne out of neccessity rather than choice at home. My internet 'router' is a Win2000 Pro box running Winroute, and my three Linux boxen (running 7.2/8.0/9.0 RH) are networked thru to the Win box using SNAT on the box I work on (don't ask why - it's just pratical, and I can't afford a hub/switch to do this). In any case, it's a

Hi, I have a problem fitting the following Weibull Model to a set of data. The model is this one: a-b*exp(-c*x^d) If I fitted the model with CurveExpert I can find a very nice set of coefficients which create a curve very close to my data, but when I use the nls.lm function in R I can't obtain the same result. My data are these: X Y 15 13 50 13 75 9 90 4 With the commercial

This data are kilojoules of energy that are consumed in starving fish over a time period (Days). The KJ reach a lower asymptote and level off and I would like to use a non-linear plot to show this leveling off. The data are noisy and the sample sizes not the largest. I have tried selfstarting weibull curves and tried the following, both end with errors. Days<-c(12, 12, 12, 12, 22, 22, 22,

I have a RHEL5 Server and some dual-boot XP/CentOS 5 systems (Linux systems all 64-bit). All Linux is out-of-box, with all packages, minus international languages, installed. No patching has been done. On the server, I selected system-config-authentication and enabled LDAP for User Information, Kerberos, LDAP, and SMB for Authentication, and Shadow and MD5 Passwords, along with

Dear R specialists, I'm working with time series and want to investigate the relationship between two time series by correlation analysis or by fitting a gen. lineair model to the plot of x(timeserie1) and y(timeserie2). Lin1 <- data.frame( Nr = c(1:lengte), NDII = window(ts.mNDII,c(1998,10),c(2003,11)), InvERC = window(Inv.ERC,c(1998,10),c(2003,11)) )

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

I am using Microsoft .NET (2.0) PrintPreviewDialog to preview reports, and all works well when the Windows default printer is set to a local printer or a network printer on a Windows server, but if I set it to a Samba network printer, it slows down to a crawl. With a Windows server network printer, it takes about 10 seconds to render 100 pages. With a Samba network printer, it takes over 5

hi! when i have busy signal on analog line (zap card) it doesn't detect that line is busy ? is it possible to change detected sequence (frequency) of busy tones on line (zapata.conf ??) tnx, Thomas my zapata.conf [channels] language=en context=lin1 signalling=fxs_ks channel => 1 group=1 echocancel=yes echocancelwhenbridged=yes rxgain=3.0 txgain=3.0 busydetect=yes

Dear R users, When running the program below I receive the following error message: fit <- optim(parm, objective, yt = tyield, hessian = TRUE) Error in as.vector(data) : no method for coercing this S4 class to a vector I can't figure out what the problem is exactly. I imagine that it has something to do with "tyield" being a matrix. Any help on explaining what's going on

Dear useRs, A few days ago, the initial version 0.1.0 of a new package 'random' was uploaded to CRAN. The random packages provides convenient access to the non-deterministic random numbers provided by the random.org site created by Mads Haahr (http://www.random.org). While certain hardware and software solutions that provide access to non-deterministic random-numbers exist, few if any

Dear useRs, A few days ago, the initial version 0.1.0 of a new package 'random' was uploaded to CRAN. The random packages provides convenient access to the non-deterministic random numbers provided by the random.org site created by Mads Haahr (http://www.random.org). While certain hardware and software solutions that provide access to non-deterministic random-numbers exist, few if any

Guys, It seems as though the llvm system doesn't deterministically iterate over Module::iterator, or global_iterator. To make myself clearer, the iterators iterate over all the global_variables but on different llvm passes (different calls to opt -load), the iterators iterate over them in different orders. I was thinking that it has something non deterministic to do with byte code reading or

I mean SVN head, but it might have been fixed in 2.6 as well, I don't remember. On Jan 25, 2010, at 2:47 PM, Augustine Mathews wrote: > Thanks chris. > > Do you mean llvm 2.6 when you say mainline or something else? > > Augustine > > On Mon, Jan 25, 2010 at 4:27 PM, Chris Lattner <clattner at apple.com> > wrote: > > On Jan 26, 2010, at 7:23 AM,

On 12.11.2010 15:26, Erik Cederstrand wrote: > I have noticed that two consecutive builds of clang, clang++ and tblgen don't produce identical binaries (as in md5 sums) on identical source code (I'm on FreeBSD). I ran strings(1) on the two clang binaries, and I get the following: > 251862c251862 > <

Oh, I should point out that using randomness is the unfortunate part, but you can get deterministic builds by passing -frandom-seed=0 to clang and/or gcc. -Chris On Nov 12, 2010, at 9:19 AM, Chris Lattner wrote: > Yes, this is controlled by the -frandom-seed flag. Sadly, it is expected behavior. > > -Chris > > On Nov 12, 2010, at 7:10 AM, Sebastian Redl <sebastian.redl at

On Nov 12, 2010, at 7:10 AM, Sebastian Redl wrote: > On 12.11.2010 15:26, Erik Cederstrand wrote: >> I have noticed that two consecutive builds of clang, clang++ and tblgen don't produce identical binaries (as in md5 sums) on identical source code (I'm on FreeBSD). I ran strings(1) on the two clang binaries, and I get the following: >> 251862c251862 >> <

Den 12/11/2010 kl. 18.48 skrev Chris Lattner: > Oh, I should point out that using randomness is the unfortunate part, but you can get deterministic builds by passing -frandom-seed=0 to clang and/or gcc. I'll try that. I assume that the randomness has a purpose, so what are the negative implications of setting -frandom-seed=0? Will i still have a functioning compiler? :-) Thanks, Erik