similar to: change likelihood function in survreg

Hello all, I'm a new user of R and I have used the spdep package to calculate Moran's I estimates. I used the dnearneigh() function and moran.test() function and obtained what I wanted, but now I want to plot my results with the p-values. I regrouped my results for the Moran's I like >m=c(m1$e,m2$e...) #m1 contain the results of the moran.test for lag1

Dear R users, I'm trying to fit a parametric survival model using the survreg function with a Weibull distribution. I'm studying the time to death of individuals from different families and I would like to fit different shape parameters (ie 1/scale in R) for each of the families. I looked it up in the help pdf and on the internet, but I couldn't find anything. Would it be possible to

Hi all, I am a very recent user of R. Mine is probably a basic question. I would like to know how do we create a certain set of points with geographic coordinates, between two endpoints. Please explain me how to do this using R. Which are the packages to be used for doing this. Thanks in advance, Ravichandra Mondreti. Ravichandra Mondreti PhD student/ Doctorant Spatial

Hello all, I'm a new user on R, and I used the abline function to draw a line on my graph but it doesn't work with xyplot(). With the plot function, abline() is ok but with plot it doesn't make what I want. Maybe it's a little thing to change in plot() but I don't find what! With xyplot I write this : Xyplot(m~ordered(l,levels=1),

Help, I tried to install the R-package on my computer, but R-Gui seems to provoke a failure in the R.dlll. I reinstalled the windows programs that R requires, but nothing changed at all (my portable computer seems to load smthg a long time before warning R-gui failed and that I have to close R- gui). What can I do??? ========================== Fr?d?ric ARCHAUX Equipe Interactions Individus

ls is a list of character vectors created by strsplit() I want to concatenate the 1st 4 character elements of each list item as a new vector called file. I admit to being confused about list syntax even after numerous readings. Here's what I tried: ls <- list(c("Focused", "10k", "A12", "t04.tif", "+", "µm"),

Dear R users, I face to a nested pattern and despite the numerous examples in the help I am still confused. I sampled bugs in different habitats within sites which were within rivers themselves within different regions. The habitat correspond to different substrata (not systematically present in all sites). For rivers and sites, I have environemental variables (e.g. altitude and slope of

Hi, I'm trying to make a regression of the form : formula <- y ~ Asym_inf + Asym_sup * ( (1 / (1 + (n1 * (exp( (tmid1-x) / scal1) )^(1/n1) ) ) ) - (1 / (1 + (n2 * (exp( (tmid2-x) / scal2) )^(1/n2) ) ) ) ) which is a sum of the generalized logistic model proposed by richards. with data such as these: x <- c(88,113,128,143,157,172,184,198,210,226,240,249,263,284,302,340) y <-

Colin, If your propose is to create a ternary plot with points and vectors, I think easier do this with graphics based plots instead of trellis based plots. Although, with little work you can do with trellis too. I gave you a reproducible code to put an arrow in a ternary plot. I use the function locator() to extract coordinates. The code is the following

I am trying to fit regression models type II with R, but it seems to me that most (all) of the linear model functions are for type I regressions. Does anyone knows whether type II regressions functions exist in R. Benjamin -- Benjamin Planque IFREMER Laboratoire d'Ecologie Halieutique BP 21105 44311 Nantes Cedex 03 France Tel: +33 (0)2 40 37 41 17 Fax: +33 (0)2 40 37 40 75 e-mail:

Please could you help me to find a package to apply a depth-constrained cluster analysis on palaeoecological data (in order to zone subfossil diagram)? Thanks, Dr. Emmanuel GANDOUIN Université d'Aix-Marseille 3 Institut Méditerranéen d'Ecologie et de Paléoécologie UMR CNRS 6116 Europôle Méditerranéen de l'Arbois, Bat. Villemin, BP 80 13545 Aix-en-Provence cedex 4 04 42 90 84 40

Bonjour, we wanted to fit a manova as descripted in Marieta /et al./ 2003, "convergent habitat segregation of /Aedes aegyptii/ and /Aedes albopictus/ (Diptera: /Culicidae/) in Southern Brazil and Florida", /J. Med. Entomol./, *40* (6), 785-794. They did their analysis with SAS software. We have the same kind of data with two insect species and city, habitat and season factors. We

Dear all, I'm puzzled by the definition of r-squared used by 'lm' when the model is without intercept. The help for summary says: >r.squared: R^2, the ``fraction of variance explained by the model'', > > R^2 = 1 - Sum(R[i]^2) / Sum((y[i]- y*)^2), > > where y* is the mean of y[i] if there is an intercept and > zero

Hello, I tried to find a way to having the setup of cups printing working for windows client. I end up now with "cupsaddsmb -v -a -U root" finishing by "Printer Driver mesange successfully installed. Running command: rpcclient localhost -N -U'root%xxxxxxxx' -c 'setdriver mesange mesange' cmd = setdriver mesange mesange result was NT_STATUS_UNSUCCESSFUL" My

Hi, I have two series of data set (it's measurment of growth but under two different conditions). To model these data I use the same function which is : formula <- y ~ Asym_inf + Asym_sup * ( (1 / (1 + (n1 * (exp( (tmid1-x) / scal1) )^(1/n1) ) ) ) - (1 / (1 + (n2 * (exp( (tmid2-x) / scal2) )^(1/n2) ) ) ) ) After the estimation of the parameters thanks to "nls", I have 2

Hi, I have to deal with a huge .txt table (~485.577 rows and 469 columns, > 1.5 Go file) I used the read.table function > tmp=read.table("data.txt", header=TRUE, sep="\t", fill=TRUE, na.strings="NA", comment.char="", stringsAsFactors = FALSE) However, I encounter troubles in interpreting some "\t" separator.. A column displays \t\t\t\t

Hi all, I am wondering whether there is a small bug in the binom.test function of the ctest library (I'm using R 1.6.0 on windows 2000, but Splus 2000 seems to have the same behaviour). Or perhaps I've misunderstood something. the command binom.test(11,100,p=0.1) and binom.test(9,100,p=0.1) give different p-values (see below). As 9 and 11 are equidistant from 10, the mean of the

Hi there, A couple years ago, I had patched the "fentonups" nut driver to allow it to work with my Unitek Alpha 500 UPS. My patch had been accepted and incorporated in the standard fentonups driver. I have had no time to stay involved in the "nut" project since (and besides this I'm not a developper...). 2 days ago I upgraded from Ubuntu Gutsy to Ubuntu Hardy and got

Hi All, I tried to apply "server=domain" for a samba server domain member with samba-3.0alpha23 and samba-3.0alpha24. In each cases, " ads_connect: Connection refused .." message appeared when trying to join the domain. But "net rpc testjoin" says ok. The main point of this mail is that " smbclient -L mysamba -U username" will work on

Hi, I would like to migrate from NT PDC to a Samba PDC. My main concern is that all my users have local profile (i.e. there is no roaming profiles). It seems that users profiles on a workstation are not valide when migrating to the Samba PDC. It seams that the users cannot keep the same identity. On the other hand, I would really like to change PDC without users noticing anything. Tools