similar to: Still confused with lapply

Hello, I have a simple xyplot with rlm lines. I would like to add the a and b coefficients (y=ax+b) of the rlm calculation in each panel. I know I can do it 'outside' the xyplot command but I would like to do all at the same time. I found some posts with the same question, but no answer. Is it impossible ? Thanks in advance for your help. Ptit Bleu. x11(15,12) xyplot(df1$col2 ~

Hello, I'm trying to fit some points with a 8-degrees polynom (result of lm is stored in pfit). In most of the case, it is ok but for some others, some coefficients are "NA". I don't really understand the meaning of these "NA". And the problem is that I can't perform a derivation (pderiv<-as.function((deriv(polynomial(pfit$coefficients))))) on pfit due to the

Hello, In one of my functions, I need to extract the hour from a date. For example: as.POSIXlt("2010-03-27 02:00:00")$hour gives 2 as expected. It works for all the dates I've tested except the following one which is in my list of dates: as.POSIXlt("2010-03-28 02:00:00")$hour which gives 0. I don't understand why it doesn't work with this particular date. Any

Hello, I'm still trying to make the life of my colleagues easier. Nice, isn't it ? At the moment, I'm looking for a way to generate multiple "single report". In fact I have a .rnw file which send a query to a MySQL database (rs<-dbSendQuery(con, statement="select * from treatdata where name='Device1'") But of course my colleagues have many devices and

Hi, I recently discovered the R program and I thought it could be useful to me. I have to analyse data saved as .Px file (x between 0 and 8 - .P0 files have 18 lines at the beginning that I have to skip). New files are generated everyday. This is my strategy : In order to analyse the data, I first want to copy the new data in a database in MySQL (which already contains the previous data). So

Hello, I have a text file with one column containing long number but stored as string. I download the file with read.table (and colClass) and the first row of this column is : "095842087016731010" As I need to make some calculations with these numbers, I tried to convert them using as.numeric. But then I get as.numeric("095842087016731010") 95842087016731008 I understand

Hello, I would like to plot the following xyplot : for each date of fff (1 date per panel), bbb=f(aaa) for the two groups (ddd=1 and ddd=2) superimposed. I can do it by group (see below) but not together. I looked at http://lmdvr.r-forge.r-project.org/figures/figures.html but I haven't found what I was looking for (to be honest, I haven't understood all the examples). Any help will be

Hello, I would like to know if it exists a package including something equivalent to this page: http://www.developpez.net/forums/d740403/autres-langages/algorithmes/contribuez/image-geometrie-projective-homography/ http://www.developpez.net/forums/d740403/autres-langages/algorithmes/contribuez/image-geometrie-projective-homography/ I've searched with these keywords: R-project homography

Hi, I have a data frame of 2 columns with the following types : data$day char data$value num And I plot my data with : plot(strptime(donnees$day,format="%Y-%m-%d %H:%M:%S"),donnees$value, type="l") And I'd just like to get the numerical value of the integration of this graph. I looked at ?integrate but, as far as I understood (that is, not very much, due to my poor

Thanks to your help, I managed to link R and a Mysql Database, send queries, plot the results and put everything in a pdf document (with Sweave). My co-workers find the job not bad but they would like to have a friendly interface to send queries and see the graphs (the pdf document is not necessary). Something like this coul be a good starting point : --------------------------------------- |

Hi, I have a vector called nfichiers of 138 names of file whose extension is .P0 or P1 ... to P8. The script is not the same when the extension is P0 or P(1 to 8). Examples of file names : [128] "Output0.P0" [129] "Output0.P1" [130] "Output0.P2" [131] "Output01102007.P0" [132] "Output01102007.P1" [133]

Hello, I'm just looking for an easy way to find the positions of a complete sequence in a bigger vector. For example : c("a","z","e") in c("a","z","e","r","t","a","z","a","z","e","c") and the result should be 1 8 that is the positions of the beginning of

Hello, I have 3 columns : a, b and a*b I would like to find the pair (a,b) so that a*b is the minimum but not from the points I measured but from the fit of the curve (I have more points that the ones given below but I fit only on this part because I know that the minimum a*b is in this interval). I thought doing it this way : - to fit a*b=f(a) abfit<-lm(ab ~ poly(a,8,raw=T)) - to use the

Hello, I convinced our new Ph student to work with R and I managed to convince the guy in charge of the computer networks to install R 2.12.0 (we haven't the admin rights). I wanted to install on his computer the same packages I use (but I have R 2.11) and especially RMySQL. Problem : there is no more RMySQL package in the list pf the packages. I've tried to install my version but it is

I'm plotting a lattice bwplot from a data.frame and would like to add some horizontal lines indicating some reference values (preferably with text annotation) as well as a grid of regular log tics. The following illustrates the plot: data(ToothGrowth) require("lattice") p = bwplot(exp(len) ~ supp, ToothGrowth, panel = function(..., box.ratio) { panel.grid(h=-1,

Hello, I would like to fit data with the following formula : y=V*(1+alpha*(x-25)) where y and x are my data, V is a constant and alpha is the slope I'm looking for. How to translate this into R-language ? At the moment, I only know : lm(y ~ x) Sorry for such a basic question. I thought I could find the solution in a post but I have to confess that, up to know, I'm not able to understand

Hello, I discovered R two years ago and thanks to the R-community I managed to write some scripts to analyze my data stored in mysql databases. The only problem is that I am the only one using R in the lab. Colleagues mainly use Matlab (but not with mysql, only with text files) but regularly come to me to get data treated with R-scripts !!!. To allow the use of my scripts by other people, my

Hello, I've two data.frames (data1 and data4), dec="." and sep=";". http://r.789695.n4.nabble.com/file/n4199964/data1.txt data1.txt http://r.789695.n4.nabble.com/file/n4199964/data4.txt data4.txt When I do plot(data1$nx,data1$ny, col="red") points(data4$nx,data4$ny, col="blue") , results seem very similar (at least to me) but the R-squared of

Hello, After a (bad) first try some months ago, I'm trying again to use RODBC package instead of DBI+RMySQL packages to populate MySQL database. The main command is sqlSave(channel, data_df, tablename, append=T, rownames=FALSE, colnames=FALSE, fast=FALSE, varTypes=vartypes) where data_df is the data.frame I want to save and vartypes is obtained following the tip given here

Hello, When I plot y=f(x) from the file xy.txt ( http://www.nabble.com/file/p18773387/xy.txt xy.txt ), I can clearly see a trend. Is there a function or a package able to take the median value of y for an interval of x (x +/- a defined value) to plot nice graph (at least a better one) ? Thanks in advance, Have a nice week-end, Ptit Bleu. -- View this message in context: