similar to: normalised Voigt random numbers

Hi, I am sorry for this question, but I am trying to speed up an application.... I will need to fit many x-y data sets (input from text files) to 4-parameter Pseudo-Voigt peak function. Until now I used SigmaPlot macro to do it (enclosed just in case...) peaksign(q) = if(total(q)>q[1], 1, -1) xatymin(q,r) = xatymax(q,max(r)-r) [Parameters] a = if(peaksign(y)>0, max(y), min(y))

First, thanks to all who helped me with my question about rescaling axes on the fly. Using unlist() and range() to set the axis ranges in advance worked well. I've since plotted about 300 datasets with relative ease. Now I'm trying to fit a lossy oscillator resonance to (the square root of) a lorentzian (testframe$y is oscillator amplitude, testframe$x is drive frequency): lorentz

Dear all, I am using the nlsLM function to fit a Lorentzian function to my experimental data. The LM algorithm should allow to specify limits, but the upper limit appears not to work as expected in my code. The parameter 'w', which is peak width at half maximuim always hits the upper limit if the limit is specified. I would expect the value to be in-between the upper and lower limit with

Dear expeRts, first of all I'd like to thank you for the quick help on my last which() problem. Here is another one I could not tackle: I have data on an absorption measurement which I want to fit with an voigt profile: fn.1 <- function(p){ for (i1 in ilong){ ff <- f[i1] ex[i1] <- exp(S*n*L*voigt(u,v,ff,p[1],p[2],p[3])[[1]]) } sum((t-ex)^2) } out <-

Hi, how may I define in a callfile the CallerID presentation to be used for the requested call, eg. set it to prohibited? TIA, Bruno -------------- next part -------------- A non-text attachment was scrubbed... Name: Bruno.Voigt.vcf Type: text/x-vcard Size: 270 bytes Desc: not available Url : http://lists.digium.com/pipermail/asterisk-users/attachments/20050913/fcb5c595/Bruno.Voigt.vcf

Hi all, I'm looking for a way with any asterisk-version with TE410P (cpe EuroISDN, Q931) for sending an INFORMATION ELEMENT KeypadFacility, eg. *87, during a connected call to the PSTN switch. Are there existing functions in asterisk to generate & send such IE ? If not what existing modules would be best to derive from? TIA, Bruno -------------- next part -------------- A non-text

Dear List, Consider the following example x=c(1,2,3,4,5) y=c(2,4,6,8,10) linearmodel=lm(y~x) To predict a y-value if you know the corresponding x value is very easy with the command predict. predict(linearmodel, newdata=(x=1.5)) The other way around, to predict an x-value with a corresponding y-value, doesn't work unfortunately. Is there another function that can do that, or do I need to

Hi, I'm trying to use Robust non-linear regression to fit dose response curves. Maybe I didnt look good enough, but I dind't find robust methods for NON linear regression implemented in R. A method that looked good to me but is unfortunately not (yet) implemented in R is described in http://www.graphpad.com/articles/RobustNonlinearRegression_files/frame.htm

Dear all, I freshly installed a new smb-server that should act as file server with some smb-shares. The pdc is a Windows 2008 r2 server. My smb.conf is as follows: [global] workgroup = agroup server string = ssmbserver #passdb backend = tdbsam map to guest = Bad User usershare allow guests = No security = domain wins support = No

Hi, I have a problem with the filesize of the samba 3.0.0 package. after compiling the source of the samba-3.0.0 i made the "sh makepkg.sh". the package i got was about 901 MB. does anybody know why the package is as great as mine and what can i do to get a samba package which is smaller than the package i have. i compiled under solaris 8. i made all what samba.org says that i

>> Do you (or the ogg vorbis community) understand how SBR works? There is >> info on http://www.codingtechnologies.de/technology/sbr.htm >> bnut I am guessing this infop is not enough to explain the technology. Robert Voigt: >I haven't heard about SBR before. After reading that webpage I can say the >following: [ . . . ] I don't think SBR will give an

Dear Sir/Madam, I could succesfully normalise my microarray data using marrayNorm package. However, i have not been able to get normalised red and green channel intensities through R package. Is there a possibility to write a formula to calculate back the red and green channel intensities after normalisation of the data. Do I need to incorporate this formula in my R script? I am biologist

Dear all, I want to see if the treatment of an animal with a specific compound has an effect on the expression of certain genes. Though my question is based in biology, it really is all about how to deal with the standard deviation in normalised data. I have three groups of animals; untreated, treated with placebo, and treated with a single concentration of the compound in question. Gene

Dear List, My data consist of nine columns and about 50,000 rows. It looks like this. -9.0225 3.46464 2.80926 -0.3847 3.73735 1.1058 -2.98936 1.38901 -8.1846 -2.4315 -5.1189 1.8225 3.3798 1.7874 4.693 -3.9286 1.4266 5.7849 -3.4894 -4.0305 3.7879 3.5195 2.9186 2.8685 -6.126 4.978 4.9381 4.5282 3.62558 -3.0455 4.6518 1.39746 0.68652 3.5708 -3.6404 -4.2963 -1.3183 0.6752 -4.0382 -2.5386

Hi, I have the following R-code: pdf(file="test.pdf") hist(rcauchy(100)) dev.off() How can I set the size of the histogram in my pdf? I tried to do with par(...) but I didn't find the right parameter yet. Thanks, Max -- View this message in context: http://www.nabble.com/set-size-of-a-plot-tp22477083p22477083.html Sent from the R help mailing list archive at Nabble.com.

Dear R Users, Are there any packages in R which carries out a normalisation to variables as follows: - find the empirical distribution function, using perhaps ecdf - use the empirical distribution function to transform the variables into a series between 0 and 1 - use this series to map the variables into the normal distribution function, using qnorm - perform a regression on the transformed

Hi R-users, Recently, I use maxNR function to find maximizer. I have error appears as follows Error in maxNRCompute(fn = fn, grad = grad, hess = hess, start = start, : NA in the initial gradient My code is mu=2 s=1 n=300 library(maxLik) set.seed(1004) x<-rcauchy(n,mu,s) loglik<-function(mu) { log(prod(dcauchy(x,mu,s))) } maxNR(loglik,start=median(x))$estimate Does anyone know how

I quote from "Principles of Digital Audio" by Ken C. Pohlmann: "Because perceptual coders tailor the coded signal to the ear's acuity, they similarly tailor the required response of the playback system itself. Live music does not pass through amplifiers and loudspeakers, it goes directly to the ear. But recorded music must pass through the playback signal chain. Much of the

Hello! My Linux-server is every day attacked with brute-force password cracking attacks. I use openssh-3.9p1 (SuSE Linux 9.2) with standard setup (PAM, LoginGraceTime 2m, MaxAuthTries 6). Unfortunately, I see cracking attempts with very short delays (1 second): Jan 31 00:46:53 XXX sshd[10774]: Invalid user backup from ::ffff:66.98.176.50 Jan 31 00:46:54 XXX sshd[10776]: Invalid user server

Oooh Plate Echo was what people had before digital reverbs! (When they didn't have spring reverb!) - Long ago... They are quite rare now because only big studios had them as they took up a lot of room! (and then they threw them all away) *jealous* Do you really have a plate reverb or just a digital simulation? love Freya --- Robert Voigt <robert.voigt@gmx.de> wrote: > uh, you