similar to: Handling Negative value due to logarithm of probabilities.

Hi, In order to use the mle2-function, one has to define the likelihood function itself. As we know, the likelihood function is a sum of the logarithm of probability density functions (pdf). I have implemented myself the pdfs that I am using. My problem is, that the pdfs values are negative and I cann't take the logarithm of them in the log-likelihood function. So how can one take the

Hello I want to calculate natural logarithm of sum of combinations as follow: (R code) { com_sum=choose(2000000,482)*choose(1000000,118)+choose(2000000,483)*choose(1000000,117)+...+choose(2000000,i)*choose(1000000,600-i)+...+choose(2000000,600)*choose(1000000,0) #calculate the sum result=log(com_sum) #calculate the log of the sum } But

Hi, I am quite new to R and I have two questions regarding ROCR. 1. I have tried to understand how to extract area-under-curve value by looking at the ROCR document and googling. Still I am not sure if I am doing the right thing. Here is my code, is "auc1" the auc value? " pred1 <- prediction(resp1,label1) perf1 <- performance(pred1,"tpr","fpr") plot(

Dear all, i want to know if in R exist a method that give a logarithm interpolation of a two columns of dates. I try to do the lm of logarithm of the dates: a<-lm(log(vettore)~log(as.numeric(pippo2[3,]))) but this is a linear correlation and not that i want..... thanks giuseppe

Is it the tick labels that you want to change? -----Original Message----- From: "Elisabeth Bjerke Rastad" <ebr024 at post.uit.no> To: "r-help at r-project.org" <r-help at r-project.org> Sent: 5/10/10 11:20 AM Subject: [R] [Fwd: Re: Plotting log-axis with the exponential base to a plot with the default logarithm base 10] Hello! Thank you for answering! What I am

Hello! Thank you for answering! What I am trying to do is to plot my raw values (biomass of different species) on a logaritmic y-axis with the base of e. When I type "log="y"", the axis transforms into a logaritmic axis with the base of 10. Best regards, Elisabeth > Dear Elisabeth, > > I'm not sure if I have understood your question -- are you trying > to

Hello! I have a problem which I have tried to solve for several days now.. I have plottet a lineplot.CI in the library "sciplot", and I am trying to plot it with a logaritmic y-axis (with exponential base). The problem is that; when I type "log "y"", the axis transforms into the logaritmic of base 10. I wonder if someeone could tell me how to specify that I would

Hi, I have implemented initial prototype of the Xapian::Weight subclass for Unigram Language Modelling to support UnigramLM weighing in xapian.Other changes include adding collection_frequency to TermFreqs struct to store collection frequency of terms and some changes to support it xapian Framework,Changing simplesearch.cc to search using UnigramLMWeight class. Following issues have not being

I am working on "#699 Better tokenisation of mixed CJK numbers", and have implemented a partial patch of Chinese for this ticket. Current code works well with special test cases and all tests in xapian-core could still pass. But I'm confused with exact requirements of the question, for how much we could pay with performance on enabling more cases, and if there are better methods to

Dear friends, hope I could be able to explain my problem through following example. Please consider this: > set.seed(1) > input <- rnorm(10) > input [1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078 -0.8204684 0.4874291 0.7383247 0.5757814 -0.3053884 > tag <- vector(length=10) for(i in 1:10) # if there is any ****error**** in evaluating

Hi, I was reading a paper published in JCO "Prediction of risk of distant recurrence using 21-gene recurrence score in node-negative and node-positive postmenopausal patients with breast cancer treated with anastrozole or tamoxifen: a TransATAC study" (ICO 2010 28: 1829). The author uses a method to estimate the 9-year risk of distant recurrence as a function of continuous recurrence

Hi! I'm using a large panel data, and now I have faced some difficulties with my analysis. The predictors are not normally distributed and there are quite many outliers (some of them are influential though). I have tried to change the logarythm, but i'm not sure, how to do that. I want also draw a plot picture in which logarythms of predictors x and y are changed. How could I do that?

I estimated a negative binomial model using zelig. z.out<- zelig(NEWBHC~ PW80 + CHNGBLK + XBLK,data=data, model="negbin") How do I calculate predicted probabilities for this model? Is it the same process as a poisson regression? Thanks in advance Joe [[alternative HTML version deleted]]

Dear all, I would like to (i) produce boxplot graphs with axis in logarithm in base 10 and (ii) showing the values on the axis in 10^exponent format rather than 10E+exponent. To illustrate with an example, I have some widely spread data that I chart plot using boxplot() [figure on the left]; the log="y" option of boxplot() I obtained the natural logarithm conversion of the data and

Greetings all, So, I''ve finally found a place to play around with the Ajax.PeriodicalUpdater. In looking at the API, I''m liking the decay option -- not necessarily for my current purpose, but just to keep in mind -- and I have a question: can the decay be a function which returns an integer? Basically, why I''m looking for is a logarithmic decay (where the system

I am using SVM to classify categorical data and I would like the probabilities instead of the classification. ?predict.svm says that its only enabled when you train the model with it enabled, so I did that, but it didn't work. I can't even get it to work with iris. The help file shows that probability = TRUE when training the model, but doesn't show an example. Then I try to

[Hope this is the right list where to send...] An attempt to update package 'mnormt' involves the addition of a small new function called 'pd.solve'. When I come to the package checking stage, an error occurs in parsing pd.solve.Rd. The full transcript of the outcome is copied below (it includes details on my installation) but the critical point is where the \examples{} section

Dear All, I am using the plotrix library to plot some matrices. I have a problem: some of my data are outliers, hence using a linear color scale does not work very well (you would see too many cells having a similar, indistinguishable color). See the code snipped at the end of the email. Plotting the logarithm of the data gets the job done, but my problem is that I would like to write in every

Hi All, This is surely an easy question but somehow I am not being able to get it. I am using R 2.13.2 and have a data set where variable names like this appear: pci1990, pci1991, ... , pci2009. "pci1990" has data on per capita income for 1990, "pci1991" has data on per capita income for 1991, and so on. I would like to create the logarithm of per capita for each of the

I have some data where the frequency is heavily weighted on the lower end. So I have lots of low values with very few higher values. I would like to find breakpoints that cover the data with as much detail as possible. I find that if I use hist() to automatically find the breaks for me it finds breaks that are too coarse for the low values. I have tried the other algorithms (like 'Scott'