Hi,
Your idea reminds me of an example in the documentation of the brew package,
featuring the generation of a template. You might want to check it out.
baptiste
On Wed, Nov 12, 2008 at 2:01 PM, Kem Phillips <[EMAIL PROTECTED]>wrote:
> Dirk,
>
>
>
> I came upon your message below in searching f
Dear R help team,
I am sorry to put urgent in the subject but to be honest I am in big
troubles if I cannot make R to work in my computer.
I am a PhD student in statistics and most of my homeworks and projects are
in R.
I have been using R in my mac book pro, mac os X 10.5.5 without problems
unti
Daniel Fernandez wrote:
> Dear R help team,
> I am sorry to put urgent in the subject but to be honest I am in big
> troubles if I cannot make R to work in my computer.
>
> I am a PhD student in statistics and most of my homeworks and projects are
> in R.
>
> I have been using R in my mac book pr
Dear R help team,
I am sorry to put urgent in the subject but to be honest I am in big
troubles if I cannot make R to work in my computer.
I am a PhD student in statistics and most of my homeworks and projects are
in R.
I have been using R in my mac book pro, mac os X 10.5.5 without problems
unti
Gam in mgcv package is predicting negative values which should not be
the case despite all the predictors and response variables are positive.
Tried to use log link function but it did not help. Please help
S
> > Could I suggest that citation() be modified to include
> > the URL automatically?
> >
> I second this suggestion. I experienced similar case once as well.
Thanks for pointing that function out to me. When I run it, the URL is
included:
citation()
To cite R in publications use:
R Devel
> Hi,
>
> I am trying to run a code in Excel using VBA from R, using
> Rexcel.but everytime I am getting this error.
Error - 2147220502 in module Recel.Rserver
Error running expression
Eval(parse(text="setwd(\"H:\\
> I am using R-2.8.0 and R.2.7
> Please help.
>
>
>
> << OLE Object
Hi,
Browse[1]> d4
EVDO_Rev Session_Setup FCA bin counts
50 NA 0 5 1 1
51 NA 0 5 2 1
52 NA 0 5 3 1
53 NA 0 5 4 1
54 NA 0 5 5 1
55 NA
I want to form a 3x3 crosstabulation for the signs of two vectors (i.e.
Negative, Zero, Positive). The problem is that I am simulating the data so
for some iterations one of the categories is absent. Thus the resulting
table shrinks to 3x2. I want it to be 3x3 with zero column corresponding to
the
yes thank you! it is perfect.
I was using lmrob in package robustbase and it did not have that option in
the summary.
Laura
2008/11/13 Mark Difford <[EMAIL PROTECTED]>
>
> Hi Laura,
>
> >> I was searching for a way to compute robust R-square in R in order to
> get
> >> an
> >> information simila
Goodevening to everyone,
I am trying to create a clustered barplot for the following 2 variables who
have more cases but as an example i am giving you 10 items:
sex socio-economic status
1Â Â Â Â 1
2Â Â Â Â 2
2Â Â Â Â 2
2Â Â Â Â 3
2Â Â Â Â 2
1Â Â Â Â 4
1Â Â Â Â 1
1Â Â Â Â 1
2Â Â Â Â
Achim Zeileis wrote:
On Wed, 12 Nov 2008, Michael Friendly wrote:
In two recent in press articles where I cited R packages,
I was asked by the copy-editor to add a location or url to the
references. Could I suggest that citation() be modified to include
the URL automatically?
That's a good
This was AFAICS a bug introduced in 2.7.1 by
o plot(, which=5) uses more correct Cook's distance contours;
(fix to fix to PR#9316).
Thanks to Greg for the diagnosis.
Will be fixed in R-patched later today.
On Wed, 12 Nov 2008, Greg Snow wrote:
Just a clarification on one of my
> "LP" == Laura Poggio <[EMAIL PROTECTED]>
> on Thu, 13 Nov 2008 10:43:14 + writes:
LP> yes thank you! it is perfect.
LP> I was using lmrob in package robustbase and it did not have that option
in
LP> the summary.
Yes
lmRob() from "robust" is from a company which
I was searching for a way to compute robust R-square in R in order to get an
information similar to the "Proportion of variation in response(s) explained
by model(s)" computed by S-Plus. This post is dealing with that. Would be
possible to have some hints on how to calculate this parameter within R
THere is a separate list for R and the mac:
https://stat.ethz.ch/mailman/listinfo/r-sig-mac
On Thu, Nov 13, 2008 at 3:38 AM, Daniel Fernandez <[EMAIL PROTECTED]> wrote:
> Dear R help team,
> I am sorry to put urgent in the subject but to be honest I am in big
> troubles if I cannot make R to work
Hi Laura,
>> I was searching for a way to compute robust R-square in R in order to get
>> an
>> information similar to the "Proportion of variation in response(s)
>> explained
>> by model(s)" computed by S-Plus.
There are several options. I have had good results using wle.lm() in package
wle and
"Kurapati, Ravichandra \(Ravichandra\)"
<[EMAIL PROTECTED]> napsal dne 13.11.2008 11:37:58:
> But I observed that it was throwing an error.
Really. If what you send is all your data then the problem is that
EVDO_Rev is all NA. The problem is in your data, change your data and you
will get resu
Hi
[EMAIL PROTECTED] napsal dne 13.11.2008 10:51:39:
> Hi,
>
>
>
> Browse[1]> d4
>
>EVDO_Rev Session_Setup FCA bin counts
>
> 50 NA 0 5 1 1
>
> 51 NA 0 5 2 1
>
> 52 NA 0 5 3 1
>
> 53 NA
Thank you DAVID, GABOR, PHIL, WACEK, PETER and GREG
I tried:
Data2 <- data.frame(x=DATA$x, y = gsub('[^[:digit:]]', '', DATA$x), z=
gsub('[^[:alpha:]]', '', DATA$x))
and everything is ok
Wacek Kusnierczyk wrote:
>
> CE.KA wrote:
>>
>> x
>> 1 12F
>> 2 13 AD
>> 3 356PO
>> 4 1D
>> 5
Hello...
Which code are you using to perform the breakup into the three classes?
Can you be more specific on that?
Best regards,
Carlos J. Gil Bellosta
http://www.datanalytics.com
On Thu, 2008-11-13 at 09:57 +, Sohail wrote:
> I want to form a 3x3 crosstabulation for the signs of two vecto
Argyro Antaraki yahoo.gr> writes:
> I am trying to create a clustered barplot for the following 2 variables
> who have more cases but as an example i am giving you 10 items:
>
> sex socio-economic status
> 1Â Â Â Â 1
> 2Â Â Â Â 2
> 2Â Â Â Â 2
> 2Â Â Â Â 3
..
> where sex: 1:man, 2:woman
>
I am aware of the limits of the parameter R^2 in this case. However often it
is required for many different reasons. And it is helpful to have a function
that does it. The most important is to know the drawback of the"number", I
think.
Laura
2008/11/13 Martin Maechler <[EMAIL PROTECTED]>
>
Hi,
I would like to plot a matrix which comprises 3 columns.
So first column should be the x-axis.
Second the y-axis and the third one should be the z-axis.
I know that I can plot such data using the scatterplot3d, but I would like
to
have a surface plot like you can do with wireframe and persp.
Dear all,
I believe there is a problem with the corfdrci and the kendallfdrci
function from the GeneNT package. The 2nd screening, that is taking what
is not in the intersect of the CI and the user-defined MAS criterium
does not work properly on my dataset. It seems to work properly on the
da
Hi list,
I'm new to R and I'm currently using ROCR package.
Data in input look like this:
DIAGNOSIS SCORE
1 0.387945
1 0.50405
1 0.435667
1 0.358057
1 0.583512
1 0.387945
1 0.531795
1 0.527148
0 0.526397
0 0.372935
1 0.8610
On 11/13/2008 8:51 AM, mentor_ wrote:
Hi,
I would like to plot a matrix which comprises 3 columns.
So first column should be the x-axis.
Second the y-axis and the third one should be the z-axis.
I know that I can plot such data using the scatterplot3d, but I would like
to
have a surface plot li
Hi All,
I'm pretty new to R, so would really appreciate it if someone could point me
in the right direction on this problem.
I am trying to "unflatten" a table in R, and can't seem to find a function
or method to complete this task, (hopefully efficiently).
My data table is full of historical st
Kaliss wrote:
Hi list,
I'm new to R and I'm currently using ROCR package.
Data in input look like this:
DIAGNOSIS SCORE
1 0.387945
1 0.50405
1 0.435667
1 0.358057
1 0.583512
1 0.387945
1 0.531795
1 0.527148
0 0.526397
0 0.372935
"Deepayan Sarkar" <[EMAIL PROTECTED]> wrote on 13/11/2008
02:23:25:
> I will try to track down the reason for the odd results that we
> currently get if I have time, but it's not going to be a priority.
>
> -Deepayan
>
In any case, thanks.
Vincenzo
Ive got problems with preparing my own experiment data. I need to see what is
the format, so there is a need of extracting rda file to tab delimited txt.
After that I will probably have to get txt again to rda. How to do that ? I
can use load() option but no clue how to save it (either txt or rda)
Dear R-Help,
I need to convert a set of '.pdf' files into an equivalent set of
'.txt' files. This is so that i can do some text mining on the
content.
In the latest R-News letter (http://cran.r-project.org/doc/Rnews/
Rnews_2008-2.pdf), the package 'tm' for text mining is mentioned. In
that lovely
Hi,
two days ago I have asked the list if there is a better way to plot a matrix
using
the following plot command:
plot(matrix[1,], type="l")
for (i in 2:dim(matrix)[1]) {
lines(matrix[i,], type="l")
}
I have been told to use the matplot function, but unfortunately I matplot
does not plo
Dear all,
I am trying to make boxplots in one diagram but it is not working. Actually I
have a large dataset, around 35 water quality variables with different units
and of course in different ranges (some are orders of magnitude bigger than
others). I want to produce one diagram with all the b
We need to have some idea of what the data looks like and what code you have
been trying to make any sensible suggestions.
Can you supply a small example ?
--- On Thu, 11/13/08, Lathouri, Maria <[EMAIL PROTECTED]> wrote:
> From: Lathouri, Maria <[EMAIL PROTECTED]>
> Subject: [R] Boxplots w
Dear all,
I am using the nlme package and trying to modelate the variance-
covariance matrix of random effects.
However, I do not understand why do I get this error message:
> M4CompSym <- update(M4,random=A~pdCompSymm(~CED))
Error in pdMatrix.pdCompSymm(x) :
Cannot extract the matrix from
on 11/13/2008 09:49 AM Brigid Mooney wrote:
> Hi All,
>
> I'm pretty new to R, so would really appreciate it if someone could point me
> in the right direction on this problem.
>
> I am trying to "unflatten" a table in R, and can't seem to find a function
> or method to complete this task, (hopef
?write.table
?read.table
--- On Thu, 11/13/08, memento1984 <[EMAIL PROTECTED]> wrote:
> From: memento1984 <[EMAIL PROTECTED]>
> Subject: [R] How to get txt file out of .rda file and vice versa ?
> To: r-help@r-project.org
> Received: Thursday, November 13, 2008, 10:32 AM
> Ive got problems wi
on 11/13/2008 09:39 AM mentor_ wrote:
> Hi,
>
> two days ago I have asked the list if there is a better way to plot a matrix
> using
> the following plot command:
>
> plot(matrix[1,], type="l")
> for (i in 2:dim(matrix)[1]) {
>lines(matrix[i,], type="l")
> }
>
> I have been told to use t
Hi List,
If I have a date format as:
d <- "2001/1/1"
I can easily convert it to number by using as.Date(d).
But if I have d<-"1/1/2001", it does not work. Does anyone know how I can
convert it using pre-written function in R?
Regards,
Alireza
[[alternative HTML version deleted]]
___
Hi Frank,
Thank you for your answer.
In fact, I don't use this for clinical research practice.
I am currently testing several scoring methods and I'd like
to know which one is the most effective and which threshold
value I should apply to discriminate positives and negatives.
So, any idea for my
Hi,
I have a vector where an integer is followed by a randomly varying number of NA
elements. I want to replace the NA's with the preceding integer. Let me explain
with a simple example :
avec<- c(1,NA,NA,5,NA,3,4,NA,NA,NA,NA,9,NA,0,NA,NA)
I want to get the following vector from a :
bvec<-c(1,1,1
Dr. Alireza Zolfaghari wrote:
> Hi List,
> If I have a date format as:
> d <- "2001/1/1"
> I can easily convert it to number by using as.Date(d).
Yes. If it means "January the 1st" and not "1st of January", that is...
> But if I have d<-"1/1/2001", it does not work. Does anyone know how I can
>
See na.locf in the zoo package. It does exactly that.
On Thu, Nov 13, 2008 at 12:07 PM, ravi <[EMAIL PROTECTED]> wrote:
> Hi,
> I have a vector where an integer is followed by a randomly varying number of
> NA elements. I want to replace the NA's with the preceding integer. Let me
> explain wit
Hi Laura,
The fact that you had copied Prof. Ripley's response suggested to me that
this was so. Nevertheless, I think Martin Maechler was wise to emphasize the
problem, just in case.
Bye, Mark.
Laura POggio wrote:
>
> I am aware of the limits of the parameter R^2 in this case. However often
Hi Pierre-Jean,
Sensitivity (Se) and specificity (Sp) are calculated for cutoffs stored in
the "performance" x.values of your prediction for Se and Sp:
For example, let's generate the performance for Se and Sp:
sens <- performance(pred,"sens")
spec <- performance(pred,"spec")
Now, you can have a
Dear All,
Am I correct in believing that logical operators such as && and ||
cannot be used as part of
a 'subset' logical expression?
Example:
>foo = as.data.frame(matrix(c(1:9), nrow=3, ncol=3,
dimnames=list(c("r1","r2","r3"),c("c1","c2","c3"
>foo
c1 c2 c3
r1 1 4 7
r2 2 5 8
r3
2008/11/13 Reitsma, Rene - COB <[EMAIL PROTECTED]>:
>>subset(foo, c1 > 1 && c2 > 5)
> [1] c1 c2 c3
> <0 rows> (or 0-length row.names)
Too many ampersands:
> subset(foo, c1 > 1 & c2 > 5)
c1 c2 c3
r3 3 6 9
See help("&")
Barry
__
R-help@r-projec
I am having a problem installing and updating packages. It doesn't happen
with every package I try to install or update, but most of them.
I get this warning message:
Warning: unable to move temporary installation 'C:\Program
Files\R\R-2.8.0\library\file5f906952\rpart' to '[EMAIL PROTECTED] $2'
On 12/11/2008 7:59 AM, David Croll wrote:
>
> Hello dear R people!
>
>
> Several times it occurred to me that a function that uses a variable
name as a parameter would be helpful, but I did not find out how to
write such a function. I have experience in several programming
language, but I did no
[EMAIL PROTECTED] wrote:
Hi Frank,
Thank you for your answer.
In fact, I don't use this for clinical research practice.
I am currently testing several scoring methods and I'd like
to know which one is the most effective and which threshold
value I should apply to discriminate positives and neg
Hello,
I am working with version 2.7.2 on a PC and have 2.8.0 available, but have
not upgraded completely yet.
I have 6 species for which I've run 6 unique environmental management
scenarios against for comparison purposes. Each scenario is run for 36
years (1965 - 2000). Some of the species h
N. Lapidus wrote:
Hi Pierre-Jean,
Sensitivity (Se) and specificity (Sp) are calculated for cutoffs stored in
the "performance" x.values of your prediction for Se and Sp:
For example, let's generate the performance for Se and Sp:
sens <- performance(pred,"sens")
spec <- performance(pred,"spec")
On Thu, 13 Nov 2008, Sohail wrote:
I want to form a 3x3 crosstabulation for the signs of two vectors (i.e.
Negative, Zero, Positive). The problem is that I am simulating the data so
for some iterations one of the categories is absent. Thus the resulting
table shrinks to 3x2. I want it to be 3x3
Dear,
Im trying to run the R program, provided by Wang (2007), for the analysis of
nested frailty models. While running this program I get the following error:
Error in inherits(x, "ordered") : subscript out of bounds
Could someone tell me what the problem is in this case?
Thanks in advance.
If anyone can assist with this problem you have my great thanks:
I am trying to establish and plot confidence intervals on a bootstrapped
function. I have a more complicated function that has no problems with
determining the confidence intervals using the quantile command. This is
outside the b
Yes you can automate it. The exact way would depend on how your data
is structured and what you want the graphs to look like. It is hard
to say without an example of the data and at least the command that
your are using to create a single plot.
On Thu, Nov 13, 2008 at 2:16 PM, <[EMAIL PROTECTED
Hi again,
Yes you are right I should have included an example of the file
It looks like this after reading in the data and using tapply I arrive at
the following
> Alt5rimc_mean
A B C D E F
1965 8.423645 0.208 1.6296296 0.000 0.000
Dear R gurus:
Here is the following from Montgomery's Design and Analysis of
Experiments, 5th edition.
> str(rout1.df)
'data.frame': 16 obs. of 3 variables:
$ resp: num 18.2 18.9 12.9 14.4 27.2 24 22.4 22.5 15.9 14.5 ...
$ A : Factor w/ 2 levels "-1","1": 1 1 1 1 2 2 2 2 1 1 ...
$ B :
hello,
i have something like:
out<-list()
for(i in 1:n){
data<-gen(...) #fixed size data
out[[i]]<- fun(data)
}
> object.size(out[[1]])
6824
In principle 1 GB should allow
n = 1024^3/6824 = 157347?
i have about 2GB are not taken by other processes. however, I can see the
memory shrinking q
library(elliptic)
f <- function (x) {myintegrate(function(t){Im(as.complex(1+(x*t)i))},
1,2)-1}
uniroot(f,lower = -2, upper = 2)
> f <- function (x) {myintegrate(function(t){Im(as.complex(1+(x*t)i))},
> 1,2)-1}
Error: unexpected symbol in "f <- function (x)
{myintegrate(function(t){Im(as.comp
Edna Bell wrote:
Dear R gurus:
Here is the following from Montgomery's Design and Analysis of
Experiments, 5th edition.
str(rout1.df)
'data.frame': 16 obs. of 3 variables:
$ resp: num 18.2 18.9 12.9 14.4 27.2 24 22.4 22.5 15.9 14.5 ...
$ A : Factor w/ 2 levels "-1","1": 1 1 1 1 2 2 2
Please show us what you get when you execute:
str(spead)
typeof(spread)
class(spread)
It is doubtful that either of your theories is correct. It's more
likely that spread is a more complex object than you realize. We need
to know what spread looks like to the R interpreter in the environment
Please read and understand the response that I already sent you
on this issue.
Rolf Turner
On 14/11/2008, at 10:44 AM, Edna Bell wrote:
Dear R gurus:
Here is the following from Montgomery's Design and Analysis of
Experiments, 5th edition.
str(rout1.df)
'data.frame': 16 obs. of
on 11/13/2008 03:30 PM [EMAIL PROTECTED] wrote:
> Hi again,
>
> Yes you are right I should have included an example of the file
>
> It looks like this after reading in the data and using tapply I arrive at
> the following
>
>> Alt5rimc_mean
> A B C D
Convert it to long form and plot like this (or look at ggplot2 for
another approach):
DF <- structure(list(A = c(8.423645, 30.881773, 29.783251, 17.827586,
27.507389, 17.46798, 8.167488, 32.502463, 44.610837), B = c(0.208,
27.344, 15.824, 31.776, 30.52, 27.232, 0, 19.736, 11.208), C = c(1.6296296,
You need to load the foreign package first.
library(foreign)
- Original Message -
From: "b g" <[EMAIL PROTECTED]>
To:
Sent: Monday, November 10, 2008 12:28 PM
Subject: [R] read SAS file
Looked at a lot of documentation and listserv postings and still can't solve
this problem. I ne
Dear list,
When I try to use lines.formula with subset and another argument I
get an error.
e.g.
x<-1:5
y<-c(1,3,NA,2,5)
plot(y~x, type="n") # set up frame
lines(y~x, subset=!is.na(y)) # works OK
lines(y~x, type="o", col="blue") # works OK
# but
lines(y~x, subset=!is.na(y), col="red") # gi
Hello,
I have a stripplot with 200+ labels i.e i have network connections.
The y-axis are the server port numbers and am graphing the number of
packets in a connection. Roughly
stripplot(totpacks~portnum,data=network)
Due to the large number of server ports in my dataset
I am running R.exe/Rterm.exe 2.7.1 under Emacs 22.1 under the DOS
shell under Windows XP. (My Cygwin installation is broken!) I am not
currently using ESS:
> R.exe --vanilla
R version 2.7.1 ...
...
> 1+1
2
This works fine until I hit an R error -- pretty much any error --
when th
Is it possible to plot multiple polygons on radial.plot(plotrix)? The
new=FALSE argument (as used to add additional points or lines with the plot
function to an existing plot) doesn't appear to work. I would like to
overlay the outlines of multiple polygons of different colors on the same
radial pl
On Thu, Nov 13, 2008 at 4:57 PM, Saptarshi Guha
<[EMAIL PROTECTED]> wrote:
> Hello,
>I have a stripplot with 200+ labels i.e i have network connections.
> The y-axis are the server port numbers and am graphing the number of
> packets in a connection. Roughly
>stripplot(totpacks~port
Please pardon an extremely naive question. I see related earlier
posts, but no responses which answer my particular question. In
general, I'm very confused about how to do variance decomposition with
random and mixed effects. Pointers to good tutorials or texts would
be greatly appreciated.
To
On Monday 10 November 2008 01:23:30 pm Prof Brian Ripley wrote:
> See ?Sys.timezone, which help.search("timezone") points you to.
>
> On Mon, 10 Nov 2008, stephen sefick wrote:
> > I have looked at
> > ?as.POSIXct
>
> That also links to Sys.timezone from the 'tz' argument.
>
> > ?POSIXct
> > and ma
Rather than estimating the variance components via the
method-of-moments estimators, have a look at the 'nlme' and 'lme4'
packages, which provide likelihood-based tools for estimating random
and mixed models (and in the case of nlme, gls models too).
Advantages of the likelihood-based approaches in
The notion that you can just add or subtract 0.03 from estimate is
obviously incorrect.
Presuming you meant to call you lower bound q05 and the upper bound
q95, the numbers I get are in your 10,000 iteration loop are 4.97 and
5.18 (around a mean of 5.08). So roughly a .1 swing on each side
> -Original Message-
> From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
> On Behalf Of John Field
> Sent: Thursday, November 13, 2008 4:41 PM
> To: r-help@r-project.org
> Subject: [R] lines.formula with subset
>
> Dear list,
>
> When I try to use lines.formula with subset and anot
On Thu, Nov 13, 2008 at 9:10 PM, J Dougherty <[EMAIL PROTECTED]> wrote:
> On Monday 10 November 2008 01:23:30 pm Prof Brian Ripley wrote:
>> See ?Sys.timezone, which help.search("timezone") points you to.
>>
>> On Mon, 10 Nov 2008, stephen sefick wrote:
>> > I have looked at
>> > ?as.POSIXct
>>
>>
Goodmorning to everyone,
I am trying to create a clustered
barplot for the following 2 variables who have more cases but as an
example i am giving you 10 items:
sex socio-economic status
1Â Â Â Â 1
2Â Â Â Â 2
2Â Â Â Â 2
2Â Â Â Â 3
2Â Â Â Â 2
1Â Â Â Â 4
1Â Â Â Â 1
1Â Â Â Â 1
2Â Â Â Â 3
Hi dear list,
I'm scrambling with this:
i have a p-dimensional data frame for which i want to have a pairs() plot, with
ellipse added in each lower.panel.
I want the ellipse's parameter (mu, sigma) to be extracted from the appropriate
entries in a pre-computed matrix ("dab" in the code below)
Hello all,
I'm trying to use text symbols in R, such are greek letters and
mathematical operators. Do you know where I can find a list telling me
how to use all those symbols? I'd like to have them displayed in my
graphs.
Best,
Guillaume
__
R-help
Thanks. I did the search before I posted and found those threads.
However, it does not seem to do what I want. All I want to do is
estimate the sample size for a point estimate, not do a GLM. I just
want the mean within a margin of error, and to a given CI.
I've tried writing some code to do
This is strange the first function that I am doing this for is saying in the
typeof() is a double for some reason this is coming up as a list.
> str(spread)
'data.frame': 18 obs. of 1 variable:
$ bootdataframe.age: num 1.46e-09 1.67e-09 1.46e-09 1.45e-09 1.37e-09 ...
> typeof(spread)
[1] "
Hi
Gam in mgcv package is predicting negative values which should not be
the case despite all the predictors and response variables are positive.
Tried to use log link function but it did not help. Please help
sunil
--
View this message in context:
http://www.nabble.com/negative-prediction
Thank you all for your help about the encoding of my yahoo.gr.
I am sending you now the original data. Hope tha you can read them.
I am trying to create a clustered barplot for the following 2 variables:
> h$rocbiopsy1
 [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0
On Thu, 2008-11-13 at 21:32 -0800, sunil_POUDEL wrote:
> Hi
>
> Gam in mgcv package is predicting negative values which should not be
> the case despite all the predictors and response variables are positive.
> Tried to use log link function but it did not help. Please help
>
>
> sunil
Pleas
86 matches
Mail list logo