similar to: R package for visualizing/analyzing accelerometry data?

Hello - I am attempting to add interaction terms to a model using the gamm function in the mgcv package. I stripped it down to only two terms to simplify it. I have no trouble with individual terms, continuous or factors, nor with interacting factor terms; I cannot, however, figure out how to add interactions between continuous*factor terms. I applied a smoother to my continuous term because it

For a paper dealing with generalized ellipsoids, I want to illustrate in 3D an ellipsoid that is unbounded in one dimension, having the shape of an infinite cylinder along, say, z, but whose cross-section in (x,y) is an ellipse, say, given by the 2x2 matrix cov(x,y). I've looked at rgl:::cylinder3d, but don't see any way to make it accomplish this. Does anyone have any ideas? thx,

Not sure if the question is appropos, but I have multiple csv's which are read into an xts object, corresponding to telemetry data (accelerometer, magnetometer/compass, and gyroscope). For examination, it would be ?ber useful if plot.zoo (or something similar) allowed me to zoom in and out of a subset of the time axes. -- Sent from my mobile device Envoyait de mon portable

I have a vibration signal coming accelerometer. I converted this signal from* m/s^2* to *mm/s*. Now I am supposed to predict this vibration signal in R using historical data. (Please see the attached picture of vibration signal). Can I use this vibration signal just like that to the prediction using predictive Models or Do i need to do some processing technique before doing prediction? Are there

I have an irregular time series in a Zoo object, and I've been unable to find any way to do an FFT on it. More precisely, I'd like to do an NFFT (non-equispaced / non-uniform time FFT) on the data. The data is timestamped samples from a cheap self-logging accelerometer. The data is weakly regular, with the following characteristics: - short gaps every ~20ms - large gaps every ~200ms

Dear List, I'm trying to draw vector in XYZ with rgl under use of a cylinder3d. Therefore I scale and rotate a basis-cylinder). However, somehow the rotation is wrong as verified by overplotting arrow3d(). Where is my mistake? library(heplots) library(rgl) # ... 2 vectors data=data.frame(row.names=c('X','Y','Z'), x1=c(2,1,5),y=c(4,3,2))

Hello everyone, I wanted to copy a function from one file to another. The file that I wanted to copy the function into contains a function with the same name and same number of instructions. I decided to just replace the instructions with those of the other function. I am doing all of this from within a function pass. *1.* I created a function pass in which I extract a function using

Hi, Fitting all possible models (GLM) with 10 predictors will result in loads of (2^10 - 1) models. I want to do that in order to get the importance of variables (having an unbalanced variable design) by summing the up the AIC-weights of models including the same variable, for every variable separately. It's time consuming and annoying to define all possible models by hand. Is there a

Since 2008, Microsoft staff and guests have written about R at the Revolutions blog (http://blog.revolutionanalytics.com) and every month I post a summary of articles from the previous month of particular interest to readers of r-help. In case you missed them, here are some articles related to R from the month of February: The R Consortium opens a new round of grant applications for R-related

Dear useRs, Is there a way to draw helix spring in R of desired diameter and number of turns? regards, Eliza [[alternative HTML version deleted]]

Hi!! I need to draw a cylinder/tube joining two points. I am trying to make something presentable, I have been able to do it using lines3d. But is it possible to increase the thickness of the lines? The size parameter does not work. Does any one have any suggestions? Thanks ../Murli library(rgl) pts<-structure(list(x = c(-0.975688, -0.975688), y = c(9.258795, -9.258795 ), z = c(-1.8, 1.8)),

Dear r users, I have two kinds of data I'm trying to represent in Rgl. The first are measurements of soils: heights, diameters, and centerpoints of each. I would typically represent these in 3d space as cylinders or cones. The second are radar data of the area under the solids , xyz plus amplitude that I would like to represent as a point cloud. Is there a currently implemented method to

Hello I want to batch job the calculation of many GLM-models, extract some values and store them in a file. Almost everything in the script below works (read file, extract values and write them to file) except I fail in indexing the GLM with the modelstructure it should run. Running GLM's conventionally is no problem. Conventionally a GLM is calculated as:

Buenos días. Tengo dos mapas que hago con ggplot (abajo) y necesito superponer el 2º al 1º. En el 1º se representa la variable PA, binomial (0, 1), en gris, y en el 2º la variable ExtTG7085, que es numérica y va de 0 a 200. La 2ª es un gradiente en rojo. Lo ideal sería que donde coinciden el rojo de la 2ª se viera más oscuro, pero supongo que eso no es fácil. Podría también ser en otro color, pero

Hello Instead of testing against 0 i would like to test regression slopes against -1. Any idea if there's an R script (package?) available. Thanks for any hint. Cheers Lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220

heya, i am fitting linear mixed effect model to a response Y. Y shows an s-shaped distribution when using QQ-plots (some zero values and some very high values). hence, which transformation should i apply that Y follows a normal distribution? any r-function/package available to do this? thanks for any hint, regards, lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of

Hi, I would like to evaluate the frequency of the variables within the best selected model by AIC among a set of 12 competing models (I fit them with GLM) with a bootstrap procedure to get unbiased results. So I would ike to do the ranking of the 12-model-set 10'000 times separately and calculate the frequency of variables of the 10'000 best ranked models. I wrote a script doing the model

Dear all, I am a new user of R. I would like to know how to get fund's returns in percentage (%). For example, I use: R <- ts(read.xls("FundData"), frequency = 12, start = c(1996, 1)) Whith this program, the returns are like 0.0152699. But, I would like to have 1.52%. Please advise me about the function. Thanks! Fabrice

display, accelerometer/motion sensor being the first 4 for release (though 81 have been mocked up so far). The long term concept is if you want a 'weather station with live video feeds and gps location control you can add various modules together to deliver what you are looking to achieve. I have high hopes for the concepts, and wish the guys well as it seems their hearts are in the right

Buenos días. Hago mis mapas con el código que os copio abajo, pero me gustaría que me dibujase las líneas del mapa después de poner los puntos para que se vean. Puedo hacer los puntos más pequeños, pero entonces quedan separados y yo quiero que rellenen el mapa. Gracias, como siempre, Manuel print(ggplot(legend=FALSE)+geom_path(data=map_data('world'), aes(x=long,