similar to: Sorting data from one column with strings

Hi I have a small problem, I'm new in using R, so I hope you can help me... I'm running a logistic regression model and want to do a nice plot. The plot I have made is done a plot with the command jitter: plot(jitter(overto$age[overto$sex=="F"]),jitter(overto$neg.pos[overto$sex=="F"]),xlab="age",ylab="neg and pos") my responsvariable is positive

R-help, I wish to replace NULL elements(or missing) in the following list : > z2 $cod mean sd 62.56190 12.65452 $haddock mean sd 36.61490 11.50365 $ling mean sd 86.17949 20.43587 $saithe mean sd 50.275847 5.453606 $whiting NULL $"norway pout" mean sd 13.739623 1.393104 $"great silver smelt" mean

I have 5 GAMs ( model1, model2, model3, model4 and model5) Before I use some data X(predictor -January to June data) to form a equation and calculate the expected value of Y (predictand -January to June). After variable selection, GAMs (Model 1)were bulit up! R-square :0.40 NOW, I want to use new X'( predictor -July - December data) and put into Model 1, then get the expected value of Y'

Dear list, for a data structure like in df: set.seed(100) Treatment<-rep(c("Nitrogen","Carbon", "Sulfur"),each=9) week<-rep(c(1,5,9),3,each=3) genes<-rep(c("18s", "EF1b", "NR"),9) copies<-rnorm(27, 1000000,400000) df<-data.frame(Treatment,week,genes,copies) i wrote this code for a xyplot: library(lattice)

Hi, I have a data frame with columns as follows: tree species (independent variable) and several dependent variables (e.g. carbon, nitrogen, phosphorus). Each row represents one tree sample. Some tree samples are unique species in the data frame while other species were replicated (i.e. rows 1,2,3 may be identical for the "tree species" column but have different values for the

Hi, Is there a function similar to excel's hlookup in R ? Thanks, Silje

I don't know if there is a way to use the scale relation free argument in ggplot2 like in lattice. I have a feeling that there is not, but I would like to make a plea for this feature. It would be nice to be able to plot Total Inorganic Nitrogen Total Phosphorus and the ratio of the two- the numbers on the axis are not related, but the previous two are surely related to the last (this ratio

Hi, I am using the mgcv package (version 1.7-22.) running the model works fine, but when I want to have a plot with residuals I get an error. fit29<-gam(IV~s(G3)+s(V3)+factor(AAR)+s(D3)+s(RUTE,bs="re"),data=subsf,gamma=1.4,method="ML") plot(fit29,residuals=T) Error in X[, first:last] %*% object$coefficients[first:last] : non-conformable arguments does some one know

Hi! I try to explain the efffect of (1) forest where i took samples's soils (* Lugar*: categorical variable with three levels), (2) nitrogen addition treatments (*Tra*: categorical variable with two levels) on total carbon concentration's soil samples (*C: *continue* *variable) during four months of sampling (*Time:* categorical and ordered variable with four levels). I fitted the

Hi, I am using the gam function in the mgcv package, I have random effects in my model (bs="re") this has worked fine, but after I updated the mgcv package to version 1.7-19 I recive an error message when I run the model. > fit1<-gam(IV~s(RUTE,bs="re")+s(T13)+s(H40)+factor(AAR)+s(V3)+s(G1)+s(H1)+s(V1)+factor(LEDD),data=data5,method="ML") > summary.gam(fit1)

hi, I have a question regarding cp values in rpart(). When I use plotcp() I get a figure with cp values on the x-axsis, but then I use printcp() the cp values in that list are different from the values in the figure by plotcp(). Does someone know why? Silje [[alternative HTML version deleted]]

We have upgraded a box from 6.2 to 6.3-RELEASE. Afterwards the box does not recognize its ST Lab I-160 serial card with Netmos 9845 Chipset. It worked flawlessly on 6.2 with puc(4) driver. >From dmesg: pci1: <simple comms, UART> at device 8.0 (no driver attached) ~# pciconf -l -v | grep -B 4 UART none2@pci1:8:0: class=0x070002 card=0x00041000 chip=0x98459710 rev=0x01 hdr=0x00

I am trying to figure out the correct syntax for the CLI command "SIP SHOW CHANNELS". I have tried SIP SHOW CHANNELS SIP/200 and I've even tried to do this when a call is connected such as: -- Zap/15-1 is ringing -- Zap/15-1 answered SIP/206-4299 asterisk*CLI> sip show channel SIP/206-4299 No such SIP Call ID 'SIP/206-4299' I always get the "No such SIP

What test do I use to determine if there is a correlation between a discrete variable and a continuous variable? For example - I have water quality ratings for streams (excellent, good, fair, poor) and a corresponding nitrogen concentration for each rating. I want to know if the the ratings correlate with the concentration of nitrogen in the stream. Help?

Hi all, I am trying to convert a data set of physician death codings (each individual's cause of death is coded by multiple physicians) from long to wide format, but the "reshape" function doesn't seem to work because it requires a "time" variable to identify the sequence among the repeated observations within individuals. My data set has no order, and different

I have no problem setting up my mixed model, or performing anova or lsmeans on my model?s outputs. However, performing lsd mean separation is giving me fits. So I do not have a problem when using two-way anova model. When using the code: fit.yield.add <- lm(data = ryzup, Yield ~ Rep + Nitrogen + Treatment) LSD.test(fit.yield.add, trt = "Nitrogen", alpha = 0.1, console = TRUE)

I have a question regarding NA in randomForest (in R). I have a dataset which include both numerical and non-numerical variables, and the data includes some NA. I tried to use na.roughfix but then i get an error message "na.roughfix only works for numeric or factor". I also tried rfImpute but this does not work either because I have some NA in my response variable. Does anyone have som

> On 30 Sep 2017, at 14:22 , Robert Baer <rbaer at atsu.edu> wrote: > > > > On 9/29/2017 3:37 PM, Rolf Turner wrote: >> On 30/09/17 07:45, JLucke at ria.buffalo.edu wrote: >> >> <SNIP> >> >>> >>> The conceptual paradigm for R is only marginally commensurate with >>> that of >>> standard statistical

Hello, I have two different dataset, and wanted to join the two. For example I have a table of codes of cities, and other with with the codes of travels, [source for the destine]. what he wanted was to have a new data.frame with all the information city<-data.frame(city="Barcelona",cod=1) city<-rbind(city,data.frame(city="Madrid",cod=2))

Hi, I am using the: pcf.aov<-aov(meas~op+part, data=pcf.ex2), command to perform a two way ANOVA. When I save the: sumpcf<- summary.aov(pcf.aov), result of the summary.aov command in a variable I need to access the individual pieces of information in the summary. The summary appears to be a list and I am having a hard time finding a way to get at the information in the summary.