similar to: integrate fails with errors

Thks for your answer, here is an exemple of what i do with the errors in french... > tmp [1] 200 150 245 125 134 345 320 450 678 > beta18 Erreur : Objet "beta18" not found //NORMAL just to show it > eta [1] 500 > func1<-function(beta18) dweibull(tmp[1],beta18,eta) > func1<-func1(beta18) * function(beta18) dweibull(tmp[2],beta18,eta) Erreur dans dweibull(tmp[1],

hi all I wish to thanks every body on the R mailing list for answering very fast, directly in my mail box ;). I've finish my work with R and i can say that it is very difficult at the beginning, and when you succeed you are stopped by a stack overflow when you call your nice recursive function (which was working with a tab of 100 element) with a tab of 900 elements, but R just do what you

hi all, I wish to draw on the same graphic device 3 functions. But i don't want them to be on different graph, i want to compare them on the same I don't need mfrow or mfcol, I need something else... 1 graph on 1 device inside this graph 3 ploted function. I saw something unsing data.frame, but i think it's overkill, and something less complicated must exist, if not why? why not

hi all, I need to generate a function inside a loop: tmp is an array for (i in 1:10) { func<- func * function(beta1) dweibull(tmp[i],beta1,eta) } because then i need to integrate this function on beta. I could have written this : func<-function(beta1) prod(dweibull(tmp,beta1,eta)) (with eta and beta1 set) but it is unplottable and no integrable... i could make it a bit different but

hi all, I wish to draw on the same graphic device 3 functions. But i don't want them to be on different graph, i want to compare them on the same I don't need mfrow or mfcol, I need something else... 1 graph on 1 device inside this graph 3 ploted function. I saw something unsing data.frame, but i think it's overkill, and something less complicated must exist, if not why? why not

Hi all i need to put on one graph 2 functions who's x axis is the same and y not. I mean on horizontal the time, and on vertical left: pressure, on vertical right: rpm of a motor, is R able to do that? i've found this that i could adapt maybe (i don't need time series really?) :/ : (http://tolstoy.newcastle.edu.au/R/help/04/03/1456.html) ## ## Description: A simple function which

hy all, I wish to switch in a part in my code that use "read.table" to "scan", actually i use this: tab<-scan("data.dat",what=integer(),skip=1) dim(tab)<-c(75,length(tab)/75) tab<-t(tab) It gives me a new tab with 75 columns, but i when i was using read.table with headers then attach i could use the columns names to access the data values, now how can i

Hy all, I wish to use the date function to draw againt the lifetime of a motor. This lifetime is given to me in Hours (it can go over 5000 hours) I'm unable to find how to convert this lifetime value to something like %H:%M:%S because when R see 24H it says 1 day, i don't want that, i just want %H:%M:%S with a value of %H higher than 24... for example: i've got this value in hours:

hi all Here is a small part of my code: tab_tmp<-tab[1:(no[off_set[i-1]+1]+(no[off_set[i]+1]-no[off_set[i-1]+1])),length(tab)]; tab_tmp1<-tab[(no[off_set[i-1]+1]+(no[off_set[i]+1]-no[off_set[i-1]+1])):length(TotalFillTimeHours),length(tab)]; tab<-c(tab_tmp,tab_tmp1); attach(tab); Here is the output: Error in attach(tab) : attach only works for lists and data frames Execution halted

Hi all I wish to draw an histogram... with dates but the following append, i don't know where is the problem, help(hist.Date) works and i don't see any usefull information on what i'm doing wrong... > hist.Date(dt_cycles) Error: couldn't find function "hist.Date" > hist.date(dt_cycles) Error: couldn't find function "hist.date" > cycles [1] 7 1

Greetings all- I'd like to install Oreka on a Centos 5.x server for monitoring my Asterisk systems(using port mirroring) but I find I'm having problems with the version of libpcap installed. The latest is 0.9.4 but the orkaudio RPM (built for Centos 4.2) requires libpcap 0.8.3. I've tried making symlinks to overcome the issue but with no success. Without having to build Oreka from

Hello folks, I wanted to setup Oreka to monitor calls on a trixbox box I have setup. Oreka doesn't seem to be catching all of the calls though.... I have port mirroring setup on the port that trixbox is connected to, mirrored to the port Oreka is connected to. I have read that Asterisk doesn't work as a SIP Proxy, so I wondered if this meant that some phones, after checking in with

Hi, I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and \beta_1, this can be achieved by solving the following three equations: n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) - \sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} 1/(psihat(i)) - (\alpha+1) \sum\limits_{i=1}^{n} ( 1 / (x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} (

Hi all-- We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the "power" code is just returning the coefficients, as something isn't working quite right. The

I use "while" loop but it produces an errro. I have no idea about this. Error in "[<-"(`*tmp*`, i, value = numeric(0)) : nothing to replace with The problem description is The likelihood includes two parameters to be estimated: lambda (=beta0+beta1*x) and alpha. The algorithm for the estimation is as following: 1) with alpha=0, estimate lambda (estimate beta0

Dear Users, I am new to both of things, so do not blame me too much... I am busy with semiparametric regression and use WinBUGS to sample posteriors. The code to call Winbugs is as follows: data <- list("y","X","n","m") #My variables inits.beta <- rep(0,K) inits.beta0 <- 0 inits <-

hy all I wish to draw a graph against a date, I have a set of date like this DD/MM/YYYY, corresponding to it a set of integer , i wish to draw on x side the dates (the space between the dates have to be constant, not based on the time between 2 dates but on the number of dates) and on y side the integers. Do i have to make a tricky function under R ? I've search the help for graphical

Hy all I use barplot to draw frequencies by dates. On the x axis it shows only 1 date for 2 bars, i've understand why, because R cant put so much date on the same axis, it will get out of the graph. that's why i tought about reducing the size of the font used to stamp the x axis. anyone knows the right way? i've tryed par(cin) and barplot(cin=3) but it says i cannot do it at this

hi all, I load a table with headers that enable me to acces it by the column names: tab<-read.table("blob/data.dat",h=T) attach(tab) everythings are OK, but i try to sort this table against one of his column like this: tab<-tab[order(tab$IndexUI),]; It is still ok, the table is sorted, if i type "tab" i see a sorted table. but, when i call the column by their names,

Hi all, Before the amount of data given has grown i was initially using read.table to load the values inside R. It was feeding my needs because i could tell read.table h=T, then use attach to access the values by columns names. Now it takes 20 seconds to load the data's and the first enhancement i could do is to win some time on the load of the data's... How could i use the scan