similar to: make R program faster

Hello, I've come across a problem in developing a set of custom functions to calculate the number of hours of daylight at a given latitude, and the number of days a date precedes or secedes the summer solstice. I discovered an inconsistency concerning leap years between my derived values and those from the US naval databases. It seems as far as I can figure that my inconsistency arises either

Dear List, I am having problems running stepAIC with a negative binomial regression model.  I am working with data on manta ray abundance, using 20 predictor variables.  Predictors include variables for location (site), time (year, cos and sin of calendar day, length of day, percent lunar illumination), oceanography (sea surface temp mean and std, sea surface height mean and std), weather (cos

Hello all, I have a (hopefully) simple question. If I have a data table (named "x") such as this: daylengthleavesflowering.datefruits LD1.90055581.64353250.76668986 LD-0.9399410-1.2592968-0.60141837 LD1.0222946-0.9697459 -0.60141837 SD-0.2668132-0.9697459-0.47363322 SD -1.5892393-0.70540660.66469785 SD -0.5222961-0.7054066-0.10701187 How do create a new data table that contains,

Hello all, I'm fitting an allometric equation that looks like a really clean fit in the log-log space, but when I back transform the fit of the curve appears to need an adjustment - the fitted curve appears to predict values a good deal higher than those from the data. I included a bias correction per Newman, M. (1993). Regression analysis of log‐transformed data: Statistical bias and its

Dear R community, I am trying to write my own user defined split function for rpart. I read the example in the tests directory and I understand the general idea of the how to implement user defined splitting functions. However, I am having troubles with addressing the data frame used in calling rpart in my split functions. For example, in the evaluation function that is called once per node,

hi, i'd like to calculate a temperatursum, adding the value of each element. let's say the data looks like this: x<-c(1,2,3,4,5) what i want to do, is ploting not the sum in the end but all the subresults, too, so my vector holds: x[i] [1] 1,3,6,10,15 here is what i tried, which seems to be right to me, bu doesn't work out: x<-c(1,2,3,4,5) i<-1 j<-1

Dear all, To modelize the abundance of fish (4 classes) with a set of environmental variables, I used the polr and predict.polr functions. I would like to know how to bring the cumulated probabilities back to a discrete ordinal scale. For the moment I used the predict.polr function with the argument "class". Is there an other way? polrf <- polrf <- polr_mod(formula =

Ben, I tried playing with the "na.action" but that didn't seem to help.  I didn't know of the na.omit function.  It made it much easier to get rid of all the missing values in the file, and now the regression model is running. I hadn't put interactions in the model yet since it wasn't even running.  Cosday, sinday, and daylength are strongly collinear.  The full model

Hello, i want to estimate a complex GARCH-model (see below). http://r.789695.n4.nabble.com/file/n4112396/GJR_Garch.png W stands for the Day of the Week Dummies. r stands for returns of stock market indices. I stands for the GJR-term. I need some help with three problems: 1.) implementation of the GJR-term in the variance equation 2.) compute robust covariance matrix

Dear help community, I'm a R-beginner and use it for my master thesis. I've got a mixed model and want to analyse it with lme. There are a lot Cofactors that coult be relevant. To extract the important ones I want to do the stepAIC, but always get an error warning. Structure of my data: data.frame': 72 obs. of 54 variables: $ Block : Factor w/ 3 levels

Dear all I am a BEGINNER and have R on my Mac. I saved my excel file as .txt file, I have just one column with first row as the column name. My file when read by R looks like this. After reading the table I try to make a histogram by hist(dbh), it says object dbh not found. What am I doing wrong? thanks dbh 1 11.53 2 16.05 3 7.36 4 16.05 5 8.66 6 12.74 7 22.93 8 7.55 9

Hello I am trying to import an Excel file but I am loosing the headers, My headers are in the first to rows of the EXCEL file. In the following R output, the NA are supposed to be the second item in the Header. Is there any way to Import more than one row as headers?. Thank you Felipe Parra > Data<-odbcConnectExcel('Curva LIBOR.xlsx',readOnly=T,) >

Hi All , Can anybody tell me if there's any way to get the summarized anova values.Now i will explain what i mean , when i say "*summarized*". Below you can see the anova table of recmeanC1 with rest* all* i.e from recmeanC2 to i15(predictors),as shown in table. Df Sum Squares Mean Square F value Significance [Pr(>F)] recmeanC2 1 89.272 89.272

Hello Scott, On Wed, 2007-02-21 at 18:57 -0500, me22 wrote: > I'm trying to follow docs/Projects.html to use llvm in a project of my own. > > If I just copy the sample directory, running configure in the new one > works fine: > ~/programming/llvm-1.9/projects $ cp -a sample/ zhe > ~/programming/llvm-1.9/projects $ cd zhe/ > ~/programming/llvm-1.9/projects/zhe $

On 21/06/2012, Douglas Gregor <dgregor at apple.com> wrote: > > On Jun 20, 2012, at 5:13 PM, Nick Lewycky <nlewycky at google.com> wrote: > >> Is there anybody who is certain that our autoconf dependency needs to stay >> around? Are there developers stuck on systems that don't have a recent >> enough cmake in their most recent release, or maybe are using

Hi everybody! I make and install autoconf 2.61, after I wget samba-3.0.25, untar and do ./autogen ./configure without any problem!!! But now I do "make" and after many minutes I get: : mush more printing/print_cups.c etc etc etc : printing/print_cups.c:1238: `IPP_GET_PRINTER_ATTRIBUTES' undeclared (first use in this function) printing/print_cups.c:1243: `IPP_TAG_OPERATION'

I'm trying to follow docs/Projects.html to use llvm in a project of my own. If I just copy the sample directory, running configure in the new one works fine: ~/programming/llvm-1.9/projects $ cp -a sample/ zhe ~/programming/llvm-1.9/projects $ cd zhe/ ~/programming/llvm-1.9/projects/zhe $ ./configure configure: creating ./config.status config.status: creating Makefile.common config.status:

On 2/23/07, Reid Spencer <rspencer at reidspencer.com> wrote: > > However, if I run the AutoRegen.sh script (after modifying it to > > accept autoconf 2.61, which I think is fine as gentoo slots it the > > same as the 2.5[0-9] versions), the configure no longer works: > > ~/programming/llvm-1.9/projects/zhe $ cd autoconf/ > >

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message is always estimated breakpoints are out of range. Since my problem is time related problem. So I

Hi, I am using 2.6.36.3 kernel with btrfs, 512MB memory and a very slow disk, no special options for mounting btrfs except noatime. Now I found it very slow. When I rm a 5GB movie, it took 20 secs. -- 竹密岂妨流水过 山高哪阻野云飞 -- To unsubscribe from this list: send the line "unsubscribe linux-btrfs" in the body of a message to majordomo@vger.kernel.org More majordomo info at