search for: wavelengths

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

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

...nother variable PE (normal distribution). Each analysed separately. - WA has 6 different levels: 400, 430, 470, 540, 570, 620 nm. - the light intensity has 4 different levels: 1, 2, 7, 50. - the same animal was used in several tests. I suspect the response varies quadratically with WA (intermediate wavelengths more efficient to evoke a response than extreme wavelengths) but I would like to know if 470 is different from 540 nm, if 430 is different from 570 nm... I analyse the response with WA as an ordered factor, either with lmer (for PE) or with lme (for LA). All models tested converge but the by-defau...

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

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 Wavelength 1 (in samples): 18.18182 Minimum samples (1): 200 (11.000000.3 wavelengths) Wavelength 1 (in samples): 16.66667 Minimum samples (1): 50 (...

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

...ax2-provision.o [LD] chan_iax2.o iax2-parser.o iax2-provision.o -> chan_iax2.so [CC] chan_local.c -> chan_local.o [LD] chan_local.o -> chan_local.so [LD] gentone.c -> gentone ./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 Wavelength 1 (in samples): 18.18182 Minimum samples (1): 200 (11.000000.3 wavelengths) Wavelength 1 (in samples): 16.66667 Minimum samples (1): 50 (...

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

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

My tuneR sawtooth wave function generator is broken. When I use the sine function, I get exactly what I expect: a sine wave whose frequency is defined by the freq parameter. In particular, higher frequencies have shorter wavelengths (more cycles per second means shorter waves). When I create a sawtooth wave, the opposite seems to occur: higher frequencies result in longer waves. But that?s not all: as frequencies increase, it appears that wavelengths increase to infinite length, then get shorter again as the wave reverses, t...

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

...the code to process my data. My current data structure consists of a .csv file read into R containing 15 columns (a charcoal dilution series going from 100% to 0%) and 1050 rows of absorbance data from 400 nm to 2500 nm at 2 nm interval. I think I need to transpose the data such that the specific wavelengths become my columns and dilutions are defined in rows, but after that point I am lost. Should I (and how do I) make my absorbance data into individual matrices that read into a data frame with only two columns (dilution vs spectra or something of the sort that both the yarn and gasoline data use)? Th...

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

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

...L''. The final x momentum in each case must equal the initial x momentum, since momentum is never lost (it is conserved). Therefore, the final x momenta are equal to each other: P'x + (h sinA) / L' = P''x - (h sinA) / L'' (3) If A is small, then the wavelengths are approximately the same, L' ~ L" ~ L. So we have P''x - P'x = dPx ~ 2h sinA / L (4) Since dx = L/(2 sinA), we obtain a reciprocal relationship between the minimum uncertainty in the measured position,dx, of the electron along the x axis and the uncer...

...L''. The final x momentum in each case must equal the initial x momentum, since momentum is never lost (it is conserved). Therefore, the final x momenta are equal to each other: P'x + (h sinA) / L' = P''x - (h sinA) / L'' (3) If A is small, then the wavelengths are approximately the same, L' ~ L" ~ L. So we have P''x - P'x = dPx ~ 2h sinA / L (4) Since dx = L/(2 sinA), we obtain a reciprocal relationship between the minimum uncertainty in the measured position,dx, of the electron along the x axis and the uncer...

Hello everyone, I have an overlay plot it's nice but you can't see all the data. I would like to know if there is a way to get a plot that gives a side by side plot so that each plot would be next to each other. The two plots have the same data are of different species. At the moment this is the code I'm using: exp<-cbind(abs(round(rnorm(10),2)*10), seq(100, 200, by=10))

...T, Michael Tiemann <mdtiemann at gmail.com> wrote: >My tuneR sawtooth wave function generator is broken. > >When I use the sine function, I get exactly what I expect: a sine wave >whose frequency is defined by the freq parameter. In particular, >higher frequencies have shorter wavelengths (more cycles per second >means shorter waves). When I create a sawtooth wave, the opposite >seems to occur: higher frequencies result in longer waves. But that?s >not all: as frequencies increase, it appears that wavelengths increase >to infinite length, then get shorter again as the...

...email, but i just showed part of my data so you could better understand my problem): ______________ X1242 X1242.5 X1243 X1243.5 X1244 X1244.5 X1245 X1245.5 X1246 X1246.5 X1247 X1247.5 X1248 X1248.5 X1249 => names of each data set ( corresponds to wavelengths, ranges from 400 to 2500 with 0.5 steps => a lot of points!) 18.14 0.9860316 0.9860272 0.9860203 0.9860121 0.9860044 0.9859994 0.9859971 0.9859976 0.9859999 0.9860035 0.9860069 0.9860103 0.9860128 0.9860151 0.9860178 15.8 0.9857134 0.9857106 0.9857063 0.9857011 0.9856958 0.9856917 0.9856887 0.9...