similar to: Multiple dependent variables

Dear R-helpers, a S-code for piecewise regressions was provided by Toms & Lesperance (2003) Ecology, 84, 2034-2041 (paper can be found on the web).The code is quite complete with different types of transitions around breakpoints and model selection fonctions. It doesn't work directly under R due to some "translation" problems I guess. However I reckon that it would be a

Hello everybody!!!! Quick question, if you'd like to throw a little tip: does anyone knows a function that runs piecewise regression models with coefficients estimation and inferences ? Thank you [[alternative HTML version deleted]]

I typed the following command *install.packages('bayesQR')'* to install bayesQR(my R version is 2.14) i am encountered the following error. Installing package(s) into ?C:/Users/knreddy.IDRBTVM/Documents/R/win-library/2.14? (as ?lib? is unspecified) Warning: unable to access index for repository http://essrc.hyogo-u.ac.jp/cran/bin/windows/contrib/2.14 Warning: unable to access index

Dear all, I did a non-linear least square model fit y ~ a * x^b (a) > nls(y ~ a * x^b, start=list(a=1,b=1)) to obtain the coefficients a & b. I did the same with the linearized formula, including a linear model log(y) ~ log(a) + b * log(x) (b) > nls(log10(y) ~ log10(a) + b*log10(x), start=list(a=1,b=1)) (c) > lm(log10(y) ~ log10(x)) I expected coefficient b to be identical

Dear all, sorry if this is a bit on the sidetrack for R-help. As a regular R user I have developed quite a lot of custom R functions, to the point of not always remembering what I have already programmed, where the file is and so on. I was wondering what other people do in this regards. A basic file with all your functions, or a custom R package, or directly integrated into a profile file

Hi, I am an old hand at chemistry but a complete beginner at statistics including R computations. My question is whether you can carry out nonlinear multivariate regression analysis in R using neural networks, where the output variable can range from -Inf to + Inf., unlike discriminant analysis where the output is confined to one or zero. The library nnet seems to work only in the latter

Dear R-helpers, Happy New Year to all the helpful members of the list. Here is the behavior I'm looking for: > v1 <- c("a","b","c") > n1 <- c(0, 1, 2) > v2 <- c("c", "a", "b") > n2 <- c(0, 1 , 2) > (f1 <- data.frame(v1, n1)) v1 n1 1 a 0 2 b 1 3 c 2 > (f2 <- data.frame(v2, n2))

Hello, I would like to fit a linear regression and when I use summary(), I got the following result: Call: lm(formula = weight ~ group - 1) Residuals: Min 1Q Median 3Q Max -1.0710 -0.4938 0.0685 0.2462 1.3690 Coefficients: Estimate Std. Error t value Pr(>|t|) groupCtl 5.0320 0.2202 22.85 9.55e-15 *** groupTrt 4.6610 0.2202 21.16 3.62e-14

Hello, I am using package MuMIn to calculate AIC for a full model with 10 explanatory variables. Thanks in advance in sharing your experience. Q1 In the AIC list of all models, each model is differentiated by model number. Please kindly advise if it is possible to find the corresponding explanatory variable(s) for the model number. Q2 error message I tried to display sub-model with only

I have fitted the faults.data to glm.nb and to the function negbin from the package aod. The output of both is the following: summary(glm.nb(n~ll, data=faults)) Call: glm.nb(formula = n ~ ll, data = faults, init.theta = 8.667407437, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.0470 -0.7815 -0.1723 0.4275 2.0896 Coefficients:

I'm using stepAIC to fit a model. Then I'm trying to use that model to predict future happenings. My first few variables are labeled as their column. (Is this a problem?) The dataframe that I use to build the model is the same as the data I'm using to predict with. Here is a portion of what is happening.. This is the value it is predicting = > [1] 9.482975 Summary of the

Hi I have a panel dataset with large number of groups and differing number of observations for each group. I want to randomly select say, 20% of the groups or 200 groups, but along with all observations from the selcted groups (with the corresponding data). I guess it is possible to generate a random sample from the groups ids and then match that with the entire dataset to have the intended

Hi, i was searching in the source code, but i didn't find any MAX_MPX or something like that. How can i change MAX_MPX_REQUESTS that smbd can handle?

I have empirical Bayes estimates for slopes and intercepts for a number of subjects and I would like to plot the slopes and intercepts with confidence ellipses. These ellipses would be based on the confidence intervals for the slope and intercepts (forming the major and minor axis of each ellipse), and the correlation in the slope and intercepts. The ellipse function in the car library

<select name="media" id="media" onchange="<%= remote_function( :url => { :action => :show_performer_categories, :rate_group_id => @rate_group }, :with => "''media_id=''+this.value''" ) %>"> I''m guessing that''s wrong since WebRick does nothign when I

I''m stumped... I need to return a number of rows with the top & bottom x percentage of rows removed in mysql 4.x e.g. SELECT COUNT(*) value_count, SUM(deals.value) / COUNT(*) average_value FROM deals ORDER BY deals.value Now say that returns the sum of 100 deals, what I actually need to do is return the count & sum of rows 10-90 (i.e. the top & bottom 10% rows removed)

Hi, I have read your post about plotting multiple CI ellipses and I would like too, but I don?t know how. I would like to plot an ellipse for each group and I have slope and intercepts, but I don?t know what to put in the arguments shape and radius. CAn you help me please?

I am using R0.64.1 under Windows 95. The help for `pmatch' states that: If `duplicates.ok' is false multiple matches will result in the value of `nomatch' being returned, and if it is true, the index of the first matching value will be returned. The help for pmatch gives the following examples: pmatch("m", c("mean",

Dear all, I have about 50 pages of R codes and ran it in both OS X and Windows. It takes at least haft a day to have the results. The running time is not very different in both Systems. I found that R does not use all cores of CPU by default. Can anybody help me to use all cores of CPU in my programming from the beginning or through programming in both OS X and Windows? Many thanks your

I am trying to make a graph with 4 scatter matrixes plots and couldn't do it. While trying to find a solution for this I also came across the idea of giving different values to the same argument for each of the lower and upper function but couldn't do it. (Examplified below with the col argument). The first problem of plotting 4 scatter matrixes in a graph is a problem of real