similar to: vardiag Package and nlregb

Hi out there, I am trying to fit a species accumulation curve (increase in number of species known vs. sampling effort) for multiple regions and several bootstrap samples. The bootstrap samples represent different arrangements of the actual sample sequence. I fitted a series of nlme-models and everything seems OK, but since the observations are correlated I tried to include some correlation

I have an interesting lme - problem. The data is part of the Master Thesis of my friend, and she had some problems analysing this data, until one of her Jurors proposed to use linear mixed-effect models. I'm trying to help her since she has no experience with R. I'm very used to R but have very few experience with lme. The group calls of one species of parrot were recorded at many

Hi, First of all a big THANK YOU to all who answered me yesterday. I am back to my problems with outliers. I did a qqnorm on my data (as somebody suggested), but I will like to compare results doing other plots as well. I used ?cook.distance? on lm and glm objects, but still I am not happy with that. So now I am loading the ?vardiag? package. First thing is to transform my data

Hi, I've implemented parallelization in one of my packages using the 'parallel' package -- many thanks for providing it! In my package I'm importing 'parallel' and so added it to the DESCRIPTION file's 'Import:' tag and also added a 'importFrom("parallel", ...)' statement in the NAMESPACE file. Parallelization works nicely, but my package

My data looks like this: > data name G_hat_0_0 G_hat_1_0 G_hat_2_0 G_0 G_hat_0_1 G_hat_1_1 G_hat_2_1 G_1 1 rs0 0.488000 0.448625 0.063375 1 0.480875 0.454500 0.064625 1 2 rs1 0.002375 0.955375 0.042250 1 0.000000 0.062875 0.937125 2 3 rs2 0.050375 0.835875 0.113750 1 0.877250 0.115875 0.006875 0 4 rs3 0.000000 0.074750 0.925250 2 0.897750 0.102000

Hello, I have a small statistics question, and as I'm quite new to statistics and R, I'm not sure if I'm doing things correctly. I am looking at two quantitative variables (x,y) that are correlated. When I divide the data set according to a categorical variable z, then x and y are more poorly correlated when z = A than when z = B (see attached figure). In fact x and y are two

Hello, I am trying to get the P values from the output of a summary for lm. lm <- lm(y ~ age + sex) s <- summary(lm) I thought that I might be able to get them using a combination of scan, grep and sub. But I got stuck on the first step - being able to process "s" as a text string. I could perhaps write it to file than scan it back but there is probably an easier way to do

I would like to fit a model of the type Y = Xf + e where f are cover fractions (or percent cover). Therefore, I must constrain the fit to 1. sum(f) = 1 2. 0<= f <=1 How can I introduce these 2 constraints in a ML fit? Currently, I'm using nlregb and I'm using the "upper and low" keywords to set the bounds and use the following to ensure that sum(f) = 1: let's say

Hello everyone, I'm trying to test the accurracy of R on the Eckerle4 dataset from NIST and I don't understand how the control option of the nls function works. I tought nls(...) was equivalent to nls(...control=nls.control()) i.e nls.control() was the default value of control, but here is the error I get : > n2=nls(V1~(b1/b2) *

Hi, I've been using "nlregb" routine in splus3.4, but cannot find it here in R. Is it hidden in some library? Or there is a routine equivalent to it but with a different name? I'm planning to move part of my from S to R so that I can work comfortablely at home. Any suggestion is highly appreciated. regards, Tony -- \|||/ Q o o Q ==========================oooQ= ( -

Hi All Please see the Rprofile file which i have modified as follows and after that when I start R then I see that R says to me "TRUE" for all the packages implying that all loaded at once. But when i try to use commands as simple as help("lm"), it doesnt work nor any of the menu "Packages" is not working. Although the regression using lm ( Y ~ X ) is working

Hello, I've had no luck finding an R package that has the ability to estimate a Tobit model allowing for heteroskedasticity (multiplicative, for example). Am I missing something in survReg? Is there another package that I'm unaware of? Is there an add-on package that will test for heteroskedasticity? Thanks for your help. Cheers, Alan Spearot -- Alan Spearot Department of Economics

Hello, I would like to inlude the Q-Q plot by "qqmath" into a panel with other plots, say, using par(mfrow=c(1,2)). How can this be done given that "qqmath" refreshes the plotting window and there seems to be no series coming out of it? Thanks Serguei [[alternative HTML version deleted]]

Hi, I'm trying to aggregate date values using the aggregate function. For example: aggregate(data,by=list(weekdays(LM),months(LM)),FUN=length) I would also like to aggregate by year but there seems to be no years() function. Should there be one? Is there any alternative choice? Also, a hours() function would be great. Any tip on this? Thanks in advance! S?rgio Nunes

> I believe the problem is related to the change in the > installation directory for the updated JVM from Sun. > > If you use the jre1.5.0_02 from java.sun.com, then the plugin > is installed in: > > /usr/java/jre1.5.0_02/plugin/i386/ns7/libjavaplugin_oji.so Ok, the above pathname is the exact one. Yes, this is the one I used. > > I just downloaded the updated

Hi there, I'm writing a function that calculates the probability of different outcomes of dice rolls (e.g., the sum of the highest three rolls of five six-sided dice). I'm using the "combinations" function from the "gtools" package, which is great: it gives me a matrix with all of the possible combinations (with repetitions allowed). Now I want to count the number

User reports: 0% of 640MB -- which cannot be true. I logged in as masteruser and was able to confirm this. I'm getting "0.0% von 640.0MB" == "0.0% of 640.0MB" His maildirsize file is attached. His mailbox is nearly full. Manual login: # telnet localhost 143 Trying 127.0.0.1... Connected to localhost.localdomain. Escape character is '^]'. * OK [CAPABILITY

CRAN (and crantastic) updates this week New packages ------------ * allan (1.0) Alan Lee http://crantastic.org/packages/allan Automates Large Linear Analysis Model Fitting * andrews (1.0) Jaroslav Myslivec http://crantastic.org/packages/andrews Andrews curves for visualization of multidimensional data * anesrake (0.3) Josh Pasek http://crantastic.org/packages/anesrake This

The Wine development release 1.6-rc4 is now available. What's new in this release (see below for details): - Bug fixes only, we are in code freeze. The source is available from the following locations: http://prdownloads.sourceforge.net/wine/wine-1.6-rc4.tar.bz2 http://mirrors.ibiblio.org/wine/source/1.6/wine-1.6-rc4.tar.bz2 Binary packages for various distributions will be available

Comunidad superpower, tengo la siguiente pregunta-problema. Tengo una matriz de datos como la siguiente (adjunto archivos de datos): sol con dia rep nca nin iso are 1 con 0 1 1 16 0 s 0.3866520976 2 con 0 1 1 16 1 c 0.0008842527 3 con 0 1 1 18 0 s 0.2409465396 4 con 0 1 1 18 1 c 0.0021092822 5 con 0 3 1 16 0 s 0.6005446999 6 con