search for: stimulation

...reatment groups in the parameter of the geometric distribution, but with a cut-off (i.e. right censored). In my experiment I stimulated animals to see whether I got a response, and stopped stimulating if the animal responded OR if I had stimulated 10 times. Since the response could only be to a stimulation, the distribution of response times is discrete (so I can't use survival analysis). I am interested in assessing at what stimulation number the animals responded, and whether this differed between groups. Assuming that the animals have a fixed probability (p) of responding to each stimulat...

Hi all, I am currently doing a consistency analysis using an lmer model and trying to use parametric bootstrapping for the confidence intervals. My model is like this: model<-lmer(y~A+B+(1|C/D)+(1|E),binomial) where E is the individual level for consistency analysis, A-D are other fixed and random effects that I have to control for. Following Nakagawa and Scheilzeth I can work out the

Hi, I have tried to use uniroot to solve a value (value a in my function) that gives f=0, and I repeat this process for 10000 times(stimulations). However error occures from the 4625th stimulation - Error in uniroot(f, c(0, 2), maxiter = 1000, tol = 0.001) : f() values at end points not of opposite sign I have also tried interval of (lower=min(U), upper=max(U)) and it won't work as well. Can anyone help me as I have struggled for...

The exceptional Discovery of Human Phero mones Our perfume has real human sexual attractant in it. Direct from the laboratory to you! After I saw that story on Dateline NBC about your perfume I ordered some and I am "extremely" satisfied. Women are now coming up and talking to ME! I don't even have to use pick-up lines anymore. Thank you. - Robert Y. in Tuscon, AZ. To Whom

Often, when exploring a dataset, I'd like to plot several very different Y variables against the same X variable, in panels stacked one over the other. Is there an easy way to do this? I'd like to achieve an elegant look similar to the look achieved by lattice in conditioned plots--for instance no space between panels. But unlike in straightforward conditioned plot, each panel may be on a

The (*windows only*) xlsReadWrite package has been updated and the new version is available on CRAN. xlsReadWrite natively reads and writes Excel files (v97 - 2003). o Changes from v1.1.1 to v1.3.3 (versions between never really published) - ROWNAMES-support also for matrices (not only for data.frame) - explicit rowNames argument in interface. (stimulated by email Adrian Dusa) -

Greg, I'm going to join issue with your here! Not that I'll go near advocating "Excel-style" graphics (abominable, and the Patrick Burns URL which you cite is remarkable in its restraint). Also, I'm aware that this is potential flame-war territory -- again, I want to avoid that too. However, this is the second time you have intervened on this theme (previously Mon 6

...2,1)/10) x=runif(n) if(x<=1/10){y[i]=1} else if(x<=2/10){y[i]=2} else if(x<=4/10){y[i]=3} else if(x<=7/10){y[i]=4} else if(x<=9/10){y[i]=5} else{y[i]=6} } bar<-barplot(table(y)) table<-table(y) return(list(bar,table)) } mydie2() but I also want to try to create the same dice stimulation using ifelse() statement, but all hell break loose when I attempt to use the command below (the barplot only shows the results to be 6), can anyone tell me what went wrong please: mydie2<-function(n=1000,y=NULL,...){ for(i in 1:n){ x<-sample(1:6,n,replace=TRUE,prob=c(1,1,2,3,2,1)/10) x=runif(...

Gordon Henriksen wrote: > On 2007-09-15, at 23:55, Talin wrote: > >> Gordon Henriksen wrote: >> >>> Can you be more specific the algorithm for which you need type >>> metadata in a write barrier? No algorithms I am aware of perform any >>> tracing from a write barrier. >> >> This one does: >> >>

These are some points stimulated by reading about C history (and related in their implementation). 1) On some platforms > as.integer("0xA") [1] 10 but not all (not on Solaris nor Windows). We do not define what is allowed, and rely on the OS's implementation of strtod (yes, not strtol). It seems that glibc does allow hex: C99 mandates it but C89 seems not to allow it. I

...al success, and hopefully, when the wells are dewatered and fully evaluated, it will be shown to be an economic success as well. Other CBM areas have seen similar success. For example, Barrett Resources experienced a 1.5 fold increase in gas production from 14 CBM wells in the Raton Basin following stimulation with coiled tubing versus conventional treatments. "Fracturing is probably the area with the greatest growth potential," says Rike, "and coalbed methane is about a third of that market." lepromata banishers NEW COIL TUBING MARKETS / FIELD APPLICATIONS CT first established it'...

After attending my first useR! conference I want to thank the organizers for doing a wonderful job and the presenters for their high quality presentations and stimulating ideas. The conference venue was excellent and of course Vienna is one of the greatest cities in the world to visit. useR! is one of the most fun conferences I've attended. Thanks again! -- Frank E Harrell Jr

Hi list, I use a .xml file for a function's demo in the R package I'm creating. Since it doesn't belong to any of the 'standard' folders, i.e. those mentioned in the 'Writing R Extension', I put it in a folder call "myXML", much like the 'iris.xl' file in 'xls' folder from 'gdata' package, for example. After running R CMD build,

...nt when we try to convert a .wav file to an .opus file. What settings in Opus that can accomplish this? Sincerely, Vinson Go cid:image001.png at 01C94EE1.4DA09EF0 TBSI - SIMPLICITY in Neuro-Solutions <http://www.tbsi.biz/> Tethered and Wireless Solutions for Neural Recording and Stimulation Vinson Go Director of Research & Development EMail: <mailto:dfan at trianglebiosystems.com> vgo at trianglebiosystems.com Website: <http://www.trianglebiosystems.com/> www.trianglebiosystems.com Triangle BioSystems International a division of Harvard Bioscience &l...

Dear R list, I would like to perform an ANOVA in a set of measurements, but I have problems formatting the data. The data is a two dimensional array containing two columns: - "Stim" : the type of stimulation (string) - "Ratio" : a ratio of two numeric values Now, because some values are missing in the data (defaulting to zero), part of this array will be populated with NA ratios. Maybe this is important later. In order to make the ANOVA analysis, I need to turn my vector into a data.frame....

Breaking News - Powell resigns as FCC Chair Federal Communications Commission (FCC) Chairman Michael Powell announced his resignation today (January 21, 2005) and will step down in March. In a brief statement on the FCC's website (www.fcc.gov) Powell said, "Having completed a bold and aggressive agenda, it is time for me to pursue other opportunities and let someone else take the reins

CENTER FOR VETERINARY BIOLOGICS UNITED STATES DEPARTMENT OF AGRICULTURE Join our team of statisticians and data scientists in 2018 We bring a broad perspective to our consulting and analytical work with scientists at the USDA National Centers for Animal Health. Our Statistics Section is looking for a good communicator who is comfortable dealing with the analytical aspects of practical

Dear R-Listers, to measure the psychometric curve of pitch discrimination, one sequentially presents two tones of slightly different pitch to an observer (animal will do), and asks "which is higher". The pschometric curve is the fraction of correct responses plotted against the pitch difference. It passes through 50% (pure guessing) at zero and normally approaches 100% at large

...Typically bursts occur in 60-80% of the heart cycles * sensory stimuli can modify these burst probabilities. I give 48 stimuli-trials at random intervals and count the number of bursts associated with the stimuli. For example, a person with 75% burst probability at rest (e.g., 36/48) may have an stimulation induced increase to 87.5% (42 bursts in 48 trials). There are 14 subjects in each of 3 different patient groups. Simple enough. But what if the stimulus reactions are modified over time? The surprise of the stimulus (electric shock) soon wears off and the responses (e.g., increased burst probab...

Folks: In answer to a query, Andy Liaw recently submitted some code to construct a block diagonal matrix. For what seemed a fairly straightforward task, the code seemed a little "overweight" to me (that's an American stock analyst's term, btw), so I came up with a slightly cleaner version (with help from Andy): bdiag<-function(...){ mlist<-list(...) ## handle case in