[R] Accounting for clustered data in Rpart or Mvpart
I am working on data that is in clusters (eg events from persons in the family).So I have a clustered data set up and would like to build a tree(regression /classification tree) taking into account this clustered nature of the data. Do anybody know how to do this or maybe the code to take into account the clustered data in rpart __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Soil texture triangle in R?
Dear R users, has anybody made an attempt to create the soil texture triangle graph in R? For an example see here: http://www.teachingkate.org/images/soiltria.gif I would like to get the lines in black and texture labels in gray to allow for plotting my texture results on top. Any examples or suggestions are very welcome! Thanks in advance, Sander. -- Dr Sander P. Oom Animal, Plant and Environmental Sciences, University of the Witwatersrand Private Bag 3, Wits 2050, South Africa Tel (work) +27 (0)11 717 64 04 Tel (home) +27 (0)18 297 44 51 Fax +27 (0)18 299 24 64 Email [EMAIL PROTECTED] Web www.oomvanlieshout.net/sander __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] mixed-integer optimisation
Does anybody know of a R package to solve mixed-integer quadratic optimisations with constraints? Alternatively is anybody working with an R interface to ILOG CPLEX which seems to be the 'best' commercial product in this area? Finally, is there any relationship between the R-project and AMPL - a modelling language for mathematical programming www.ampl.com - which also originated in Bell Labs? Regards, John Marsland [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Contingency tables from data.frames
The final version with the help of Gabor Grotendieck (thanks Gabor, very much!)
###
# EasieR - Package #
###
# Common function
er.make.table - function(x,
start,
end,
h,
right) {
# Absolut frequency
f - table(cut(x, br=seq(start, end, h), right=right))
# Relative frequency
fr - f/length(x)
# Relative frequency, %
frP - 100*(f/length(x))
# Cumulative frequency
fac - cumsum(f)
# Cumulative frequency, %
facP - 100*(cumsum(f/length(x)))
fi - round(f, 2)
fr - round(as.numeric(fr), 2)
frP - round(as.numeric(frP), 2)
fac - round(as.numeric(fac), 2)
facP - round(as.numeric(facP),2)
# Make final table
res - data.frame(fi, fr, frP, fac, facP)
names(res) - c('Class limits', 'fi', 'fr', 'fr(%)', 'fac', 'fac(%)')
return(res)
}
#With Gabor Grotendieck suggestions (thanks Gabor, very much!)
er.table - function(x, ...) UseMethod(er.table)
er.table.default - function(x,
k,
start,
end,
h,
breaks=c('Sturges', 'Scott', 'FD'),
right=FALSE) {
#User define nothing or not 'x' isn't numeric - stop
stopifnot(is.numeric(x))
#User define only 'x'
#(x, {k, start, end, h}, [breaks, right])
if (missing(k) missing(start) missing(end) missing(h) ){
x - na.omit(x)
brk - match.arg(breaks)
switch(brk,
Sturges = k - nclass.Sturges(x),
Scott = k - nclass.scott(x),
FD = k - nclass.FD(x))
tmp - range(x)
start - tmp[1] - abs(tmp[2])/100
end - tmp[2] + abs(tmp[2])/100
R - end-start
h - R/k
}
#User define 'x' and 'k'
#(x, k, {start, end, h}, [breaks, right])
else if (missing(start) missing(end) missing(h)) {
stopifnot(length(k) = 1)
x - na.omit(x)
tmp - range(x)
start - tmp[1] - abs(tmp[2])/100
end - tmp[2] + abs(tmp[2])/100
R - end-start
h - R/abs(k)
}
#User define 'x', 'start' and 'end'
#(x, {k,} start, end, {h,} [breaks, right])
else if (missing(k) missing(h)) {
stopifnot(length(start) = 1, length(end) =1)
x - na.omit(x)
tmp - range(x)
R - end-start
k - sqrt(abs(R))
if (k 5) k - 5 #min value of k
h - R/k
}
#User define 'x', 'start', 'end' and 'h'
#(x, {k,} start, end, h, [breaks, right])
else if (missing(k)) {
stopifnot(length(start) = 1, length(end) = 1, length(h) = 1)
x - na.omit(x)
}
else stop('Error, please, see the function sintax!')
tbl - er.make.table(x, start, end, h, right)
return(tbl)
}
er.table.data.frame - function(df,
k,
breaks=c('Sturges', 'Scott', 'FD'),
right=FALSE) {
stopifnot(is.data.frame(df))
tmpList - list()
logCol - sapply(df, is.numeric)
for (i in 1:ncol(df)) {
if (logCol[i]) {
x - as.matrix(df[ ,i])
x - na.omit(x)
#User define only x and/or 'breaks'
#(x, {k,}[breaks, right])
if (missing(k)) {
brk - match.arg(breaks)
switch(brk,
Sturges = k - nclass.Sturges(x),
Scott = k - nclass.scott(x),
FD = k - nclass.FD(x))
tmp - range(x)
start - tmp[1] - abs(tmp[2])/100
end - tmp[2] + abs(tmp[2])/100
R - end-start
h - R/k
}
#User define 'x' and 'k'
#(x, k,[breaks, right])
else {
tmp - range(x)
start - tmp[1] - abs(tmp[2])/100
end - tmp[2] + abs(tmp[2])/100
R - end-start
h - R/abs(k)
}
tbl - er.make.table(x, start, end, h, right)
tmpList - c(tmpList, list(tbl))
}
}
valCol - logCol[logCol]
names(tmpList) - names(valCol)
return(tmpList)
}
Best,
--
Jose Claudio Faria
Brasil/Bahia/UESC/DCET
Estatistica Experimental/Prof. Adjunto
mails:
[EMAIL PROTECTED]
[EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: 73-3634.2779
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Simplify formula for heterogeneity
Thank you very much Ted! I have been looking at your simplification for more then an hour, but I don't see how you did it. Could you perhaps, if it is not to much work, explain me how you reduced H? It would help me to understand what I am realy doing. Looking at the result, it seems indeed that H does add more information than sd already did. Intuitively I thought the square of the sum of all possible differences would not be related to the standard deviation. Looking at your result it seems it is related by a factor sqrt(2*(n-1)) so there is no special point in calculating H and I know I cannot trust my intuition anymore. Thanks again! Kind regards, Stef __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Testing Nonlinear Restrictions
Dear all, I'm interested in testing 2 nonlinear restrictions on coefficients of a nls object. Is there a package for doing this? Something in the lines of `test(nls object, res=c(res 1,res 2),...)' I only found the function delta.method in the alr3 library that calculates the se of a singleton nonlinear restriction of a nls object using the delta method. Thanks in advanced for your help and suggestions. David __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] cancor
Hello, I'm a beginner to use R. When use cancor to analysis two data set X,Y, for example: can-cancor(X, Y) then 5 components are returned, that is cor, xcoef, ycoef, xcentre, ycentre; The explained variance can be calculated by X %*% can$xcoef or Y %*% can$ycoef, but how to alculate the canonical maps of field X and Y? Thanks! yana __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] R commandline editor question
I am using R 2.1 on Apple OS X.
When I get the prompt, I find it works well with emacs commandline
editing. Keys like M-f C-k etc. work fine.
The one thing that I really yearn for, which is missing, is bracket
matching When I am doing something which ends in it is really
useful to have emacs or vi-style bracket matching, so as to be able
to visually keep track of whether I have the correct matching
brackets, whether ( or { or [.
I'm sure this is possible. I will be most grateful if someone will
show the way :-) Thanks,
--
Ajay Shah Consultant
[EMAIL PROTECTED] Department of Economic Affairs
http://www.mayin.org/ajayshah Ministry of Finance, New Delhi
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] R commandline editor question
Hi Ajay
well ESS has such a facility.
However, I think Mathematica has a super scheme: unbalanced brackets
show up
in red, making them obvious.
This is particularly good for spotting wrongly interleaved brackets, as
in
([ blah di blah )]
note bracket closure is out of order
in which case both opening braces are highlighted in red: and the
system won't
accept a newline until the closures are all correctly matched.
Would anyone else find such a thing useful?
Could the ESS team make something like this happen?
On May 27, 2005, at 12:11 pm, Ajay Narottam Shah wrote:
I am using R 2.1 on Apple OS X.
When I get the prompt, I find it works well with emacs commandline
editing. Keys like M-f C-k etc. work fine.
The one thing that I really yearn for, which is missing, is bracket
matching When I am doing something which ends in it is really
useful to have emacs or vi-style bracket matching, so as to be able
to visually keep track of whether I have the correct matching
brackets, whether ( or { or [.
I'm sure this is possible. I will be most grateful if someone will
show the way :-) Thanks,
--
Ajay Shah Consultant
[EMAIL PROTECTED] Department of Economic Affairs
http://www.mayin.org/ajayshah Ministry of Finance, New Delhi
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
--
Robin Hankin
Uncertainty Analyst
National Oceanography Centre, Southampton
European Way, Southampton SO14 3ZH, UK
tel 023-8059-7743
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Power set
Laura Holt wrote: Hi again! I have a data.frame with the columns y, x1, x2, x3. I would like to fit linear models with one variable at a time, then 2 variables at a time, and then 3. Makes me think of a power set. Makes me think of irreproducible results if you use the output to select a single model :-) Frank Harrell Anyhow, is there a function to produce the right hand side of the formulas, please? thanks, Laura Holt mailto: [EMAIL PROTECTED] R Version 2.1.0 Windows -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Function environments lm() weights
I am writing a function of weighted regression, as a procedure for
heteroskedasticity.
The function runs an auxiliary regression whose fitted values I assign
to fit, and then I go:
w - 1/(exp(fit/2))
## Rerun the old regression ##
if(gls) {
wtd.model - glm(model, weights=w)
}
if(!gls) {
wtd.model - lm(model, weights=w, x=TRUE)
}
In this version, R complains that it can't find w. How can I tell it to
look for w in the function's environment, rather than in environment 1
or whatever?
An easy workaround, of course, is to superassign w and remove it
afterwards, but that's a little messy, in case the user already has a
variable called w in his environment.
Thanks,
Tobias Muhlhofer
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] R commandline editor question
well ESS has such a facility. However, I think Mathematica has a super scheme: unbalanced brackets show up in red, making them obvious. This is particularly good for spotting wrongly interleaved brackets, as in ([ blah di blah )] note bracket closure is out of order in which case both opening braces are highlighted in red: and the system won't accept a newline until the closures are all correctly matched. Would anyone else find such a thing useful? Could the ESS team make something like this happen? ess is great, but I was asking about the R commandline. I tend to write a lot of stuff on the fly at the R commandline. Yes, colours are a great way to deal with this, and this feature should ideally be in ESS. -- Ajay Shah Consultant [EMAIL PROTECTED] Department of Economic Affairs http://www.mayin.org/ajayshah Ministry of Finance, New Delhi __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Power set
Spencer Graves wrote: Hi, Laura: Have you considered regsubsets in library(leaps)? Also, have you done an R site search www.r-project.org - search - R site search for something like all subsets regression? Or much simpler if you are running R-2.1.0, use RSiteSearch(all subsets regression) Best, Jim -- James W. MacDonald Affymetrix and cDNA Microarray Facility University of Michigan Comprehensive Cancer Center 1500 E. Medical Center Drive Ann Arbor MI 48109 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Postscript
Philippe Lamy wrote: Hi, I would like to create a multi-page postscript file. How can I do that in R ? Is it possible ? Yes, simply draw more than one plot. See ?postscript and its argument onefile, which already defaults to TRUE. Uwe Ligges Thanks for help. Philippe __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Testing Nonlinear Restrictions
What kind of nonlinear restriction? Can you solve for one or more of the parameters in terms of the other(s) [either directly or implicitly]? If yes, then let fit1 - nls(... full model ... ) fit2 - nls(... restricted model ...) anova(fit1, fit2) If my memory is correct, Doug Bates, in his PhD dissertation ~25 years ago, decomposed the nonlinearity in nonlinear least squares into intrinsic curvature and parameter effects curvature. The Wald test is distorted by both sources types of nonlinearity, but the standard likelihood ratio anova is affected only by intrinsic curvature, and not parameter effects (provided the algorithm actually converges appropriately). Moreover, by reanalyzing a fair number of published data sets, Doug demostrated that in a nearly all practical application, the parameter effects curvature was much larger than the intrinsic curvature, and the latter was close to negligible in nearly all cases, while the parameter effects curvature was often of sufficient magnitude to substantively distort the answers. For more information, see Bates Watts (1988) Nonlinear Regression Analysis and Its Applications (Wiley). hope this helps. spencer graves Jacho-Chavez,DT (pgr) wrote: Dear all, I'm interested in testing 2 nonlinear restrictions on coefficients of a nls object. Is there a package for doing this? Something in the lines of `test(nls object, res=c(res 1,res 2),...)' I only found the function delta.method in the alr3 library that calculates the se of a singleton nonlinear restriction of a nls object using the delta method. Thanks in advanced for your help and suggestions. David __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Round a line
Luis Ridao Cruz wrote: R-help, I have lloked in the archives found no answer to how to round the line joint. I have usedthe arguments lnd, ljoin in par but I get no differences in the plotting. x=1:10 par(ljoin=round,lend=round) plot(x,sin(x),type=l,lwd=2) Any suggestions? Well, round is the default! You have to zoom in or make even thicker lines. Hence you might want to try out the folowing to see differences: x - 1:10 par(mfrow = c(1, 2)) plot(x, sin(x), type=l, lwd=10) par(ljoin=mitre, lend=butt) plot(x, sin(x), type=l, lwd=10) Uwe Ligges I run on a Windows XP machine. version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major2 minor1.0 year 2005 month04 day 18 language R Thank you in advance __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] r: LEXICOGRAPHIC ORDERING
HI all i have a seemingly simple question. given a sequence of numbers say, 1,2,3,4,5. i would like to get all of the possible two number arrangments (combinations), all 3 number arrangents ... 5 number arrangements (combinations). i.e. in the 2 number case: 12,13,14,15,23,24,25,34,35,45 any ideas?__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] how to get this kind of binomial distribution simulation number?
hai, I want to perform a simulation like this Suppose that I have one population ,it's size is 500, is composed of x,y and z. The probability of x, y and z is respectively is 0.3, 0.5, 0.2. I wan to simulate a new same size population based ratio of x, y and z, how can I get and assess the number of x, y and z. ? luan Key Laboratory for Sustainable Utilization of Marine Fisheries Resources, Ministry of Agriculture ,Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071,China __ G __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] longitudinal survey data
On 5/26/05, Thomas Lumley [EMAIL PROTECTED] wrote: If you *want* to fit mixed models (eg because you are interested in estimating variance components, or perhaps to gain efficiency) then it's quite a bit trickier. You can't just use the sampling weights in lme(). You can correct for the biased sampling if you put the variables that affect the weights in as predictors in the model. Cluster sampling could perhaps then be modelled as another level of random effect. I've been struggeling with case weights (in the case of unequal selection probabilities) in mixed effects models. Those are not possible in lme(). Isn't it, however, possible to use case weights in glmmPQL from MASS? Koen Pelleriaux Sociologist University of Antwerp __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] how to get this kind of binomial distribution simulation number?
rmultinom(n=1,size=500,prob=c(0.3,0.5,0.2)) to get n samples each of size 500 just use the n= argument. Reid Huntsinger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of luan_sheng Sent: Friday, May 27, 2005 9:26 AM To: r-help@stat.math.ethz.ch; R-help@stat.math.ethz.ch Subject: [R] how to get this kind of binomial distribution simulation number? hai, I want to perform a simulation like this Suppose that I have one population ,it's size is 500, is composed of x,y and z. The probability of x, y and z is respectively is 0.3, 0.5, 0.2. I wan to simulate a new same size population based ratio of x, y and z, how can I get and assess the number of x, y and z. ? luan Key Laboratory for Sustainable Utilization of Marine Fisheries Resources, Ministry of Agriculture ,Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071,China __ G __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] 3D density estimation with library sm - no estimate returned
Dear List, I have been trying to use library sm to do density estimation on a 3D dataset. I am using the current MacOS X binary of sm from CRAN. If I do this on a 2D dataset, sm.density returns a list including the component estimate which contains the density estimate over a uniform grid. When doing this with 3D data, although I get a nice plot (even when I don't ask for one), the returned list only contains the original data (see below). It would appear that the internal function sm.density.3d is returning NULL. Have I misunderstood what should be returned? Is this a platform specific problem? Can someone suggest a fix or an alternative library for my application? Very many thanks for your help, Greg Jefferis. R.version _ platform powerpc-apple-darwin7.9.0 arch powerpc os darwin7.9.0 system powerpc, darwin7.9.0 status Patched major2 minor1.0 year 2005 month05 day 12 language R library(sm) str(sm.density(matrix(rnorm(300),ncol=3),display=none)) List of 2 $ data:List of 3 ..$ x: num [1:100, 1:3] 0.9470 -1.5112 1.0589 -0.0884 -0.1900 ... .. ..- attr(*, dimnames)=List of 2 .. .. ..$ : NULL .. .. ..$ : chr [1:3] ..$ nbins: num 0 ..$ freq : num [1:100] 1 1 1 1 1 1 1 1 1 1 ... $ call: language sm.density(x = matrix(rnorm(300), ncol = 3), display = none) str(sm.density(matrix(rnorm(200),ncol=2),display=none)) List of 10 $ eval.points: num [1:50, 1:2] -2.67 -2.57 -2.47 -2.38 -2.28 ... ..- attr(*, dimnames)=List of 2 .. ..$ : NULL .. ..$ : chr [1:2] xnew ynew $ estimate : num [1:50, 1:50] 3.76e-05 5.10e-05 7.57e-05 1.22e-04 2.05e-04 ... $ h : Named num [1:2] 0.414 0.512 ..- attr(*, names)= chr [1:2] $ h.weights : num [1:100] 1 1 1 1 1 1 1 1 1 1 ... $ weights: num [1:100] 1 1 1 1 1 1 1 1 1 1 ... $ se : num [1:50, 1:50] 0.0306 0.0306 0.0306 0.0306 0.0306 ... $ upper : num [1:50, 1:50] 0.00454 0.00468 0.00489 0.00523 0.00571 ... $ lower : num [1:50, 1:50] 0 0 0 0 0 0 0 0 0 0 ... $ data :List of 3 ..$ x: num [1:100, 1:2] -0.694 -0.192 0.149 -0.718 -0.357 ... .. ..- attr(*, dimnames)=List of 2 .. .. ..$ : NULL .. .. ..$ : chr [1:2] ..$ nbins: num 0 ..$ freq : num [1:100] 1 1 1 1 1 1 1 1 1 1 ... $ call : language sm.density(x = matrix(rnorm(200), ncol = 2), display = none) -- Gregory Jefferis, PhD and: Research Fellow Department of Zoology St John's College Downing Street Cambridge Cambridge, CB2 3EJ CB2 1TP Tel: +44 (0)1223 336683 +44 (0)1223 339899 Fax: +44 (0)1223 336676 +44 (0)1223 337720 [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] longitudinal survey data
Thank you for your reply. Does that mean that in order to take in account the repeated measures I denote these as another cluster in R? Dassy Quoting Thomas Lumley [EMAIL PROTECTED]: On Thu, 26 May 2005 [EMAIL PROTECTED] wrote: Dear R-Users! Is there a possibility in R to do analyze longitudinal survey data (repeated measures in a survey)? I know that for longitudinal data I can use lme() to incorporate the correlation structure within individual and I know that there is the package survey for analyzing survey data. How can I combine both? I am trying to calculate design-based estimates. However, if I use svyglm() from the survey package I would ignore the correlation structure of the repeated measures. You *can* fit regression models to these data with svyglm(). Remember that from a design-based point of view there is no such thing as a correlation structure of repeated measures -- only the sampling is random, not the population data. If you *want* to fit mixed models (eg because you are interested in estimating variance components, or perhaps to gain efficiency) then it's quite a bit trickier. You can't just use the sampling weights in lme(). You can correct for the biased sampling if you put the variables that affect the weights in as predictors in the model. Cluster sampling could perhaps then be modelled as another level of random effect. -thomas Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] xmlAttrs and problems with reading node attributes of XML file (b ug?)
Hi,
Consider the following code:
require(XML)
xmlFile = paste( ?xml version=\1.0\
encoding=\ISO-8859-1\?\n,
mzXML xmlns=\a\ xmlns:xsi=\b\
xsi:schemaLocation=\c\\n,
parentFile a=\a\ b=\b\ /\n,
/mzXML\n)
cat(xmlFile)
a = function(x,...){
cat(Attributes of , xmlName(x), : );
print(xmlAttrs(x));
cat(\n);
NULL
}
xmlTreeParse(file=xmlFile, asText=TRUE,
handlers=list(startElement=a) )
And its output:
?xml version=1.0 encoding=ISO-8859-1?
mzXML xmlns=a xmlns:xsi=b xsi:schemaLocation=c
parentFile a=a b=b /
/mzXML
Attributes of parentFile : a b
a b
Attributes of mzXML : schemaLocation
c
It seems to me that XML parser was able to correctly list all the attributes
of the parentFile node but it failed on mzXML node.
Am I missusing the functions somehow or is it a bug?
Attribute names in mzXML node are part of mzXML file format and can not be
changed (removing all ':' would fix the problem) and I would like to store
than so whan I write mzXML file I have the same header.
Jarek
\===
Jarek Tuszynski, PhD. o / \
Science Applications International Corporation \__,|
(703) 676-4192 \
[EMAIL PROTECTED] `\
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] plotting box plots on same x
I am trying to construct a graph of 6 box plots of blood pressures. I want them to be on a single set of axis and I want the SBP to be ontop of the DBP. I have an array bp with the data in it and I tried a[1,]-c(145,60,147,62,140,57) a[2,]-c(160,75,160,74,160,70) a[3,]-c(140,55,140,65,142,55) boxplot(data.frame(a), main = Blood Pressures, at=c(1,1,2,2,3,3), names=c(sit,,lie,,stand,)) which is close to what I want, but it gives me a bunch of empty space at the end. is there a better way to do this to avoid this? As always, Thank You. ~Erithid __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Round a line
Luis Ridao Cruz [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] R-help, I have lloked in the archives found no answer to how to round the line joint. I have usedthe arguments lnd, ljoin in par but I get no differences in the plotting. x=1:10 par(ljoin=round,lend=round) plot(x,sin(x),type=l,lwd=2) On my Windows 2000 machine using R 2.1.0, par()$ljoin and par()$lend are already round by default. The par()$lmitre parameter is 10, Paul Murrell's article Fonts, lines, and transparency ... in R News 4/2 (Sept 2004) gives some clues under The end of the line: http://cran.stat.auckland.ac.nz/doc/Rnews/Rnews_2004-2.pdf All lines are drawn using a particular style for line ends and joins, though the difference only becomes obvious when lines become thick. Is it possible that with par()$lmitre at 10, and a lwd=2, you won't see any difference? efg __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] box plots without whiskers
I searched the archives, but couldnt find any way to do this with boxplot. Is there a way? Thanks again ~BJ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Chars as numbers
Josef,
Not sure if this is exactly what you mean, but there is a generic function
as()
to which you can then specify numeric as an argument and which then
coerces stuff into numeric format.
Tobias
Josef Eschgfaeller wrote:
Is there a proper function for transforming
a character to a number instead of using
i=match('c',letters)
# 3
Thanks.
Josef Eschgfäller
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
--
**
When Thomas Edison invented the light bulb he tried over 2000
experiments before he got it to work. A young reporter asked
him how it felt to have failed so many times. He said
I never failed once. I invented the light bulb.
It just happened to be a 2000-step process.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Chars as numbers
as.numeric() hih __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Chars as numbers
Josef Eschgfaeller wrote:
Is there a proper function for transforming
a character to a number instead of using
i=match('c',letters)
# 3
I'd suggest to use the above if you really mean it. Note that
transforming a character to a number is not well defined, because you
have to think about encodings and characters at first. See the article
by Brian Ripley in the most recent issue of R News, for example.
Uwe Ligges
Thanks.
Josef Eschgfäller
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] box plots without whiskers
BJ wrote: I searched the archives, but couldnt find any way to do this with boxplot. Is there a way? Thanks again ~BJ See ?boxplot and its argument par which points you to ?bxp: Now, you can do stuff like boxplot(1:10, pars=list(staplelty=0, whisklty=0)) Uwe Ligges __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Using R for classifying new samples
Hello, I do not have any statistical background, So I shall apologise if I am asking trivial question or help. I am trying to work with R. The problem I have on hand is: I have 2 sets of data(means SD for each sample in the group of both sets). The sample size is massive(2000+ in each grp). I have a new set of experimental data and I like to classify this to either of the grps based on statistical evidence. I thought using Linear discriminant analysis of R(MASS package) I can solve this problem. But unfortunately I am not quite sure its the right way and further with given my background I am not even sure how to prepare the data to test. To explain my problem further, I have prepared a small set of data to be testedso that it might be easy to suggest the right way to go ahead.The data is in the formObjectQuerymean_grp1SD_grp1mean_grp2SD_grp2 A3.8903.3151.1051.3950.345 B0.9151.1930.3340.6380.034 C2.0592.1550.6141.0420.159 D1.3720.9010.3140.3840.174 What I would like to test is if the query belongs to grp 1 or 2(for all my samples,2000+)and if it can read out to me teh query belongs to grp_1 or grp_2. Is there some good examples that I can refer to which can teach me step wise right from preparing the data to testing???or would there be someone who can help me with this problem? Thanks a lot for your time...look forward to hear some help and suggestions Manav __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] box plots without whiskers
Hmm ok, that seems to work with that example, but when I add an at= option i get graphs that are only verticle lines. Sorry to be difficult. ~BJ Uwe Ligges wrote: BJ wrote: I searched the archives, but couldnt find any way to do this with boxplot. Is there a way? Thanks again ~BJ See ?boxplot and its argument par which points you to ?bxp: Now, you can do stuff like boxplot(1:10, pars=list(staplelty=0, whisklty=0)) Uwe Ligges __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] box plots without whiskers
BJ wrote: Hmm ok, that seems to work with that example, but when I add an at= option i get graphs that are only verticle lines. Sorry to be difficult. He? So you need to pe much more precisely in this message as well! boxplot(data.frame(1:10, 2:11), at = 1:2, pars=list(staplelty=0, whisklty=0), at=1:2) works perfectly for me! People are not that happy to invest their time for helping you on bad specified problems... *REALLY*: PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Uwe Ligges ~BJ Uwe Ligges wrote: BJ wrote: I searched the archives, but couldnt find any way to do this with boxplot. Is there a way? Thanks again ~BJ See ?boxplot and its argument par which points you to ?bxp: Now, you can do stuff like boxplot(1:10, pars=list(staplelty=0, whisklty=0)) Uwe Ligges __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] plotting box plots on same x
On Fri, 2005-05-27 at 10:17 -0400, BJ wrote: I am trying to construct a graph of 6 box plots of blood pressures. I want them to be on a single set of axis and I want the SBP to be ontop of the DBP. I have an array bp with the data in it and I tried a[1,]-c(145,60,147,62,140,57) a[2,]-c(160,75,160,74,160,70) a[3,]-c(140,55,140,65,142,55) boxplot(data.frame(a), main = Blood Pressures, at=c(1,1,2,2,3,3), names=c(sit,,lie,,stand,)) which is close to what I want, but it gives me a bunch of empty space at the end. is there a better way to do this to avoid this? As always, Thank You. ~Erithid To answer both of your posts, use the following: # Review how your data is structured above. Your code for creating a # is not replicable. For those lacking clinical insight, explaining # your acronyms would also be helpful... :-) SBP - matrix(c(145, 160, 140, 147, 160, 140, 140, 160, 142), ncol = 3) DBP - matrix(c(60, 75, 55, 62, 74, 65, 57, 70, 55), ncol = 3) colnames(SBP) - colnames(DBP) - c(sit,lie,stand) # The key here is to only plot three at a time, lest boxplot() # default to a 'xlim' of 0.5 to 6.5 (1:# of groups +/- 0.5) # Then use 'add = TRUE' to plot the second group of 3 # Note also that I set the 'ylim' to the range of the combined # values in the first plot. boxplot(data.frame(SBP), main = Blood Pressures, ylim = range(c(SBP, DBP)), whisklty = 0, staplelty = 0) boxplot(data.frame(DBP), add = TRUE, whisklty = 0, staplelty = 0) Note the final two arguments, which result in the whiskers being drawn with an invisible line. See ?boxplot and ?bxp for more information. HTH, Marc Schwartz __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] Using R for classifying new samples
Read this book, Multivariate Statistical Analysis: A Conceptual Introduction by Sam Kash Kachigan. I think it's *great*, and perfect for someone without any statistical background. -Original Message- From: manav ram [mailto:[EMAIL PROTECTED] Sent: Friday, May 27, 2005 10:56 AM To: r-help@stat.math.ethz.ch Subject: [R] Using R for classifying new samples Hello, I do not have any statistical background, So I shall apologise if I am asking trivial question or help. I am trying to work with R. The problem I have on hand is: I have 2 sets of data(means SD for each sample in the group of both sets). The sample size is massive(2000+ in each grp). I have a new set of experimental data and I like to classify this to either of the grps based on statistical evidence. I thought using Linear discriminant analysis of R(MASS package) I can solve this problem. But unfortunately I am not quite sure its the right way and further with given my background I am not even sure how to prepare the data to test. To explain my problem further, I have prepared a small set of data to be testedso that it might be easy to suggest the right way to go ahead.The data is in the formObjectQuerymean_grp1SD_grp1mean_grp2SD_grp2 A3.8903.3151.1051.3950.345 B0.9151.1930.3340.6380.034 C2.0592.1550.6141.0420.159 D1.3720.9010.3140.3840.174 What I would like to test is if the query belongs to grp 1 or 2(for all my samples,2000+)and if it can read out to me teh query belongs to grp_1 or grp_2. Is there some good examples that I can refer to which can teach me step wise right from preparing the data to testing???or would there be someone who can help me with this problem? Thanks a lot for your time...look forward to hear some help and suggestions Manav __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] box plots without whiskers]
---BeginMessage--- boxplot(data.frame(1:10, 2:11), at = 1:2, + pars=list(staplelty=0, whisklty=0), at=1:2) Error in boxplot.default(data.frame(1:10, 2:11), at = 1:2, pars = list(staplelty = 0, : formal argument at matched by multiple actual arguments a-array(dim=c(3,8)) a[1,]-c(140,60,145,60,147,62,140,57) a[2,]-c(160,70,160,75,160,74,160,70) a[3,]-c(140,60,140,55,140,65,142,55) plot.new() boxplot(data.frame(a),pars=list(staplelty=0, whisklty=0), main = Boxplots, at=c(1,1,4,4,6,6,8,8), names=c(Base Line,,Sit,,Stand,,Lie,),ylab=DBP SBP) This does not work. The output is a collection of verticle lines that are as long as the boxes should be. boxplot(data.frame(a),staplelty=0, whisklty=0, main = Boxplots, at=c(1,1,4,4,6,6,8,8), names=c(Base Line,,Sit,,Stand,,Lie,),ylab=DBP SBP) does though. So I thank you for your help. Didnt mean to get yellled at over what seemed like a decent question. The mailing list archives say you have to use the lattice package to achieve the effect I wanted.. Anyhow. Sorry to everyone else for the excessive trafffic. ~Erithid Uwe Ligges wrote: BJ wrote: Hmm ok, that seems to work with that example, but when I add an at= option i get graphs that are only verticle lines. Sorry to be difficult. He? So you need to pe much more precisely in this message as well! boxplot(data.frame(1:10, 2:11), at = 1:2, pars=list(staplelty=0, whisklty=0), at=1:2) works perfectly for me! People are not that happy to invest their time for helping you on bad specified problems... *REALLY*: PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Uwe Ligges ~BJ Uwe Ligges wrote: BJ wrote: I searched the archives, but couldnt find any way to do this with boxplot. Is there a way? Thanks again ~BJ See ?boxplot and its argument par which points you to ?bxp: Now, you can do stuff like boxplot(1:10, pars=list(staplelty=0, whisklty=0)) Uwe Ligges __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html ---End Message--- __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] plotting box plots on same x
BJ-
For the record as well I could not make heads or tails of your question.
The dedicated folks who day in, day out field poorly specified questions have
no obligation to do so, and have every right to get frustrated when it seems
that their time is being taken advantage of.
On Fri, 27 May 2005, BJ wrote:
Date: Fri, 27 May 2005 11:01:39 -0400
From: BJ [EMAIL PROTECTED]
To: Uwe Ligges [EMAIL PROTECTED]
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] plotting box plots on same x
Does it matter what they are? they are just names. The six box plots are
from the array I created with columns. I forgot to add teh dim
statement, my mistake. But I thought it was obvious that i was using a
matrix from the data frame call and the assignments I provided. My
question was simply about setting the x axis so that it stopped after
x=3, even though I had 6 plots. For teh record, before I get jumped on
for the statistics, I am just using a 3 x 6 matrix to test the code
before applying it to actual data. Also, I was not being scarcastic. I
have recieved a lot of help from this mailing list. The R documentation
is hard to hunt down and not complete. Without the help of actual people
I would be dead in the water. I am sorry for any hostility that I have
incurred. ~Erithid
Uwe Ligges wrote:
BJ wrote:
I am trying to construct a graph of 6 box plots of blood pressures. I
want them to be on a single set of axis and I want the SBP to be
ontop of the DBP. I have an array bp with the data in it and I tried
Folks, please invest 1 minute of time to rethink whether other people
will understand your question!
What is SBP, DBP, and where are the 6 boxplots from?
a[1,]-c(145,60,147,62,140,57)
a must already be defined here, or we cannot replace a column!
a[2,]-c(160,75,160,74,160,70)
a[3,]-c(140,55,140,65,142,55)
boxplot(data.frame(a), main = Blood Pressures, at=c(1,1,2,2,3,3),
names=c(sit,,lie,,stand,))
which is close to what I want, but it gives me a bunch of empty space
at the end. is there a better way to do this to avoid this?
Well, the first point is that you should write questions that you
would understand yourself. In a next step you might find someone who
is able to answer it
As always, Thank You. ~Erithid
As always? Does not sound very enthusiastic ...
Uwe Ligges
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
J.R. Lockwood
412-683-2300 x4941
[EMAIL PROTECTED]
http://www.rand.org/statistics/bios/
This email message is for the sole use of the intended recip...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] I never made any assumption athat anyone had any obligation to do anything
I just asked a question. If I was too vague, then i am sorry. I dont expect anyone to help me, but I thought that it was ok to put the question out there in case someone wanted to help me. I didnt expect abject hostility for it. Human decency was the only thing I did expect. If my question was a pain,badly formatted, or too juvinile, then i have no problem being ignored.I asked a question that the archives did not satisfactorily answer. I appreciate the help of the poeple on this list as they are an invaluable resource, especially since the R documentation is sketchy at times. I am sorry for wasting everyones time. I just dont like being yelled at first thing in the morning for asking for help on a help list. See you all around ~Erithid __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] longitudinal survey data
On Fri, 27 May 2005 [EMAIL PROTECTED] wrote: Thank you for your reply. Does that mean that in order to take in account the repeated measures I denote these as another cluster in R? Yes, but unless you have multistage finite population corrections to put in the design object only the first stage of clustering affects the results, so you may not need to bother. -thomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Soil texture triangle in R?
Right,
Got the data points plotted on top of the soil texture background,
thanks to Jim and ternaryplot{vcd}! See code below.
Now there is some fine tuning to do, as it should really look like this
graph:
http://soil.scijournals.org/content/vol65/issue4/images/large/1038f2.jpeg
Things to do:
- rotate axis labels;
- correct small errors in class divisions;
- correct the partial covering of the bottom tick labels;
- rotate ticks in order to simplify viewing the graph.
Any help still appreciated!
Cheers,
Sander.
soil.triangle - function() {
oldpar - par(no.readonly=TRUE)
## now the bottom internal ticks
x1-seq(0.1,0.9,by=0.1)
x2-x1
y1-rep(0,9)
y2-rep(-0.02,9)
segments(x1,y1,x2,y2)
text(x1,y1-0.03,as.character(rev(seq(10,90,by=10 #, cex=0.8)
## now the left internal ticks
y1-x1*sin60
x1-x1*0.5
x2-x1+0.02*sin60
y2-y1-0.02*0.5
segments(x1,y1,x2,y2)
text(x1-0.03,y1+0.015,as.character(seq(10,90,by=10)))
## now the right internal ticks
x1-rev(x1+0.5-0.02*sin60)
x2-x1+0.02*sin60
segments(x1,y2,x2,y1)
text(x2+0.03,y1+0.015,as.character(rev(seq(10,90,by=10
## the labels at the corners
par(xpd=TRUE)
# text(0.5,0.9,100% clay)
# text(-0.1,0,100% sand)
# text(1.1,0,100% loam)
text(0.09,0.43,% Clay)
text(0.90,0.43,% Silt)
text(0.5,-0.1,% Sand)
# boundary of clay with extensions
x1-c(0.275,0.35,0.6)
x2-c(0.4,0.79,0.7)
y1-c(0.55*sin60,0.41*sin60,0.41*sin60)
y2-c(0.285*sin60,0.41*sin60,0.6*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.5,0.57,Clay, col=grey)
# lower bound of clay loam silty divider
x1-c(0.4,0.68)
x2-c(0.86,0.6)
y1-c(0.285*sin60,0.285*sin60)
y2-c(0.285*sin60,0.41*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.7,0.49*sin60,Silty, col=grey)
text(0.7,0.44*sin60,clay, col=grey)
text(0.73,0.37*sin60,Silty clay, col=grey)
text(0.73,0.33*sin60,loam, col=grey)
text(0.5,0.35*sin60,Clay loam, col=grey)
x1-c(0.185,0.1,0.37)
x2-c(0.36,0.37,0.4)
y1-c(0.37*sin60,0.2*sin60,0.2*sin60)
y2-c(0.37*sin60,0.2*sin60,0.285*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.27,0.43*sin60,Sandy, col=grey)
text(0.27,0.39*sin60,clay, col=grey)
text(0.27,0.3*sin60,Sandy clay, col=grey)
text(0.27,0.26*sin60,loam, col=grey)
# sand corner
x1-c(0.05,0.075)
x2-c(0.12,0.3)
y1-c(0.1*sin60,0.15*sin60)
y2-c(0,0)
segments(x1,y1,x2,y2, col=grey)
text(0.25,0.13*sin60,Sandy loam, col=grey)
text(0.13,0.075*sin60,Loamy, col=grey)
text(0.15,0.035*sin60,sand, col=grey)
text(0.055,0.021,Sand, col=grey)
x1-c(0.37,0.42,0.5,0.8,0.86)
x2-c(0.42,0.54,0.65,0.86,0.94)
y1-c(0.2*sin60,0.08*sin60,0,0,0.12*sin60)
y2-c(0.08*sin60,0.08*sin60,0.285*sin60,0.12*sin60,0.12*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.49,0.18*sin60,Loam, col=grey)
text(0.72,0.15*sin60,Silt loam, col=grey)
text(0.9,0.06*sin60,Silt, col=grey)
par(oldpar)
}
tmp - array(dim=c(10,3))
tmp[,1] - abs(rnorm(10)*20)
tmp[,2] - abs(rnorm(10)*10)
tmp[,3] - 100-tmp[,1]-tmp[,2]
tmp
library(vcd)
## Mark groups
ternaryplot(tmp,
grid=FALSE,
dimnames.position = none,
pch=1, col=black,
scale=1, main=NULL,
prop.size=FALSE,
)
soil.triangle()
Sander Oom wrote:
Hi Jim,
This looks impressive! It gives me the 'background' graph. However, I'm
not sure how I can use this function to plot my soil texture values! Can
you explain?
I would like to be able to plot my soil texture samples in the same
graph as the one your function plots.
Maybe I should try to figure out how to replicate your code inside a
ternaryplot{vcd} call.
Cheers,
Sander.
Jim Lemon wrote:
Sander Oom wrote:
Dear R users,
has anybody made an attempt to create the soil texture triangle graph
in R? For an example see here:
http://www.teachingkate.org/images/soiltria.gif
I would like to get the lines in black and texture labels in gray to
allow for plotting my texture results on top.
Any examples or suggestions are very welcome!
It's not too hard to write a plot function to do this, but I'm not sure
whether this will be what you want. Anyway, try it out.
Jim
soil.triangle-function() {
oldpar-par(no.readonly=TRUE)
plot(0:1,type=n,axes=FALSE,xlim=c(0,1.1),ylim=c(0,1),
main=Soil Triangle,xlab=,ylab=)
# first draw the triangle
x1-c(0,0,0.5)
sin60-sin(pi/3)
x2-c(1,0.5,1)
y1-c(0,0,sin60)
y2-c(0,sin60,0)
segments(x1,y1,x2,y2)
# now the bottom internal ticks
x1-seq(0.1,0.9,by=0.1)
x2-x1
y1-rep(0,9)
y2-rep(0.02,9)
segments(x1,y1,x2,y2)
text(x1,y1-0.03,as.character(rev(seq(10,90,by=10
# now the left internal ticks
y1-x1*sin60
x1-x1*0.5
x2-x1+0.02*sin60
y2-y1-0.02*0.5
segments(x1,y1,x2,y2)
text(x1-0.03,y1+0.015,as.character(seq(10,90,by=10)))
x1-rev(x1+0.5-0.02*sin60)
x2-x1+0.02*sin60
segments(x1,y2,x2,y1)
text(x2+0.03,y1+0.015,as.character(rev(seq(10,90,by=10
[R] nlminb to optmin
Hi! I want to convert S-Plus 6.2 code to R 2.1.0. Instead of the function nlminb I use the function optmin optmin(start,fn,gr,method=L-BFGS-B, lower, upper, hess,...) But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn Then I used optmin(start,fn,gr,method=BFGS, hess, ...) But then I get the Error in optmin ...: initial value in vmmin is not finite I know the final parameter estimates from S-Plus which I use as starting values in R. The upper and lower bounds are close around the final estimates. So there is not much to maximize. What can I do? Thank you for help, Peter [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] R commandline editor question
Of course it should be in ESS. In fact, it is for me, anyway
(different colors depending on whether they match or not).
Check out the paren-hilit or paren-match (or something like that)
customize options.
I.e.
M-x customize-groups ret paren spc
(which ought to complete on names of groups/options starting with
paren) or blink-paren.
(I'm not on an Emacs-enabled computer right now).
On 5/27/05, Robin Hankin [EMAIL PROTECTED] wrote:
Hi Ajay
well ESS has such a facility.
However, I think Mathematica has a super scheme: unbalanced brackets
show up
in red, making them obvious.
This is particularly good for spotting wrongly interleaved brackets, as
in
([ blah di blah )]
note bracket closure is out of order
in which case both opening braces are highlighted in red: and the
system won't
accept a newline until the closures are all correctly matched.
Would anyone else find such a thing useful?
Could the ESS team make something like this happen?
On May 27, 2005, at 12:11 pm, Ajay Narottam Shah wrote:
I am using R 2.1 on Apple OS X.
When I get the prompt, I find it works well with emacs commandline
editing. Keys like M-f C-k etc. work fine.
The one thing that I really yearn for, which is missing, is bracket
matching When I am doing something which ends in it is really
useful to have emacs or vi-style bracket matching, so as to be able
to visually keep track of whether I have the correct matching
brackets, whether ( or { or [.
I'm sure this is possible. I will be most grateful if someone will
show the way :-) Thanks,
--
Ajay Shah Consultant
[EMAIL PROTECTED] Department of Economic Affairs
http://www.mayin.org/ajayshah Ministry of Finance, New Delhi
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
--
Robin Hankin
Uncertainty Analyst
National Oceanography Centre, Southampton
European Way, Southampton SO14 3ZH, UK
tel 023-8059-7743
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
--
best,
-tony
Commit early,commit often, and commit in a repository from which we can easily
roll-back your mistakes (AJR, 4Jan05).
A.J. Rossini
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] nlminb to optmin
Stefan Pohl wrote: Hi! I want to convert S-Plus 6.2 code to R 2.1.0. Instead of the function nlminb I use the function optmin optmin(start,fn,gr,method=L-BFGS-B, lower, upper, hess,...) But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn Then I used optmin(start,fn,gr,method=BFGS, hess, ...) But then I get the Error in optmin ...: initial value in vmmin is not finite I know the final parameter estimates from S-Plus which I use as starting values in R. The upper and lower bounds are close around the final estimates. So there is not much to maximize. What can I do? Thank you for help, Peter What is optmin? Do you mean optim? Either way, you can always try your function at the initial values outside of optim. If it returns Inf, you have a problem with your objective function. Have you considered that your objective function in S-PLUS didn't port the way you expected to R? I.e. if you call your objective function in S-PLUS and then in R, using the same inputs, do you get identical outputs? --sundar __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] nlminb to optmin
Stefan Pohl wrote: Hi! I want to convert S-Plus 6.2 code to R 2.1.0. Instead of the function nlminb I use the function optmin optmin(start,fn,gr,method=L-BFGS-B, lower, upper, hess,...) But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn Then I used optmin(start,fn,gr,method=BFGS, hess, ...) But then I get the Error in optmin ...: initial value in vmmin is not finite I know the final parameter estimates from S-Plus which I use as starting values in R. The upper and lower bounds are close around the final estimates. So there is not much to maximize. What can I do? Thank you for help, I have a test version of a package available as http://www.stat.wisc.edu/~bates/port_0.1-1.tar.gz that provides nlminb for R. If you can install packages from source code then you may want to try that. If you are running under Windows and only install binary packages then we will need to ask for a volunteer to create a Windows binary from the source package. I do not plan to upload this package to CRAN. Instead I plan to incorporate nlminb into r-devel in time to have it become part of R-2.2.0 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] longitudinal survey data
Sorry, still confused. If I dont have fpc's ready in my dataset (calculate myself?) that means that R will use the weight of an individual for each of his repeated observations. But is that then still correct? The cluster individual is ignored and each observation of an individual has the same weight. Thanks a lot. Dassy Quoting Thomas Lumley [EMAIL PROTECTED]: On Fri, 27 May 2005 [EMAIL PROTECTED] wrote: Thank you for your reply. Does that mean that in order to take in account the repeated measures I denote these as another cluster in R? Yes, but unless you have multistage finite population corrections to put in the design object only the first stage of clustering affects the results, so you may not need to bother. -thomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Soil texture triangle in R?
Cleaned up the class divisions and created a full function.
Still to do:
- rotate axis labels;
- correct the partial covering of the bottom tick labels;
- rotate ticks in order to simplify viewing the graph.
See:
http://soil.scijournals.org/content/vol65/issue4/images/large/1038f2.jpeg
Wonder whether triangle.plot{ade4} will give more flexibility!?
Anyway, hopefully the result so far is useful for other people.
Cheers,
Sander.
plot.soiltexture - function(x,pch,col) {
## triangle plots:
## triangle.plot {ade4}
## triplot{klaR}
## ternaryplot {vcd}
require(vcd)
require(Zelig)
ternaryplot(x,
grid=FALSE,
dimnames.position = none,
pch=pch, col=col,
scale=1, main=NULL,
prop.size=FALSE
)
oldpar - par(no.readonly=TRUE)
ticlength - 0.01
## now the bottom internal ticks
x1-seq(0.1,0.9,by=0.1)
x2-x1
y1-rep(0,9)
y2-rep(ticlength,9)
segments(x1,y1,x2,y2)
text(x1,y1-0.03,as.character(rev(seq(10,90,by=10 #, cex=0.8)
## now the left internal ticks
y1-x1*sin60
x1-x1*0.5
x2-x1+ticlength*sin60
y2-y1-ticlength*0.5
segments(x1,y1,x2,y2)
text(x1-0.03,y1+0.015,as.character(seq(10,90,by=10)))
## now the right internal ticks
x1-rev(x1+0.5-ticlength*sin60)
x2-x1+ticlength*sin60
segments(x1,y2,x2,y1)
text(x2+0.03,y1+0.015,as.character(rev(seq(10,90,by=10
## the labels at the corners
par(xpd=TRUE)
# text(0.5,0.9,100% clay)
# text(-0.1,0,100% sand)
# text(1.1,0,100% loam)
## the axis labels
text(0.09,0.43,% Clay)
text(0.90,0.43,% Silt)
text(0.5,-0.1,% Sand)
# boundary of clay with extensions
x1-c(0.275,0.355,0.6)
x2-c(0.415,0.8,0.7)
y1-c(0.55*sin60,0.4*sin60,0.4*sin60)
y2-c(0.285*sin60,0.4*sin60,0.6*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.5,0.57,Clay, col=grey)
# lower bound of clay loam silty divider
x1-c(0.415,0.66)
x2-c(0.856,0.6)
y1-c(0.285*sin60,0.285*sin60)
y2-c(0.285*sin60,0.40*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.7,0.49*sin60,Silty, col=grey)
text(0.7,0.44*sin60,clay, col=grey)
text(0.72,0.36*sin60,Silty clay, col=grey)
text(0.73,0.32*sin60,loam, col=grey)
text(0.5,0.35*sin60,Clay loam, col=grey)
x1-c(0.185,0.1,0.37)
x2-c(0.37,0.37,0.415)
y1-c(0.37*sin60,0.2*sin60,0.2*sin60)
y2-c(0.37*sin60,0.2*sin60,0.285*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.28,0.43*sin60,Sandy, col=grey)
text(0.27,0.39*sin60,clay, col=grey)
text(0.27,0.3*sin60,Sandy clay, col=grey)
text(0.27,0.26*sin60,loam, col=grey)
# sand corner
x1-c(0.05,0.075)
x2-c(0.15,0.3)
y1-c(0.1*sin60,0.15*sin60)
y2-c(0,0)
segments(x1,y1,x2,y2, col=grey)
text(0.25,0.13*sin60,Sandy loam, col=grey)
text(0.14,0.07*sin60,Loamy, col=grey)
text(0.18,0.03*sin60,sand, col=grey)
text(0.06,0.021,Sand, col=grey)
x1-c(0.37,0.435,0.5,0.8,0.86)
x2-c(0.435,0.537,0.64,0.86,0.94)
y1-c(0.2*sin60,0.08*sin60,0,0,0.12*sin60)
y2-c(0.08*sin60,0.08*sin60,0.285*sin60,0.12*sin60,0.12*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.49,0.18*sin60,Loam, col=grey)
text(0.72,0.15*sin60,Silt loam, col=grey)
text(0.9,0.06*sin60,Silt, col=grey)
ternarypoints(x, pch = pch, col = col)
par(oldpar)
}
tmp - array(dim=c(10,3))
tmp[,2] - abs(rnorm(10)*20)
tmp[,3] - abs(rnorm(10)*10)
tmp[,1] - 100-tmp[,2]-tmp[,3]
col - rep(black,10)
pch - rep(1, 10)
plot.soiltexture(tmp,pch,col=black)
Sander Oom wrote:
Right,
Got the data points plotted on top of the soil texture background,
thanks to Jim and ternaryplot{vcd}! See code below.
Now there is some fine tuning to do, as it should really look like this
graph:
http://soil.scijournals.org/content/vol65/issue4/images/large/1038f2.jpeg
Things to do:
- rotate axis labels;
- correct small errors in class divisions;
- correct the partial covering of the bottom tick labels;
- rotate ticks in order to simplify viewing the graph.
Any help still appreciated!
Cheers,
Sander.
soil.triangle - function() {
oldpar - par(no.readonly=TRUE)
## now the bottom internal ticks
x1-seq(0.1,0.9,by=0.1)
x2-x1
y1-rep(0,9)
y2-rep(-0.02,9)
segments(x1,y1,x2,y2)
text(x1,y1-0.03,as.character(rev(seq(10,90,by=10 #, cex=0.8)
## now the left internal ticks
y1-x1*sin60
x1-x1*0.5
x2-x1+0.02*sin60
y2-y1-0.02*0.5
segments(x1,y1,x2,y2)
text(x1-0.03,y1+0.015,as.character(seq(10,90,by=10)))
## now the right internal ticks
x1-rev(x1+0.5-0.02*sin60)
x2-x1+0.02*sin60
segments(x1,y2,x2,y1)
text(x2+0.03,y1+0.015,as.character(rev(seq(10,90,by=10
## the labels at the corners
par(xpd=TRUE)
# text(0.5,0.9,100% clay)
# text(-0.1,0,100% sand)
# text(1.1,0,100% loam)
text(0.09,0.43,% Clay)
text(0.90,0.43,% Silt)
text(0.5,-0.1,% Sand)
# boundary of clay with extensions
x1-c(0.275,0.35,0.6)
x2-c(0.4,0.79,0.7)
y1-c(0.55*sin60,0.41*sin60,0.41*sin60)
y2-c(0.285*sin60,0.41*sin60,0.6*sin60)
segments(x1,y1,x2,y2, col=grey)
text(0.5,0.57,Clay, col=grey)
# lower bound of clay loam silty divider
[R] Windows binary version of port_0.1-1 available
In an earlier message today I mentioned that a source package with a version of nlminb for R was available as http://www.stat.wisc.edu/~bates/port_0.1-1.tar.gz Thanks to Kjetil Halvorsen there is now a Windows binary version available as http://www.stat.wisc.edu/~bates/port_0.1-1.zip __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] I never made any assumption athat anyone had any obligation to do anything
Hello, Erithid: I'm very sorry you feel you got yelled at first thing in the morning. Please try not to take it personally, though I know that may be difficult. After one superficially insulting reply, I just laughed and told my manager, Look what I learned for exposing myself to ridicule! My philosophy on this was expressed last December and subsequently added to the fortunes pacakage, from which I will quote for you now: library(fortunes) fortune(Graves) Our great-great grandchilren as yet unborn may read some of the stupid questions and/or answers that I and perhaps others give from time to time. I'd rather get flamed for saying something stupid in public on this list than to continue to provide substandard service to the people with whom I work because I perpetrated the same mistake in an environment in which no one questioned so effectively my errors. -- Spencer Graves (in a discussion on whether answers on R-help should be more polite) R-help (December 2004) I've heard from some of the best on this listserve that it happens to everyone. That doesn't make it right. If you can learn to accept this as a different subculture operating by different rules of diplomacy, you might get more out of this list, do better work with whatever you are attempting to do, and enjoy life more. (My wife's middle name is Rose. When she gives me a hard time about something, I just try to listen and console myself with the thought that, If I want to sleep in a bed of roses, I must get used to the thorns.) Best Wishes, spencer graves BJ wrote: I just asked a question. If I was too vague, then i am sorry. I dont expect anyone to help me, but I thought that it was ok to put the question out there in case someone wanted to help me. I didnt expect abject hostility for it. Human decency was the only thing I did expect. If my question was a pain,badly formatted, or too juvinile, then i have no problem being ignored.I asked a question that the archives did not satisfactorily answer. I appreciate the help of the poeple on this list as they are an invaluable resource, especially since the R documentation is sketchy at times. I am sorry for wasting everyones time. I just dont like being yelled at first thing in the morning for asking for help on a help list. See you all around ~Erithid __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] images and maps in R
Hi:
I have a question arising from my project.
A sample of the data is below. The first row stands
for the names of state in USA. The second row stand
for some numeric value in that state. Some of them are
NA. I can use the commands data(stateMapEnv) and
map('state', fill = F) in library maps to make a
plot of USA states. What I want to do is: 1. put the
corresponding state name on the Map 2. give different
state different colors which is related to their
value. For example, red for values ranging from
0-30, green for values from 80-90, etc. 3. if
possible, put the value of each state within the state
on the map. here, take the numeric value as some text.
Thank you very much!
AB AK AL AR AZ CT CA
91 80 NA NA 17 33 20
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] installing spatstat in OSX
Friends, I am trying to install the current version of spatstat on Mac OS 10.3.9, but the compilation fails with the following messages at the end: ld: warning -L: directory name (/usr/local/lib/gcc/powerpc-apple-darwin6.8/3.4.2) does not exist ld: can't locate file for: -lg2c make: *** [spatstat.so] Error 1 ERROR: compilation failed for package 'spatstat' The current version of sm (which spatstat needs) also fails to compile when I try that separately, with the same error message. I remember that binaries of these packages were formerly available for OSX, but not anymore. I would appreciate your help with installing these. Thanks, -Robert __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] images and maps in R
This is is not difficult from the online help for the map() function.
Here is an example.
map('state', region = c('new york', 'new jersey', 'penn'),fill=TRUE,col=1:4)
There is also an example there in how to add text to the map.
Another way uses the maptools package.
require(maptools)
?plot.Map
and then follow the examples given in the help page for plot.Map().
At 2:42 PM -0700 5/27/05, yyan liu wrote:
Hi:
I have a question arising from my project.
A sample of the data is below. The first row stands
for the names of state in USA. The second row stand
for some numeric value in that state. Some of them are
NA. I can use the commands data(stateMapEnv) and
map('state', fill = F) in library maps to make a
plot of USA states. What I want to do is: 1. put the
corresponding state name on the Map 2. give different
state different colors which is related to their
value. For example, red for values ranging from
0-30, green for values from 80-90, etc. 3. if
possible, put the value of each state within the state
on the map. here, take the numeric value as some text.
Thank you very much!
AB AK AL AR AZ CT CA
91 80 NA NA 17 33 20
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
--
--
Don MacQueen
Environmental Protection Department
Lawrence Livermore National Laboratory
Livermore, CA, USA
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] longitudinal survey data
On Fri, 27 May 2005 [EMAIL PROTECTED] wrote: Sorry, still confused. If I dont have fpc's ready in my dataset (calculate myself?) that means that R will use the weight of an individual for each of his repeated observations. But is that then still correct? The cluster individual is ignored and each observation of an individual has the same weight. Well, it depends to some extent on what inferences you are making, but yes, you probably do want each observation to have the same weight. Suppose you have 4 measurements on each person, and you are working with a simple random sample of 1000 people from a population of 1,000,000. If you had done these 4 measurements on the whole population you would have 4,000,000 measurements, so the 4000 measurements you have are 1/1000 of the population. This is the same weighting as if you had a single measurement person person, giving 1000 measurements in the sample and 1,000,000 in the population. If different individuals have different numbers of measurements then things get a bit trickier. It depends then on why there are different numbers of measurements.If they are the result of non-response you might want to rescale the weights at later time points to give the right population totals. If they are part of the sampling design then the design will specify what to do with them. -thomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] I never made any assumption athat anyone had any obligation to do anything
To add to Spencer's reply, I have two things. First, if you look back at your earlier list entries, Uwe asked you politely at least twice to read the posting guide. The reason for this is that the posting guide gives guidance to allow you to post questions that are clear enough to answer, and points you to other sources of information that you might not already know about. It was only after you appeared to ignore his advice that he 'yelled' at you. Second, you ultimately learned several things. As far as I can tell, Uwe did answer your questions, but more importantly, you learned how to interact on a list where you get free advice on how to use free software (the operative term here being *free*). Isn't that better than being ignored and not learning anything? Best, Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 Spencer Graves [EMAIL PROTECTED] 05/27/05 5:34 PM Hello, Erithid: I'm very sorry you feel you got yelled at first thing in the morning. Please try not to take it personally, though I know that may be difficult. After one superficially insulting reply, I just laughed and told my manager, Look what I learned for exposing myself to ridicule! My philosophy on this was expressed last December and subsequently added to the fortunes pacakage, from which I will quote for you now: library(fortunes) fortune(Graves) Our great-great grandchilren as yet unborn may read some of the stupid questions and/or answers that I and perhaps others give from time to time. I'd rather get flamed for saying something stupid in public on this list than to continue to provide substandard service to the people with whom I work because I perpetrated the same mistake in an environment in which no one questioned so effectively my errors. -- Spencer Graves (in a discussion on whether answers on R-help should be more polite) R-help (December 2004) I've heard from some of the best on this listserve that it happens to everyone. That doesn't make it right. If you can learn to accept this as a different subculture operating by different rules of diplomacy, you might get more out of this list, do better work with whatever you are attempting to do, and enjoy life more. (My wife's middle name is Rose. When she gives me a hard time about something, I just try to listen and console myself with the thought that, If I want to sleep in a bed of roses, I must get used to the thorns.) Best Wishes, spencer graves BJ wrote: I just asked a question. If I was too vague, then i am sorry. I dont expect anyone to help me, but I thought that it was ok to put the question out there in case someone wanted to help me. I didnt expect abject hostility for it. Human decency was the only thing I did expect. If my question was a pain,badly formatted, or too juvinile, then i have no problem being ignored.I asked a question that the archives did not satisfactorily answer. I appreciate the help of the poeple on this list as they are an invaluable resource, especially since the R documentation is sketchy at times. I am sorry for wasting everyones time. I just dont like being yelled at first thing in the morning for asking for help on a help list. See you all around ~Erithid __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html ** Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] I never made any assumption athat anyone had any obligation to do anything
Perhaps the list would work better if posters who do not follow the posting guide are nicely told that they will more likely get helpful replies if they read and follow the posting guide and, in particular, provide a reproducible example. If the poster does not follow that advice then the poster can just be ignored. This thread pointed out that even though the original poster did not follow the guide and ignored direct suggestions that he still got his answer. The problem is that he felt insulted in the process. If the above procedure were followed he might not have received an answer to his problem but there would have been less animosity all around and maybe that's more important. Perhaps the posting guide could even suggest such an approach. On 5/27/05, James MacDonald [EMAIL PROTECTED] wrote: To add to Spencer's reply, I have two things. First, if you look back at your earlier list entries, Uwe asked you politely at least twice to read the posting guide. The reason for this is that the posting guide gives guidance to allow you to post questions that are clear enough to answer, and points you to other sources of information that you might not already know about. It was only after you appeared to ignore his advice that he 'yelled' at you. Second, you ultimately learned several things. As far as I can tell, Uwe did answer your questions, but more importantly, you learned how to interact on a list where you get free advice on how to use free software (the operative term here being *free*). Isn't that better than being ignored and not learning anything? Best, Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 Spencer Graves [EMAIL PROTECTED] 05/27/05 5:34 PM Hello, Erithid: I'm very sorry you feel you got yelled at first thing in the morning. Please try not to take it personally, though I know that may be difficult. After one superficially insulting reply, I just laughed and told my manager, Look what I learned for exposing myself to ridicule! My philosophy on this was expressed last December and subsequently added to the fortunes pacakage, from which I will quote for you now: library(fortunes) fortune(Graves) Our great-great grandchilren as yet unborn may read some of the stupid questions and/or answers that I and perhaps others give from time to time. I'd rather get flamed for saying something stupid in public on this list than to continue to provide substandard service to the people with whom I work because I perpetrated the same mistake in an environment in which no one questioned so effectively my errors. -- Spencer Graves (in a discussion on whether answers on R-help should be more polite) R-help (December 2004) I've heard from some of the best on this listserve that it happens to everyone. That doesn't make it right. If you can learn to accept this as a different subculture operating by different rules of diplomacy, you might get more out of this list, do better work with whatever you are attempting to do, and enjoy life more. (My wife's middle name is Rose. When she gives me a hard time about something, I just try to listen and console myself with the thought that, If I want to sleep in a bed of roses, I must get used to the thorns.) Best Wishes, spencer graves BJ wrote: I just asked a question. If I was too vague, then i am sorry. I dont expect anyone to help me, but I thought that it was ok to put the question out there in case someone wanted to help me. I didnt expect abject hostility for it. Human decency was the only thing I did expect. If my question was a pain,badly formatted, or too juvinile, then i have no problem being ignored.I asked a question that the archives did not satisfactorily answer. I appreciate the help of the poeple on this list as they are an invaluable resource, especially since the R documentation is sketchy at times. I am sorry for wasting everyones time. I just dont like being yelled at first thing in the morning for asking for help on a help list. See you all around ~Erithid __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html ** Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues. __
