search for: stim

I am trying to fit a nonlinear model using nlm(). My application of nlm() is a bit complicated. Here is the story behind the model being fit: The observer is trying to detect a signal corrupted by noise. On each trial, the observer gets stim=signal+rnorm(). In the simulation below I have 500 trials. Each row of stim is a new trial. On each trial, if the cross-correlation between the stim and the signal is above some criterion level (crit=.5 here), the observer says "signal" (resp=1), else he says "no signal" (resp=...

1)I am wandering how the following SQL statement can be written in R language w/o using sqldf: create table detail2 as select a.* from detail a, pdetail b where a.TDATE=b.TDATE and (a.STIM >= b.STIM and a.STIM <=b.MAXTIM) 2) when try if then in R it only applies to the 1st row & not to whole dataset like in SAS. How do you get round that? in SAS: data summary; set all1; if entry='a:prop' then pctexec=stkful/stocks*100; run; Thanks in advance for your help....

...e may or may not get an answer on this list, I decided to try it anyway, hoping not to incur in somebody's ire for misusing the tool. For the sake of clarity, I will explain the problem. We conducted an experiment measuring average response times in a cognitive task. The task has two types of stimuli (stim1, stim2) and was performed in a 6-run session, where subjects performed the task under condition A on odd-numbered runs and under condition B on even-numbered runs. Thus, the temporal sequence of the runs was the following: - run 1: cond A - run 2: cond B - run 3: cond A - run 4: cond B...

Hi I'm reading a CSV file using read.table, and it keeps importing a text variable as a factor. To overcome this, I've used the as.is command referring to the variable in question (called "stim") data<-read.table(file.choose(), header=T, sep=",", as.is = "stim") However, "stim" is still imported as a factor. I notice there are other read.table options related to this issue, but the help files don't appear to contain much detail on these. Any he...

Hi, I would have thought that these two constructions would produce the same result but they do not. Resp <- rbinom(10, 1, 0.5) Stim <- rep(0:1, 5) mm <- model.matrix(~ Stim) Xb <- mm %*% c(0, 1) ifelse(Resp, log(pnorm(Xb)), log(1 - pnorm(Xb))) pnorm(as.vector(Xb), lower.tail = Resp, log.p = TRUE) > ifelse(Resp, log(pnorm(Xb)), log(1 - pnorm(Xb))) [1] -0.6931472 -1.8410216 -0.6931472 -0.1727538 -0.6931472 [6] -0...

Hi I'm an R newbie. I can't seem to add new variables to data frames. Here are the stages (1) I import the data using read.csv. (2) I fix it using fix(data) (3) I create a new variable using spos<-tagPOS(stim,language="en",model=NULL,tagdict=NULL). (tagPOS is a function in the OpenNLP toolkit, which tags a string for part of speech. "stim" is a variable in the current data frame) What happens is that the new variable, "spos" is saved as an object, but not added to the dat...

Dear all, I am estimating a mixed-model in Ubuntu Raring (13.04¸ amd64), with the code: fm0 <- lme(rt ~ run + group * stim * cond, random=list( subj=pdSymm(~ 1 + run), subj=pdSymm(~ 0 + stim)), data=mydat1) When I check the approximate varianc...

...use neither dist() nor as.dist() seem to be able to work with data in this format. I could try brute-force, but maybe there is an easier way. Can someone help? As a distance matrix, the column helm.raw$CD1 for subject CD1 should appear something like shown below (without the Obs column, where stim is the rowname) --------------------------------- Subject=CD1 ---------------------------------- Obs stim RPur Red Yel Gy1 Gy2 Green Blue BlP Pur1 Pur2 1 RPur . . . . . . . . . . 2 Red 11.5 . . . ....

...the beta diversity index value between two grids. Now I would like to get the average value of each GridCell. Please kindly advise how to make the calculation. Thank you. Elaine The matrix looks like (cited from Michael Friendly) I would like to get the average value of each color. Obs stim RPur Red Yel Gy1 Gy2 Green Blue BlP Pur1 Pur2 1 RPur . . . . . . . . . . 2 Red 11.5 . . . . . . . . . 3 Yel 13.1 6.0 . . ....

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

I will have data in the following form: Time resp type stim type 300 a A 200 b A 155 a B 250 b B 80 c A 1000 d B ... c is left censored observation; d is right censored This sort of problem is discussed in Chap 9 of Cox & Oakes Analysis of Survi...

...ter:ALLOWED/0/FULL when the file or directory is copied to Samba the ACLs become ACL:FESFOO\root:ALLOWED/0/RW ACL:FESFOO\InfoCenter:ALLOWED/0/RW ACL:\Everyone:ALLOWED/0/ I can manually go in and fix these via a windows GUI but when I try to use smbcacls smbcacls -d 3 //spiderman-new/infocenter stim.zip -Uadministrator -M ACL:FESFOO/InfoCenter:ALLOWED/0/FULL I get Connecting to 192.168.242.129 at port 445 Doing spnego session setup (blob length=58) got OID=1 3 6 1 4 1 311 2 2 10 got principal=NONE Got challenge flags: Got NTLMSSP neg_flags=0x60890215 NTLMSSP: Set final flags: Got NTLMSSP ne...

Here is a programming question. The code I am using is quite slow and I was wondering if there is a way to get rid of the for loop. I am dealing with "interaction" in 2x2 table, and am using Edwards's G_I (Likelihood, p. 194). I label the cells in the table as follows stim response "y" "n" total -------------------------------- y hit miss nsignal -------------------------------- n false correct nnoise alarm rejection --------------------------------- Gsens<-function(nhit,nmiss,nfa,ncr) { # Edwards's G_I, measure of interaction, i...

ok as everyone knows, (and to be frank i didnt think it would in my lifetime ;P) fallout 3 gets released in 3 days or 6 for those who live in europe. now for those who where playing games since the 3rd grade the fallout universe has a special place in their hearts since literally i grew with it! so in mho fallout 3 fells in the category of games that someone must simply have, and thats why im

Hi, I am running a Cox Mixed Effects Hazard model using the library coxme. I am trying to model time to onset (age_sym1) of thought problems (e.g. hearing voices) (sym1). As I have siblings in my dataset, I have decided to account for this by including a random effect for family (famid). My covariate of interest is Mother's diagnosis where a 0 is bipolar, 1 is control, and 2 is major

Hi I am running Mac OX tiger with the installed ruby 1.8.2, and mysql5 (I checked that it was running in preference). I tried to run the various fixes recommended in agileweb. I am stuck on page 57, script/generate scaffold... > Could you paste in your database.yml file, from the > "config" folder inside your Rails app? If you have any > passwords entered there, you should