search for: nimh

Dear Dr. Sheth, This spring, Dr. Cottler and the EPRG submitted a renewal for the T32 NIMH Training Program in Psychiatric Epidemiology and Biostatistics. We are submitting a letter to Council for the competing continuation for NIMH which better showcases the accomplishments of our past trainees. It is vital that you provide us with an update that includes grants received since the comp...

I want to run deconvolution of a time series by an impulse or point-spread function through Wiener filter, regularized filter, Lucy-Richardson method, or any other approaches. I searched the CRAN website and the mailing list archive, but could not find any package for such a deconvolution analysis. Does anybody know an existing R function for deconvolution? TIA, Gang

Hello, I'm simply trying to import a .txt file that has a column of about 30,000 pts. I found the file, but I'm getting an error: > m <- read.table(choose.files()) Error in "row.names<-.data.frame"(`*tmp*`, value = row.names) : duplicate row.names are not allowed Any help with getting around this? I really appreciate all the help. Thanks dave

Hi I am struggling with nested random effects and hope someone can help. I have individuals (ID) who are nested within families (FAM). I want to model an outcome variable, and take account of the intercorrelation of individuals within each family. I think this amounts to two random effects, one nested within the other. How can I model this in R? So far I have tried using the

Hello, I'm new to R, and I'm having trouble importing a text file (I'm on Windows XP) > m <- read.table("/Desktop/work/128_L") Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `/Desktop/work/128_L' do you know why this isn't working? All I have is a bunch of text files, each with a

...ar. In this plot, each boxplot is labeled by the corresponding value of var, but the absissas are evenly spaced. Is it possible to generate a boxplot with the absissas of each boxplot equal to the corresponding value of var? -- Alan Barnett, PhD Imaging Physicist National Institutes of Health NIMH/CBDB 301 402 3507

Hi, Suppose I have a multidimensional array: tmp <- array(1:8, c(2,2,2)) is there a function out there that, given a one-dimensional array index, will return the separate indices for each array dimension? for instance, tmp[8] is equivalent to tmp[2,2,2]. I'd like to derive the vector (2,2,2) from the index 8. thanks, Brad Buchsbaum

...;ve also looked at the Aikake and Bayes information criteria for this, which give different results from one another. (Assuming that what I should look for is just a decrease to mean better prediction.) Any advice and/or references appreciated. Thanks, Matt Wiener Laboratory of Neuropsychology NIMH/NIH Bethesda, MD 20892 301-496-5625 x254 mcw at ln.nimh.nih.gov -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in th...

Hello, I'm trying to do a cluster analysis on a large data set. I tried it out with a smaller one first, but I got this error: > hc<-hclust(dist(x),"ave") Error: cannot allocate vector of size 4129151 Kb The data sample used (i.e. "x") is a numerical data set of size 32513 by 31 Does anyone know how I can do this analysis? Is R capable of this

...st. of Mental Health Computer Scientist National Institutes of Health Laboratory of Neuropsychology Bldg. 49, Room 1B-80 Neural Coding and Computation Bethesda, MD USA 20892 tel: 301-496-5625 ext. 251 http://neuron.nimh.nih.gov <><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>

The following commands illustrate a problem with graphing relatively large data sets using the "line" option to plot: temp <- runif(200000) # or whatever other numbers you like; # same thing happens with a sine wave plot(1:200000, temp) # everything fine plot(1:200000, temp, type="l") # data gets cut off The number of points isn't crucial. The cutoff point

I have several newbie questions - any help much appreciated: 1. I run a combination Windows 98 and XP home network - have been looking at Linux on two boxes at home - one desktop - one laptop - RH 7.3 2 I can get samba to run okay and in each of the two Win98SE machines can see them as localhost 3. I cannot access any resources on the Linux machines - Win98 asks for a password - \\localhost\IPC$

Hello, I'm attempting to run a PCA on an example data set. I ran it just fine, but I don't know how to few the output? I listed what the variable got stored in it, but I don't know how I can get anything else out of it. Are there other ways to view the results? Also, I'm confused about the last line "6 variables and 8 observations" Aren't the rows the

...;ve also looked at the Aikake and Bayes information criteria for this, which give different results from one another. (Assuming that what I should look for is just a decrease to mean better prediction.) Any advice and/or references appreciated. Thanks, Matt Wiener Laboratory of Neuropsychology NIMH/NIH Bethesda, MD 20892 301-496-5625 x254 mcw at ln.nimh.nih.gov -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in...

...mbership model." Right now, I would like to keep this simple, using only a continuous dependent variable, but would like to also extend this to a repeated measures design. This doesn't seem to be possible with the lme package. Is there something else I could consider? Thanks, Brian NIMH Training Fellow GWB School of Social Work, PhD Program Washington University in St. Louis One Brookings Drive St. Louis, MO 63130 [[alternative HTML version deleted]]

Bradley Taylor shot me this *goldmine* of ugly ugly ugly dogs: http://www.sonoma-marinfair.org/uglydogvote.shtml I *have* to use one of these for the Mongrel BUGS Mascot. He''ll go on our bug list page and replace the little beetles on the left. Pick the dog you think best represents a lovable but defective pooch and reply to this with your +1. Let the voting begin! -- Zed A. Shaw

Hello All First - a warning. I'm not very R or programming savvy. I am trying to do something without much luck, and have scoured help-pages, but nothing has come up. Here it is: I have a matrix (m) of approx 40,000 rows and 3 columns, filled with numbers. I would like to convert the contents of this matrix into another matrix (m_p), where the numbers of (m) have been coerced into a

...e 0 or 1 Directory: Documentation * PR# 323 * Subject: { and } in examples get confused in help From: Prof Brian Ripley <ripley@stats.ox.ac.uk> Date: Sun, 14 Nov 1999 22:51:43 GMT Directory: Graphics * PR# 120 * Subject: graphics bug: type="l" From: Matthew Wiener <mcw@ln.nimh.nih.gov> Date: Tue, 16 Feb 1999 15:30:25 -0500 (EST) --X server problem drawing > 64K lines at once. Workaround? * PR# 149 * Subject: plot.formula and pch= From: Peter Dalgaard BSA <p.dalgaard@biostat.ku.dk> Date: 25 Mar 1999 21:22:44 +0100 --plot(y~x,pch=z) may subset x and y but leave...

Matthew Wiener <mcw@ln.nimh.nih.gov> writes: > t1 <- data.frame(matrix(rnorm(16), nc=4)) > t1$X1 <- 1 > t1$X2 <- 2 > print(t1) > Error: dim<- length of dims do not match the length of object Well, it is prototype-compatible. Splus 5.3 does likewise. A way out is t1<-data.frame(unclass(t1)...

Hi, all. Using the patched version of R-0.99.0, I cannot reproduce the following example from the strsplit documentation: (Example) unlist(strsplit("a.b.c", ".")) ## [1] "" "" "" "" "" ## Note that `split' is a regexp! ## If you really want to split on `.', use unlist(strsplit("a.b.c", "\."))