similar to: advice on using this C code in R

Hi there, I try to plot visible light spectrum (380nm~780nm) with color corresponding to the specific wavelength. However, I don't find a function that could do this. Another question, it's possible to plot a color space chromaticity diagram like this: http://upload.wikimedia.org/wikipedia/commons/thumb/0/02/CIExy1931.svg/300px-CIExy1931.svg.png Thanks in advance! Jinsong

Dear all, I wrote a wrapper to a FORTRAN program using R. The main program uses a text file (~200 lines) as an input describing the simulation to be run. I typically generate the file once with the right parameters using a combination of file(), paste(), cat(). This is fine, and it works well, however I then need to update only a few values in the file many times (~200 times,

Hi I'm trying to write a generic script for processing some data which finishes off with some plots. Given Im never sure how many columns will be in my dataframe I wanted to using the following plot(spectra.wavelength, cormat, type = "l", ylim=c(-1,1), xlab="Wavelength (nm)", ylab="Correlation") however even if I specify as type="l" it appears plot

Hello Im trying to install an old version of Asterisk. But it isnt working: when I run "make install": gcc -o gentone gentone.c -lm ./gentone busy 480 620 Wavelength 1 (in samples): 16.66667 Minimum samples (1): 50 (3.000000.3 wavelengths) Wavelength 1 (in samples): 12.90323 Minimum samples (1): 400 (31.000000.3 wavelengths) Need 400 samples Wrote busy.h ./gentone ringtone 440 480

Hi, I have some difficulty in figuring out whether I am doing correct or not. A brief introduction about the work: It is a light/dark choice test conducted in insect larvae.  The response is binary (0- present in dark area, 1-present in light area) and the experiment is run for 15 min, so there are 15 repeated measurements per individual larva at 1 min intervals.  The factors which affect

Hi, am building the latest version of Asterisk (1.4.13) on a self-build Linux host (based on LFS-6.3). Version 1.4.12 built installed and worked fine. Last night I upgraded the kernel to 2.6.23 and rebuilt the zaptel driver package 1.4.5.1 against it. That seemed to build and install O.K. too. I dl'd the 1.4.13 source tarball and tried to build that: ./configure ran O.k. 'make

Dear R users, I am a newbie. Just switched from MATLAB. So thanks a lot for your patience. I have 50000 spectra collected in field. Each spectra has two columns : Wavelength (56) and the actual measurement. Each measurement came in a different .txt file on disk (50000 files in total). I wrote a script that reads every spectra in a for loop and constructs two variables : Wavelength (56) and

First time ever that I try to call subroutines in a Win DLL using R. Have done this before using VBA and Python. The C code's function argument list contains only double pointers (double *x). The function is declared void, the output value is one of the arguments. C calling sequence is used. When inspecting the exports of the DLL with dumpbin I can see that the function I need is

I have a spectrophotometric dataset with repeated measures of a value at 200 wavelengths for each of 150 individuals. I would like to use the repeated samples to at each wavelength to look at measurement/observer error, compared to difference between individuals error I have looked at doing this with icc{irr} or using an anova approach, but I am unclear how to acheive this given that there

please reply to hdgbyi@public.guangzhou.gd.cn or bgpgong@hotmail.com, thank you. UNCERTAINTY PRINCIPLE IS UNTENABLE By reanalysing the experiment of Heisenberg Gamma-Ray Microscope and one of ideal experiment from which uncertainty principle is derived , it is found that actually uncertainty principle can not be obtained from these two ideal experiments . And it is found that

please reply to hdgbyi@public.guangzhou.gd.cn or bcpgong@hotmail.com, thank you. UNCERTAINTY PRINCIPLE IS UNTENABLE By reanalysing the experiment of Heisenberg Gamma-Ray Microscope and one of ideal experiment from which uncertainty principle is derived , it is found that actually uncertainty principle can not be obtained from these two ideal experiments . And it is found that

Apologies if this is already answered elsewhere, but I could not find the answer searching the archives. I am trying to bring my data in from an .xlsx file. The data of interest is an n x m array, but it is preceded with ~2 columns and 30 rows of instrumental "information." An abbreviated example (cols denoted with "," for this example: Subject , 01 wavelength, 400 (30

Hi, I am currently analysis Raman spectroscopic data with the hyperSpec package. I consulted the documentation on this package and I found an example work-flow dedicated to Raman spectroscopy (see the address : http://hyperspec.r-forge.r-project.org/chondro.pdf) I am currently trying to remove outliers thanks to PCA just as they did in the documentation, but I get a message error I can't

I see about 7 different R packages for multi-process parallel programming. Which do you think is the best, most complete, and most robust to pick for general purpose Erlang-style message-passing programming in R, and why? First here's my use case, and then my analysis so far. I often have code whose basic organization looks something like this: 1. Fetch step: For each date, gather up or

Dear lme and lmer -ers, I have some problems using "home-made" contrast matrix in lme and lmer. I did an experiment to investigate the relationship between the response of an animal and some factors, namely the light wavelength (WA), the light intensity to which this animal was exposed and the sex of the animal tested. - The response can be a variable LA (normal distribution) or

I am not sure if I should ask this question in this list. But I'll try. Currently I am trying to analyze images using EBImage and biOps. One of the features that I need to extract from various images is the color spectrum, namely, which colors each image consists of. So, each image hopefully can be converted into some sort of color histogram so that color ingredients are easily comparable

Dear R users, I have spectral data (say, wavelength vs. extinction coefficient) for which I??d like to calculate an integral (i.e. the area underneath the curve). Suppose the (artificial) dataset is lambda E 1 2 2 4 3 5 4 8 5 1 6 5 7 4 8 9 9 8 10 2 How can I calculate an integral for these values using R? Many thanks for any help! Regards Christoph

I am currently work on a research project which requires parallel technique. I have successfully adapted an Rmpi code from the following website. http://math.acadiau.ca/ACMMaC/Rmpi/task_pull.R and be able to run my own code now. However, since I am developing my own code which involves many calculations, each time I had an error with my Rmpi code, I got an error message like the

Hi, I'm new to R and looking for a way to use a variable to reference to an object. I'm plotting several graphs on top of each other and want to do this by a for loop. The field I want to graph have names like a1,a2,a3, ... I can't figure out how to get this working: It should look like this: lines(spectral$Wavelength,paste("spectral$a",j,sep=""),col=j) but

Hi All, I'm trying to parallelize some code using Rmpi and I've started with a sample 'hello world' program that's available at http://math.acadiau.ca/ACMMaC/Rmpi/sample.html. The code is as follows; # Load the R MPI package if it is not already loaded. if (!is.loaded("mpi_initialize")) { library("Rmpi") } # Spawn as many slaves as possible